{"id":"9e92a05b-d714-457e-aa2e-f0c947d452fc","arxiv_id":"2501.13210","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"DHOST inflation with shift-symmetry breaking in the F2 term yields a nearly scale-invariant tensor spectrum and a constant Ω_GW h^2 ≈ 5.2×10^-17, strongly violating the single-field consistency relation.","lead":"Cosmologists show that a class of modified-gravity inflation models (DHOST) can produce gravitational waves with a nearly flat spectrum, in strong violation of the standard single-field inflation consistency relation. The result makes a testable prediction: a constant relic gravitational-wave background in the mHz range that future space interferometers such as μARES could distinguish from ordinary inflation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tensor flatness and perturbativity are not verified at mHz scales: Fig. 1 covers only k = 10^-4–10^-1 Mpc^-1, while mHz GWs correspond to k ~ 10^11–10^13 Mpc^-1, where φ(η) at horizon crossing is orders of magnitude larger and the shift-symmetry-breaking corrections grow.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the detection claim in Section 3 depends on the tensor spectrum remaining flat at frequencies many orders of magnitude above the computed range. This is indeed the most critical gap. The paper's internal consistency is not at issue; the model is coherent, and the parameters are tuned to satisfy Planck constraints at the displayed scales. But the mHz prediction is a direct extrapolation, and the field φ evolves by ~10^6 over the additional ~30 e-folds to mHz horizon crossing, so the perturbative coefficients and their derivatives—which set n_T—are not bounded by the analysis. Fine-tuning of λ and m^2 is an aesthetic concern, not a logical flaw; the unexplained end of inflation is acknowledged but does not directly threaten the tensor prediction. The missing high-k verification is a concrete, addressable omission: extending the Mukhanov-Sasaki computation to mHz scales would settle whether the constant Ω_GW signature survives. Since the paper explicitly computes only CMB scales and offers no argument for the far-infrared extrapolation, the conditional verdict is appropriate. We recommend UNCHANGED: the reader's conditional verdict correctly captures that the central claim is plausible but not yet demonstrated.","tokens_in":10114,"tokens_out":17337,"duration_ms":171261,"concrete_test":"Solve the tensor Mukhanov-Sasaki equation for the model parameters (α_H = 1.04, α_B = 1, β_K = 3.97343, c = 1, h_dS = 5e-5, f2 = 4.34, λ = -1.1e-30, m^2 = -1e-21) over k from 0.05 to 10^12 Mpc^-1, using the full time-dependent φ(η) = c + (1/h_dS) ln(-η) + C (with C fixed so that the CMB spectra match Fig. 1) and the explicit A_tensor, B_tensor from App. A. Compute P_T(k) and n_T(k) at each k* = -1/k at horizon crossing; if n_T stays below ~10^-4 and the perturbation corrections stay below ~10^-2 across the band, the flat-spectrum prediction survives. Alternatively, derive an analytic estimate of n_T(k) from the φ-dependence of the perturbed coefficients and check whether the correction exceeds the 4×10^-8 value at k = 10^12 Mpc^-1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—constant Ω_GW h^2 = 5.18×10^-17 in the mHz band—rests on the tensor spectrum remaining nearly scale-invariant from CMB scales (k ~ 0.05 Mpc^-1) to k ~ 10^12 Mpc^-1 (f ~ mHz). The paper computes spectra only in Fig. 1 over k ∈ [10^-4, 10^-1] Mpc^-1 and does not evaluate the model at the much larger k. The shift-symmetry breaking enters through φ^2 and φ^4 terms (Eq. 4, App. A7–A8), where φ(η) = c - t(η) with t = -(1/h_dS) ln(-η) + const, so at horizon crossing (η* = -1/k), φ(k) = const - (1/h_dS) ln k. With h_dS = 5×10^-5, |φ| changes by ≈ 2×10^4 per e-fold. From k = 0.05 to k = 6×10^11 Mpc^-1, ln(k/k*) ≈ 30, so φ shifts by ≈ 6×10^5. The perturbed coefficients A_tensor_pert, B_tensor_pert scale as φ^2 and φ^4, so their magnitude and, more importantly, their logarithmic derivatives (which set n_T) change substantially; the slow-roll-like treatment used to solve the Mukhanov-Sasaki equation is not justified when φ' = -k/h_dS becomes enormous at late times. The paper neither extends the computation nor bounds n_T at high k; the claim of stability from 'large values of the coefficients' is only argued for the CMB range. Hence the constant-Ω_GW prediction and the μARES detectability claim are extrapolations, not demonstrated results.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies DHOST inflationary models whose background is a de Sitter phase generated by shift-symmetric operators, with the shift symmetry broken in the F2 sector by small φ^2 and φ^4 terms. The authors derive the scalar and tensor power spectra using a Mukhanov-Sasaki quantization procedure, fit the model parameters to the Planck 2018 scalar spectral observables, and report a tensor-to-scalar ratio r(k*)=0.0147879 with a nearly scale-invariant tensor spectrum (n_T ~ 4×10^-8), thereby strongly violating the single-field consistency relation. They then use the constant tensor spectrum to predict a gravitational-wave density parameter Ω_GW h^2 = 5.18×10^-17 in the mHz band and argue that this is within reach of the proposed μARES interferometer.","tokens_in":10479,"tokens_out":13417,"duration_ms":137207,"significance":"If the mHz extrapolation is valid, the paper provides a concrete counterexample to the single-field consistency relation in a well-defined DHOST framework, with explicit spectra that are compatible with Planck scalar constraints. The computations are substantial, the fiducial model is clearly specified, and the appendix gives the perturbed scalar and tensor coefficients. The main caveat is that the central detectability claim rests on an uncomputed extrapolation of the tensor spectrum from CMB scales to mHz scales, and the extremely small couplings λ ~ 10^-30 and m^2 ~ 10^-21 are tuned without a dynamical explanation. The paper is therefore interesting and publishable in principle, but the mHz prediction needs additional support.","major_comments":[{"comment":"The constant Ω_GW h^2 = 5.18×10^-17 is obtained by taking the tensor power spectrum to remain essentially scale-invariant from CMB scales (k ~ 0.05 Mpc^-1) to mHz frequencies (k ~ 6×10^11 Mpc^-1). However, Fig. 1 only displays spectra over k ∈ [10^-4, 10^-1] Mpc^-1, and the paper does not compute P_T(k) at the high-k, late-time horizon crossings relevant for mHz detectors. Since the shift-symmetry-breaking corrections in Eqs. (A7) and (A8) depend on φ at horizon crossing, which evolves as φ(k) = const - (1/h_dS) ln k, the size of these corrections changes by many orders of magnitude over the extrapolated range. The authors should either extend the numerical computation of P_T and n_T to mHz scales or provide an analytic bound on n_T at high k derived from the explicit expressions in Appendix A. This is load-bearing for the μARES detectability statement in the conclusions.","section":"Section 3, Fig. 2 and Eq. (30)"},{"comment":"The perturbativity argument states that \"Due to the large values of the coefficients K^2_λ and K^2_m2, λ≲O(10^-30) and m^2≲O(10^-18) ensure that the perturbative regime is indeed valid,\" but the coefficients K^2_λ and K^2_m2 are never displayed, and the statement appears to be checked only over the CMB k-range shown in Fig. 1. At mHz horizon crossings, φ is roughly 6×10^5 times larger than at CMB scales for the fiducial h_dS = 5×10^-5, so the same coefficients must be re-evaluated there. The paper should state the range of k for which the perturbative treatment is claimed to hold and verify that the high-k part of Fig. 2 is within that range.","section":"Section 2, after Eq. (29)"}],"minor_comments":[{"comment":"There is a stray closing parenthesis in the definition of α_s(k): \"d log k)\" should read \"d log k.\"","section":"Section 2, Eq. (20)"},{"comment":"The value N* = 59.78 is quoted without a derivation of the end of inflation or a reheating mechanism; the model is presented as a truncated effective theory, and it would be helpful to state the assumptions used to estimate N*.","section":"Section 2, Eq. (10) and surrounding text"},{"comment":"The formulas for c̄_s^2 and r are cited from earlier work, but the notation (e.g., β_H, c̄_s) is introduced only in passing; a short definition or reference to the specific equations in Brax & Lazanu (2021) would improve readability.","section":"Section 2, Eq. (26)-(27)"},{"comment":"The panel labels for α_s and β_s are difficult to read, and the captions abbreviate \"the running of its running\" without naming α_s and β_s in the same order as the panels; please make the caption-to-panel correspondence explicit.","section":"Figure 1"},{"comment":"The statement that the model is \"within the detectable range\" of μARES is stronger than the body text, where the amplitude is described as \"very close to the sensitivity threshold\"; please use consistent wording and quantify the signal-to-noise ratio or sensitivity margin.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for MNRAS and the central idea is interesting. The main technical gap is the unverified extrapolation of the tensor spectrum to mHz scales; this can likely be fixed by a direct computation or an analytic estimate using the Appendix A expressions, so I do not recommend rejection. The tuning of λ and m^2 is acknowledged by the authors and is a common feature of such model-building papers; it should not by itself block publication, but the detectability claim should be stated with the appropriate caveat."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe new thing here is breaking the shift symmetry in F2 (the Ricci coefficient) with tiny phi^2 and phi^4 terms rather than in F0, and showing the resulting tensor spectrum is nearly flat while scalars match Planck. The computation is straightforward but honest, with explicit appendix expressions, and the fiducial model gives r(k*) ~ 0.015 with n_T ~ 4e-8, a genuine strong violation of the consistency relation. The constant Omega_GW h^2 ~ 5e-17 follows from that flatness. The paper also notes its own caveats—parameter tuning unexplained, end of inflation not modeled, only a subset of DHOST theories. The core model-building result is solid.\n\nThe soft spot is the mHz extrapolation. Fig. 1 only spans k = 10^-4 to 10^-1 Mpc^-1. mHz GWs are k ~ 10^11 to 10^13 Mpc^-1, thirty e-folds away. The background phi = c - t shifts by ~5e5 over that range, and the shift-breaking terms are phi^2 and phi^4. The paper does not compute the spectra, n_T, or perturbative stability there. The preprint asserts stability only for the CMB range. Without a high-k computation or at least a bound, the constant Omega_GW and the muARES detectability claim are extrapolations, not demonstrated predictions. That does not kill the paper—the consistency-relation violation remains a CMB-scale result—but it means the headline detection prospect is unproven.\n\nThe stress-test note from the reader is correct on this. I'd also note the tuning: lambda ~ 1e-30 and m^2 ~ 1e-21 are put in by hand, and the paper says so. That is a weakness, but not fatal in effective-field-theory model building.\n\nWho this is for: people working on DHOST inflation and alternative inflationary signatures. It deserves a serious referee—the math is standard, the result is new, and the gap can be addressed with a high-k extension. I'd recommend engage, but require the high-k analysis before publication claims about mHz detection.\n\nBest,\n[Your name]","headline":"Solid DHOST inflation extension, but the mHz detection claim is an uncomputed extrapolation.","tokens_in":11069,"tokens_out":2698,"would_cite":true,"duration_ms":27073,"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":"Tiny shift-symmetry-breaking terms in DHOST inflation flatten the tensor spectrum, violating the single-field consistency relation and putting a constant gravitational-wave background within reach of a future mHz detector.","keywords":["DHOST inflation","primordial gravitational waves","tensor power spectrum","consistency relation","shift symmetry breaking","mHz interferometers","inflationary perturbations","Planck constraints"],"falsifier":"Solve the model's Mukhanov-Sasaki equation at the high-$k$ modes corresponding to mHz gravitational waves ($k \\sim 10^6$ to $10^{15}\\,{\\rm Mpc}^{-1}$) and check whether $\\Omega_{\\rm GW} h^2$ remains constant at $5.18 \\times 10^{-17}$; if the tensor spectrum tilts or the small couplings stop being perturbative at those scales, the claimed detectability fails.","tokens_in":9823,"feed_emoji":"🌊","tokens_out":17391,"duration_ms":141975,"temperature":0.7,"pith_summary":"This paper claims that a class of DHOST theories, gravitational theories with a scalar field that avoid ghost instabilities, can inflate the early universe with an exactly de Sitter background while generating small deviations from scale invariance through tiny shift-symmetry-breaking terms added to the coefficient of the Ricci scalar. The result is a tensor power spectrum that is almost perfectly flat, with a tensor spectral index $n_T \\approx 4 \\times 10^{-8}$ and a tensor-to-scalar ratio $r(k_*) = 0.0148$, which strongly violates the single-field consistency relation $r = -8 n_T$. If the claim is right, the gravitational-wave density parameter $\\Omega_{\\rm GW} h^2$ remains constant across frequencies, at $5.18 \\times 10^{-17}$, and a future space-based mHz interferometer could detect it and tell this model apart from slow-roll inflation. The paper also finds that the running of the scalar spectral index is positive, opposite to what most single-field models predict.","feed_headline":"Inflationary gravitational waves stay flat, breaking a key relation","feed_subtitle":"A nearly flat gravitational-wave background would let a mHz interferometer tell DHOST inflation from slow-roll models.","key_machinery":"The central object is the shift-symmetry-breaking perturbation in $F_2(\\phi,X)$, the coefficient of the Ricci scalar, taking the form $F_2(X) \\to F_2(X) - \\left(\\frac{m_{\\rm phys}^2}{2}\\phi^2 + \\frac{\\lambda_{\\rm phys}}{4!}\\phi^4\\right)$. With the background solution $\\phi = c - t$ (so $X = -1$), the functions $f_i$ and their derivatives are constant, and the model's parameters can be expressed through the $\\alpha_i$ and $\\beta_i$ coefficients. The shift-breaking terms modify the coefficients $A$ and $B$ in the second-order action for both scalar and tensor perturbations, and the power spectra are obtained by quantising the Mukhanov-Sasaki variable from a Bunch-Davies vacuum. This machinery is what decouples the background expansion (set by the shift-invariant part of the action) from the generation of perturbations (set by the tiny shift-breaking parameters), which is the root of the consistency-relation violation.","core_discovery":"The central discovery is that breaking the shift symmetry of a DHOST action in the $F_2$ function, by adding small $\\phi^2$ and $\\phi^4$ pieces to the coefficient of the Ricci scalar, leaves the de Sitter background intact while making the scalar and tensor perturbations acquire their observed tilts. For the fiducial parameters the model matches Planck's scalar spectral index, running, and amplitude constraints and gives $r(k_*) = 0.0147879$ with $n_T \\approx 4 \\times 10^{-8}$, so that $r \\gg 8|n_T|$, a strong violation of the consistency relation of single-field inflation. Because the tensor spectrum is so flat, the gravitational-wave density parameter is constant, $\\Omega_{\\rm GW} h^2 = 5.18 \\times 10^{-17}$, and the paper argues this lies within the detectable range of the proposed $\\mu$ARES space-based interferometer in the mHz band.","pith_inferences":["If the flat tensor spectrum persists at the much higher wavenumbers ($k \\sim 10^6$ to $10^{15}\\,{\\rm Mpc}^{-1}$) that source mHz gravitational waves, the model predicts a nearly white stochastic background in $\\Omega_{\\rm GW}$; but the paper only computes the spectrum on CMB scales, so this is an untested extrapolation.","The same mechanism of decoupling the background from perturbation generation could be realised in other shift-symmetric scalar-tensor theories, making a flat tensor spectrum a more general signature than this particular DHOST construction; the paper does not explore that generality.","A future detection of a positive scalar spectral-index running would be a strong discriminator against single-field slow-roll, but current Planck uncertainties are too large to distinguish the sign; high-resolution CMB experiments could settle it.","The extreme smallness of the symmetry-breaking couplings resembles the cosmological-constant fine-tuning problem, and explaining it might require a symmetry or anthropic argument; the paper leaves this open."],"forward_implications":["A future measurement that finds $r \\approx -8 n_T$ would contradict this model; the paper predicts $r(k_*) = 0.0148$ with $n_T \\approx 4 \\times 10^{-8}$.","The gravitational-wave density parameter is predicted to be frequency-independent at $\\Omega_{\\rm GW} h^2 \\approx 5.2 \\times 10^{-17}$ in the mHz band, in contrast to slow-roll models where it decays with frequency, so the shape of the background alone can discriminate the two.","The model predicts a positive running of the scalar spectral index, $\\alpha_s \\approx 0.0011$, opposite to the negative running expected in most single-field slow-roll models, testable with future CMB data.","The tiny shift-symmetry-breaking couplings ($\\lambda \\approx -1.1 \\times 10^{-30}$, $m^2 \\approx -1 \\times 10^{-21}$) keep the perturbative treatment valid at the scales where the paper checks it, but their smallness is not explained.","The model yields $N_* \\approx 59.8$ e-foldings between pivot-scale horizon exit and the end of inflation, assuming inflation ends by a sudden change of physics."],"supporting_citations":[{"why":"Supplies the DHOST perturbation formalism and the $\\alpha_i$/$\\beta_i$ parametrisation that this paper extends to shift-symmetry breaking in $F_2$.","marker":"Brax & Lazanu 2021"},{"why":"Provides the Planck 2018 measurements of $n_s$, $\\alpha_s$, $\\beta_s$, and $A_s$ that the fiducial model is tuned to satisfy.","marker":"Akrami et al. 2020"},{"why":"Sets the current upper bound $r < 0.032$ (reduced to $0.028$) that the model's $r(k_*) = 0.0148$ must beat.","marker":"Tristram et al. 2021, 2022"},{"why":"Tightens the tensor-to-scalar ratio constraint by adding LIGO-Virgo-KAGRA data; the model remains compatible.","marker":"Galloni et al. 2023"},{"why":"Gives the formula converting the tensor power spectrum into the gravitational-wave density parameter $\\Omega_{\\rm GW} h^2$.","marker":"Clarke et al. 2020"},{"why":"Provides the gravitational-wave sensitivity curves used to assess whether the predicted $\\Omega_{\\rm GW}$ is detectable.","marker":"Roshan & White 2025"},{"why":"Proposes the $\\mu$ARES interferometer whose sensitivity the detectability claim is compared against.","marker":"Sesana et al. 2021"},{"why":"Supplies the no-ghost and gravitational-wave constraints that reduce the general DHOST action to the simplified form used.","marker":"Crisostomi & Koyama 2018"},{"why":"Underlies the field quantisation procedure used to derive the scalar and tensor power spectra.","marker":"Gorji et al. 2021"}],"fun_headline_variants":["Flat gravitational waves break inflationary consistency rule","DHOST inflation flattens tensor spectrum, defying standard link","Nearly flat gravitational waves herald mHz detections","Shift-symmetry breaking yields flat tensor spectrum in inflation","Inflation model with flat gravitational waves poised for mHz tests"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The detection prospect rests on the tensor spectrum staying almost perfectly flat at much higher frequencies than the paper actually computes, with the tiny symmetry-breaking couplings remaining perturbative even as the background field grows with time.","fun_headline_variants_meta":{"raw":{"variants":["Flat gravitational waves break inflationary consistency rule","DHOST inflation flattens tensor spectrum, defying standard link","Nearly flat gravitational waves herald mHz detections","Shift-symmetry breaking yields flat tensor spectrum in inflation","Inflation model with flat gravitational waves poised for mHz tests"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000155,"raw_usage":{"total_tokens":1143,"prompt_tokens":804,"completion_tokens":339,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":420,"completion_tokens_details":{"reasoning_tokens":270}},"tokens_in":420,"tokens_out":339,"duration_ms":4209,"temperature":1.0,"reasoning_tokens":270,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:22:24.225692+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the model's Mukhanov-Sasaki equation at the high-$k$ modes corresponding to mHz gravitational waves ($k \\sim 10^6$ to $10^{15}\\,{\\rm Mpc}^{-1}$) and check whether $\\Omega_{\\rm GW} h^2$ remains constant at $5.18 \\times 10^{-17}$; if the tensor spectrum tilts or the small couplings stop being perturbative at those scales, the claimed detectability fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the DHOST perturbation formalism and the $\\alpha_i$/$\\beta_i$ parametrisation that this paper extends to shift-symmetry breaking in $F_2$."},{"cited_title":"J., Copeland E","cited_arxiv_id":null,"evidence_quote":"Gives the formula converting the tensor power spectrum into the gravitational-wave density parameter $\\Omega_{\\rm GW} h^2$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the gravitational-wave sensitivity curves used to assess whether the predicted $\\Omega_{\\rm GW}$ is detectable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the no-ghost and gravitational-wave constraints that reduce the general DHOST action to the simplified form used."},{"cited_title":"A., Motohashi H., Mukohyama S., 2021, @doi [JCAP] 10.1088/1475-7516/2021/03/081 , 03, 081","cited_arxiv_id":null,"evidence_quote":"Underlies the field quantisation procedure used to derive the scalar and tensor power spectra."}],"review_version":1}