{"id":"29cc7b8c-25b4-4c4b-803c-39dfa2bbac5f","arxiv_id":"2507.18934","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A binary black hole moving through scalar dark matter experiences a drag force that is not the simple sum of two individual black hole drags, because the two gravitational wakes interfere.","lead":"This paper simulates two black holes orbiting inside a cloud of lightweight scalar dark matter and measures the drag and torque the dark matter exerts on them. It provides the first relativistic numerical benchmarks for how such dark matter could slightly alter gravitational waves from black hole mergers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline claim that binary drag is not a simple superposition of isolated black holes is asserted but never directly tested against single-BH runs.","rationale":"The reader's weakest assumption was the 2D vertical integration, which is a legitimate concern about quantitative accuracy and is acknowledged in the conclusion. However, the most load-bearing issue for the advertised central claim is that the non-superposition statement is never compared against a single-BH baseline. Even a fully 3D simulation would not establish that the binary drag is not a simple superposition unless the single-BH sum were computed. The appendix validation mismatch is also real and worrisome: the convergence and extraction-radius tests appear to use different configurations than the production runs, so confidence in the numerical values is weaker than the text suggests. I still agree with the CONDITIONAL verdict: the work is a useful first exploration, the conservation checks in Appendix C for Case A provide some independent support, and the qualitative scenario is plausible. But the requested condition should include an explicit single-BH comparison, not only a 3D cross-check or clarification of the appendices.","tokens_in":18375,"tokens_out":9410,"duration_ms":106996,"concrete_test":"Run the same wind-tunnel setup for two isolated single BHs of masses M1 and M2, using identical v, mu, grid resolution, box size, and extraction radii as Case A; compute F_x, F_y, and tau_z for each, sum the two single-BH results with the appropriate orbital phases and center-of-mass lever arms, and compare the sum to the orbit-averaged binary Case A values in Fig. 5. If the summed single-BH result agrees with the binary result within the reported error bars, the central non-superposition claim is not supported; if it differs by more than the error bars, the claim is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that the binary's drag is not a simple superposition of two isolated black holes—is stated in the abstract and conclusion but is not actually tested in Sec. III. No comparison is made between the binary force/torque and a sum of two single-BH wind-tunnel results at the same v, mu, resolution, and extraction prescription, nor against the single-BH results of Refs. [21,56,57]. The evidence offered instead consists of orbital-period flux modulation, a dual beat frequency attributed to q=0.8, and separation/velocity trends in Fig. 5. None of these discriminate between true companion-modified wakes and a superposition model, since even independent single-BH drag forces would yield separation-dependent, orbit-modulated forces and a nonzero net torque for an unequal-mass binary. Thus the phrase 'significant nonlinearities' is an interpretation rather than a measured result. A related internal-consistency issue compounds this: Appendix A describes a convergence test on an equal-mass, mu_M_tot=0.05 binary, and Appendix B uses a spinning BH with a/M=0.7, while the production runs are q=0.8 and non-spinning. Those tests therefore do not directly validate the configurations from which the headline results are drawn. If the non-superposition claim is the paper's main contribution, it currently lacks the direct comparison needed to support it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents two-dimensional general-relativistic simulations of a complex scalar field on a fixed, constraint-satisfying binary black hole spacetime, with the binary placed in a 'wind tunnel' of asymptotically homogeneous scalar dark matter. The authors extract the drag force, torque, mass accretion, and scalar charge accretion as functions of the background velocity, scalar mass, and binary separation, covering four parameter regimes. Their central qualitative claim is that the drag on a binary is not the simple superposition of the drags on two isolated black holes, since the companion modifies the gravitational wake and produces nonlinear interference effects; they further argue that this additional force and torque could, in principle, dephase gravitational wave signals. The results include a parameter survey (Fig. 5), a Magnus-force measurement (Fig. 6), and conservation checks of linear and angular momentum (Appendix C).","tokens_in":18593,"tokens_out":3183,"duration_ms":33579,"significance":"If the central claim holds, this is the first relativistic treatment of scalar dark matter dynamical friction on a binary black hole, and it provides a concrete starting point for environmental corrections to waveform models in ultralight scalar dark matter scenarios. The paper's strengths include explicit conservation checks (Appendix C) showing momentum and angular momentum balance to numerical precision, a resolution study (Appendix A) and an extraction-radius study (Appendix B), and a systematic parameter survey across velocity, scalar mass, and separation. The simulations are also run for many orbital periods, giving a well-defined steady state for the orbital averages. However, as detailed below, the headline non-superposition claim is not directly tested, the 2D reduction is not validated against 3D for scalar fields, and the validation tests are run on configurations different from the production runs. These issues are load-bearing for the quantitative and qualitative conclusions, so the paper needs major revision.","major_comments":[{"comment":"The central claim that 'the binary's drag is not a simple superposition of two isolated black holes' is asserted in the abstract and conclusion but is never directly tested. No comparison is made between the binary force/torque and the sum of two single-black-hole wind-tunnel results at the same masses, velocity, scalar mass, resolution, and extraction prescription, nor against the single-BH results of Refs. [21,56,57]. The evidence offered—orbital-period flux modulation, a dual beat frequency attributed to q=0.8, and separation/velocity trends in Fig. 5—is consistent with a superposition of independent single-BH wakes, since even two separated, unequal-mass single-BH drag forces would produce orbit-modulated, separation-dependent total forces and a nonzero net torque. I recommend adding a direct control: a single-BH run with the same code and extraction at mass M1 and M2, summing the two forces, and comparing to the binary run. Unless such a test is performed, the phrase 'significant nonlinearities' remains an interpretation rather than a measured result.","section":"Abstract and Sec. III (Figs. 3-5)"},{"comment":"The quantitative results rest on a two-dimensional (vertically integrated) reduction of the scalar field problem, as stated in Sec. II.F: 'we follow similar works for hydrodynamic torque calculations... using vertically integrated accretion flows.' This approximation is reasonable for thin hydrodynamic disks, but the scalar wake and its interference around a relativistic binary are genuinely three-dimensional phenomena, and the manuscript provides no 3D scalar-field comparison or estimate of the resulting error. The Conclusion acknowledges this only as a future direction. Since the drag and torque values in Fig. 5 are the main quantitative output, the 2D reduction should either be validated with at least one representative 3D run or be accompanied by a clear statement of the expected magnitude of 3D corrections.","section":"Sec. II.F and Conclusion"},{"comment":"The numerical validation does not cover the production configurations from which the headline results are drawn. Appendix A describes a convergence test on an equal-mass, non-spinning binary with mu M_tot = 0.05, while the production runs use q = 0.8 and mu M = 0.2, 0.8, and 0.05. Appendix B tests extraction-radius sensitivity using a spinning BH with a/M = 0.7, although the production binaries are non-spinning. Moreover, Sec. II.F states that the reference convergence simulation has mu M = 0.2, d_BBH = 26M, v = 0.5c, which contradicts the equal-mass mu M_tot = 0.05 setup described in Appendix A. As a result, there is no convergence or inner-boundary study at the actual mass ratio, spin, and scalar mass of the runs that support the central claims. I request a resolution study (at least two additional resolutions) and an extraction-radius test for a representative production configuration, e.g., q = 0.8, mu M = 0.2, d_BBH = 26M, v = 0.5c.","section":"Appendix A, Appendix B, and Sec. II.F"},{"comment":"The velocity dependence of the drag force is described as 'approximately linear' near v ~ 0.5 and 'saturating' above v = 0.7c, but no fit or quantitative criterion is provided beyond a blue dashed line in the top-left panel of Fig. 5. Since the paper later invokes 'a modified dynamical friction law,' the functional form of the velocity dependence should be stated explicitly (e.g., the best-fit power-law or the analytic expression from Ref. [56]) so that the claim can be checked. The current presentation leaves the velocity scaling ambiguous and the 'saturation' claim unsupported without error bars on the fit.","section":"Sec. III (Fig. 5 and text)"}],"minor_comments":[{"comment":"In Sec. II.A, 'helliptical Killing vector' should be 'elliptical Killing vector.'","section":"Sec. II.A"},{"comment":"In Sec. II.F, 'Table II F' should be 'Table I' (or a proper table number), and the table caption should be self-contained.","section":"Sec. II.F"},{"comment":"In Sec. III, 'T he accretion is generally small' contains a spacing typo and should read 'The accretion is generally small.'","section":"Sec. III"},{"comment":"In the Conclusion, 'It would be interested to systematically compare' should be 'It would be interesting to systematically compare.'","section":"Sec. IV (Conclusion)"},{"comment":"The 'dual beat frequency' attributed to the mass ratio q = 0.8 is not defined quantitatively; the text should state the expected beat frequency or the period of the amplitude modulation so that the reader can verify the claim from the plotted time series.","section":"Sec. III (Fig. 4)"},{"comment":"Reference [12] appears with odd spacing ('M uffled murmurs') in the bibliography; this is a formatting artifact that should be cleaned before publication.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable first step and the numerical infrastructure appears solid, particularly the conservation checks in Appendix C. However, the central non-superposition claim lacks a direct test, and the validation runs do not match the production parameters. These are fixable within the scope of a revision, but they are load-bearing for the paper's main message. I would encourage the editor to request the single-BH control run and a convergence/extraction study at the production mass ratio, in addition to a representative 3D benchmark or a clearly quantified estimate of the 2D error. The mismatch between Sec. II.F's stated reference simulation and Appendix A's actual setup is a red flag that needs clarification, not just a change of wording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first relativistic simulation of scalar-field dark matter drag on a binary black hole, with a clean parameter sweep and credible conservation checks, but the advertised result that the binary drag is not a superposition of single-black-hole drags is never directly tested.\n\nWhat's new: previous GR simulations of scalar DM drag treated a single BH (Traykova et al. 2021, 2023; Wang et al. 2024). This work extends to a binary, in a wind-tunnel setup on a fixed metric from FUKA initial data. It maps drag force, torque, mass and charge accretion against velocity, separation, and scalar mass, and identifies interesting qualitative features: orbital-period modulation of the fluxes, a dual beat frequency tied to q=0.8, a sign-changing transverse force (Magnus-like), and a regime classification by de Broglie wavelength vs separation. The numerics look solid: AMReX/GPU evolution, convergence and extraction-radius tests, and momentum/angular-momentum conservation checks in App. C that close to the expected level. The paper is also honest about its fixed-metric approximation and about the small astrophysical effect under typical densities.\n\nSoft spots, in order of importance. First, the paper's central claim—'not a simple superposition of two isolated black holes'—is not directly tested. There is no comparison between the binary force/torque and the sum of two single-BH wind-tunnel results at the same v, mu, resolution, and extraction radius. The evidence offered (modulated fluxes, beat frequency, separation/velocity trends) would also arise from independent single-BH drags on an unequal-mass binary, so it does not discriminate. The abstract should either soften the claim or include a direct superposition test. This is a real gap, but not fatal to the paper's other content. Second, the validation tests do not match production runs: App. A uses an equal-mass, muMtot=0.05 binary, and App. B uses a spinning BH with a/M=0.7, while all production runs are q=0.8 and non-spinning. The authors need to clarify whether these are code checks or physical validation of the actual configurations. Third, the 2D vertically-integrated reduction is not cross-checked in 3D; the authors acknowledge this, and it undermines the quantitative torque values, though not necessarily the qualitative behavior.\n\nVerdict: this is a useful first exploration and deserves peer review. I would send it back for a revision that either adds the single-BH superposition test or moderates the claim, and that clarifies which validation runs apply to which production setups. The 2D issue can remain as a stated limitation for now.","headline":"First relativistic binary scalar-DM drag simulations with solid numerics, but the headline non-superposition claim is not directly tested against single-BH runs.","tokens_in":19133,"tokens_out":2918,"would_cite":true,"duration_ms":30307,"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":"The drag that scalar dark matter exerts on a binary black hole is not the sum of the drags on each black hole alone, and the extra torque can dephase the binary's gravitational-wave signal.","keywords":["scalar dark matter","dynamical friction","binary black hole","gravitational wave dephasing","Klein-Gordon equation","general-relativistic simulation","accretion drag","ultralight bosons"],"falsifier":"Run the same fiducial case (scalar mass times total mass 0.2, binary separation 26M, wind speed 0.5c) in a full three-dimensional general-relativistic scalar-field simulation with identical extraction surfaces, and compare the orbit-averaged drag force and torque; if the three-dimensional values fall outside the quoted error bars of the two-dimensional run, the quantitative claim fails.","tokens_in":18131,"feed_emoji":"🕳️","tokens_out":9720,"duration_ms":96568,"temperature":0.7,"pith_summary":"This paper argues that when a binary black hole moves through a homogeneous background of relativistic scalar dark matter, the dynamical friction and torque it experiences cannot be obtained by adding the single-black-hole results: the companion perturbs the gravitational wake, and the resulting interference produces a binary-specific drag. To show this, the authors run two-dimensional general-relativistic simulations of a binary in a wind tunnel, with a complex massive scalar field streaming past a fixed binary spacetime, and extract the drag force, torque, and mass and charge accretion for a range of binary separations, scalar masses, and wind speeds. They find an additional force and a spin-down torque that, in principle, accelerate the inspiral and imprint a dephasing in the gravitational-wave signal. The point of the exercise is that future searches for dark matter around stellar-mass binaries need binary-specific environmental corrections, not a superposition of isolated-hole formulas.","feed_headline":"Binary black holes feel a drag no single hole can produce","feed_subtitle":"Scalar dark matter drag on a binary is nonlinear, so gravitational-wave dephasing needs binary-specific models.","key_machinery":"The load-bearing object is the relativistic momentum and angular-momentum balance law built from the scalar stress-energy tensor. For a chosen approximate Killing vector, the paper casts conservation as a continuity equation and decomposes the force on the binary into a volume-integrated drag sourced by the curved spacetime, a surface accretion flux through the two excision spheres, and an outer-boundary flux. Together with the Noether current of the complex scalar, this bookkeeping turns the raw field evolution into the forces, torque, and accretion rates that are compared across parameters. The wind-tunnel setup, an extended conformal thin-sandwich binary spacetime Lie-dragged so the black holes stay on a circular orbit while a plane-wave scalar field streams past, is the stage on which this bookkeeping operates, and the ratio of de Broglie wavelength to binary separation organizes the results into Fraunhofer-like, Fresnel-like, and unimpeded regimes.","core_discovery":"The central claim is that the drag force on a binary black hole in a scalar dark matter wind is nonlinear in the binary: the companion changes the structure of the wake, so the total drag and torque differ from the sum of two isolated black holes. The paper demonstrates this by computing, in a fixed binary spacetime, the momentum and angular-momentum exchange between a complex massive scalar field and the binary, separated into a volume-integrated drag, a horizon accretion flux, and an outer-boundary flux. Over a parameter survey of separation, scalar mass, and wind speed (about 0.2c to 0.7c), the orbit-averaged drag is roughly linear in velocity near 0.5c, saturates for larger velocities, and is suppressed when the scalar de Broglie wavelength is much smaller than the binary separation. The torque always acts to spin the binary down, has a non-monotonic velocity dependence, and there is a transverse anti-Magnus force whose sign can flip in some parameter ranges. Accretion of scalar particles is confirmed by a positive Noether charge flux, and mass accretion becomes important for heavy scalars at low speeds.","pith_inferences":["Beyond the paper, the 3D question is the first test: if a fully three-dimensional wake changes the orbit-averaged numbers, the magnitudes here are estimates even if the qualitative non-additivity survives.","Beyond the paper, the diffraction analogy points to a searchable parameter window: binaries with separation within roughly one de Broglie wavelength of the scalar should show the largest dephasing, which future space-based detectors could target.","Beyond the paper, the half-orbital-period alternation tied to the mass asymmetry of the binary suggests that a measured beat in the environmental dephasing could carry information about both the dark matter and the binary's mass ratio.","Beyond the paper, since the authors use a complex scalar to avoid Compton-frequency oscillations, a real scalar field would add a periodic modulation of the drag; this could be tested with a straightforward modification of their setup."],"forward_implications":["Gravitational-wave models for binaries embedded in scalar dark matter must include a binary-specific drag and torque; summing single-black-hole forces will not reproduce the dephasing.","The environmental correction is strongest when the scalar de Broglie wavelength is comparable to the binary separation and for wind speeds near 0.5c, and suppressed when the wavelength is much shorter than the separation.","The scalar field consistently removes orbital energy and angular momentum, so it accelerates the inspiral and acts as a spin-down torque on the binary.","Under typical galactic dark matter densities the effect is tiny compared with gravitational-wave energy loss, but in dense scalar clouds or spikes the induced dephasing could be detectable by future space-based detectors.","The transverse Magnus-like force and its possible sign change mean scalar environments can also deflect the binary's motion, not just slow it down."],"supporting_citations":[{"why":"The single-black-hole relativistic scalar dynamical friction simulations that define the baseline the binary results are compared against.","marker":"[56]"},{"why":"Extends the single-hole drag to clouds of all sizes and supplies the analytical accretion-efficiency expectation used to interpret the binary's mass accretion.","marker":"[57]"},{"why":"The single-Kerr gravitational Magnus force calculation whose sign and velocity dependence the binary's transverse force is compared with.","marker":"[21]"},{"why":"Derives the continuity-equation formalism for general matter that the paper uses to split the force into drag and accretion terms.","marker":"[18]"},{"why":"The spectral initial-data solver used to construct the constraint-satisfying extended conformal thin-sandwich binary black hole spacetime.","marker":"[72]"},{"why":"The adaptive mesh refinement infrastructure used to run the high-resolution scalar-field evolutions.","marker":"[77]"},{"why":"The vertically integrated circumbinary accretion simulation approach invoked to justify the two-dimensional setup.","marker":"[78]"},{"why":"Another vertically integrated hydrodynamic torque calculation used to motivate the 2D treatment.","marker":"[79]"}],"fun_headline_variants":["Nonlinear dark matter drag alters binary black hole orbits","Dark matter wake makes binary drag non-additive","Scalar dark matter torque spins down black hole binaries","Gravitational wave dephasing from dark matter drag on binaries","Binary black holes feel dark matter drag nonlinearity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative drag and torque numbers rest on the assumption that a two-dimensional, vertically integrated scalar-field simulation captures the same wake physics a fully three-dimensional simulation would, so a large three-dimensional discrepancy would change the magnitudes.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinear dark matter drag alters binary black hole orbits","Dark matter wake makes binary drag non-additive","Scalar dark matter torque spins down black hole binaries","Gravitational wave dephasing from dark matter drag on binaries","Binary black holes feel dark matter drag nonlinearity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000949,"raw_usage":{"total_tokens":4047,"prompt_tokens":938,"completion_tokens":3109,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":3032}},"tokens_in":554,"tokens_out":3109,"duration_ms":24217,"temperature":1.0,"reasoning_tokens":3032,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:04:35.257642+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same fiducial case (scalar mass times total mass 0.2, binary separation 26M, wind speed 0.5c) in a full three-dimensional general-relativistic scalar-field simulation with identical extraction surfaces, and compare the orbit-averaged drag force and torque; if the three-dimensional values fall outside the quoted error bars of the two-dimensional run, the quantitative claim fails.","supporting_citations":[{"cited_title":"Ferreira, and Lam Hui","cited_arxiv_id":null,"evidence_quote":"The single-black-hole relativistic scalar dynamical friction simulations that define the baseline the binary results are compared against."},{"cited_title":"Annulli, V","cited_arxiv_id":null,"evidence_quote":"Extends the single-hole drag to clouds of all sizes and supplies the analytical accretion-efficiency expectation used to interpret the binary's mass accretion."},{"cited_title":"Ferreira, and Lam Hui","cited_arxiv_id":null,"evidence_quote":"The single-Kerr gravitational Magnus force calculation whose sign and velocity dependence the binary's transverse force is compared with."},{"cited_title":"Bondi and F","cited_arxiv_id":null,"evidence_quote":"Derives the continuity-equation formalism for general matter that the paper uses to split the force into drag and accretion terms."},{"cited_title":"Jens Papenfort, Samuel D","cited_arxiv_id":null,"evidence_quote":"The spectral initial-data solver used to construct the constraint-satisfying extended conformal thin-sandwich binary black hole spacetime."},{"cited_title":"KADATH: A spectral solver for theo- retical physics","cited_arxiv_id":null,"evidence_quote":"The adaptive mesh refinement infrastructure used to run the high-resolution scalar-field evolutions."},{"cited_title":"Most and Alexander A","cited_arxiv_id":null,"evidence_quote":"Another vertically integrated hydrodynamic torque calculation used to motivate the 2D treatment."}],"review_version":2}