{"id":"1f1766ed-1f09-4dc0-ab6e-e690a06e812f","arxiv_id":"2505.24093","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Very wide-orbit planets form naturally when planet-planet scattering flings a planet outward and a nearby stellar flyby lifts its perihelion, with an estimated occurrence of about 10^-3 per star.","lead":"Planets scattered onto wide orbits by instabilities in their own planetary systems can be rescued by a chance stellar flyby while the host star is still in its birth cluster, landing on stable, eccentric orbits hundreds of au out. The authors estimate such wide-orbit planets form at a rate near one in a thousand stars, and give odds for a Planet Nine-like body in our Solar System.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Trapped-orbit definition may include planets that do not survive Gyr-timescale evolution; if the long-term survival fraction is lower than asserted, the quoted efficiencies and the 10^-3 occurrence claim lose support.","rationale":"The reader's weakest assumption is the long-term stability of orbits counted as successes, and my independent reading arrives at the same load-bearing concern. The mechanism itself is plausible and the simulation campaign is large and carefully reported, so the paper deserves a conditional verdict rather than rejection. The decisive issue is whether the q > 150 au inclusion criterion, applied at 10-20 Myr, is a reliable proxy for Gyr-scale survival. The paper's own exclusion of q < 35 au and q = 35-150 au high-e orbits shows that the authors recognize a stability boundary, but they do not test the boundary they adopt. The claimed occurrence rate of 10^-3 per star is built directly on efficiencies defined with this untested boundary, and the Solar System-specific 40% probability additionally depends on both instability events occurring before cluster dispersal, which the cited 10-100 Myr constraint does not guarantee. A single targeted long-term integration study would settle whether the boundary is adequate; until then the conditional verdict is appropriate, and no change to the reader's recommendation is needed.","tokens_in":27474,"tokens_out":3941,"duration_ms":46487,"concrete_test":"Take the final snapshots of all trapped orbits counted as successes in the nominal Solar System ice-giant formation and gas-giant instability runs (Rc = 40 kau, Nstar = 200) and continue them for 4.5 Gyr after removing cluster members, adding the Galactic tide (rho0 = 0.15 Msun pc^-3) and a stochastic field-star encounter model. Recompute epsilon counting only planets that remain bound with q > 150 au, a < 10^4 au, and e < 0.9 at 4.5 Gyr. If the survival fraction is below 50%, or if the surviving q-distribution is depleted below ~300 au, the quoted efficiencies and the 10^-3 occurrence claim require downward revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper defines trapping efficiency as epsilon = Nwide/(Nwide+Nejec), with Nwide counting planets at simulation end (10-20 Myr) that satisfy q > 150 au, a < 10^4 au, e < 0.9 (Main text, trapping-efficiency definition). The central claim that these are durable wide-orbit planets rests on the sentence 'most trapped wide-orbit planets survive for tens of billions of years,' but no simulation in the paper runs beyond 20 Myr. The q > 150 au boundary is motivated only by excluding orbits that 'may not maintain long-term stability'; the same logic applies to the thousands of counted planets with q between 150 and ~550 au visible in Fig. 3 panels b and c. Many of these have a approaching 10^4 au and e near 0.9, where Galactic tides, passing field stars, and secular forcing from inner giant planets can alter q on Gyr timescales. If the true stable-perihelion boundary is larger, the quoted 5-10% Solar System efficiencies and the 1-5% exoplanet efficiencies are overestimates, and the abstract's 'at least 10^-3 per star' would need to be treated as an optimistic value rather than a lower bound. The Solar System probability also assumes both instability events occurred during the ~3-10 Myr cluster phase; the final giant-planet instability is constrained to 10-100 Myr, and if it happened after cluster dispersal the 40% combined figure should be reduced to the single-event estimate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that very wide-orbit planets (semimajor axes of a few hundred to 10^4 au) can be produced when a planet scattered during a dynamical instability receives a gravitational kick from a passing star or from the cluster potential, raising its pericenter and decoupling it from the inner planetary system. The authors report N-body simulations of several classes of planetary systems embedded in model stellar clusters: the Solar System's ice-giant formation phase, the final giant-planet instability, extrasolar gas-giant and ice-giant instability systems, and circumbinary gas-giant systems. They define a trapping efficiency and measure values of about 5-10% for Solar System-like cases, 1-5% for extrasolar gas-giant systems, and lower values for ice-giant-only and circumbinary cases. From this they estimate that wide-orbit planets occur at least 10^-3 per star and discuss implications for Planet Nine and for the connection between wide-orbit and free-floating planets.","tokens_in":27658,"tokens_out":3921,"duration_ms":41391,"significance":"If the result holds, the paper offers a concrete and physically plausible formation channel for a population that is otherwise hard to explain, and it makes falsifiable predictions: wide-orbit planets should be preferentially found around stars that host gas giants, hence around metal-rich, somewhat massive stars. The study is quantitatively transparent: it reports at least 100-1000 realizations per configuration, uses binomial error bars, explicitly states exclusion rules for orbits that may not be stable, and flags exploratory runs (ice-giant-only and circumbinary cases) as observationally unconstrained. The trapping efficiencies are measured rather than fitted to the target occurrence rate. The main weaknesses are the absence of long-term integrations supporting the claim that the counted orbits survive for Gyr, and the timing assumptions that underpin the combined Solar System probability.","major_comments":[{"comment":"The central claim that the counted wide-orbit planets are durable rests on the sentence \"most trapped wide-orbit planets survive for tens of billions of years,\" but no simulation in the paper runs beyond 10-20 Myr. The efficiency epsilon = Nwide/(Nwide+Nejec) counts orbits at simulation end with q > 150 au, a < 10^4 au, and e < 0.9. Figure 3, panels b and c, shows large numbers of counted planets with q between 150 and roughly 550 au, many with a approaching 10^4 au and e near 0.9. These occupy the same dynamical regime as the excluded q = 35-150 au, e ~ 0.9-1 orbits that the paper itself says \"may not maintain long-term stability.\" Galactic tides, passing field stars, and secular forcing from inner giants can alter q on Gyr timescales. Unless the authors supply long-term integrations of the trapped population (or an empirically calibrated stability boundary), the quoted efficiencies and the \"at least 10^-3 per star\" claim should be treated as upper-envelope estimates rather than established lower bounds. This is a load-bearing issue because it directly affects every quantitative result in the abstract.","section":"Main text, trapping-efficiency definition and final paragraph"},{"comment":"The combined 40% probability assumes that both Solar System scattering events—the ice-giant formation phase and the final giant-planet instability—occurred while the Sun was still in its natal embedded cluster (3-10 Myr). The paper itself cites constraints that the final instability happened within the first 10-100 Myr (refs. 33, 34), and only about 1% of clusters remain bound at ~100 Myr. If the final instability occurred after cluster dispersal, only the single-event 5-10% efficiency applies and the abstract's \"rising to 40% if both were\" is not supported. The manuscript should either quantify what fraction of allowed final-instability timings falls within the cluster phase or state explicitly that the 40% figure is conditional on both events occurring before cluster dispersal.","section":"Main text, Solar System probability paragraph"},{"comment":"The \"at least 10^-3 per star\" claim is built from a chain of literature and simulation factors: giant-planet occurrence of 1-10%, an instability fraction of 75-90%, an assumed lower limit of one ejected planet per instability, and a trapping efficiency of 1-10%. The text itself gives a plausible range of 7.5 x 10^-5 to 9 x 10^-3; the central value is near the geometric middle of this range, not obviously a lower bound. Each factor could be smaller, and the trapping-efficiency denominator Nwide+Nejec excludes planets that remain bound but not wide, so the status of \"at least\" should be justified or replaced by a central estimate with a stated uncertainty.","section":"Main text, occurrence-rate estimate"}],"minor_comments":[{"comment":"The caption reads \"he eccentricity distribution of radial velocity exoplanets\" and should read \"The eccentricity distribution.\"","section":"Extended Data Figure 1 caption"},{"comment":"\"Burlisch-Stoer\" should be \"Bulirsch-Stoer\" in the two places where it appears.","section":"Methods, Circumbinary Star Systems"},{"comment":"The verb \"refereed\" is used in place of \"referred\" twice (\"hereafter refereed to as\" and \"usually refereed to as open-clusters\"); please correct these.","section":"Main text and Methods"},{"comment":"Some data points in Figure 4 appear to exceed 10% trapping efficiency, while the text summarizes the efficiency as 1-10%; please clarify in the caption or text whether these points correspond to the exploratory Solar-System-like (J-S-U) runs or to particular cluster configurations, and ensure the legend distinguishes the no-cluster crosses consistently.","section":"Figure 4"},{"comment":"The Data availability statement says source data are \"available at this link\" but no URL appears in the manuscript text; please include the actual repository link.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid simulation study with large statistics and transparent methodology. The main concern is not circularity or fabrication but the gap between the simulated 10-20 Myr timescale and the asserted Gyr-timescale survival, which directly affects the headline efficiencies and occurrence rate. I would support publication after the authors either add long-term stability integrations or soften the claims to explicitly acknowledge the boundary uncertainty. The paper is well within scope for a planetary-dynamics journal; my recommendation of major_revision is driven by the load-bearing nature of the survival assumption."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a serious simulation study, and the core mechanism is not new — Bailey & Fabrycky and the Raymond, Izidoro & Kaib Oort-planet paper already showed that a stellar flyby can lift the perihelion of a planet scattered to a wide orbit. What this paper adds is the systematic demographics: ~30,000 runs across five planetary architectures and a wide range of cluster densities, with clear definitions, binomial error bars, and an explicit statement of which setups are exploratory. The eccentricity-distribution calibration for the gas-giant systems is a good methodological choice. The trapping efficiencies are measured, not fitted to the target result.\n\nThe real soft spot is the one your stress-test identifies: the trapping-efficiency definition counts orbits with q > 150 au, a < 10^4 au, e < 0.9 as successes at simulation end (10-20 Myr), and the claim that most of these survive tens of Gyr is asserted, not demonstrated. Since the sample includes orbits just above the q threshold with e near 0.9, long-term survival is not guaranteed. The fix is a set of Gyr integrations of representative trapped orbits with the Galactic tide and field stars. Until then, I would treat the absolute efficiencies and the occurrence rate as provisional, though the mechanism itself is solid.\n\nThe Solar System timing issue is real but more minor. The 5-10% and 40% figures are conditional on the instabilities happening during the cluster phase; the final instability is constrained to 10-100 Myr, so the marginal probability is lower than 40%. The paper states the conditionality, but the abstract's \"rising to 40%\" will likely be read as an actual probability rather than an upper bound. Also, the abstract's \"at least 10^-3 per star\" is not exactly consistent with the body's quoted range of 7.5e-5 to 9e-3; the lower end is an order of magnitude below 10^-3. The range is honest, but the abstract overstates it.\n\nOne more distinction worth keeping straight: the 5-10% \"wide-orbit planet\" efficiency is a broad category (a < 10^4 au, q > 150 au, e < 0.9), not the much narrower Planet Nine box (a ~ 250-1000 au, e < 0.6). The paper is transparent that P9-specific matches occur in ~1% or fewer of ice-giant simulations, so the headline numbers should not be quoted as P9 probabilities.\n\nI'd send this to a serious referee. The authors know the field, the stats are well reported, and the long-term stability question is answerable. The paper will be a useful baseline for wide-orbit planet formation even if the survival corrections shift the numbers. I'd also bring it to our reading group; it's a good example of how to turn an existing mechanism into a quantitative prediction.","headline":"A large, honestly reported simulation campaign that turns a known flyby-trapping mechanism into quantitative demographics; the headline occurrence rate rests on a long-term stability assertion that isn't yet tested.","tokens_in":28369,"tokens_out":5320,"would_cite":true,"duration_ms":49026,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Wide-orbit planets are a natural byproduct of dynamical instabilities in stellar birth clusters.","keywords":["wide-orbit planets","Planet Nine","dynamical instabilities","stellar birth clusters","planet-planet scattering","stellar flybys","trapping efficiency","solar system formation"],"falsifier":"Take the full set of simulated trapped orbits with perihelia between 150 and a few hundred au and eccentricities up to 0.9, then integrate them for 4.5 Gyr including Neptune's perturbations, galactic tides, and passing field stars; if the survival fraction is substantially below one, the quoted efficiencies and the $10^{-3}$ per star rate must be reduced by that factor. A survey measuring wide-orbit planet occurrence around metal-poor versus metal-rich gas-giant hosts would also test the predicted metallicity enhancement directly.","tokens_in":27136,"feed_emoji":"🪐","tokens_out":7064,"duration_ms":63384,"temperature":0.7,"pith_summary":"Planets on very wide, eccentric orbits are hard to explain with planet formation alone because most protoplanetary disks are far smaller than the orbits observed. This paper argues that such planets are a natural byproduct of dynamical instabilities happening while a planetary system is still embedded in its natal stellar cluster. In this picture, planet-planet scattering first flings a planet onto an eccentric orbit, and a stellar flyby or a change in the cluster potential raises its perihelion, detaching it from further scattering. The simulations give trapping efficiencies of 5-10 percent for Solar System-like instabilities, 1-5 percent for gas-giant exoplanet instabilities, and a resulting population of at least $10^{-3}$ wide-orbit planets per star. A sympathetic reader should care because this connects widely separated exoplanets and the hypothesized Planet Nine to a common early phase of every star's life.","feed_headline":"Birth-cluster flybys trap planets on very wide orbits","feed_subtitle":"Simulations put the trapping odds at 1-10 percent, implying at least one wide-orbit planet per thousand stars.","key_machinery":"The load-bearing mechanism is two-step orbital trapping. First, planet-planet scattering during a dynamical instability drives a planet onto a highly eccentric orbit with an apoastron of several hundred au, too far to be circularized by the disk and vulnerable to ejection. Second, an external perturbation from a star passing within roughly 1000 au, a change in the cluster's gravitational potential near the cluster center, or an encounter with a dense gas filament kicks the planet so that its periastron rises; once the periastron is large enough, with $q>150$ au, $a<10^4$ au, and $e<0.9$, the planet is dynamically decoupled from the inner system and survives after cluster dispersal. The paper quantifies this with the trapping efficiency $\\epsilon = N_{\\rm wide}/(N_{\\rm wide}+N_{\\rm eject})$, which links the number of wide-orbit planets to the number of free-floating planets produced by the same instability.","core_discovery":"The paper's central claim is that very wide-orbit planets with semimajor axes between roughly 100 and 10,000 au are not rare anomalies but a predictable outcome of planetary dynamical instabilities that occur before the natal cluster disperses. A planet scattered by other planets onto an eccentric orbit with a large apoastron can be trapped if a nearby stellar flyby or a rapid change in the cluster potential delivers a kick that raises its periastron, decoupling it from the inner planets. Using about 30,000 N-body simulations of five classes of systems, the authors find that Solar System-like configurations trap 5-10 percent of scattered planets, extrasolar gas-giant instabilities trap 1-5 percent, and ice-giant-only or circumbinary systems trap under 1 percent. Applied to the Solar System, the model gives a 5-10 percent chance that either the ice-giant growth phase or the final giant-planet instability produced a Planet Nine-like object, rising to roughly 40 percent if both happened during the cluster phase. Combining the efficiencies with the known frequency of giant-planet instabilities, the paper concludes that wide, eccentric planets occur at least $10^{-3}$ per star.","pith_inferences":["One testable extension is to search for wide-orbit planets in very young clusters a few million years old, before dispersal; if trapping is the main route, these clusters should already show a measurable population at 100-1000 au, possibly with randomized orbital inclinations.","If the metallicity correlation holds, wide-orbit planet searches around metal-rich, gas-giant-hosting stars could raise the detection rate by a factor of several relative to unbiased surveys, because the underlying instability rate is higher there.","The same trapping logic may apply to lower-mass ejected bodies: planetesimals scattered by giant planets could be parked at the inner edge of the Oort cloud by the cluster, which the paper notes as a byproduct consistent with Oort cloud formation.","A direct long-term stability check, using gigayear integrations of the trapped orbits with galactic tides and passing stars, would sharpen the quoted efficiencies; the paper does not perform those integrations, so the survival assumption remains the main uncertainty."],"forward_implications":["If the mechanism is correct, the Solar System's two early instabilities together give a meaningful chance, up to about 40 percent, that a Planet Nine-like planet was trapped and still exists.","Wide-orbit, eccentric planets should be most common around stars that already host gas giants and are metal-rich, giving a concrete observational target for direct-imaging surveys.","The same instabilities that create free-floating planets should also create a comparable population of bound wide-orbit planets, with the ratio set by the trapping efficiency.","Current observational upper limits of a few percent at 100-5000 au are consistent with the predicted roughly $10^{-3}$ per star occurrence, so deeper surveys can test the prediction.","Trapping efficiency depends on cluster density and lifetime: compact or long-lived clusters trap more planets, while systems that eject planets too quickly, such as circumbinary gas giants, or too slowly, such as pure ice giants, trap almost none."],"supporting_citations":[{"why":"Provides the setup for accretion of Uranus and Neptune from inward-migrating planetary embryos, the basis of the Solar System ice-giant formation simulations.","marker":"[25]"},{"why":"Statistical study of the early Solar System instability with four, five, and six giant planets, used for the final giant-planet instability scenario and ejected ice-giant counts.","marker":"[27]"},{"why":"Planet-planet scattering in planetesimal disks, used to set gas-giant exoplanet instability initial conditions that match observed eccentricity distributions.","marker":"[32]"},{"why":"Supports an early Solar System instability triggered by dispersal of the gaseous disk, constraining the instability to the cluster phase.","marker":"[34]"},{"why":"Establishes the embedded-cluster Plummer potential model for external perturbations on planetary systems, foundational to the cluster simulations.","marker":"[66]"},{"why":"Supplies the algorithm for initializing gas-giant planetary systems in the extrasolar scattering simulations.","marker":"[68]"},{"why":"Provides the modified MERCURY integrator with cluster stellar objects that is used for all simulations.","marker":"[98]"},{"why":"Defines the five-giant-planet Solar System scenario used as the starting point for the early instability simulations.","marker":"[112]"}],"fun_headline_variants":["Cluster flybys lock in scattered planets' wide orbits","Wide orbits from stellar birth chaos, simulations show","Planet Nine odds rise to 40% with birth-cluster events","Birth-cluster scattering makes wide-orbit planets common","How stellar flybys stabilize far-flung planets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The efficiencies assume that planets counted as trapped, with perihelia beyond 150 au and semimajor axes below 10,000 au, actually remain on those orbits for gigayears after the cluster disperses, even though the simulations only run 10-20 million years.","fun_headline_variants_meta":{"raw":{"variants":["Cluster flybys lock in scattered planets' wide orbits","Wide orbits from stellar birth chaos, simulations show","Planet Nine odds rise to 40% with birth-cluster events","Birth-cluster scattering makes wide-orbit planets common","How stellar flybys stabilize far-flung planets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000408,"raw_usage":{"total_tokens":2159,"prompt_tokens":1029,"completion_tokens":1130,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":1051}},"tokens_in":645,"tokens_out":1130,"duration_ms":10165,"temperature":1.0,"reasoning_tokens":1051,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:37:34.215656+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the full set of simulated trapped orbits with perihelia between 150 and a few hundred au and eccentricities up to 0.9, then integrate them for 4.5 Gyr including Neptune's perturbations, galactic tides, and passing field stars; if the survival fraction is substantially below one, the quoted efficiencies and the $10^{-3}$ per star rate must be reduced by that factor. A survey measuring wide-orbit planet occurrence around metal-poor versus metal-rich gas-giant hosts would also test the predicted metallicity enhancement directly.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the embedded-cluster Plummer potential model for external perturbations on planetary systems, foundational to the cluster simulations."},{"cited_title":"& Nesvorný, D","cited_arxiv_id":null,"evidence_quote":"Supplies the algorithm for initializing gas-giant planetary systems in the extrasolar scattering simulations."},{"cited_title":"A., White, E","cited_arxiv_id":null,"evidence_quote":"Provides the modified MERCURY integrator with cluster stellar objects that is used for all simulations."},{"cited_title":"Young Solar System’s Fifth Giant Planet? Astrophys","cited_arxiv_id":null,"evidence_quote":"Defines the five-giant-planet Solar System scenario used as the starting point for the early instability simulations."}],"review_version":1}