{"id":"d71546a0-d6d9-41c5-bc59-940a7da8fb1f","arxiv_id":"2607.17887","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The scalar, electromagnetic, and tensor contributions to the angular momentum flux from noncircular EMd binaries are derived through 1PN order.","lead":"This paper derives the angular momentum lost by a black-hole binary in a modified gravity theory (Einstein-Maxwell-dilaton) as it emits radiation, for elliptical orbits, to first post-Newtonian order. The result supplies a missing ingredient for future gravitational-wave models of such binaries.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (II.23) and (II.24) appear to have the γ and c²/γ factors interchanged relative to action (II.2), giving wrong rest-frame scalings for T^{00} and T_φ; later equations use the corrected forms, but the manuscript should fix this.","rationale":"The reader's weakest_assumption identifies the skeletonized point-particle action (II.2) and its expansion as the key physical input. I agree that this is the right sector to scrutinize, but I find a more concrete and checkable issue inside that sector: the printed matter stress-energy and scalar source equations (II.23) and (II.24) do not follow from the action (II.2). In the rest frame, the printed T^{00} is m c⁴ instead of m c², and the printed T_φ is m α instead of m α c². The later equations in the paper appear to use the correct forms, so this is probably an error in Eq. (II.23)–(II.24) rather than in the central derivation. Nevertheless, because the source densities feed every multipole moment and therefore the angular momentum flux, the discrepancy must be resolved publicly. If the corrected factors match the later equations, the central claim survives; if not, the flux formulas would need re-derivation. This is consistent with a conditional-acceptance verdict: the result is probably correct, but the manuscript should be amended and the omitted algebra shown. I therefore do not move the reader's verdict.","tokens_in":33914,"tokens_out":22284,"duration_ms":209300,"concrete_test":"Recompute δS_m/δg^{μν} and δS_m/δφ from Eq. (II.2) at rest, and verify that (II.23) should read T^{μν}=ρ_g c²/(γ√-g)v^μv^ν and (II.24) should read T_φ=ρ_g γ/√-g [α+(α²+β)φ]. Then trace these corrected densities through (III.14), (III.35), and the moments (III.15)–(III.18), (III.42)–(III.45), re-evaluating the flux (V.13)–(V.25) if any factor of c² is found to have propagated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Starting from the point-particle action (II.2), a direct variational derivative gives T^{μν} = ρ_g c²/(γ√-g) v^μ v^ν and T_φ = ρ_g γ/√-g α, with γ = c√(-g_{μν}v^μv^ν). In the rest frame γ=c², so the physical results are T^{00}=m c² δ³ and T_φ=m α c² δ³. As printed, Eq. (II.23) has γ in the numerator of T^{μν}, yielding T^{00}=m c⁴ δ³; Eq. (II.24) has c²/γ, yielding T_φ=m α δ³. Both printed equations are thus off by a factor c² and are mutually inconsistent with the action. The subsequent source densities — for example the scalar source (III.14) and the gravitational source μ00 in (III.35) — use the corrected factors (γ in T_φ, c²/γ in T^{μν}), so the final flux is likely unaffected if this is a typesetting or transcription error. But because all source multipole moments in Sec. III and hence the flux in Sec. V depend on these matter densities, the manuscript must correct (II.23)–(II.24) and confirm that no c² factor was propagated into later equations. The omitted functional-derivative algebra makes this a concrete, checkable weak point in the load-bearing input, more immediate than the general skeletonization concern. This is not an attack on the physical framework, which is standard, but on the internal consistency of the printed equations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the instantaneous angular-momentum flux for nonspinning, noncircular black-hole binaries in Einstein-Maxwell-dilaton (EMd) theory through relative first post-Newtonian (1PN) order. Working in the Einstein frame and using the direct integration of the relaxed field equations (DIRE) approach, the authors construct the scalar, electromagnetic, and tensor fields in the wave zone, including next-to-leading 1/R² finite-distance terms, and evaluate the resulting multipole moments in the center-of-mass frame. The main results are the separate scalar, electromagnetic, and tensor contributions to the angular-momentum flux, Eqs. (V.13)-(V.25), plus a total flux. The scalar and electromagnetic channels begin at dipole order O(c^{-3}) (with the expected dipole-suppression limits when α1=α2 or q1/m1=q2/m2), while the tensor channel begins at quadrupole order O(c^{-5}). The paper reports three consistency checks: the quasicircular limit satisfies ΩJ_z=F, the GR limit recovers the known generic-orbit result, and the dipole-suppression limits are recovered.","tokens_in":34286,"tokens_out":10124,"duration_ms":110348,"significance":"If the final results are correct, this is a useful and nontrivial extension of post-Newtonian radiation-reaction calculations to a three-channel theory with scalar, electromagnetic, and tensor radiation. The paper is careful about PN counting, provides explicit multipole moments and waveforms, and separately documents the 1/R² wave-zone terms, which are usually omitted. The consistency checks, once verified, would give strong support to the long calculation. The main value is that the flux provides the second dissipative balance equation needed for future eccentric-orbit evolution and phasing in EMd theory. The manuscript is not machine-checked, and the differentiability of the moment algebra is not fully demonstrated, but the level of detail is appropriate for the intended audience.","major_comments":[{"comment":"There is an internal inconsistency in the printed matter densities. From the point-particle action (II.2), a direct variational derivative gives T^{μν}=ρ_g c²/(γ√-g) v^μ v^ν and T_φ=ρ_g γ/√-g α, with γ=c√(-g_{μν}v^μv^ν); in the rest frame γ=c². As printed, Eq. (II.23) has γ in the numerator of T^{μν}, yielding T^{00}~m c⁴δ³, and Eq. (II.24) has c²/γ, yielding T_φ~m α δ³. Both are off by a factor c² and mutually inconsistent with the action. The later source densities (e.g., μ in Eq. (III.14) and μ00 in Eq. (III.35)) use the corrected normalizations, so the final flux is probably unaffected if this is a transcription error. However, because every multipole moment in Sec. III depends on these densities, the manuscript must correct (II.23)-(II.24) and explicitly confirm that no c² factor has propagated into the derived moments and flux formulas.","section":"II.B, Eqs. (II.23)-(II.24)"},{"comment":"The three claimed consistency checks are asserted without showing the reduction. In particular, after Eq. (V.26) the paper states that the quasicircular limit of the computed flux satisfies ΩJ_z=F and agrees with Ref. [25] through 1PN order, but no comparison of the individual terms is shown. Likewise, in the GR limit the statement that the tensor flux reduces to the quadrupole and octupole results of Refs. [52,55] would benefit from explicit matching of the coefficients, especially the 1PN correction jLL_quad in Eq. (V.24). Since these checks are the primary validation of a long multipole calculation, the authors should provide enough algebra in an appendix or as a supplementary derivation for the reader to verify the reductions.","section":"V.D, Sec. V.D consistency checks"}],"minor_comments":[{"comment":"The generalized Lorentz factor γ has dimensions of c² in this convention. This is unusual and should be stated explicitly, as the later equations (II.23)-(II.24) are only dimensionally consistent once this convention is understood.","section":"II.B, Eq. (II.22)"},{"comment":"The quantity δ is defined in Eq. (IV.15) and used extensively in the center-of-mass reduced moments. It would aid the reader to add a sentence explaining that δ is the 1PN center-of-mass correction and that it does not carry any explicit scalar-charge dependence because the dipole condition is imposed on the total effective density μ00.","section":"IV.B, Eq. (IV.15)"},{"comment":"The derivation of the angular-momentum flux formulas is compact. The definition of the transverse derivative ð_T^j and the treatment of the retarded-time implicit dependence should be stated in one or two sentences before Eq. (C.3), so the reader can follow the order-counting that drops R^{-2} terms in the flux.","section":"Appendix C, Eqs. (C.8)-(C.10)"},{"comment":"The abstract and introduction refer to 'post-Newtonian order' and 'relative first post-Newtonian order' almost interchangeably. Since the scalar dipole starts at formally negative PN order relative to the tensor quadrupole, the paper should consistently use 'relative 1PN order' and state this convention near Eq. (III.4).","section":"General"},{"comment":"A few references are incomplete (e.g., Refs. [17,20,26,28,35] lack article numbers or page counts). This is a presentation issue but should be cleaned up before publication.","section":"References"}],"recommendation":"minor_revision","confidential_remarks":"I found the calculation plausible and the final flux formulas likely correct, but the typo in Eqs. (II.23)-(II.24) is concerning because it sits at the very beginning of the derivation. I traced that the later source densities use the corrected factors, so I do not think the main results are invalidated. The other substantive request, to expand the consistency-check algebra, is more about reproducibility than correctness. This is a competent technical paper that would be acceptable after a focused revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something genuinely new: it computes the instantaneous angular momentum flux for generic noncircular, nonspinning binaries in Einstein-Maxwell-dilaton theory through 1PN order, split into scalar, electromagnetic, and tensor channels. That was a real gap—Khalil et al. gave the energy flux and conservative dynamics, and Juli'e handled circular orbits, but nobody had the noncircular angular momentum flux. The DIRE derivation is long, explicit, and structurally sound: multipole moments, waveforms, and final flux formulas are all written out, and the three consistency checks (circular balance, GR limit, dipole suppression) are the right checks to run. The result itself is an enabling step for eccentric phasing and orbital evolution studies in this theory.\n\nThe soft spots are real but manageable. First, Eqs. (II.23) and (II.24) as printed have the γ and c^2/γ factors interchanged relative to what you get from varying the point-particle action (II.2). Direct variation gives T^{μν} ∝ c^2/γ v^μv^ν and T_φ ∝ γ, not the other way around—in the rest frame the printed versions are off by a factor c^2 in each. The later equations (III.14, III.35, etc.) use the correct factors, so this looks like a transcription error rather than a load-bearing flaw, but it needs to be fixed and the subsequent moments re-checked for any propagated factor. Second, the quasicircular and GR-limit comparisons are asserted without showing the algebra. An independent reader cannot verify those limits from the text alone. That is a minor-to-moderate weakness—the result is probably right, but the missing algebra makes the checks less convincing. Third, the point-particle skeletonization with the expansion (II.3) is standard and grounded in Juli'e's matching, so I am not worried about it at this PN order.\n\nWho is this for? People computing waveform phasing or eccentric evolution in modified-gravity theories. They should have it. It deserves a serious referee, not a desk rejection. My recommendation: send it to peer review, with the request that the authors fix the c^2/γ typo and either show the consistency-check algebra or point to a source where it is displayed.","headline":"First computation of the noncircular EMd angular momentum flux through 1PN; the physics looks right, but two printed source equations carry a c^2/γ typo that should be fixed before this goes out.","tokens_in":34816,"tokens_out":2962,"would_cite":true,"duration_ms":29218,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C25","83C35","83C57"],"pacs":["04.30.-w","04.25.Nx","04.50.Kd"],"model":"deepseek-v4-flash","headline":"For noncircular black-hole binaries in Einstein–Maxwell–dilaton theory, the angular momentum flux is now known through first post-Newtonian order, with scalar and electromagnetic dipole radiation appearing a full order before tensor quadrup","keywords":["Einstein-Maxwell-dilaton theory","angular momentum flux","post-Newtonian approximation","noncircular orbits","black-hole binaries","scalar dipole radiation","electromagnetic dipole radiation","gravitational waves"],"falsifier":"Compute the orbit-averaged angular momentum loss of two eccentric EMd black holes using a fully nonlinear numerical relativity simulation and compare with Eqs. (V.13)–(V.25); if the scalar/electromagnetic dipole channels do not appear at O(c^-3), or if α1 = α2 does not kill the scalar dipole flux through 1PN order, the central claim fails.","tokens_in":33740,"feed_emoji":"🌀","tokens_out":3251,"duration_ms":36504,"temperature":0.7,"pith_summary":"This paper computes the instantaneous angular momentum loss from a nonspinning, eccentric black-hole binary in Einstein–Maxwell–dilaton theory, through relative first post-Newtonian order. It separately gives the scalar, electromagnetic, and tensor gravitational contributions, showing that the scalar and electromagnetic channels start with dipole radiation at O(c^-3), while the tensor channel starts at quadrupole order at O(c^-5). These results are needed because, for noncircular orbits, the energy flux alone cannot determine the orbital evolution; two balance equations require both energy and angular momentum flux. The paper checks that its flux satisfies the quasicircular balance relation, reduces to the known general-relativistic result when scalar and electric charges vanish, and vanishes appropriately when scalar or electromagnetic dipole radiation is suppressed. A sympathetic reader would care because this provides the dissipative input for modeling eccentric inspirals and gravitational-wave phasing in a theory with scalar and vector radiation.","feed_headline":"Eccentric EMd binaries: angular-momentum loss now at 1PN","feed_subtitle":"Scalar and electromagnetic dipole radiation appear a full order before the tensor quadrupole channel.","key_machinery":"The calculation rests on the skeletonized point-particle action with a scalar-dependent mass expanded as m(φ) = m[1 + α δφ + ½(α²+β)δφ² + O(c^-6)], where α and β are body-dependent scalar-charge parameters. From this, the paper constructs the relaxed field equations for gravity, scalar, and electromagnetic fields, evaluates the scalar, electric, magnetic, and Epstein–Wagoner gravitational source multipole moments in the center-of-mass frame, and expands the wave-zone fields in powers of 1/R, keeping both the radiative 1/R terms and the finite-distance 1/R² terms. The angular momentum flux is then obtained by inserting the radiative fields into the scalar, electromagnetic, and Landau–Lifshitz","core_discovery":"The paper's central result is the instantaneous angular momentum flux for generic noncircular orbits up to relative 1PN order, expressed separately for the scalar, electromagnetic, and tensor channels. The scalar dipole flux scales as O(c^-3) and is proportional to (α1−α2), the electromagnetic dipole flux also scales as O(c^-3) and is proportional to (q1/m1 − q2/m2), and the tensor quadrupole flux scales as O(c^-5); all contributions are proportional to n×v, so radiation changes the magnitude of the orbital angular momentum without tilting the orbital plane at this order. The paper demonstrates consistency with the quasicircular balance relation ΩJ_z = F, with the generic-orbit GR limit, and","pith_inferences":["If the scalar charges α and β are not truly constant during an eccentric inspiral—because the scalar field varies along the orbit or because the black-hole charges are not conserved—the flux formulas would need additional terms beyond those computed here.","The retained 1/R² finite-distance waveform terms could be tested against numerical-relativity waveforms extracted at finite radius, where such corrections matter before taking the null-infinity limit.","For eccentric binaries, the dipole channels at O(c^-3) may drive a qualitatively different eccentricity evolution than in GR; averaging these fluxes over an orbit is a natural next step that the paper leaves to future work.","The α and β parameters are taken from an earlier matching calculation for EMd black holes; a direct comparison of the predicted angular-momentum loss with a fully relativistic simulation of two charged, scalar-hairy black holes would test whether the point-particle truncation is sufficient at 1PN order."],"forward_implications":["An orbit-averaged eccentric binary in EMd theory can be evolved using the balance equation ⟨dL/dt⟩ = −⟨J⟩, with J now known to relative 1PN order.","Scalar and electromagnetic dipole radiation dominate angular momentum loss at a lower PN order than the gravitational quadrupole channel, making them potentially important for eccentric inspiral rates.","The quasicircular limit of the flux satisfies ΩJ_z = F, so the new result is consistent with the previously known energy flux for circular orbits.","In the GR limit α = β = q = 0, the scalar and electromagnetic contributions vanish and the tensor flux reduces to the known generic-orbit GR result.","The dipole-suppression conditions α1 = α2 and q1/m1 = q2/m2 make the respective dipole contributions vanish through the order considered, confirming the expected structure."],"fun_headline_variants":["1PN angular momentum flux for noncircular EMd binaries","Dipole radiation leads EMd angular momentum loss at 1PN","Eccentric EMd binaries: angular momentum flux to 1PN","Scalar and EM dipoles precede tensor quadrupole in EMd flux"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper assumes each black hole is described by a point-particle action with a scalar-dependent mass expanded to second order and uses constant scalar-charge parameters α and β taken from an earlier matching calculation, without re-deriving them or checking whether higher-order sensitivity terms are negligible for strong-field EMd black holes; if these parameters vary during inspiral or the truncation misses important terms, the flux needs modification.","fun_headline_variants_meta":{"raw":{"variants":["1PN angular momentum flux for noncircular EMd binaries","Dipole radiation leads EMd angular momentum loss at 1PN","Eccentric EMd binaries: angular momentum flux to 1PN","Scalar and EM dipoles precede tensor quadrupole in EMd flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000251,"raw_usage":{"total_tokens":1410,"prompt_tokens":773,"completion_tokens":637,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":560}},"tokens_in":517,"tokens_out":637,"duration_ms":6956,"temperature":1.0,"reasoning_tokens":560,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T16:42:23.233620+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the orbit-averaged angular momentum loss of two eccentric EMd black holes using a fully nonlinear numerical relativity simulation and compare with Eqs. (V.13)–(V.25); if the scalar/electromagnetic dipole channels do not appear at O(c^-3), or if α1 = α2 does not kill the scalar dipole flux through 1PN order, the central claim fails.","supporting_citations":[],"review_version":1}