{"id":"a3e7b746-462e-4fe7-a3bc-b1c795971ab1","arxiv_id":"2505.04398","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Double hot Jupiters can form via simultaneous von Zeipel-Lidov-Kozai migration of two cold Jupiters in a stellar binary, with up to ~9% of hot-Jupiter-hosting binaries predicted to host a second.","lead":"This paper uses N-body simulations to show that two cold giant planets, one around each star in a binary system, can migrate inward together and become two hot Jupiters through the same gravitational mechanism. It predicts that up to about 9% of close and moderate separation binary systems that already host one hot Jupiter could host a second, a targetable prediction for future searches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ~9% rate in Section 4 fixes the companion planet's inclination at the optimal I2=93° and marginalizes only over I1; an independent isotropic I2 would lower the expected rate to ~1.4%, so the abstract's assumptions do not fully specify the headline number.","rationale":"The reader's verdict of CONDITIONAL is well founded. The paper's N-body machinery is standard (REBOUND/IAS15 with GR and equilibrium tides), and the proof-of-concept that two cold Jupiters can undergo simultaneous ZLK migration is convincing. The central prediction is presented as an upper limit, and the paper explicitly labels assumptions (1) and (2). My stress-test focuses on a narrower but numerically important point that the reader flagged only in passing: Section 4 fixes I2=93° and marginalizes only over I1. The 9% is therefore the probability of a second hot Jupiter conditional on the companion planet having a near-optimal inclination, not the probability averaged over the orientation distribution that the paper otherwise treats as isotropic. If I2 is drawn independently from an isotropic distribution, the rate drops to ~1.4%; if I2 is perfectly aligned with I1, the 9% is recovered. The paper does not state in its abstract which orientation model the headline corresponds to, making the number easy to over-read. This is load-bearing because it directly affects the numerical claim and is testable with the existing simulation setup. I concur with the reader that the conditional rate should not be interpreted as a firm occurrence prediction, and I recommend the paper remain CONDITIONAL with a request to clarify or marginalize the companion inclination distribution.","tokens_in":20875,"tokens_out":18414,"duration_ms":179405,"concrete_test":"Re-run the Section 4 population synthesis with the same binary sample and I1 grid, but sample I2 from (a) an isotropic distribution independent of I1 and (b) a correlated distribution such as I2 = I1 plus Gaussian scatter with σ = 5°. Compare the marginalized double-hot-Jupiter rate with the 9% headline. If the independent-isotropic rate falls below ~2%, the headline is a strong upper limit rather than an expected rate under the stated assumptions, and the paper should present the marginalized rate as the primary prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline ~9% rate in Section 4 is computed by fixing the companion planet's initial mutual inclination at I2=93° (Figure 5, left) and varying only I1 over a grid. The paper then multiplies the width of the successful band in cos I1 (0.306, or 15.3% of the isotropic range of length 2) by the 59.9% in-band success probability. That product is P(HJ2 | HJ1, I2=93°, I1 isotropic), not the marginal rate over the companion's orientation. If I2 is instead drawn independently from the same isotropic distribution used for I1, the joint probability that both planets fall in the favorable inclination window is 0.153^2 ≈ 2.3%, and with ~60% in-band success the expected rate falls to roughly 1.4%—about a factor of 6 below 9%. Even under the more physical assumption of perfectly aligned planetary orbits (I2=I1, as might occur for co-aligned protoplanetary disks), the rate is 9%, but the paper does not state that its synthesis assumes this alignment. The abstract's assumption (2) concerns cold Jupiter presence, not orbital geometry, so a reader interpreting 9% as the expectation under the stated assumptions would overestimate the occurrence rate substantially unless I2 is also assumed optimal. This is a load-bearing modeling choice that directly affects the central numerical claim and is readily testable.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the formation of double hot Jupiter (HJ) systems in stellar binaries through simultaneous von Zeipel-Lidov-Kozai (ZLK) migration of two cold Jupiters, one around each star. Using four-body N-body simulations with REBOUND/IAS15 including general relativity and equilibrium tides, the authors first demonstrate a proof of concept in a perfectly mirrored fiducial system, then examine robustness to asymmetries in stellar mass, planetary mass, and initial inclination. They also compare four-body and three-body secular behavior in Appendix A. The paper's main quantitative claim is a population synthesis based on Gaia-resolved HJ-hosting binaries with projected separation s ≤ 2000 AU, from which they derive an upper limit of ~9% for the occurrence of double HJ systems among such binaries, under the stated assumptions that most HJs in binaries form via ZLK migration and that a cold Jupiter forms around both components. The paper recommends that searches for ZLK-migrated double HJs prioritize twin binaries with pericenter approaches of a few hundred AU.","tokens_in":21200,"tokens_out":11615,"duration_ms":117327,"significance":"The core proof of concept is a genuine and useful contribution: the four-body simulations directly demonstrate that mirrored cold Jupiters can both migrate to form double HJs, and they appropriately model the planet-planet interaction that is absent from the standard three-body treatments. The robustness study and the concrete target selection (twin binaries, q ≈ 200 AU) are valuable for future observing programs, and the prediction is falsifiable in principle. The main weakness is that the headline ~9% rate is not the marginal expectation under the assumptions stated in the abstract: it is conditional on a fixed near-orthogonal inclination for the second planet and on a favorable binary pericenter, and it also assumes a cold Jupiter around both stars. These conditionals materially change the rate that a comprehensive, unselected survey would see, and the manuscript should either marginalize over the second planet's orientation or present the number transparently as a conditional upper bound rather than a projected occurrence rate.","major_comments":[{"comment":"The headline ~9% rate is computed with the second planet's initial inclination fixed at I2 = 93° (Section 4: 'We kept the initial inclination I2 fixed at 93°'), while the introduction states that 'Assuming isotropic initial orbital orientations' the authors find a ~9% rate. For an independent isotropic distribution of I2, the joint probability that both planets fall in the successful band in cos I1 is 0.153^2 ≈ 2.3%, and with the reported 59.9% in-band success fraction the expected rate drops to ~1.4%. Even under perfectly aligned planetary orbits (I2 = I1) the rate is 9%, but that is not one of the two assumptions listed in the abstract. Because the abstract and introduction present 9% as the outcome of a comprehensive search under the stated assumptions, a reader would overestimate the expected occurrence rate by roughly a factor of six. The authors should either run simulations that marginalize over I2 with an isotropic prior, or explicitly and prominently state that the 9% is conditional on the second cold Jupiter's orbit being near-orthogonal (I2 ≈ 90°) and revise the abstract accordingly.","section":"Section 4 and Abstract"},{"comment":"The sentence 'For an isotropic set of inclinations for the second planet, the overall success rate of producing double hot Jupiters in the lifetime of the universe is ~46%' is ambiguous and potentially misleading in light of the simulation description that immediately precedes it. The text and Figure 3 state that the simulations are 'initiated uniformly in cos I1', but do not state how I2 was sampled. If I2 was set to the mirror value 180° − I1, then the second planet's orientations are not drawn independently from an isotropic distribution, and the ~46% is a rate for perfectly mirrored systems, not a marginal rate over independent I2. If instead I2 was varied independently, that should be stated explicitly along with the sampling procedure. This ambiguity also affects how the rate connects to Equation (4), where the probability for the second planet is meant to be a function of its own geometry only.","section":"Section 3.1"},{"comment":"The population synthesis assigns every binary a fixed pericenter a*(1−e*) = 200 AU and adopts a* = s (the projected separation), so the resulting ~9% is a simultaneous best-case evaluation over several favorable choices: the second planet's inclination (I2 = 93°), the binary pericenter, and the assumption of a cold Jupiter around both stars. The text does call the result an 'upper limit,' but the abstract and conclusions do not carry the same hedging, and the word 'up to' does not communicate how many separate favorably chosen conditions are being stacked. For the paper's central forecast to be usable, the abstract and the Section 4 discussion should explicitly state that the 9% is conditional on a near-orthogonal second planet, a few-hundred-AU binary pericenter, and the presence of a primordial cold Jupiter around both stars, none of which are guaranteed in an unselected sample of known HJ-hosting binaries.","section":"Section 4, population synthesis setup"}],"minor_comments":[{"comment":"The phrase 'to producedoublehot Jupiter systems' is missing spaces between 'produce', 'double', and 'hot'; this appears to be a LaTeX rendering artifact and should be fixed in the source.","section":"Abstract"},{"comment":"The caption says 'The shaded region spans a ±1 ZLK cycle error bar as defined in Equation 2,' but Equation (2) is a proportionality relation for a timescale and does not define an error bar; please state how the ±1-cycle uncertainty is computed from the simulation scatter.","section":"Figure 4 caption"},{"comment":"The sentence 'The rest of the system parameters, including masses, were set to those shown in Table 1' is confusing because Table 1 includes the fiducial inclinations (I1 = 83°, I2 = 97°) that are explicitly overridden in this section; please state which parameters are retained and which are replaced by the grid or the fixed I2 = 93°.","section":"Section 4"},{"comment":"The phrase 'isochrones are so closely separated' is awkward; consider 'isochrones are so closely spaced in age' or similar.","section":"Section 5.4"},{"comment":"The manuscript would benefit from a data and code availability statement, since the simulations are central to the results and the public packages are named but the initial conditions, analysis scripts, and reproduction details are not deposited.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and the core mechanism result is credible and likely of interest to the exoplanet dynamics community. The main obstacle is the presentation of the 9% rate as a marginal expectation under the abstract's assumptions; this is a load-bearing issue for the central quantitative claim, but it is fixable by rerunning or correctly framing the I2 marginalization. The authors should be asked to address the ambiguity in Section 3.1 as well. I see no citation or novelty concerns: the self-citations support background and methods, and the four-body treatment and observational targeting are original."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper does two things: it shows, with N-body integrations that include tides and GR, that two cold Jupiters on mirrored, mutually inclined orbits around the two stars in a binary can both be driven into hot Jupiter orbits by ZLK migration; and it turns that proof of concept into a predicted occurrence rate for double hot Jupiters among Gaia-resolved hot-Jupiter-hosting binaries. The first half is solid. The fiducial simulation is clean, the asymmetry tests are sensible, and the appendix correctly shows that the standard three-body approximation breaks down in exactly the high-inclination, close-pericenter region where the mechanism operates, which justifies the four-body approach. The qualitative recommendation to search twin binaries with pericenters of a few hundred AU is reasonable.\n\nThe soft spot is the headline number. Section 4 fixes I2 at 93 degrees — the optimal inclination for the second planet — and marginalizes only over I1. The width of the successful band in cos I1 is 0.306, and with a 59.9% in-band success rate the paper reports ~9%. But under the isotropic orientation assumption stated in the abstract, I2 should also be drawn from the same distribution. The joint probability that both planets land in the favorable window is then (0.153)^2, and with the same success rate the expected rate falls to about 1.4%. The 9% is therefore an upper limit, not the expectation under the stated assumptions. The abstract never tells the reader that the companion planet's orbit is assumed to be optimally aligned. That is a load-bearing omission.\n\nThere are smaller issues. The population synthesis adopts a*=s and fixes q=200 AU for every binary, which biases the eccentricity distribution. The rate also assumes a cold Jupiter around both stars; the paper acknowledges cold Jupiter occurrence is only ~10-15% and that correlated formation is unconstrained, so the absolute rate is an upper limit on that count as well. The release has no code or data, which makes the population calculation harder to audit. And the removal of ten 'outlier' points from Figure 4 is described but not motivated in detail.\n\nNone of this kills the central claim. The mirrored double-ZLK pathway is new and worth taking seriously, and the paper is honest about many of its own limitations. I'd send it to a referee. The main request should be to either marginalize over I2 or explicitly relabel the 9% as an upper limit conditional on optimal geometry. As it stands, the mechanism and the target-selection advice are the valuable parts; the absolute rate should not be quoted without the caveat.","headline":"A genuinely new channel for double hot Jupiters via mirrored ZLK migration, backed by real N-body work, but the headline ~9% rate is a conditional upper limit that drops to ~1.4% when the second planet's inclination is treated as isotropic.","tokens_in":21720,"tokens_out":2787,"would_cite":true,"duration_ms":27521,"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":"Two hot Jupiters can form around both stars of a binary, simulations find","keywords":["hot Jupiters","von Zeipel-Lidov-Kozai migration","double hot Jupiters","stellar binaries","tidal friction","high-eccentricity migration","exoplanet dynamics","N-body simulations"],"falsifier":"A survey of the secondary stars of known hot-Jupiter-hosting binaries with separations 100-2000 AU that finds zero hot Jupiters around the secondaries in ~50 systems would rule out the ~9% prediction at high confidence; equivalently, measuring the conditional probability P(cold Jupiter around the companion | cold Jupiter around the primary) and finding it close to the field rate would show the predicted rate should be reduced by roughly an order of magnitude.","tokens_in":20702,"feed_emoji":"🪐","tokens_out":9005,"duration_ms":80848,"temperature":0.7,"pith_summary":"This paper argues that a binary star system can produce a hot Jupiter around each star at once, through mirrored von Zeipel-Lidov-Kozai (ZLK) migration: two cold Jupiters, one around each star, are driven by the companion star into high-eccentricity orbits that tidal friction shrinks and circularizes into close-in hot Jupiters. Using N-body simulations with tides and general relativity, the authors show this happens naturally in symmetric binaries and remains robust to modest asymmetries in stellar mass, planet mass, and mutual inclination. They predict that, if most hot Jupiters in binaries form this way and if a cold Jupiter around one star implies one around the other, then up to about 9% of known hot-Jupiter-hosting binaries with stellar separations out to 2000 AU should also host a second hot Jupiter. A curious reader should care because this is a concrete, testable prediction: the companion stars of known hot Jupiter hosts are largely unmonitored, and finding or failing to find these second hot Jupiters would test whether this formation channel operates.","feed_headline":"Two hot Jupiters can form in one binary, simulations show","feed_subtitle":"A second hot Jupiter may hide around the companion star in up to 9% of known systems.","key_machinery":"The mechanism is the von Zeipel-Lidov-Kozai (ZLK) effect with tidal friction, applied in a mirrored four-body configuration: an inclined stellar companion drives coupled oscillations in each planet's eccentricity and inclination, and tides at pericenter drain orbital energy, shrink the semimajor axis, and circularize the orbit into a hot Jupiter. The mirrored setup—two identical planets on orbits whose mutual inclinations with respect to the binary plane sum to 180°—makes the two evolutions identical, so success for one planet implies success for the other. The paper's population-level estimate rests on a Bayesian reduction: for twin binaries, P(both planets become hot Jupiters | one does) equals P(the second planet becomes a hot Jupiter), turning the four-body problem into a single-planet geometry problem whose 15.3%-of-orientations window, multiplied by the 59.9% simulation success rate within that window, yields the ~9% upper limit.","core_discovery":"The paper's central claim is that double hot Jupiter systems can arise through simultaneous ZLK migration in stellar binaries, and that this channel is efficient enough to be observable. In a perfectly mirrored four-body configuration—two equal-mass stars, two identical cold Jupiters, with planetary orbital planes inclined by 83° and 97° to the binary plane—both planets undergo the same secular eccentricity-inclination cycles and both circularize into hot Jupiters within a few hundred million years. Relaxing the symmetry, the authors find that inclination asymmetries dominate the outcome while stellar and planetary mass asymmetries mainly shift formation timescales; in equal-mass 'twin' binaries the formation time is minimized. Combining the simulation success rate with a Bayesian argument that the joint formation probability factors into independent per-planet probabilities, they predict that up to ~9% of known hot-Jupiter-hosting binaries with projected separations s ≤ 2000 AU could host a second hot Jupiter, and that the most favorable targets are twin binaries whose orbits bring the stars to pericenter distances of a few hundred AU.","pith_inferences":["The ~9% figure is best read as an upper limit rather than an expectation: it assumes a cold Jupiter exists around both stars, whereas the field rate of cold Jupiters around FGK stars is only ~10-15%, so the actual yield of a survey could be several times lower if planet formation is not strongly correlated across twin components.","The paper's assumption of isotropic orbital orientations may overestimate the rate, because several recent studies find an excess of low mutual inclinations in planet-hosting binaries; folding in that measured inclination distribution would be a direct test of the channel's contribution to the observed population.","The same mirrored-migration logic should apply to other close-in planet populations formed by high-eccentricity migration, such as hot Neptunes or super-Earths, so the 'double close-in planet' prediction is not restricted to Jupiter-mass planets.","Because the mechanism predicts that second hot Jupiters preferentially orbit the lower-mass star in unequal binaries, measuring the mass ratio of the stellar binary in a survey design directly concentrates the search; this is a testable prioritization strategy the paper gestures toward but does not quantify."],"forward_implications":["A deliberate search of the companion stars of known hot Jupiter hosts in binaries with separations up to 2000 AU should uncover second hot Jupiters at a rate bounded by about 9%.","Blind surveys for ZLK-migrated double hot Jupiters should prioritize twin, equal-mass binaries with stellar pericenter approaches of a few hundred AU, where formation is fastest and the success rate is highest.","In unequal-mass binaries, if a hot Jupiter formed around the more massive star, the second hot Jupiter is likely to have already formed around the less massive companion, making secondary stars the best targets.","For a close, eccentric binary with a*=200 AU and e*=0.7, roughly 46% of isotropically drawn orientations produce double hot Jupiters within 13.8 Gyr, and the rate stays near 33% within 2.3 Gyr.","If the mechanism operates, the occurrence of double hot Jupiters should correlate with binary properties—closer pericenters and equal masses—so the predicted systems should be clustered in that region of binary parameter space."],"supporting_citations":[{"why":"Supplies the fiducial physical parameters and the equilibrium-tide treatment of ZLK migration that the simulations build on.","marker":"Wu & Murray 2003"},{"why":"Establishes ZLK-mediated hot Jupiter formation and the inclination and tidal-disruption framework used to define success.","marker":"Naoz et al. 2012"},{"why":"Provides the REBOUND N-body integrator used for all simulations in the paper.","marker":"Rein & Liu 2012"},{"why":"Provides the IAS15 adaptive-timestep integrator that resolves close pericenter passages during high-eccentricity phases.","marker":"Rein & Spiegel 2015"},{"why":"Defines the q<0.1 AU pericenter criterion used to identify a hot Jupiter in the simulations.","marker":"Rice et al. 2022"},{"why":"Supplies the Gaia cross-match method and catalogue used to identify the companion-star sample of known hot Jupiter hosts.","marker":"El-Badry et al. 2021"},{"why":"Supplies the equilibrium tide prescription for tidal dissipation that drives orbital circularization.","marker":"Eggleton et al. 1998"},{"why":"Cited as evidence that the correlation of planet formation across twin binary components is unconstrained, which the paper acknowledges as a caveat.","marker":"Hand et al. 2025"}],"fun_headline_variants":["Mirrored ZLK migration yields double hot Jupiters in binaries","Double hot Jupiters: a likely outcome in twin-star binaries","Simulations predict double hot Jupiters in 9% of binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline ~9% rate rests on assuming that a cold Jupiter forms around both stars in a binary whenever it forms around one, and that the companion planet's orbit has a favorable inclination; observed cold-Jupiter occurrence around FGK stars is only ~10-15%, and the correlation of planet formation across twin binary components is not yet measured.","fun_headline_variants_meta":{"raw":{"variants":["Mirrored ZLK migration yields double hot Jupiters in binaries","Double hot Jupiters: a likely outcome in twin-star binaries","Simulations predict double hot Jupiters in 9% of binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001015,"raw_usage":{"total_tokens":4340,"prompt_tokens":1056,"completion_tokens":3284,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":672,"completion_tokens_details":{"reasoning_tokens":3224}},"tokens_in":672,"tokens_out":3284,"duration_ms":20001,"temperature":1.0,"reasoning_tokens":3224,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:30:54.204831+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A survey of the secondary stars of known hot-Jupiter-hosting binaries with separations 100-2000 AU that finds zero hot Jupiters around the secondaries in ~50 systems would rule out the ~9% prediction at high confidence; equivalently, measuring the conditional probability P(cold Jupiter around the companion | cold Jupiter around the primary) and finding it close to the field rate would show the predicted rate should be reduced by roughly an order of magnitude.","supporting_citations":[{"cited_title":"M., & Rasio, F","cited_arxiv_id":null,"evidence_quote":"Establishes ZLK-mediated hot Jupiter formation and the inclination and tidal-disruption framework used to define success."},{"cited_title":"2012, , 537, A128","cited_arxiv_id":null,"evidence_quote":"Provides the REBOUND N-body integrator used for all simulations in the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the IAS15 adaptive-timestep integrator that resolves close pericenter passages during high-eccentricity phases."},{"cited_title":"2022, , 926, L17","cited_arxiv_id":null,"evidence_quote":"Defines the q<0.1 AU pericenter criterion used to identify a hot Jupiter in the simulations."},{"cited_title":"The Case for Edge-On Binaries: An Avenue Toward Comparative Exoplanet Demographics","cited_arxiv_id":"2503.08583","evidence_quote":"Cited as evidence that the correlation of planet formation across twin binary components is unconstrained, which the paper acknowledges as a caveat."}],"review_version":1}