{"id":"827f3082-57da-4f56-919b-980428acbc8a","arxiv_id":"1908.06298","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In N-body simulations of S-stars near Sgr A*, inner planets of solar-like systems stay bound more often than outer planets, and some escaped stars may swap their planetary systems.","lead":"This short paper uses N-body simulations with post-Newtonian corrections to ask what happens to hypothetical solar-like planetary systems when the S-stars near the Milky Way's central black hole interact with it. It reports preliminary survival rates, finding inner planets stay bound to their host stars more often than outer ones and that some escaping S-stars may swap their planets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 2 survival rates rest on a single unperturbed near-circular solar twin; the 'inner planets survive better' claim needs an ensemble check.","rationale":"The reader's weakest assumption focuses on whether real S-stars actually host planetary systems, which is a legitimate external premise but is transparently labeled 'putative' and is not needed for the paper's stated aim of studying the destiny of hypothetical systems. The stronger, internal risk is that the single deterministic realization per star and inclination leaves the headline percentages underspecified: no error bars, no phase averaging, and no convergence check over initial conditions. That is precisely the kind of missing support the rubric asks to flag; the paper itself calls the results 'preliminary', which supports CONDITIONAL rather than ACCEPT. The concern is not a fatal flaw because the ordering by semi-major axis is physically plausible and would likely survive modest averaging, but the reported 100% Mercury retention could be sensitive to initial phases and to the adopted circular mock solar systems. The suggested random-phase ensemble is cheap relative to the stated 48-hour runs and would settle whether the central claim is a dynamical property or an initialization artefact. The proposed verdict is CONDITIONAL with the ensemble check as the explicit condition, matching the reader's verdict even though the load-bearing concern differs.","tokens_in":3584,"tokens_out":1277,"duration_ms":12549,"concrete_test":"Repeat the simulations for a subset of, say, 10 S-stars and two inclinations with 20 random realizations each, varying the mean anomaly of every planet and the longitude of the ascending node while keeping the same semi-major axes, masses, and the host star orbits. If the average bound fraction for Mercury–Mars drops by more than ~10 percentage points from Table 2 or the inner/outer ordering reverses, then the reported survival rates are phase-locked artifacts; otherwise the central claim is robust to this concern.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that Mercury through Mars remain bound at 61–100% versus 16–49% for outer giants is computed from only one realization per inclination for each S-star: the paper gives no multiplicity over random planetary phases, eccentricities, or mutual inclinations. Each of the approximately 40 stars is assigned a copy of the Solar System, and the planet orbits are initialized as nearly circular with all planets coplanar and with mutual separations identical to ours. Since strong encounters with Sgr A* act on orbital phase through resonant and impulsive kicks, a single phase choice per star per inclination could systematically misestimate the bound fraction by selecting configurations that happen to be protected. The stated spread between inner and outer bound rates is then not yet demonstrated to be a property of the dynamics rather than of the chosen initial conditions. This is an internal robustness concern, not a question of the migration premise: the conclusion of Section 3 and Table 2 would need to survive averaging over randomized planetary phases and eccentricities to be the physical statement the abstract makes.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper hypothesizes that the S-stars in the Galactic Center may host planetary systems and investigates their dynamical fate under repeated close encounters with Sgr A*. Using a regularized N-body code with post-Newtonian corrections to order 2.5 and the observed orbital elements of about 40 S-stars, the authors attach a solar-like planetary system to each star (7 planets for late-type stars, 3 outer planets for early-type stars) and vary the inclination of the planetary orbital plane relative to the stellar orbit from 0 to 180 degrees. They report preliminary bound fractions in Table 2, concluding that inner planets (Mercury through Mars) remain bound at higher rates than outer planets (Jupiter through Uranus) for late-type hosts, and that some early-type stars escape the Galactic Center with their planetary systems swapped by a nearby star. The paper is explicitly a preliminary proceedings contribution.","tokens_in":3865,"tokens_out":4854,"duration_ms":47858,"significance":"If the central claim holds, this is a first estimate of the survival probability of putative planetary systems around S-stars, with direct relevance to the possibility of planets in the Galactic Center. The strengths of the work are that the simulations use a well-tested regularized N-body code with PN corrections, adopt observed S-star orbits as initial conditions, and contain no parameter fitted to the survival fractions—the bound rates emerge from the assumed initial conditions. However, the result is conditional on a speculative planetary architecture and on the migration premise, and the reported rates lack ensemble statistics and sensitivity tests. The paper therefore constitutes an interesting preliminary finding rather than a fully established result.","major_comments":[{"comment":"The central claim that inner planets remain bound at higher rates (61–100% vs 16–49%) rests on a single realization per star per inclination: each of the ~40 S-stars is seeded with an identical copy of the Solar System with fixed orbital phases, eccentricities, and mutual inclinations, and the inclination of the planetary plane is the only varied parameter. Because the tidal perturbation from Sgr A* depends on planetary phase at pericenter, the reported dichotomy could be an artifact of the chosen initial phases rather than a generic dynamical property. Please rerun the simulations for an ensemble of random mean anomalies and small eccentricities for each star and inclination, and report mean bound fractions with standard deviations or confidence intervals. Without such averaging, the numbers in Table 2 overstate the precision (e.g., '61.85%' from a small sample) and the qualitative conclusion is not yet robust.","section":"Section 3, Table 2"},{"comment":"The simulated planetary systems are all exact solar twins, with planet masses, semi-major axes, and eccentricities fixed to Solar System values, and for early-type stars an ad-hoc choice to include only Jupiter, Saturn, and Uranus ('Since it is less likely that massive early-type stars could harbor planets so close'). This architecture directly determines the columns of Table 2. The paper should test sensitivity to plausible alternative architectures (e.g., different planet multiplicities, semi-major axis distributions, or an extra ice-giant belt) and to the uncertain stellar masses (8–14 M⊙ and 0.5–2 M⊙) before claiming a general result for S-star planetary systems.","section":"Section 2"},{"comment":"The premise that S-stars 'probably still keep their planetary systems throughout their voyage' is asserted without support. Given that the preferred migration scenario is binary disruption (Hills 1988), the retained companion may be ejected with a kick that could strip or destabilize its planets; the simulations instead start with the planets already in place around their host on the observed S-star orbits. Please either add a justification or test of planetary survival during the migration event, or explicitly reframe the conclusion as applying to a hypothetical population of stars that arrived with their planetary systems intact. This also affects the physical interpretation of the bound fractions in Table 2.","section":"Abstract and Section 1"}],"minor_comments":[{"comment":"The sentence 'Our simulations shows that the innermost planets i.e, Mercury, Venus, Earth and Mars...' contains a subject-verb agreement error ('simulations shows') and a punctuation error ('i.e,' should be 'i.e.,').","section":"Section 3"},{"comment":"The caption states 'the inclinations in the range of 0° − 90°' while the text (Section 2) says the inclination is varied over 0° − 180°; clarify whether Table 2 is a subset of the runs and why the 90°–180° results are omitted.","section":"Table 2 caption"},{"comment":"Figure 2 is described only by the caption; the text should state what each plotted curve represents and how the statement 'six out of seven planets remain bound to the star S89' is inferred from the plotted orbital elements. The current description is too vague for the reader to verify the result.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings contribution with a modest but interesting novel result. The recommendation of major revision is intended for a full-length version: the core simulation campaign is appropriate, but the authors need to add ensemble statistics and sensitivity runs before the central claim can be considered robust. There is no issue with circularity or improper citation; the main concern is the lack of averaging over planetary initial conditions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you care about planetary survival near Sgr A*. This proceedings paper reports N-body simulations that assign Solar-System-like planetary systems to 40 S-stars, integrate for 4000 yr with post-Newtonian terms, and report bound fractions. That is new: prior work (Trani et al. 2016) treated single planets, not full systems. The qualitative result—inner planets survive better than outer giants—comes straight out of Table 2 and is physically reasonable: tighter orbits are harder to unbind in a tidal encounter. The use of actual S-star orbital elements from Gillessen et al. is a plus, and the ARGdfcode regularization gives confidence the close encounters are handled properly.\n\nThe soft spots are real but not fatal. First, the migration premise (S-stars carry planets into the GC) is stated without evidence. The paper is upfront that these are \"putative\" systems, so it is framed as a hypothetical scenario rather than a claim about actual detections. That is acceptable for a first exploration, though it limits the astrophysical interpretation.\n\nSecond, and more important: Table 2 presents one realization per inclination per star. The planetary orbits are initialized as near-circular, coplanar copies of the Solar System, with no randomization of phases, eccentricities, or mutual inclinations. Since the impulsive kick from a close Sgr A* passage depends on orbital phase, a single phase choice can systematically skew the bound fraction. The 61–100% inner versus 16–49% outer spread could be a property of the dynamics, or it could be an artifact of the chosen initial conditions. The paper needs an ensemble average over random phases and eccentricities before the central claim is solid. This is an internal robustness issue, not a question of the migration premise.\n\nThere are also no error bars, no convergence tests, and the masses of the S-stars are set to representative values (12 M☉, 1.5 M☉) rather than drawn from the measured distribution. For a conference proceedings, this is understandable, but it means the results should be treated as preliminary, as the authors themselves say. The citation pattern is fine—they engage the relevant S-star and planetary dynamics literature without padding.\n\nBottom line: a modest but real contribution, with a clear and testable claim that currently rests on a single-per-configuration numerical experiment. The paper deserves a serious referee if expanded into a full journal article with ensemble runs and sensitivity tests. I would bring it to the reading group only if someone is actively working on GC planetary dynamics; otherwise it's a maybe. I wouldn't cite it in my own work until the ensemble check is done. Recommendation: engage with it as a legitimate preliminary study, send it to peer review in a venue that allows short contributions, but ask for the ensemble robustness analysis before accepting.","headline":"A modest but genuine first step: simulated multi-planet survival around S-stars, whose main trend is plausible but rests on a single realization per configuration.","tokens_in":4254,"tokens_out":2267,"would_cite":false,"duration_ms":20986,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations of S-stars with planetary systems show that inner planets remain bound to their hosts during close encounters with Sgr A*, while outer planets are preferentially stripped.","keywords":["Sgr A*","S-stars","planetary systems","N-body simulation","post-Newtonian dynamics","Galactic Center","planet survival","orbital stability"],"falsifier":"A transit or radial-velocity survey sensitive to close-in planets around S-stars that finds none would undercut the premise that planetary systems are present, though it would not by itself disprove the dynamics. A direct numerical check is to rerun the same 40-star simulation with randomly sampled planetary mean anomalies and mutual inclinations; if the Mercury-type bound fraction drops well below the reported 100%, the result depends on the chosen initial conditions rather than on the encounter dynamics.","tokens_in":3371,"feed_emoji":"🪐","tokens_out":6496,"duration_ms":59987,"temperature":0.7,"pith_summary":"This paper asks whether planets can survive when their host stars make close passes by Sagittarius A*, the supermassive black hole at the center of the Milky Way. It runs N-body simulations of the 40 innermost S-stars, each given a miniature Solar System, with relativistic corrections in the equations of motion. The central result is that the innermost planets—Mercury, Venus, Earth, and Mars analogues—remain bound to their host stars far more often than the outer planets, and Mercury-like planets around late-type stars stay bound in 100% of the runs. If true, planetary systems can persist in the immediate neighborhood of a supermassive black hole, at least in their inner parts.","feed_headline":"Inner planets survive close passes by Sgr A*","feed_subtitle":"N-body runs find Mercury-like planets stay bound to their S-stars; outer planets are stripped away.","key_machinery":"The load-bearing mechanism is a regularized N-body integrator with post-Newtonian corrections up to order 2.5, which lets the simulation follow close encounters with the black hole without numerical singularities. The 40 S-stars are initialized on their observed orbits, and each host star is assigned a planetary system patterned on the Solar System: seven planets for late-type stars and the three outer planets for early-type stars. The planetary system is inclined relative to the stellar orbit by angles from 0 to 180 degrees, and the quantity that carries the argument is the fraction of planets still bound to their host stars at the end of the simulation.","core_discovery":"The paper claims that during close encounters with Sgr A*, the innermost planets of a solar-like system are preferentially retained. In its simulations, Mercury-type planets around late-type S-stars remain bound in every case, Venus- and Earth-type planets in 72.5% of cases, and Mars-type in 61.85%; for early-type stars, which start with only Jupiter, Saturn, and Uranus analogues, the bound fractions are 75.62%, 47.5%, and 35.31%, respectively. The paper also reports that some early-type stars escape the Galactic Center and that their planetary systems can be swapped to a nearby star during the encounter.","pith_inferences":["The simulations do not include planet-planet scattering, stellar evolution, or the mechanism that originally placed planets around these stars; any of these could change the retention fractions, so the 100% Mercury survival should be read as a dynamical statement, not a formation statement.","A natural testable extension is to vary the planetary architecture, for example using hot Jupiters on eccentric orbits or multiple planets with mutual inclinations, and check whether inner-planet survival remains the robust pattern.","If the swap of planetary systems during encounters is common, the Galactic Center may be a place where planets change host stars, producing systems with no analog in the field.","The survival hierarchy inner > outer suggests that any planetary system discovered around an S-star is likely a compact, multi-planet system with its outermost planets missing."],"forward_implications":["Planetary systems around S-stars, if present, should arrive at the Galactic Center heavily truncated: inner planets survive, outer planets are stripped.","The stripped outer planets become free-floating planets wandering in the Galactic Center region.","Early-type S-stars are expected to lose most or all of their planetary systems, and some may exchange planets with other stars during close encounters.","Observational searches for planets around S-stars should concentrate on close-in orbits, where the survival probability is highest.","The simulated bound fractions give a quantitative prior for interpreting future detections or non-detections of planets in the S-star cluster."],"supporting_citations":[{"why":"Supplies the observed S-star orbits used as initial conditions for the simulations.","marker":"Gillessen et al. 2009"},{"why":"Provides the updated S-star sample and orbital data.","marker":"Gillessen et al. 2017"},{"why":"Gives masses of early-type S-stars used to assign host-star masses.","marker":"Habibi et al. 2017"},{"why":"Provides ages and initial masses for the late-type S-stars.","marker":"Habibi et al. 2019"},{"why":"Presents the N-body code with dynamical friction and external potentials used for the runs.","marker":"Arca-Sedda & Capuzzo-Dolcetta 2019"},{"why":"Provides the regularized N-body integration method underlying the code.","marker":"Mikkola & Tanikawa 1999"},{"why":"Extends the integrator with post-Newtonian corrections up to order 2.5.","marker":"Mikkola & Merritt 2008"},{"why":"Proposes the binary-disruption migration model that motivates why S-stars would carry planets.","marker":"Hills 1988"},{"why":"Supports the migration scenario by which S-stars arrived at the Galactic Center.","marker":"Antonini & Merritt 2013"}],"fun_headline_variants":["Mercury-like planets stick with S-stars near Sgr A*","Sgr A* close passes spare inner planets, strip outer ones","Outer planets lost, inner planets survive Sgr A* pass","Simulations: Sgr A* strips outer planets, spares inner ones"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole exercise assumes that the S-stars were born elsewhere, migrated to the Galactic Center, and kept their planetary systems during the trip; if S-stars arrived planet-free, the survival rates describe a hypothetical rather than an actual population.","fun_headline_variants_meta":{"raw":{"variants":["Mercury-like planets stick with S-stars near Sgr A*","Sgr A* close passes spare inner planets, strip outer ones","Outer planets lost, inner planets survive Sgr A* pass","Simulations: Sgr A* strips outer planets, spares inner ones"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001143,"raw_usage":{"total_tokens":4633,"prompt_tokens":723,"completion_tokens":3910,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":339,"completion_tokens_details":{"reasoning_tokens":3832}},"tokens_in":339,"tokens_out":3910,"duration_ms":27862,"temperature":1.0,"reasoning_tokens":3832,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:49:29.296134+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A transit or radial-velocity survey sensitive to close-in planets around S-stars that finds none would undercut the premise that planetary systems are present, though it would not by itself disprove the dynamics. A direct numerical check is to rerun the same 40-star simulation with randomly sampled planetary mean anomalies and mutual inclinations; if the Mercury-type bound fraction drops well below the reported 100%, the result depends on the chosen initial conditions rather than on the encounter dynamics.","supporting_citations":[{"cited_title":"2019, ApJL, 872, L15","cited_arxiv_id":null,"evidence_quote":"Provides ages and initial masses for the late-type S-stars."},{"cited_title":"2019, MNRAS, 483, 152","cited_arxiv_id":null,"evidence_quote":"Presents the N-body code with dynamical friction and external potentials used for the runs."},{"cited_title":"1999, MNRAS, 310, 745","cited_arxiv_id":null,"evidence_quote":"Provides the regularized N-body integration method underlying the code."},{"cited_title":"2008, AJ, 135, 2398","cited_arxiv_id":null,"evidence_quote":"Extends the integrator with post-Newtonian corrections up to order 2.5."},{"cited_title":"2013, ApJL, 763, L10","cited_arxiv_id":null,"evidence_quote":"Supports the migration scenario by which S-stars arrived at the Galactic Center."}],"review_version":1}