{"id":"af5597f1-4caf-4b79-bc0e-8a2e81453742","arxiv_id":"2505.00762","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"N-body simulations show that 50 to 90 percent of wide planet-planet binaries in dense star-forming regions are destroyed within a few million years, so many more than the observed 42 JuMBOs must have formed.","lead":"This paper uses N-body simulations to show that 50 to 90 percent of wide Jupiter-mass binary objects are destroyed within a few million years in dense star-forming regions. If the Orion JuMBOs are real, many more of them formed than the 42 we see, and the original pairs may have been even wider.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inference that 'many more must form' is not checked against the total planetary-mass object inventory that the same JWST surveys should have seen after binaries are ionized.","rationale":"The reader's weakest assumption focuses on the reality of the JuMBOs and the hand-selected initial conditions. I agree those are genuine caveats, but the most load-bearing gap in the central inference is internal: the paper converts a decline in binary fraction into a required primordial binary population without checking whether the ionized components, which remain as single planetary-mass objects in the simulations, are consistent with the observed total planetary-mass object counts. This is a concrete, testable omission rather than a generic parameter-sensitivity worry. It does not overturn the N-body result itself, but it determines whether the headline conclusion about the required primordial population is quantitatively sound. Because the reader already returned a conditional verdict and this concern strengthens rather than replaces that conditionality, the verdict should remain unchanged.","tokens_in":8317,"tokens_out":14177,"duration_ms":169824,"concrete_test":"Use the saved particle data, or rerun with particle tracking, to count at 1 Myr and at 4 Myr the number of single planetary-mass objects retained within the observed ONC field and in the same magnitude range as the JWST survey. Scale the simulations so that the surviving binary count equals the 42 observed JuMBOs, then compare the predicted single count with the number of free-floating planetary-mass objects reported by Pearson & McCaughrean (2023), Langeveld et al. (2024), and Luhman et al. (2024) after applying incompleteness corrections. If the predicted singles exceed the observed total, the primordial-population inference needs revision; if they are compatible, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 reads the drop in fbin (Fig. 1) as destruction, and Section 4(ii) concludes that a 'significantly higher number of primordial systems must have been present' to explain the 42 JuMBOs. But the components of ionized binaries are not lost: they become single planetary-mass objects that are included in the denominator of fbin and should be present in the same JWST field. Taking the dense-case survival fraction of about 0.1 at face value, the observed 42 binaries imply about 420 primordial binaries, i.e. roughly 840 planetary-mass members, of which about 756 would be singles after 1 Myr. The paper does not compare this predicted population of liberated singles with the observed number of free-floating planetary-mass objects in the ONC, nor does it quote a retention fraction within the observed field. This omission matters because the JWST surveys that found the JuMBOs also counted the single planetary-mass objects; if the predicted retained singles already exceed the observed total, or force the total planetary-mass inventory well above what the same surveys allow, then the 'many more must form' conclusion is not self-consistently supported. The N-body survival fractions may themselves be right; the missing object-budget check is the weakest link in the leap from fbin to the required primordial population.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents N-body simulations of subvirial, fractal star-forming regions in which 10 per cent of systems are initially planet--planet binaries whose component masses and separations are drawn from the 42 JuMBO candidates in the Orion Nebula Cluster. From the evolution of the substellar binary fraction, the authors find that 50--90 per cent of these binaries are destroyed within 1 Myr for initial median densities of 10^2 and 10^4 M_sun pc^-3 (simulations A--C), and that the surviving separation distribution is shifted to smaller separations in the dense case. They conclude that the observed JuMBO population requires a substantially larger primordial population of planet--planet binaries, and that if the ONC was initially dense, the initial separation distribution must have been wider than observed, possibly resembling a flat 50--500 au distribution. The paper is candid about its simplifications, including zero eccentricities, no stellar binaries, and no gas potential, and it explicitly acknowledges the alternative background-contamination interpretation of Luhman (2024).","tokens_in":8540,"tokens_out":2962,"duration_ms":33321,"significance":"If the dynamical destruction fractions are correct, the paper provides the first quantitative constraint on the survival of wide planet--planet binaries in star-forming regions, and it sharpens the challenge that the JuMBO candidate population poses to formation theories. The design has genuine strengths: it uses the actual observed masses and separations where available, it brackets the plausible density range of the ONC, it includes a brown-dwarf-bearing variant (simulation C), and it checks sensitivity to the initial separation distribution (simulations D and E). The central simulation result—that wide, low-mass binaries are efficiently ionized in dense environments—is physically plausible and internally consistent. The main weakness is not the N-body calculation itself but the leap from the destruction fractions to the inference that 'many more must form', which is not tested against the observable inventory of liberated single planetary-mass objects in the same JWST fields.","major_comments":[{"comment":"The inference that 'a significantly higher number of primordial systems must have been present' is not checked against the observable single-planetary-mass-object budget. Ionized binaries do not disappear; their components become single planetary-mass objects that enter the denominator of fbin in Eq. (2) and should be detectable in the same JWST fields. Taking the dense-case survival fraction of ~0.1 at face value, the 42 observed binaries imply ~420 primordial binaries and ~756 liberated singles after 1 Myr. The paper does not compare this predicted retained single population, nor the implied total planetary-mass inventory, with the observed number of free-floating planetary-mass objects in the ONC, nor does it quote a retention fraction within the observed field. Without such a consistency check, the conclusion that many more JuMBOs must have formed is not self-consistently supported, even though the simulated survival fractions may be correct.","section":"Section 3 and Section 4(ii), Fig. 1"},{"comment":"Setting all eccentricities to zero is a strong assumption that likely biases the destruction fractions. For a given semi-major axis, an eccentric binary spends part of its orbit at larger physical separations, increasing its collisional cross-section and susceptibility to ionization; conversely, circular orbits maximize binding energy per unit angular momentum for a given semi-major axis. The observed JuMBO separations are projected separations, not semi-major axes, so the mapping from observation to initial condition is also uncertain. Since the headline 50--90 per cent destruction fractions are the paper's central quantitative output, the sensitivity to eccentricity should be tested (e.g., with a thermal or uniform eccentricity distribution) or the authors should justify why the zero-eccentricity choice is conservative.","section":"Section 2, paragraph on eccentricities"},{"comment":"The claim that dynamical processing in the dense case with a flat 50--500 au initial distribution 'would reproduce the observed JuMBO distribution' (Fig. 4a) is made by visual inspection of cumulative distributions. The paper provides no quantitative comparison, such as a Kolmogorov--Smirnov or Anderson--Darling test, between the evolved simulated distribution and the observed JuMBO sample, nor error bars on the simulated distributions across the ten realizations. Because conclusion (iii)—that the observed distribution has been dynamically sculpted and the initial distribution must contain wider systems—rests on this comparison, a statistical test is needed to establish whether the agreement is meaningful and whether the low-density case (Fig. 4b) is truly excluded.","section":"Section 3, Figs. 3 and 4"}],"minor_comments":[{"comment":"Several typos should be corrected: 'poper motion' in Section 4, 'therough' in Section 3, 'he observed' in conclusion (ii), 'showns' in the Fig. 2 caption, and 'inital' in conclusion (iii).","section":"Throughout"},{"comment":"The text states that 10 per cent of systems are randomly selected to be planet--planet binaries, but it is not stated whether this 10 per cent applies independently in simulation C, where the IMF also produces brown dwarfs down to 0.01 M_sun, or whether the reported initial fbin of 0.29 follows from the IMF sampling. A sentence clarifying the construction of the initial substellar population in C would help.","section":"Section 2, Table 1"},{"comment":"The observation that the JuMBO separation histogram 'shows an increasing trend to higher separations' is used to argue for incompleteness and an even larger primordial population. This is an interesting point, but it should be quantified or referenced to the completeness limits of the Pearson & McCaughrean (2023) catalogue, since observational incompleteness is otherwise not discussed in detail.","section":"Section 3, paragraph after Fig. 2"},{"comment":"The absence of a background gas potential is acknowledged, but its effect on the early evolution of a subvirial fractal is not discussed. Including the gas potential would deepen the potential and could increase or decrease the effective destruction rate depending on the collapse timescale; a brief justification that the two density choices bracket the relevant behavior would strengthen the presentation.","section":"Section 2, paragraph on gas potential"}],"recommendation":"major_revision","confidential_remarks":"The paper's core N-body result is plausible and the authors are appropriately cautious about the Luhman (2024) contamination alternative. My main concern is the logical gap between the simulated destruction fractions and the 'many more must form' conclusion, which is exactly the kind of inference that should be checked against the same observational data that produced the 42 JuMBOs. The missing object-budget check is fixable and should be added in revision. The eccentricity sensitivity and the statistical comparison of separation distributions are also important but secondary."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is the first explicit N-body study of whether wide planet-planet binaries like the JuMBOs can survive in star-forming regions. That alone makes it worth reading. The authors take established cluster-disruption machinery (box fractals, kira, standard IMF choices) and point it at a new mass and separation regime. The headline result is believable: in a dense 10^4 Msun/pc^3 region the substellar binary fraction collapses from 1.0 to 0.1 within a Myr, and even at 10^2 Msun/pc^3 half are destroyed. The internal simulation logic is sound, the paper is candid about its simplifications, and the qualitative claim that wider systems are preferentially destroyed holds up in the cumulative distributions.\n\nThe soft spots are real but not fatal. There are no statistical uncertainties on the 50-90 per cent destruction range; ten seeds per setup is fine for a trend, but not for a headline number without a spread. The initial conditions are hand-picked: zero eccentricity, no gas potential, no stellar binaries, fixed fractal dimension and virial ratio. The authors acknowledge these, but they directly shape the quantitative fractions. Reproducibility is also weak: data are only 'available on reasonable request' and no code is public.\n\nThe bigger gap is the one the stress-test note flags. The paper reads the drop in fbin as destruction and concludes that 'many more must form' to explain the 42 JuMBOs. But the components of ionized binaries are not lost; they become single planetary-mass objects in the same cluster, and they are included in the denominator of fbin. In the dense case, a survival fraction of 0.1 implies roughly 420 primordial binaries, so about 756 liberated singles after 1 Myr. The paper never checks whether that predicted population of free-floating planetary-mass objects is consistent with the same JWST surveys that found the JuMBOs. If those singles would already exceed the observed inventory, the 'many more must form' conclusion is not self-consistently supported. The N-body fractions themselves can be right; the missing object-budget check is the weakest link in the astrophysical leap.\n\nThe paper is also upfront that all of this is conditional on the JuMBOs being real bound binaries, given Luhman's background-contamination argument. That is the right caveat, but it does mean the inference is a conditional one.\n\nWho is this for? Anyone working on substellar binaries, star-formation dynamics, or JWST observations of free-floating planets. The dynamical constraint is useful even if the primordial-population conclusion needs more work. A serious referee should engage with this; the requested revisions are mostly about adding error bars, sharing data, and running the single-object budget check. I would bring it to reading group and would cite it in my own work on binary disruption.\n\nMy recommendation: send it to peer review, but with a referee who will insist on the missing inventory check before publication.","headline":"First explicit N-body survival fractions for JuMBO-scale binaries; the core dynamical result is credible, but the leap to 'many more must form' skips a check against the observed single planet-mass inventory.","tokens_in":9097,"tokens_out":1628,"would_cite":true,"duration_ms":18434,"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":"N-body simulations show that 50–90 per cent of planet–planet binaries with JuMBO-like properties are destroyed within a few Myr in dense star-forming regions, so the 42 observed JuMBOs require a much larger primordial population.","keywords":["JuMBOs","planet-planet binaries","star-forming regions","N-body simulations","dynamical destruction","Orion Nebula Cluster","substellar binary fraction","JWST free-floating planetary-mass objects"],"falsifier":"Take multi-epoch astrometry and spectroscopy of the 42 JuMBO candidates and compare their proper motions and colours with Orion Nebula Cluster members; if most are consistent with reddened background stars rather than young cluster members, the premise that they are primordial planet–planet binaries fails and the calculated destruction fractions no longer apply to them. A complementary test would be a deep JWST survey of another dense young cluster that finds wide planet–planet binaries surviving at high frequency, which would contradict the predicted destruction.","tokens_in":8091,"feed_emoji":"🪐","tokens_out":10201,"duration_ms":94266,"temperature":0.7,"pith_summary":"The paper asks whether the wide Jupiter-Mass Binary Objects (JuMBOs) recently found by JWST in the Orion Nebula Cluster can survive the gravitational encounters of a star-forming region. The authors place a population of planet–planet binaries with JuMBO-like masses, separations, and circular orbits into N-body simulations (computer models of many gravitationally interacting objects) of substructured, subvirial clusters and follow them for 10 Myr. They find that 50–90 per cent of these binaries are destroyed within a few Myr, depending on the cluster density. If the JuMBOs are real bound systems, the observed 42 require that far more of them formed, and if the Orion Nebula Cluster was initially very dense, the original separation distribution must have been wider than the observed one. The paper thereby turns an apparent formation puzzle into a quantitative dynamical survival problem.","feed_headline":"50–90% of wide planet-planet binaries are destroyed within a few Myr","feed_subtitle":"If the Orion Nebula was dense at birth, the 42 observed JuMBOs require a far larger primordial population.","key_machinery":"The argument is carried by the substellar binary fraction, $f_{\\rm bin}=B/(S+B)$, where $B$ is the number of binary systems and $S$ the number of singles among objects below the hydrogen-burning mass limit, tracked over time. The simulations place $N_{\\rm sys}=1500$ systems in a box-fractal (substructured, filament-like) configuration with fractal dimension $D=1.6$ and subvirial velocities (virial ratio $\\alpha_{\\rm vir}=0.3$), assigning 10 per cent of systems to be planet–planet binaries whose component masses and separations are drawn from the observed JuMBO catalogue, or in two runs from a flat 50–500 au distribution. Each setup is evolved for 10 Myr with an N-body integrator, and the binary fraction and cumulative separation distribution are recorded. Because stellar evolution, gas, and stellar binaries are not included, any change in these quantities is attributable to dynamical encounters in the cluster.","core_discovery":"The central claim is that dynamical encounters in star-forming regions destroy most planet–planet binaries with JuMBO-like properties within a few Myr. In the dense simulations (initial median density $10^4\\,M_\\odot\\,\\mathrm{pc}^{-3}$) the substellar binary fraction drops from 1.0 to 0.1 within 1 Myr; in the lower-density runs ($10^2\\,M_\\odot\\,\\mathrm{pc}^{-3}$) it drops to 0.5. The authors therefore conclude that between 50 and 90 per cent of such binaries are destroyed, and that the observed population of 42 JuMBOs must have been drawn from a much larger primordial population. They further show that, in the dense case, the surviving systems are preferentially the closer ones, so the observed separation distribution has been dynamically sculpted; the initial JuMBO separations would have extended to roughly 500 au.","pith_inferences":["Extending beyond the paper, rerunning the same setup with nonzero eccentricities, stellar binaries, and a gas potential would probably raise the destruction rate, so the quoted 50–90 per cent is more likely a lower bound than an upper bound on dynamical loss.","A decisive test is astrometric: if multi-epoch imaging shows that most JuMBO candidates are reddened background stars rather than members of the Orion Nebula Cluster, then the survival calculation applies to a population that is not actually there.","The same method could predict how the wide planet–planet binary fraction should vary across clusters of different density and age, giving JWST surveys outside Orion a quantitative target to confirm or refute the destruction rates.","If the large inferred primordial population is real, formation mechanisms that produce many wide, weakly bound pairs at once become more attractive than mechanisms that produce them rarely."],"forward_implications":["The 42 observed JuMBOs must represent only part of the primordial population: in the dense case the implied initial number is about ten times larger, and in the lower-density case at least double.","If the Orion Nebula Cluster formed at a density near $10^4$ solar masses per cubic parsec, the present-day JuMBO separation distribution is not primordial; wide systems out to about 500 au must have existed and been preferentially destroyed.","If the cluster's initial density was instead similar to its present-day value near $10^2$ solar masses per cubic parsec, the shape of the separation distribution is roughly preserved even though half the binaries still die.","Because destruction preferentially removes the widest, most weakly bound systems, the surviving separation distribution shifts toward shorter separations, and any observational incompleteness at wide separations strengthens the need for a large primordial population."],"supporting_citations":[{"why":"Supplies the observed JuMBO catalogue: the 42 systems, component masses, and 28–384 au separations used to initialise the planet–planet binaries and to compare with the simulated survivors.","marker":"Pearson & McCaughrean (2023)"},{"why":"Presents the alternative interpretation that the JuMBO candidates are reddened background sources, the premise against which the survival calculation is conditioned.","marker":"Luhman (2024)"},{"why":"Provides the box-fractal method used to set up the substructured, subvirial initial spatial and kinematic conditions of the star-forming regions.","marker":"Goodwin & Whitworth (2004)"},{"why":"Supplies the N-body integrator used to evolve the 1500-system simulations for 10 Myr.","marker":"Portegies Zwart et al. (1999, 2001)"},{"why":"Provides the present-day density, age, and stellar binary fraction of the Orion Nebula Cluster that set the low-density initial conditions and the context for the dense case.","marker":"King et al. (2012)"},{"why":"Supplies the initial mass function used to draw the masses of the single-star systems in the simulations.","marker":"Maschberger (2013)"},{"why":"Establishes the earlier result that stellar and substellar binaries are disrupted in dense star-forming regions, providing the dynamical baseline this paper extends to planet–planet binaries.","marker":"Parker & Goodwin (2011)"}],"fun_headline_variants":["Most wide planet binaries are torn apart within a few Myr","Dense clusters destroy 50–90% of JuMBO-like binaries","JuMBO survival demands a far larger initial population","Observed JuMBOs are just the tip of a destroyed iceberg","Wide planet pairs need extreme birth rates to persist"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the 42 JuMBO candidates are real, initially circular, bound planet–planet binaries with the exact masses and separations from the observed catalogue, placed in a cluster with a 10 per cent binary fraction, no stellar binaries, and no gas; if the candidates are background contaminants or the initial conditions differ, the destruction fractions change.","fun_headline_variants_meta":{"raw":{"variants":["Most wide planet binaries are torn apart within a few Myr","Dense clusters destroy 50–90% of JuMBO-like binaries","JuMBO survival demands a far larger initial population","Observed JuMBOs are just the tip of a destroyed iceberg","Wide planet pairs need extreme birth rates to persist"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1297,"prompt_tokens":940,"completion_tokens":357,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":270}},"tokens_in":556,"tokens_out":357,"duration_ms":3932,"temperature":1.0,"reasoning_tokens":270,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:36:17.949770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take multi-epoch astrometry and spectroscopy of the 42 JuMBO candidates and compare their proper motions and colours with Orion Nebula Cluster members; if most are consistent with reddened background stars rather than young cluster members, the premise that they are primordial planet–planet binaries fails and the calculated destruction fractions no longer apply to them. A complementary test would be a deep JWST survey of another dense young cluster that finds wide planet–planet binaries surviving at high frequency, which would contradict the predicted destruction.","supporting_citations":[{"cited_title":"R., Parker R","cited_arxiv_id":null,"evidence_quote":"Provides the present-day density, age, and stellar binary fraction of the Orion Nebula Cluster that set the low-density initial conditions and the context for the dense case."}],"review_version":1}