{"id":"154977df-5365-4793-983b-fc9cb02667ab","arxiv_id":"2608.06487","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The claimed -1PN correction to the binary equations of motion is an artifact of an invalid expansion and never becomes larger than the Newtonian term, per the paper's own Eq. (97).","lead":"This paper studies gravitational waves from two orbiting heavy objects in a universe with small extra curled-up dimensions, claiming a new correction term that would depend strongly on the extra dimensions' size. The paper's own equations show that correction never exceeds the Newtonian force, so the central observational claim does not survive.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed -1PN term is an artifact of expanding Eq. (93) with the series (94), whose expansion parameter gamma_p/gamma = m_d r exp(m_d r) is exponentially large in the observational regime; Eq. (97) shows the real correction is bounded by 0.74 times Newton.","rationale":"The load-bearing step is the derivation of Eq. (95)'s -1PN term. The reader's weakest assumption pinpoints exactly this: the series inversion (94) is used where its expansion parameter is exponentially large. I agree. This is an internal inconsistency, not merely a disagreement with standard PN: Eq. (95) contains an uncancelled -gamma_p/gamma term, whereas the paper's own Eq. (97) has exp(-m_d r)(1 + m_d r) <= 0.74 relative to Newton. Both cannot be true for the same acceleration. The exact Lambert-W inversion confirms the series (94) has broken down. Because the abstract and Section X use the -1PN term as the main observable signature, the central claim fails. Secondary issues noted in the paper (the factor-of-2 discrepancy with [21], the block-diagonal metric caveat, and the apparent non-reduction of Eq. (98) to the standard 1PN flux in the gamma_p -> 0 limit) reinforce the rejection but are not needed for it. No adjustment to the reader's REJECT verdict is required.","tokens_in":16435,"tokens_out":7734,"duration_ms":63905,"concrete_test":"Evaluate the exact Eq. (97) at a representative physical binary (e.g., m_total = 2.8 solar masses, r = 10^8 cm, R = 58 micrometers) and compare its leading coefficient, exp(-m_d r)(1 + m_d r), with the -gamma_p/gamma term printed in Eq. (95). Alternatively, invert Eq. (93) exactly with r = W(gamma_p/gamma)/m_d and substitute into Eq. (97); if no 1/gamma term survives, the claimed -1PN effect is not present. Either check requires no new physics and can be done analytically in a few lines.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section VII defines gamma = exp(-m_d r) Gm/(c^2 r) (Eq. 93) and then inverts it with the series (94). The expansion variable is gamma_p/gamma = m_d r exp(m_d r): Eq. (94) is obtained by expanding exp(-m_d r) in powers of m_d r and replacing r by its leading value, so it is valid only for m_d r << 1. The paper's own observational target is R <= 58 micrometers (Section X) and binary separations r of kilometers or more, so m_d r is enormous and gamma_p/gamma is exponentially large, far outside the radius of convergence. Equation (95) nevertheless keeps the uncancelled -gamma_p/gamma term and labels it -1PN. This contradicts the paper's exact leading-order acceleration, Eq. (97): relative to Newton it is exp(-m_d r)(1 + m_d r), whose maximum is 2/e ~ 0.74; it never exceeds Newton, let alone dominates as an inverse power of gamma. The exact inversion of Eq. (93) is r = W(gamma_p/gamma)/m_d; using it in Eq. (97) gives a correction that decays as (gamma/gamma_p) log(gamma_p/gamma), not one proportional to gamma_p/gamma. Hence the -1PN term in Eq. (95), and the gamma^{-2}, gamma^{-1} terms in the flux Eq. (98), are artifacts of an invalid expansion. The central claim of the paper rests on this term, so the stated result fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers gravitational wave emission from a compact binary in a d+4-dimensional spacetime with d extra dimensions compactified on a hypersphere. Using the multipolar post-Minkowskian formalism, the author solves the homogeneous and inhomogeneous wave equations, derives metric potentials, equations of motion, and energy flux at 1PN order, and claims the appearance of a new term at -1PN order in the equations of motion, with corresponding modifications to the energy flux. The central physical claim is that this -1PN term would dominate the waveform and strongly constrain the size and number of extra dimensions.","tokens_in":16648,"tokens_out":3653,"duration_ms":33116,"significance":"If the central claim were correct, it would be a dramatic and observationally consequential result: a -1PN term dominating over the Newtonian acceleration would have been seen in existing gravitational-wave data, providing sharp bounds on compactified extra dimensions. The paper also attempts a comparison with prior work on Kaluza-Klein reductions of binary systems. However, the central claim fails against the manuscript's own leading-order acceleration: the claimed -1PN term is an artifact of expanding an equation in a parameter that is exponentially large in the observational regime. The paper does contain explicit formal manipulations and a clear presentation of the multipolar framework, and it correctly identifies the decaying pseudo-massive homogeneous modes, but these strengths do not rescue the main conclusion.","major_comments":[{"comment":"The inversion of Eq. (93) via the series in Eq. (94) is the load-bearing step, and it is invalid in the regime relevant to observations. Equation (94) is a power-series expansion in the variable gamma_p/gamma = m_d r exp(m_d r), which converges only for m_d r << 1. In the paper's own observational regime (Section X: R <= 58 micrometers, binary separations r of kilometers or more), m_d r is enormous and gamma_p/gamma is exponentially large. Equation (95) nevertheless keeps the uncancelled -gamma_p/gamma term and identifies it as a -1PN contribution. That term is an artifact of expanding beyond the radius of convergence, not a physical acceleration term.","section":"Section VII, Eqs. (93)-(95)"},{"comment":"The manuscript's own leading-order acceleration, Eq. (97), is exp(-m_d r)(1 + m_d r) Gm/r^2. Relative to the Newtonian acceleration Gm/r^2, the correction factor is 1 + m_d r exp(-m_d r), and the non-Newtonian part is at most e^{-1} Gm/r^2. The total ratio exp(-m_d r)(1 + m_d r) never exceeds 1. This directly contradicts the claim in Section X that the new term 'would be much larger than the Newtonian term.' The internal inconsistency confirms that the conclusion is not supported by the paper's own equations.","section":"Section VII, Eq. (97)"},{"comment":"The flux formula (98) contains terms proportional to gamma^{-2} and gamma^{-1}, and the text states that these originate from the inverse power of gamma in the equations of motion. Since that inverse power is the artifact identified in Eq. (95), the claimed modification of the energy flux is likewise unsupported. Using the exact inversion r = W(gamma_p/gamma)/m_d in the exact leading-order acceleration gives a correction that decays as (gamma/gamma_p) log(gamma_p/gamma), not one proportional to gamma_p/gamma, so no inverse-power enhancement appears.","section":"Section VIII, Eq. (98)"}],"minor_comments":[{"comment":"There are several typographical and formatting issues, including 'P AST ST A TIONARY' in the Section II heading, 'O(c)-4' in Eq. (91) instead of O(c^{-4}), 'igmatchi' in the potentials, and 'gven' in Appendix B; these should be corrected.","section":"Section II and throughout"},{"comment":"Equation (B1) writes the exponential as exp(-m_d r/R), which is dimensionally inconsistent because m_d already has dimensions of inverse length; presumably exp(-m_d r) is intended.","section":"Appendix B, Eq. (B1)"},{"comment":"The notation in Table II labels entries as 'm/R' where m is the total mass; since gamma_p = G m m_d/c^2 is dimensionless, the table entries should be described explicitly as G m/(R c^2) times a numerical factor to avoid confusion.","section":"Table II"},{"comment":"The comparison with Ref. [21] is left unresolved: the author states that a factor of 2 in the Newtonian potential cannot be reproduced. This unresolved discrepancy should be acknowledged more prominently, since it weakens the claimed validation against prior work.","section":"Section IX"}],"recommendation":"reject","confidential_remarks":"The central claim is not defensible because the -1PN term arises from an invalid series inversion, and the paper's own exact leading-order acceleration contradicts the claimed dominance. This is not a local fixable issue; the manuscript would need a fundamentally different treatment of the large m_d r regime. The unresolved discrepancy with Ref. [21] further weakens confidence in the derivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my read of arXiv:2608.06487. The headline is that the central claim, a -1PN term in the equations of motion and flux, does not hold up. It comes from expanding r in powers of m_d r exp(m_d r), which is exponentially large in the regime they care about. The paper's own Eq. (97) gives acceleration exp(-m_d r)(1 + m_d r) Gm/r^2; relative to Newton that's at most 2/e, so it never dominates. The -1PN term is an artifact. The same problem infects the flux: Eq. (98) has inverse-gamma terms that don't reduce cleanly to the standard 1PN flux when gamma_p goes to zero, and the tail term is missing.\n\nThat said, the paper is not junk. The multipolar treatment in compactified extra dimensions is real work: solving the wave equation, getting exponentially decaying KK modes, deriving the Yukawa-like Newtonian potential, and comparing with Du et al. [21]. The author flags the block-diagonal metric caveat and the unexplained factor-2 discrepancy honestly. Those are useful signposts for anyone working in this area.\n\nBut the load-bearing result is the -1PN term, and it's wrong. Section X's claim that it would already have been seen in 400 GW events is therefore unsupported. A corrected treatment using the Lambert W inversion would show a tiny, decaying correction, not an inverse power of gamma. The flux needs rechecking too. The author's self-identified assumptions (block diagonal, no mixing) only add to the uncertainty.\n\nWho is this for? Someone wanting a catalog of technical difficulties in extra-dimensional PN calculations might get something from the intermediate formulas and the comparison with [21]. But as a research result, the stated finding fails. I would not cite the -1PN claim. If the journal sends it out, a referee should be able to catch the expansion problem quickly, so I don't object to referee time, but my recommendation is reject or major revision with the central claim rederived. Not publishable as is.","headline":"The claimed -1PN term is an artifact of an invalid series inversion; the paper's own exact acceleration shows the correction is bounded by about 0.74 times Newton, so the central result does not survive.","tokens_in":17305,"tokens_out":1775,"would_cite":false,"duration_ms":17075,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83E15"],"pacs":["04.30.-w","04.50.-h","04.25.Nx"],"model":"deepseek-v4-flash","headline":"Small extra dimensions would add a -1PN term to the gravitational-wave flux, a correction that would already show up in observed binary events.","keywords":["gravitational waves","extra dimensions","compactified dimensions","post-Newtonian expansion","compact binary inspiral","energy flux","higher-dimensional gravity","Kaluza-Klein modes"],"falsifier":"Compute the center-of-mass acceleration directly from Eq. (B1) without the series inversion (94), using physical values such as $R = 10\\,\\mu\\mathrm{m}$ and $r \\sim 10^3\\,\\mathrm{km}$; if the exact expression contains no term inversely proportional to $\\gamma$ (or if the inversion series diverges), the claimed $-1$PN term is an artifact of the expansion.","tokens_in":16063,"feed_emoji":"🌌","tokens_out":10204,"duration_ms":82272,"temperature":0.7,"pith_summary":"This paper extends the standard post-Newtonian treatment of a compact binary to a spacetime with $d$ extra spatial dimensions curled into a small hypersphere. The central claim is that the gravitational-wave emission contains an infinite tower of exponentially decaying 'pseudo-massive' modes on top of the usual $1/r$ radiation, and that these modes modify the inspiral equations of motion at minus-first post-Newtonian order ($-1$PN). That modification enters the energy flux through terms proportional to $\\gamma^{-2}$ and $\\gamma^{-1}$, where $\\gamma = Gm/(c^2 r)$ is the usual post-Newtonian parameter. The author argues that the new term would be larger than the Newtonian term for small extra dimensions, so it would already have been seen in the gravitational-wave events detected so far; the result is therefore presented as a way to bound the size and number of extra dimensions. A sympathetic reading takes the paper's goal to be a concrete, testable waveform correction that data can confirm or exclude.","feed_headline":"Small extra dimensions would add a -1PN term to GW flux","feed_subtitle":"A new -1PN term would dwarf Newtonian gravity and show up in the events already observed.","key_machinery":"Central to the construction is the reduction of the $(4+d)$-dimensional d'Alembertian on the compact sphere to a Klein-Gordon operator with pseudo-mass $m_d = \\sqrt{l_d(l_d+d-1)}/R$, using the eigenvalues of the $d$-dimensional hyperspherical Laplacian. This makes the Green's function a sum of Yukawa kernels $\\exp(-m_d|\\vec x-\\vec x'|)/|\\vec x-\\vec x'|$ rather than a single $1/r$ kernel. All subsequent effects flow from one identity: the gradient of the Yukawa potential with respect to the source position produces the factor $\\exp(-m_d r)(1+m_d r)$ in the Newtonian acceleration. Re-expressing that factor in terms of $\\gamma = (Gm/c^2 r)\\exp(-m_d r)$ yields the series inversion whose leading correction produces the claimed $-1$PN term in the equations of motion and the $\\gamma^{-2}$, $\\gamma^{-1}$ terms in the flux.","core_discovery":"The paper claims that in a $4+d$ dimensional spacetime with the $d$ extra dimensions compactified on a hypersphere of radius $R$, the homogeneous wave equation separates into a four-dimensional Klein-Gordon equation whose pseudo-mass is $m_d = \\sqrt{l_d(l_d+d-1)}/R$, with $l_d$ labeling hyperspherical harmonics. Solving both the homogeneous and sourced wave equations, the metric potentials acquire sums of Yukawa-type terms $\\exp(-m_d r)/r$. When the 1PN acceleration is computed for a circular orbit in the center-of-mass frame, a term of order $-1$PN appears, proportional to $\\gamma_p/\\gamma$ with $\\gamma_p = G m m_d / c^2$; the same inverse powers of $\\gamma$ enter the energy flux. The author traces this term to the gradient of the leading-order potential, which contains the factor $\\exp(-m_d r)(1+m_d r)$, and notes that unless $\\gamma_p$ is very small the term would dominate the waveform and already be visible in the observed events. The paper therefore concludes that observations either exclude compactified extra dimensions, force them to be larger than expected, or require many dimensions whose combined $\\gamma_p$ is small.","pith_inferences":["One immediate check is to compute the circular-orbit acceleration directly from the exact Yukawa factor without the series inversion (94); the ratio of the correction term to the Newtonian term is then $m_d r \\exp(-m_d r)$, which never exceeds $e^{-1}$, so the $-1$PN term may be an artifact of expanding outside the series' radius of convergence.","Because the inversion parameter is $\\gamma_p/\\gamma = m_d r \\exp(m_d r)$, realistic binaries with $R \\lesssim 58\\,\\mu\\mathrm{m}$ and kilometer separations make that parameter astronomically large; checking whether Eq. (95) survives a direct evaluation would settle the matter.","The unresolved factor-of-2 discrepancy with the earlier 5D treatment suggests that the block-diagonal metric ansatz, rather than the physics of extra dimensions, may be responsible for the new terms; repeating the calculation with a general metric could remove or shift the $-1$PN effect.","A natural extension is to compute the 2PN flux: if the $-1$PN term is real, the higher-order terms will contain even stronger inverse powers of $\\gamma$, making the disagreement with standard templates grow rather than shrink."],"forward_implications":["If the calculation is correct, current ground-based detections already constrain the compactification radius and the number of extra dimensions, because the $-1$PN term would have shown up in roughly 400 observed events.","The flux formula (98) provides coefficients in $\\gamma^{-2}$ and $\\gamma^{-1}$ that data analysis can compare with standard templates to place limits on $\\gamma_p$.","For one extra dimension the Newtonian potential becomes the harmonic series $(Gm/r)\\sum_{l_d}\\exp(-l_d r/R)$; summing it reproduces the $1/r$ behavior only in the limit of zero radius, and comparison with an earlier 5D calculation leaves an unexplained factor of 2.","The paper's stated conclusion is a three-way alternative: no extra dimensions, larger dimensions than current bounds allow, or a sufficiently large number of extra dimensions to make $\\gamma_p$ small.","The author identifies the block-diagonal metric assumption as the leading caveat, noting that mixing between 4D and extra-dimensional components at higher order would invalidate some results."],"supporting_citations":[{"why":"supplies the multipolar post-Minkowskian formalism and the standard 1PN metric, acceleration, and flux that the paper generalizes.","marker":"[8]"},{"why":"provides the wave-equation solution and the derivation of potentials and equations of motion from which the extra-dimensional terms are computed.","marker":"[17]"},{"why":"gives the inverse-Laplacian solutions and the regularization function used in the 1PN potentials.","marker":"[12]"},{"why":"supplies the eigenvalues and eigenfunctions of the Laplacian on the $d$-dimensional sphere used to separate extra dimensions.","marker":"[13]"},{"why":"is the earlier 5D binary treatment with Kaluza-Klein reduction whose results are compared and found to differ by a factor of 2.","marker":"[21]"},{"why":"provides the GW170817/GRB170817A observational bound that sets the maximum size of extra dimensions used in the conclusions.","marker":"[7]"},{"why":"gives the torsion-balance laboratory limit on the size of extra dimensions, used to argue the -1PN term would already have been seen.","marker":"[28]"}],"fun_headline_variants":["Extra dimensions yield -1PN term in gravitational waves","New -1PN ordering from compact extra dimensions in GWs","GW equations gain -1PN term if extra dimensions exist","Compact extra dimensions force -1PN GW flux correction","Higher-dim GWs: -1PN term could be visible in events"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the quantity $\\gamma_p/\\gamma = m_d r \\exp(m_d r)$ can be treated as a small expansion parameter; for observed binaries the compactification radius is at most tens of micrometers while the orbital separation is kilometers or more, so this parameter is astronomically large and the expansion that produces the $-1$PN term is being used where it does not converge.","fun_headline_variants_meta":{"raw":{"variants":["Extra dimensions yield -1PN term in gravitational waves","New -1PN ordering from compact extra dimensions in GWs","GW equations gain -1PN term if extra dimensions exist","Compact extra dimensions force -1PN GW flux correction","Higher-dim GWs: -1PN term could be visible in events"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000705,"raw_usage":{"total_tokens":3147,"prompt_tokens":880,"completion_tokens":2267,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":2184}},"tokens_in":496,"tokens_out":2267,"duration_ms":13531,"temperature":1.0,"reasoning_tokens":2184,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:35:17.644668+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the center-of-mass acceleration directly from Eq. (B1) without the series inversion (94), using physical values such as $R = 10\\,\\mu\\mathrm{m}$ and $r \\sim 10^3\\,\\mathrm{km}$; if the exact expression contains no term inversely proportional to $\\gamma$ (or if the inversion series diverges), the claimed $-1$PN term is an artifact of the expansion.","supporting_citations":[{"cited_title":"Gravitational waves from a binary source in higher dimensional spacetime with compactified extra dimensions","cited_arxiv_id":"2608.06487","evidence_quote":"supplies the multipolar post-Minkowskian formalism and the standard 1PN metric, acceleration, and flux that the paper generalizes."},{"cited_title":"Higuchi, Journal of Mathematical Physics43, 6385 (2002)","cited_arxiv_id":null,"evidence_quote":"provides the wave-equation solution and the derivation of potentials and equations of motion from which the extra-dimensional terms are computed."},{"cited_title":"Blanchet, Living Rev","cited_arxiv_id":null,"evidence_quote":"gives the inverse-Laplacian solutions and the regularization function used in the 1PN potentials."},{"cited_title":"Higuchi, Journal of Mathematical Physics28, 1553 (1987)","cited_arxiv_id":null,"evidence_quote":"supplies the eigenvalues and eigenfunctions of the Laplacian on the $d$-dimensional sphere used to separate extra dimensions."},{"cited_title":"The treatment performed here is somewhat different from that in [21] and their results are also different","cited_arxiv_id":null,"evidence_quote":"is the earlier 5D binary treatment with Kaluza-Klein reduction whose results are compared and found to differ by a factor of 2."}],"review_version":1}