{"id":"bcf94a70-e912-4a49-b4e5-bcf3056160c8","arxiv_id":"2507.00904","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"The paper applies periodic-orbit and quadrupole-waveform techniques to an asserted higher-derivative Einstein-Æther black hole metric and claims observable phase shifts.","lead":"This paper computes periodic orbits and gravitational waveforms around black holes in a version of Einstein-Æther gravity with added curvature-squared terms. It claims the extra terms produce measurable phase shifts in the gravitational waves, which could be seen by future space-based detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (45)-(46) replace the quadrupole second time derivative with products of coordinate accelerations; for a circular orbit the resulting h_+ scales as R^2 Omega^4 instead of R^2 Omega^2, so the computed waveforms do not follow from the paper's own quadrupole definition.","rationale":"The reader identified the missing derivation of the metric correction (16)-(18) as the weakest assumption; that is a valid concern, since the text goes from the field equations (11) to the metric ansatz with no computation. But the waveform step is more decisive because it is internally inconsistent with the paper's own Eq. (34), so it does not depend on an external derivation. The abstract's claim that higher-derivative corrections produce 'characteristic phase modulations and harmonic deformations' is exactly the content of Eqs. (45)-(46); replacing the quadrupole second derivative with squared accelerations changes the harmonic structure, so no conclusion about parameter sensitivity can be drawn without redoing the computation. I also note the absence of code or parameter values, which would be the only way to check whether the figures used the printed formula; absent that, the correct reading is that the printed derivation supports REJECT. Therefore I keep the reader's verdict unchanged.","tokens_in":8240,"tokens_out":10275,"duration_ms":114841,"concrete_test":"Recompute h_+ and h_x for the (z,w,v) = (1,2,0) orbit using the expressions that actually follow from Eqs. (34)-(36): h_+ = (2 G_ae m/D_L)(xdot^2 - ydot^2 + x xddot - y yddot) and h_x = (2 G_ae m/D_L)(x yddot + 2 xdot ydot + xddot y). Compare against Eqs. (45)-(46) for the same numerical orbit. The circular-orbit limit already distinguishes them: correct h_+ ~ -4 G_ae m R^2 Omega^2 cos(2 Omega t)/D_L, while Eq. (45) gives +G_ae m R^2 Omega^4 cos(2 Omega t)/D_L. If the figures were generated with the correct formula, the paper must state so and the printed formulas corrected; otherwise the waveform results are artifacts of the incorrect equation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing defect is not the metric ansatz (though that is also unsupported) but the transition from the stated quadrupole formula to the waveform equations. Eq. (34) defines h_+ = (G_ae/D_L) Qddot_ij e+_ij, with Q_ij = m(x_i x_j - (1/3) delta_ij r^2). For an equatorial orbit this gives h_+ = (2 G_ae m/D_L)(xdot^2 - ydot^2 + x xddot - y yddot). The paper instead states Eq. (45): h_+ = (G_ae m/D_L)(xddot^2 - yddot^2), and similarly Eq. (46) for h_x. This is not a simplification: xddot^2 is not the second derivative of x^2, and the expression has dimensions of 1/L^2 rather than dimensionless (in c=1 units), while the correct combination xdot^2 + x xddot is dimensionless. The error changes both the amplitude scaling and the spectral content. A circular orbit x = R cos(Omega t), y = R sin(Omega t) gives the correct h_+ proportional to R^2 Omega^2 cos(2 Omega t), whereas Eq. (45) gives R^2 Omega^4 cos(2 Omega t). Since all figures and parameter-sensitivity claims in Sec. V are built on Eqs. (45)-(46), the central GW-signature claim fails internally, independent of whether metric (19) is the true solution.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a static, spherically symmetric black hole metric in Einstein-Æther theory supplemented by quadratic curvature corrections, with the metric function e(r) = 1 - 2M/r + α1/r^2 + α2 log r/r^2 - λ0 (2M/r)^4 (Eq. (19)). It then computes periodic timelike geodesics, classifies them with the (z,w,v) periodic-orbit table, and uses a quadrupole approximation to compute the tensor gravitational waveforms h+ and h× (Eqs. (45)-(46)). The central claim is that higher-derivative and æther corrections imprint phase shifts and harmonic distortions in the waveforms that are sensitive to α, β, γ, and c13, providing a possible observational probe with future detectors such as LISA.","tokens_in":8637,"tokens_out":3701,"duration_ms":45913,"significance":"If the calculation were correct, the paper would extend the periodic-orbit gravitational-wave program to a combined Einstein-Æther plus higher-curvature setting, potentially offering a new test of Lorentz violation and ultraviolet gravity modifications. The paper correctly notes that for c14 = 0 only the two tensor polarizations survive, which is an important consistency check with current observations. It also makes sensible use of the established Levin periodic-orbit classification framework. However, the central metric correction is introduced without derivation, and the gravitational wave formulas are dimensionally inconsistent and mathematically incorrect. As presented, the claims are not supported by the paper's own equations, so the significance is currently prospective rather than demonstrated. The paper does not provide reproducible code, error estimates, or a detectability analysis, which further limits its present contribution.","major_comments":[{"comment":"The waveform formulas do not follow from the stated quadrupole definition. With Q_ij = m(x_i x_j - (1/3)δ_ij r^2), the second time derivative for equatorial motion is Qddot_xx - Qddot_yy = 2m(\\dot x^2 - \\dot y^2 + x\\ddot x - y\\ddot y), so h+ = (2G_ae m/D_L)(\\dot x^2 - \\dot y^2 + x\\ddot x - y\\ddot y). Equation (45) instead states h+ = (G_ae m/D_L)(\\ddot x^2 - \\ddot y^2), which is not the second derivative of x^2 - y^2, has dimensions of 1/L^2 rather than being dimensionless in c=1 units, and changes the amplitude scaling; for a circular orbit x = R cos(Ωt), y = R sin(Ωt) the correct expression is proportional to R^2Ω^2 cos(2Ωt), whereas Eq. (45) gives R^2Ω^4 cos(2Ωt). Since the figures and parameter-sensitivity claims in Sec. V are built on Eqs. (45)-(46), this error invalidates the central waveform claim of the paper.","section":"§V.B, Eqs. (34)-(36) and (45)-(46)"},{"comment":"The higher-derivative correction δe_HD(r) = α1/r^2 + α2 log r/r^2, with coefficients α1 and α2 given in Eqs. (17)-(18), is presented without derivation. The abstract states that the corrected metric is derived, but the manuscript does not show how Eqs. (16)-(18) follow from the field equations (11), nor does it cite a source for the result. This is load-bearing because the entire geodesic and waveform analysis uses metric (19); if this correction is not the actual perturbative solution, all subsequent results collapse. The authors need either to provide the derivation explicitly or to cite a verifiable derivation.","section":"§III, Eqs. (16)-(19)"},{"comment":"The paper claims that the resulting signatures provide a potential observational window for LISA, but it contains no error estimates, no signal-to-noise estimates, and no comparison of the computed waveform amplitudes with plausible detector sensitivities or astrophysical event rates. The adiabatic approximation is asserted without checking whether the radiation-reaction timescale is long compared with the orbital timescale for the parameters used in Figs. 1 and 2. This lack of quantitative support makes the observational claim premature even if the waveform formulas were corrected.","section":"§V, overall"}],"minor_comments":[{"comment":"The field equation (3) has notation that is difficult to parse, with the subscripts on the H terms misaligned; the equation appears to be missing parentheses or indices in several terms, and the overall sign conventions in Eqs. (4)-(6) should be checked against standard references.","section":"§II, Eq. (3)"},{"comment":"The subscript for the æther correction is missing in the text; δe_(r) should be written with a clear identifier, for example δe_ae(r), to avoid confusion with the total metric perturbation.","section":"§III, Eq. (15)"},{"comment":"The polynomial P(x) is said to be sixth-order in x, but the displayed expression includes α2 x^2 log(1/x), which is not a polynomial in x; the text should describe P(x) as a function involving a logarithm rather than as a polynomial.","section":"§IV, Eq. (28)"},{"comment":"The figures are not described with sufficient detail: the axes, the values of the free parameters α, β, γ, c13, M, and L, and the normalization of h+ and h× are not stated, making it impossible for the reader to reproduce or interpret the plotted waveforms.","section":"§V, Figs. 1 and 2"},{"comment":"Reference [35] lists 'T. J. E. Barausse and T. P. Sotiriou', which appears to be an incorrect author name; the reference should be checked and corrected, and the paper would benefit from citing the original Einstein-Æther black hole papers for the c14 = 0 branch used in Eq. (15).","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper applies the well-known periodic-orbit framework to a proposed higher-derivative Einstein-Æther metric, but the central metric ansatz is unproved and the gravitational wave formulas contain an elementary dimensional inconsistency. The latter error alone means that all figures and conclusions in Sec. V are unreliable. I do not see a path to acceptance without a complete reworking of the waveform calculation, an independent derivation or citation for the metric correction, and a quantitative detectability study; this exceeds the scope of a minor or even a typical major revision. I recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a well-intentioned extension of the periodic-orbit program to a combined Einstein-Æther + curvature-squared metric, but the central gravitational-wave calculation is wrong, and the metric itself is asserted rather than derived.\n\nWhat's new: the combination of these two modifications in a periodic-orbit analysis is genuinely new relative to the cited literature. The paper correctly sets up the (z,w,v) classification and the frequency ratio q, and the figures show plausible parameter sensitivity. The references are appropriate, including the recent Lu-Zhu paper that it builds on.\n\nThe soft spots, in order of seriousness. The biggest is eqs. (45)-(46). The paper starts with the standard quadrupole formula h+ = (G_ae/D_L) \\ddot Q_ij e+_ij and Q_ij = m(x_i x_j - (1/3)\\delta_ij r^2). For an equatorial orbit that gives h+ = (2 G_ae m/D_L)(\\dot{x}^2 + x\\ddot{x} - \\dot{y}^2 - y\\ddot{y}). Instead the paper writes h+ = (G_ae m/D_L)(\\ddot{x}^2 - \\ddot{y}^2). That is not a simplification; it is dimensionally wrong (\\ddot{x}^2 has units 1/L^2) and changes the circular-orbit scaling from R^2\\Omega^2 to R^2\\Omega^4. All of Sec. V, including the waveforms and the claimed phase/harmonic signatures, is built on that equation, so the paper's main claim fails internally, independent of the metric question.\n\nThe second problem is the metric. Eq. (19) is introduced as 'the additional correction is given by' with no derivation, although the abstract says the metric is derived perturbatively. The coefficients alpha1 and alpha2 in eqs. (17)-(18) are simply stated. If that ansatz is not the actual solution to the field equations, the orbit and waveform results are moot. I can't verify it from the paper, and the burden is on the authors.\n\nThird, there is no detectability estimate and no sensitivity analysis. The abstract promises 'distinguishable shifts' but nothing connects the parameters to LISA noise curves or event rates. The paper also gives no code and no explicit parameter values for the figures, so the numerics are not reproducible.\n\nThe citation pattern is fine. The paper acknowledges the prior periodic-orbit work and the Einstein-Æther waveform work it draws on.\n\nBottom line: the project has a plausible kernel, but the main computation is wrong and the input metric is unsupported. I would not send this to a referee in its current form.","headline":"The paper applies periodic-orbit techniques to a combined Einstein-Æther plus curvature-squared metric, but a dimensionally wrong quadrupole-to-waveform step invalidates the main results.","tokens_in":9077,"tokens_out":5295,"would_cite":false,"duration_ms":48654,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C10","83C35","83C57","83D05"],"pacs":["04.30.-w","04.50.Kd","04.70.Bw"],"model":"deepseek-v4-flash","headline":"Quadratic curvature corrections to Einstein-Æther black holes imprint parameter-dependent phase modulations and harmonic deformations on the gravitational waves emitted by periodic orbits.","keywords":["higher-derivative Einstein-Æther gravity","periodic orbits","zoom-whirl orbits","gravitational wave waveforms","quadrupole approximation","extreme-mass-ratio inspirals","Lorentz violation","black hole metric corrections"],"falsifier":"Substitute the ansatz $e(r) = 1 - 2M/r + \\alpha_1/r^2 + \\alpha_2 \\log r/r^2 - \\lambda_0(2M/r)^4$ into the full field equations and expand to $O(1/r^4)$; if the residual stress-energy does not vanish at that order, the metric is not the solution and the central claim collapses.","tokens_in":8087,"feed_emoji":"🌊","tokens_out":12765,"duration_ms":125016,"temperature":0.7,"pith_summary":"The paper works in Einstein-Æther gravity supplemented by quadratic curvature terms $R^2$, $R_{\\mu\\nu}R^{\\mu\\nu}$, and $R_{\\mu\\nu\\lambda\\rho}R^{\\mu\\nu\\lambda\\rho}$, solving perturbatively for a static spherically symmetric black hole metric. It then studies periodic timelike geodesics in that spacetime, organized by the rational frequency ratio $q = w + v/z$, and computes the tensor-mode gravitational waveforms $h_+$ and $h_\\times$ in the quadrupole approximation. The central claim is that even small higher-derivative corrections shift the orbital dynamics and imprint characteristic phase modulations and harmonic deformations in the waveform, with the effect controlled by $\\alpha$, $\\beta$, $\\gamma$, and the aether parameter $c_{13}$. If correct, future space-based detectors watching extreme-mass-ratio inspirals could read ultraviolet modifications of gravity and Lorentz violation off the fine structure of the signal.","feed_headline":"Periodic orbit waveforms carry fingerprints of higher-curvature gravity","feed_subtitle":"Small curvature-squared and aether corrections phase-shift the h+ and h× waveforms, giving future detectors a concrete target.","key_machinery":"The central object is the corrected metric function $e(r) = 1 - 2M/r + \\alpha_1/r^2 + \\alpha_2 \\log r/r^2 - \\lambda_0(2M/r)^4$, which enters the effective potential $V_{\\mathrm{eff}}(r) = (1 + L^2/r^2)e(r)$ for equatorial timelike geodesics. This potential turns the radial geodesic equation into a sixth-order polynomial $P(x)$ in $x = 1/r$, and periodic orbits occur when the frequency ratio $q = (1/\\pi)\\int_{x_1}^{x_2} L/\\sqrt{P(x)}\\, dx - 1$ takes a rational value $w + v/z$, encoded in the $(z,w,v)$ triplet. The waveforms are generated from the second time derivatives of the quadrupole moment along these geodesics, using the polarization tensors of the two tensor modes that survive when $c_{14} = 0$.","core_discovery":"On the paper's own terms, the discovery is a worked-out chain from action to waveform: the corrected metric $e(r) = 1 - 2M/r + \\alpha_1/r^2 + \\alpha_2 \\log r/r^2 - \\lambda_0(2M/r)^4$, with $\\alpha_1 = -(16\\pi G_{ae}/3)(3\\alpha+\\beta+2\\gamma)M$, $\\alpha_2 = (16\\pi G_{ae}/3)(\\beta+4\\gamma)M$, and $\\lambda_0 = 27 c_{13}/[256(1-c_{13})]$, changes the effective potential for timelike geodesics, which in turn shifts the orbital frequency ratio $q$ and the rational $(z,w,v)$ classification of periodic orbits. The quadrupole waveforms $h_+$ and $h_\\times$ computed along these geodesics show phase modulations and harmonic deformations relative to general relativity, with corrections proportional to $\\gamma$ singled out as especially visible. The result is presented as a potential observational window into ultraviolet deviations and Lorentz-symmetry breaking in the strong-field regime.","pith_inferences":["If the metric correction is the true solution, the same $\\alpha_1$ and $\\alpha_2$ terms will shift light deflection and perihelion precession at the same post-Newtonian order, giving independent checks on the coefficients the waveform analysis is sensitive to.","The quadrupole calculation covers only a few orbits and ignores radiation reaction, but EMRI signals accumulate over thousands of cycles; a full inspiral waveform would show whether the phase modulation grows secularly or cancels.","A parameter-estimation forecast on synthetic data from future space-based detectors would convert the claimed phase sensitivity into quantitative bounds on $\\alpha$, $\\beta$, $\\gamma$, and $c_{13}$; the paper identifies the observable but does not quantify the reach.","The unusual $\\log r$ term in the metric correction, if physical, may also alter horizon thermodynamics and quasinormal-mode spectra, providing additional signatures of the same underlying corrections."],"forward_implications":["With $c_{14} = 0$ the polarization content is purely tensor, so any deviation from general relativity in the computed $h_+$ and $h_\\times$ is a clean signal for quadratic-curvature and aether parameters.","The frequency ratio $q$ that fixes the $(z,w,v)$ orbit family depends on $\\alpha_1$, $\\alpha_2$, and $\\lambda_0$, so measured orbital harmonics can be mapped back to the coupling constants once a periodic orbit is identified.","Even small values of $\\gamma$ are claimed to produce visible harmonic deformations, making $\\gamma$ the most promising of the curvature couplings for future detectors.","Because the adiabatic approximation holds for extreme-mass-ratio inspirals, the same waveform model applies directly to the EMRI signals that future space-based detectors are designed to observe."],"supporting_citations":[{"why":"Supplies the baseline Einstein-Æther periodic-orbit gravitational radiation calculation that this paper extends by adding quadratic-curvature corrections.","marker":"[1]"},{"why":"Establishes that with c14 = 0 the theory has only the two tensor polarizations, justifying the computation of solely h+ and h×.","marker":"[12]"},{"why":"Provides the (z,w,v) rational periodic-orbit classification used to organize the geodesic families.","marker":"[15]"},{"why":"Motivates the quadratic curvature terms R^2, R_{μν}R^{μν}, and R_{μνλρ}R^{μνλρ} added to the Einstein-Æther action.","marker":"[6]"},{"why":"Defines Einstein-Æther theory and the aether field whose parameter c13 enters the metric correction λ0.","marker":"[9]"},{"why":"Relates extreme-mass-ratio inspiral waves to deviations from general relativity, the observational channel the paper targets.","marker":"[37]"},{"why":"Identifies the future space-based detector that could observe the predicted phase modulations.","marker":"[38]"}],"fun_headline_variants":["Phase-shifted waveforms from aether black hole orbits","Curvature-squared terms leave marks on gravitational wave phases","Periodic orbits near aether black holes alter h+ and h× phases","Small higher-curvature corrections shift gravitational wave spectra","Aether black hole periodic orbits imprint harmonic deformations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole chain of results rests on the form of the higher-curvature metric correction $\\delta e_{\\mathrm{HD}}(r) = \\alpha_1/r^2 + \\alpha_2 \\log r/r^2$: the paper introduces this correction without deriving it from the quadratic-curvature field equations, so if the actual solution differs, the orbital shifts and waveform deformations do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Phase-shifted waveforms from aether black hole orbits","Curvature-squared terms leave marks on gravitational wave phases","Periodic orbits near aether black holes alter h+ and h× phases","Small higher-curvature corrections shift gravitational wave spectra","Aether black hole periodic orbits imprint harmonic deformations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000537,"raw_usage":{"total_tokens":2612,"prompt_tokens":1014,"completion_tokens":1598,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":1516}},"tokens_in":630,"tokens_out":1598,"duration_ms":13595,"temperature":1.0,"reasoning_tokens":1516,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:04:32.911291+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Substitute the ansatz $e(r) = 1 - 2M/r + \\alpha_1/r^2 + \\alpha_2 \\log r/r^2 - \\lambda_0(2M/r)^4$ into the full field equations and expand to $O(1/r^4)$; if the residual stress-energy does not vanish at that order, the metric is not the solution and the central claim collapses.","supporting_citations":[{"cited_title":"Einstein-aether waves,","cited_arxiv_id":null,"evidence_quote":"Establishes that with c14 = 0 the theory has only the two tensor polarizations, justifying the computation of solely h+ and h×."},{"cited_title":"A periodic table for black hole orbits,","cited_arxiv_id":null,"evidence_quote":"Provides the (z,w,v) rational periodic-orbit classification used to organize the geodesic families."},{"cited_title":"Renormalization of higher derivative quantum gravity,","cited_arxiv_id":null,"evidence_quote":"Motivates the quadratic curvature terms R^2, R_{μν}R^{μν}, and R_{μνλρ}R^{μνλρ} added to the Einstein-Æther action."},{"cited_title":"Gravity with a dynamical preferred frame,","cited_arxiv_id":null,"evidence_quote":"Defines Einstein-Æther theory and the aether field whose parameter c13 enters the metric correction λ0."},{"cited_title":"Equivalence principle and gravitational waves from extreme mass-ratio inspirals,","cited_arxiv_id":null,"evidence_quote":"Relates extreme-mass-ratio inspiral waves to deviations from general relativity, the observational channel the paper targets."}],"review_version":1}