{"id":"4916b807-99b7-42df-a339-d052f42f2618","arxiv_id":"2507.13046","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First full merger simulations of 2:1 and 3:1 black hole binaries in scalar-Gauss-Bonnet gravity, with weak-coupling dephasing matching PN predictions but strong-coupling results limited by initial-data transients.","lead":"This paper reports the first full numerical simulations of black hole mergers with mass ratios 2:1 and 3:1 in scalar-Gauss-Bonnet gravity, a modified theory where black holes can carry scalar hair. The authors show that weaker couplings produce dephasing matching post-Newtonian predictions, while stronger couplings are contaminated by eccentricity from imperfect initial data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"All quantitative claims rest on an unvalidated spatial regulator: the sGB coupling is turned off inside each horizon (Sec. II D), with no sensitivity test for the turn-off profile, so the PN match and eccentricity attribution could be regulator artifacts.","rationale":"The paper is honest about limitations, and the existence of merger evolutions for q=2 and q=3 is credible given the convergence test in Appendix A for q=2, lambda=0.04. However, the quantitative headline—PN agreement at weak coupling and eccentricity-driven deviations at strong coupling—is only meaningful if the exterior sGB dynamics are faithfully simulated. The interior turn-off is the only place where the simulated equations differ from the physical theory, and the paper provides no evidence that this difference is causally isolated in practice. In continuum GR, a modification strictly inside the event horizon should not affect future null infinity, but the implementation here is not shown to satisfy that condition: the paper does not specify the mask, and the apparent horizon (used colloquially in the text) is a slice-dependent surface, not the event horizon. A fixed radial mask around a moving puncture can drift relative to the horizon; the transition region can also be under-resolved because the small BH horizon is only ~80 points across. The reader's weakest assumption is exactly this, and I agree. I do not think the concern warrants rejection, because the authors acknowledge it and because a low-coupling run without the mask is a clean, feasible check; the current conditional verdict is the right one unless and until that check is done. I therefore recommend no change to the verdict (UNCHANGED), with the condition made explicit: the quantitative claims require the mask-sensitivity test.","tokens_in":18515,"tokens_out":8763,"duration_ms":107642,"concrete_test":"Run the q=2, lambda/m1^2=0.04 case with the spatial turn-off suppressed entirely (no mask), which should be well-posed at this low coupling, and compare the dephasing and (2,2) waveform against the masked run; then repeat the masked run with the turn-off radius moved between, say, 0.5 r_AH and 0.9 r_AH (or with a wider transition) at both low and high coupling. If the dephasing differs by more than the LR-MR-HR truncation error quantified in Appendix A, the regulator is not innocuous; if it does not differ, record the result as the missing validation. For the high-coupling case, vary the mask radius between 0.5 and 0.9 r_AH and check whether the 10-20% eccentricity-induced spread in dephasing changes; this isolates whether the eccentricity attribution is robust to the regulator.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—first q=2 and q=3 sGB mergers, with weak-coupling dephasing matching PN and strong-coupling deviations attributed to initial-data eccentricity—depends on the assumption that the spatial regulator used to avoid elliptic regions is causally inert. In Section II D the authors state that they 'smoothly turn off the coupling within the BH horizon' and immediately warn that 'there is a risk of introducing unphysical errors if the turn-off region is not fully contained within a well-resolved horizon.' No implementation detail (profile shape, transition width, how the 'horizon' is located in practice, whether it tracks the apparent horizon or is a fixed coordinate mask around each puncture) and no convergence or sensitivity test for this regulator are reported. The concern is concrete: if the mask is tied to the initial puncture coordinates rather than to the instantaneous horizon, it can approach or overlap the exterior during the late inspiral or during common-horizon formation; if the transition is under-resolved, it can add spurious curvature and scalar gradients near the extraction region. Either failure would contaminate the dephasing shown in Figures 2-4 and the ringdown fits in Figure 8. Since the same regulator is present in every sGB run but absent in the GR reference runs, a naive GR-vs-sGB comparison cannot separate physical sGB effects from regulator artifacts. The footnote that the same fix helps the code of [54] is reassuring for stability but does not establish that the exterior waveform is unaffected. Thus the quantitative content of the paper's strongest claim is conditioned on an untested ad hoc choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports numerical relativity evolutions of unequal-mass binary black hole mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet gravity, using the GRFolres extension of GRChombo with a modified CCZ4 formulation. It presents the first claimed full simulations through merger for mass ratios 2:1 and 3:1, together with equal-mass runs, and compares the gravitational-wave dephasing between GR and sGB with post-Newtonian predictions. For the weak-coupling q=2 case the dephasing is reported to agree with PN; for stronger couplings the deviations are attributed to eccentricity introduced by the initial-data scalar hair transient rather than to physical strong-field effects. The paper also studies scalar radiation modes and the ringdown, and devotes significant discussion to technical challenges including a coupling turn-off inside horizons, a slow turn-on of the coupling, and alignment-frequency choices.","tokens_in":18875,"tokens_out":6865,"duration_ms":79645,"significance":"If the technical claims are correct, this is a noteworthy advance: it extends sGB binary evolutions to unequal masses through merger, which previous work could not do at higher mass ratios. The paper is candid about its limitations, releases the code, and includes a clean convergence study in Appendix A showing fourth-to-sixth order convergence for the q=2 weak-coupling orbital phase; the independent PN comparison in Figure 3 is an appropriate validation target. However, the central quantitative claims are currently supported only weakly by the evidence: the horizon turn-off regulator used in all sGB runs is not validated, and the dephasing results lack quantified error bars. The contribution is therefore potentially important but needs the requested sensitivity and uncertainty analysis to be accepted as definitive.","major_comments":[{"comment":"The spatial regulator that smoothly turns off the sGB coupling inside the horizon is an ad hoc modification present in every sGB run and absent in the GR reference runs, yet the paper gives no implementation details (functional form of the profile, transition width, how the horizon is identified, whether the mask tracks the instantaneous horizon) and no sensitivity or convergence test for it. The authors themselves warn that unphysical errors can arise if the turn-off region is not contained within a well-resolved horizon. Because the same regulator is used in all sGB evolutions, it could contaminate the dephasing match in Figure 3 and the eccentricity attribution in Section III B; please add a sensitivity study that varies the turn-off profile and its width, and ideally a test with the regulator active in a GR run, to demonstrate that the exterior physics is unaffected.","section":"Section II D and footnote 1"},{"comment":"The dephasing curves are presented without quantified uncertainty, so the statements of \"good agreement\" with the PN prediction and \"artificial deviations\" for large coupling are not supported by error estimates. Appendix A only establishes convergence for the q=2, lambda/m_1^2=0.04 orbital phase, and not for the strong-coupling or unequal-mass runs; in addition, the alignment frequency f0 is chosen by hand and the text states that its value can change the dephasing significantly (Section III E). Please provide uncertainty bands from resolution differences, from allowed variations of f0, and from the slow turn-on variants, and show these bands on the dephasing curves.","section":"Section III E and Figures 2-4"},{"comment":"The attribution of the strong-coupling dephasing deviations to initial-data eccentricity is not directly established. The slow turn-on runs alter both the measured eccentricity and the dephasing, but correlation does not demonstrate causation, and no comparison is made with initial data from which the eccentricity has been removed or reduced by an independent method. A concrete test, for example evolving the same binary with eccentricity-reduced initial data or estimating the dephasing contribution of the residual eccentricity using the GR runs, is needed before the conclusion that the deviations are \"artificial\" can be accepted.","section":"Section III B and Appendix B"},{"comment":"The q=3 case is presented as a full simulation through merger, but the paper explicitly says that \"the constraint violations are higher\" and that the dephasing \"is also closer to the numerical error, and so we do not analyse it in detail.\" The abstract and the introductory claims should be qualified so that the reader understands that the detailed quantitative conclusions apply only to the q=1 and q=2 runs, and that the q=3 claim is limited to a successful evolution through merger.","section":"Section III B and the abstract"}],"minor_comments":[{"comment":"The caption says \"simulations presented in the papers\"; this should be \"in this paper.\"","section":"Table I caption"},{"comment":"The abbreviations sGB and EsGB are used interchangeably; please define both or choose one and use it consistently.","section":"Throughout"},{"comment":"The caption refers to the \"average value of the weak coupling condition (11)\", but Eq. (11) is an inequality; please specify exactly which quantity is averaged and plotted.","section":"Figure 1, top panel"},{"comment":"The legend does not cleanly distinguish the line styles for immediate versus slow turn-on in the left panel; please clarify.","section":"Figure 4"},{"comment":"For the q=2 ringdown fit, the imaginary part moves away from the quoted sGB prediction (0.0801) rather than toward it, so the statement that both fitted values shift \"in the right direction\" is inaccurate for at least one component.","section":"Section III C and Figure 8"},{"comment":"The sentence \"We understand from discussions with the authors of [54]...\" is informal; if this information is retained, it should be attributed precisely or rephrased as a personal communication.","section":"Footnote 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and does not show citation problems. The main risk is the unvalidated horizon turn-off regulator, which is load-bearing for the dephasing and eccentricity claims. I recommend requesting a focused sensitivity study of that regulator and explicit error bars for the dephasing rather than a full rewrite, because the required analysis is local and testable; if the authors cannot provide it, the dephasing claims should be downgraded to qualitative statements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. First, this paper really does deliver the first full merger evolutions for 2:1 and 3:1 mass-ratio binaries in shift-symmetric scalar-Gauss-Bonnet gravity. That is a genuine capability advance, and the central existence claim is credible. Second, the quantitative dephasing results are more fragile than the headline suggests, because they rely on an ad hoc spatial regulator—turning off the sGB coupling inside each horizon—that is never given a sensitivity test.\n\nThe paper does several things well. The convergence study in Appendix A shows clean fourth-to-sixth order convergence for the q=2, weak-coupling case, which supports the dephasing numbers in that regime. The comparison against independent post-Newtonian predictions is done honestly; the low-coupling match is a real result. The slow turn-on study is also useful: the authors show it helps for q=1 but does not universally improve agreement, which is a nuanced and credible finding. And they are refreshingly transparent about the q=3 run, stating outright that constraint violations are too large for a detailed dephasing analysis. The scalar-charge scalings and ringdown fits are consistent with expectations and are presented with appropriate error caveats.\n\nThe soft spots are real but not fatal. The horizon turn-off is the biggest one: the paper warns that errors can creep in if the turn-off region is not fully inside a well-resolved horizon, but it never shows a test varying the profile or width. Since the same regulator is present in every sGB run and absent in the GR references, a skeptic cannot fully rule out contamination of the dephasing. However, the argument that the regulator is causally disconnected from the exterior is physically reasonable, and the paper's own caveats prevent it from overclaiming. The dephasing curves also lack quantified error bars, and the alignment frequency is chosen by hand, though the authors state that small variations do not affect the results. These are limitations, not fatal flaws: the paper's main contribution is showing that such evolutions are possible and identifying where the problems lie, not producing final data-analysis-ready waveforms.\n\nWho is this for? Numerical relativists working on beyond-GR theories, and waveform modelers who need to know what is feasible and what is not. It deserves a serious referee. My recommendation: send it to peer review, and ask the authors to add a sensitivity study of the turn-off profile, or at least a careful justification of why it is causally inert, plus a discussion of the numerical uncertainty in the dephasing. That would turn a solid capability paper into a fully quantitative one.","headline":"First credible q=2 and q=3 sGB merger evolutions, with honest caveats; the quantitative dephasing is provisional because the horizon turn-off regulator is untested.","tokens_in":19387,"tokens_out":1558,"would_cite":true,"duration_ms":20097,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.25.D-","04.30.-w","04.50.Kd","04.70.Bw"],"model":"deepseek-v4-flash","headline":"First unequal-mass black hole mergers simulated in modified gravity","keywords":["scalar-Gauss-Bonnet gravity","binary black hole mergers","unequal mass ratio","numerical relativity","gravitational wave dephasing","initial data eccentricity","black hole scalarization","modified CCZ4 formulation"],"falsifier":"Run the $q=2$, $\\lambda/m_1^2=0.106$ merger with the coupling turn-off boundary placed at two different radii well inside the horizon, keeping all other settings fixed; if the frequency-domain dephasing or the ringdown frequency changes by more than the measured convergence error, the horizon turn-off contaminates the exterior waveform and the central quantitative claim fails.","tokens_in":1759,"feed_emoji":"🕳️","tokens_out":2055,"duration_ms":118053,"temperature":0.7,"pith_summary":"Black holes in scalar-Gauss-Bonnet gravity carry scalar hair that changes the inspiral and merger, but numerical simulations have struggled to carry unequal-mass binaries through the merger. This paper reports the first full evolutions of 2:1 and 3:1 mass-ratio binaries through merger in shift-symmetric scalar-Gauss-Bonnet gravity, using a well-posed formulation that suppresses the coupling inside each horizon. At weak coupling the gravitational-wave dephasing relative to general relativity matches post-Newtonian predictions, supporting the analytic description in that regime. At larger coupling, close to the effective-field-theory limit, apparent deviations from post-Newtonian theory are traced to eccentricity introduced by the transient formation of scalar hair in the initial data, not to new physics. The paper argues that reliable precision waveforms will require scalarized, quasistationary, constraint-satisfying initial data with eccentricity reduction.","feed_headline":"First unequal-mass black hole mergers simulated in modified gravity","feed_subtitle":"Weak-coupling dephasing matches post-Newtonian predictions; strong-coupling deviations trace to initial-data artifacts.","key_machinery":"The load-bearing mechanism is the ad hoc horizon turn-off: the scalar-Gauss-Bonnet coupling is smoothly suppressed inside each black hole horizon so that the weakly coupled, hyperbolic regime never breaks down, avoiding the elliptic regions that previously terminated unequal-mass evolutions at merger. A companion diagnostic is the weak-coupling condition $\\sqrt{|\\lambda f'(\\phi)|}/L \\ll 1$, where $L^{-1}$ is the largest curvature or scalar-gradient scale, which monitors whether the simulation stays inside the effective-field-theory's regime of validity. The dephasing measurement is carried by aligning general-relativity and scalar-Gauss-Bonnet waveforms at a chosen frequency $f_0=0.01/M_{\\mathrm{ADM}}$ and comparing the orbital phase as a function of frequency against post-Newtonian expressions from prior work.","core_discovery":"The central claim is that unequal-mass black-hole binaries in shift-symmetric scalar-Gauss-Bonnet gravity can be evolved through merger with a modified CCZ4 scheme, provided the Gauss-Bonnet coupling is smoothly turned off inside each apparent horizon to prevent the formation of elliptic regions. For the mass ratio $q=2$ at weak coupling $\\lambda/m_1^2=0.04$, the dephasing of the (2,2) gravitational-wave mode relative to general relativity agrees well with post-Newtonian predictions across the inspiral. At the larger coupling $\\lambda/m_1^2=0.106$, near the weak-coupling limit, the measured dephasing deviates from post-Newtonian values, but the paper attributes this to eccentricity induced when the initially vanishing scalar field grows into hair and disturbs the binary trajectory; the ad hoc slow turn-on of the coupling helps in some cases and hurts in others, so it is not a reliable fix. The $q=3$ case also reaches merger, though with higher constraint violations and dephasing close to numerical error. The scalar monopole and dipole radiation scale with mass ratio as expected, and the post-merger ringdown shifts in the predicted direction but lies within numerical error.","pith_inferences":["If the horizon turn-off is truly exterior-innocuous, the same strategy should extend stable evolutions to mass ratios beyond 3:1 and to other higher-curvature theories with elliptic-region breakdowns; this is testable by varying the turn-off shell radius.","The paper's identification of initial-data eccentricity as the dominant systematic suggests that current high-coupling dephasing values may shift once scalarized quasistationary initial data are used, potentially improving or eroding post-Newtonian agreement.","Because the remnant scalar charge is approximately independent of mass ratio after merger while the inspiral amplitude scales with the product of the masses, the jump in scalar monopole amplitude could serve as a clean probe of merger dynamics in future detections.","Once the initial-data problem is solved, the same pipeline can produce mismatches and parameter-estimation biases for future gravitational-wave detectors, quantifying how well parametrized general-relativity tests capture scalar-Gauss-Bonnet physics."],"forward_implications":["2:1 and 3:1 unequal-mass scalar-Gauss-Bonnet binaries can now be followed stably through merger, so beyond-GR waveform modeling extends beyond equal-mass, nonspinning cases.","In the weakly coupled regime, nonlinear numerical dephasing confirms post-Newtonian predictions, so post-Newtonian inspiral models can be used for data analysis in that regime.","Large-coupling dephasing discrepancies should be read as initial-data artifacts rather than as evidence for novel strong-field physics, until scalarized quasistationary initial data become available.","Reliable beyond-GR waveform catalogs will need constraint-satisfying scalarized initial data with eccentricity reduction; without them, dephasing measurements are sensitive to the alignment frequency and can be biased.","Post-merger ringdown frequency shifts are in the predicted direction but smaller than current numerical error, so accurate quasinormal-mode tests in this theory require higher resolution and more refined extraction."],"supporting_citations":[{"why":"Supplies the prior nonlinear sGB binary merger study that broke down for higher mass ratios, along with the post-Newtonian dephasing expressions used for comparison.","marker":"[54]"},{"why":"Establishes the first full nonlinear evolution of Einstein-scalar-Gauss-Bonnet gravity, the simulation approach this work extends.","marker":"[44]"},{"why":"Provides the well-posed modified CCZ4 formulation in singularity-avoiding coordinates that keeps the evolution hyperbolic.","marker":"[28]"},{"why":"Provides the puncture gauge formulation for Einstein-Gauss-Bonnet and four-derivative scalar-tensor theories used in these black-hole evolutions.","marker":"[29]"},{"why":"Gives constraint-satisfying initial data for modified gravity with an arbitrary scalar profile, but not a quasistationary state, identifying the missing ingredient.","marker":"[68]"},{"why":"Gives quasistationary scalar hair for binary black hole initial data in sGB gravity, the solution the paper argues is needed for reliable waveforms.","marker":"[69]"},{"why":"Supplies the 3PN puncture momenta used to construct the initial binary orbits with roughly seven orbits and low eccentricity.","marker":"[76]"},{"why":"Provides the astrophysical constraints that set the high-coupling value near the observational limit.","marker":"[66]"}],"fun_headline_variants":["First 2:1 and 3:1 black hole mergers simulated in sGB gravity","Unequal-mass BH mergers in sGB gravity: first simulations","Weak-coupling dephasing matches PN predictions in sGB binary mergers","Strong-coupling dephasing in sGB mergers traced to initial-data eccentricity","Simulating unequal-mass black hole binaries in scalar-Gauss-Bonnet gravity"],"cache_read_input_tokens":21504,"weakest_assumption_plain":"The results depend on the assumption that suppressing the modified-gravity coupling inside each black hole's horizon has no effect on the outside spacetime and gravitational waves, even though this suppression is an ad hoc procedure with no independent check of its innocuousness.","fun_headline_variants_meta":{"raw":{"variants":["First 2:1 and 3:1 black hole mergers simulated in sGB gravity","Unequal-mass BH mergers in sGB gravity: first simulations","Weak-coupling dephasing matches PN predictions in sGB binary mergers","Strong-coupling dephasing in sGB mergers traced to initial-data eccentricity","Simulating unequal-mass black hole binaries in scalar-Gauss-Bonnet gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00076,"raw_usage":{"total_tokens":3382,"prompt_tokens":961,"completion_tokens":2421,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":2320}},"tokens_in":577,"tokens_out":2421,"duration_ms":18138,"temperature":1.0,"reasoning_tokens":2320,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:31:24.754757+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $q=2$, $\\lambda/m_1^2=0.106$ merger with the coupling turn-off boundary placed at two different radii well inside the horizon, keeping all other settings fixed; if the frequency-domain dephasing or the ringdown frequency changes by more than the measured convergence error, the horizon turn-off contaminates the exterior waveform and the central quantitative claim fails.","supporting_citations":[],"review_version":1}