Pith. sign in

REVIEW 3 major objections 4 minor 28 references

A minimally modified gravity theory propagates only two tensor modes and becomes cosmologically predictive by coupling its auxiliary constraints to spatial Laplacians of the Lagrange multipliers.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 00:47 UTC pith:YICUW56M

load-bearing objection Neat Laplacian fix for the MMG multiplier problem, but the inflationary predictions are under-determined because the scalar tilt depends on an unspecified tensor-speed parameter. the 3 major comments →

arxiv 2607.26031 v1 pith:YICUW56M submitted 2026-07-28 astro-ph.CO gr-qc

Minimally modified gravity with Laplacian auxiliary constraints and an inflationary realization

Jakkrit Sangtawee , Antonio De Felice , Khamphee Karwan This is my paper
classification astro-ph.CO gr-qc PACS 04.50.Kd98.80.Cq
keywords minimally modified gravityauxiliary constraintsLaplacian couplingtensor degrees of freedominflationtensor-to-scalar ratioconstraint reductioncosmological perturbations
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a modified gravity theory that, around a flat expanding universe, carries only the two tensor gravitational-wave degrees of freedom and a single scalar matter field. The authors first show that a previously proposed four-constraint construction leaves homogeneous Lagrange multipliers undetermined, and these undetermined functions enter the tensor action, so the theory cannot predict gravitational-wave propagation from the background alone. They repair this by coupling the auxiliary constraints to spatial Laplacians of the multipliers, which makes the homogeneous multipliers drop out while still constraining inhomogeneous modes. The result is a predictive quadratic action for tensor, vector, and scalar perturbations, with no gravitational vector or scalar modes on the generic branch. In a lapse-proportional cubic-momentum subclass, the inflationary predictions reduce to a shift in the scalar spectral index and a tensor-to-scalar ratio r = 16 epsilon_s / c_T, where c_T > 1 suppresses r but also lowers the tensor kinetic coefficient.

Core claim

The central claim is that the predictivity obstruction of the original four-constraint MMG construction can be removed by replacing the direct multiplier coupling with a Laplacian coupling: (gamma^{ij} D_i D_j mu_I)(Q^I - P^I). On FLRW, the homogeneous part of mu_I is annihilated by the spatial Laplacian, so it drops out of the background equations, while inhomogeneous multiplier perturbations still enforce Q^I = P^I for k not equal to zero under standard boundary conditions. The resulting quadratic tensor action has coefficients that depend only on the Hamiltonian and background quantities, and the generic branch propagates no vector or scalar gravitational modes. In the lapse-proportional

What carries the argument

The Laplacian auxiliary coupling: the total Hamiltonian density is modified to contain (gamma^{ij} D_i D_j mu_I)(Q^I - P^I) instead of mu_I(Q^I - P^I). Variation gives D^2(Q^I - P^I) = 0, which for non-zero Fourier modes and boundary conditions excluding harmonic functions is equivalent to the original constraints Q^I = P^I, while the homogeneous multiplier part drops out of the background. This mechanism turns an inconsistent homogeneous sector into a closed predictive system. The later subclass with Hamiltonian density proportional to the lapse and containing cubic momentum invariants is what makes the k^4 term vanish and yields the simple c_T formula.

Load-bearing premise

The construction relies on the absence of non-trivial harmonic functions on the FLRW spatial sections: if boundary conditions or the spatial topology allow a homogeneous zero-mode of the Laplacian, the Lagrange multipliers re-enter the background equations and the predictivity obstruction persists.

What would settle it

Compute the quadratic tensor action around a flat FLRW background on a compact three-torus, including the homogeneous k = 0 mode. If the homogeneous Lagrange multipliers or their time derivatives reappear in the tensor coefficients, the claim that the Laplacian coupling removes the predictivity obstruction is false for that topology.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The theory propagates exactly two tensor gravitational degrees of freedom around flat FLRW, with no additional vector or scalar modes on the generic branch.
  • The tensor dispersion relation generically contains both k^2 and k^4 terms; imposing the multimessenger luminality condition selects D_T = 0 and C_T = Q_T.
  • With minimally coupled canonical matter, the only scalar mode is the inflaton fluctuation with unit sound speed.
  • In the cubic-momentum toy model, c_T^2 = 1 + (3/2) f M_P^2 (xi_2 + xi_3) and r = 16 epsilon_s / c_T, so c_T > 1 suppresses the tensor-to-scalar ratio.
  • The observational bounds adopted in the paper require c_T > 1, but large c_T lowers the tensor kinetic coefficient and may lower the perturbative cutoff; establishing the actual strong-coupling scale requires a nonlinear analysis.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same Laplacian-auxiliary trick could be applied to other constraint-based modified gravity constructions where homogeneous multipliers are undetermined, not just the specific four-constraint MMG of this paper.
  • The c_T > 1 requirement during inflation is notable given the tight late-time luminality constraint from GW170817; one could test whether a transition from large c_T during inflation to c_T close to 1 at late times is dynamically possible within the theory.
  • A concrete testable extension is to compute the cubic and quartic actions to determine whether the inflationary predictions are under perturbative control where c_T is large; the quadratic action alone cannot settle this.
  • The zero-mode caveat suggests checking the construction on compact spatial topologies such as a three-torus; if such backgrounds admit harmonic functions, the predictivity claim may need modification.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes a modification of the four-constraint minimally modified gravity (MMG) construction in which the auxiliary constraints are coupled to spatial Laplacians of the Lagrange multipliers. It first shows that in the original construction the homogeneous values of the multipliers are not fully fixed by the FLRW background equations and yet enter the quadratic tensor action, obstructing cosmological predictivity. In the Laplacian-modified version, the homogeneous multiplier sector drops out, while inhomogeneous modes still enforce the intended constraints. The paper derives the quadratic tensor, vector, and scalar actions and finds that, on a generic nondegenerate branch, only two tensor degrees of freedom propagate around flat FLRW. It then couples a canonical scalar field, introduces a lapse-proportional cubic-momentum toy model, and computes the scalar and tensor primordial spectra, concluding that observational bounds require a superluminal tensor speed c_T>1 while noting possible strong-coupling concerns.

Significance. If correct, the Laplacian auxiliary-constraint mechanism provides a genuinely predictive MMG cosmological sector, addressing a concrete obstruction in the prior four-constraint framework. The paper is self-contained, gives explicit perturbative actions, and clearly states the nondegeneracy assumptions. The inflationary realization is a concrete testbed with falsifiable predictions (ns and r), and the authors are transparent about the strong-coupling limitation. The main value lies in the construction, which appears internally consistent on the generic branch, and in the explicit demonstration that the homogeneous multipliers can be removed without destroying the inhomogeneous constraint structure.

major comments (3)
  1. [V.D, Eq. (113), (124), (157)-(159), (162)] The inflationary predictions are not uniquely determined by the parameters (c,A) used in Figs. 1-2. The scalar effective mass Cφ in Eq. (113) depends on the combination ξ2+ξ3, both through the denominator 18ξ1+7ξ2+3ξ3 = -2c/3 + (ξ2+ξ3) and through the final factor 2/M_P^2 + 3f(ξ2+ξ3); this combination is not fixed by c. Since Cφ enters the exact mode equation (157), and Cφ/H^2 is generically O(ε) on the slow-roll background, the derivation of ns from Eq. (162), which drops the Cφ contribution, is not a consistent first-order slow-roll calculation. Equation (171) and Fig. 1 therefore do not state which value of ξ2+ξ3 (or c_T) is used. Moreover, the requirement c_T>1 for r<0.06 forces ξ2+ξ3 ≠ 0, which changes Cφ and hence ns; the ns curves and the r bound in Fig. 2 cannot be combined without recomputing ns for the c_T values adopted. The authors must specify the treatment of ξ2+ξ3 and give
  2. [III.C, Eq. (38)] The elimination of the homogeneous multiplier sector is valid only for Fourier modes with k≠0 and for boundary conditions that exclude non-trivial harmonic functions. On compact spatial manifolds, which are standard in cosmological perturbation theory, the homogeneous kernel of D^2 is unconstrained, so the global zero modes of the multipliers do not drop out. The paper acknowledges this after Eq. (38), but the central claim that the theory propagates only two tensorial degrees of freedom around a spatially flat FLRW background is then not established for k=0 or for compact spatial topology. Please state the precise mode content of the claim and either prove that the k=0 sector is non-dynamical or explicitly exclude it from the claim.
  3. [IV.B, Eq. (99)-(103), Appendix A] The scalar reduction to a single δφ mode requires the nondegeneracy conditions f≠0, P^3,N ≠ 0, and Eϕ≠0. The paper states these conditions but does not verify them on the inflationary trajectories used in Sec. V. In particular, the denominator Eϕ in Eq. (103) could vanish for some parameter values, and P^3(N) is left completely free. For the claimed predictive scalar sector, the authors should either prove these inequalities hold on the relevant phase space or check them numerically for the parameter range in Figs. 1-2, and also specify a concrete choice of P^3(N) satisfying the generic-branch condition.
minor comments (4)
  1. [Figs. 1-2] The captions should specify the value of ξ2+ξ3 (or c_T) used for each curve, since the scalar spectral index depends on it (see major comment). Without this, the numerical results are not reproducible.
  2. [Eq. (124)] The parameter c is defined through the dimensionless tilde quantities; the text could clarify the dimensions and the relation to Eq. (115)-(118) for readers tracking the tildes.
  3. [Sec. V.A, after Eq. (109)] The 'illustrative subclass' with ξ2=-ξ3 gives luminal tensors, but the later observational analysis requires c_T>1. Please clarify whether the luminal subclass is used only as an example and is not the branch on which the final constraints are placed.
  4. [Sec. V.D, Eq. (171)] The 'GR limit' c→0 does not imply ξ_i→0; the scalar sector may retain deformations for c=0 if ξ2+ξ3≠0. This should be stated to avoid the impression that c=0 reproduces the full GR perturbation theory.

Circularity Check

0 steps flagged

No significant circularity: the perturbation results follow from the defined Hamiltonian; the explicit zero-mode caveat and an under-specified inflationary parameter are limitations, not circular reductions.

full rationale

All central claims are derived from the stated total Hamiltonian (37)/(40) and, in the inflationary section, from the toy Hamiltonian (104). The background equations (42)-(44), tensor action (49)-(52), vector action (61), scalar reduction (95)-(100), and slow-roll spectra (161), (167) are obtained by standard phase-space reduction; no quantity is fitted to those observables and then called a prediction. The Laplacian replacement (37) is a construction choice designed so that homogeneous multipliers drop out, and the paper explicitly flags the required condition after Eq. (38): 'For Fourier modes with k≠0, and under boundary conditions that exclude non-trivial harmonic functions, Eq. (38) is locally equivalent to Q^I−P^I=0. Its homogeneous kernel is deliberately left unconstrained... Global zero modes can depend on the spatial topology and boundary conditions and should be treated separately.' This is a scoping caveat, not a circular identification. The observational comparison uses Planck/BK15 numbers only as external benchmarks after the derivation, as stated in Sec. V.D: 'we adopt the rounded bounds... These values serve as illustrative benchmarks rather than as an updated likelihood analysis.' The self-citations [16]-[19] are contextual for matter coupling in MMG and are not load-bearing; no uniqueness theorem from the authors' prior work is invoked. One non-circular weakness should be noted: Cϕ in Eq. (113) depends on the combination ξ2+ξ3 (equivalently c_T), while the background is parameterized by c (Eq. 124) and A; the paper's Fig. 1 and Eq. (171) do not specify ξ2+ξ3. This makes the reported ns curves under-specified, but it is an ambiguity/reproducibility issue, not a reduction of the prediction to its input.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The paper introduces no new physical entities; it introduces a class of modified gravitational Hamiltonians with arbitrary P^I and ξ parameters. The main assumptions concern boundary conditions on the Laplacian and nondegeneracy of the constraint reduction.

free parameters (3)
  • P^I(N) constraint functions = unspecified (free functions)
    The auxiliary constraints are chosen as S^I=Q^I−P^I(N); P^I are arbitrary functions, subject only to generic conditions like P^3,N ≠ 0. They determine the degenerate vs generic branch and the scalar perturbation constraints (Eq. 66).
  • ξ1, ξ2, ξ3 (deformation parameters) = constrained: c ranges over about -2 to 0.5 and c_T > 1 (Fig. 2)
    In the toy model (Eq. 104), these set the background deformation c=−3(ξ̃3+9ξ̃1+3ξ̃2) and the tensor speed c_T²=1+3/2 f̃(ξ̃2+ξ̃3). They are free parameters of the model, constrained by the adopted n_s and r bounds.
  • A (quadratic potential coefficient) = scanned over 10^-5–10^-2
    The potential V=A φ² is used; A is scanned in the numerical analysis, not fixed by the observed scalar amplitude.
axioms (5)
  • standard math ADM Hamiltonian formulation with phase space (γ_ij, π^ij, N, π_N, N^i, π_i) and the generalized Cayley-Hamilton basis for scalars built from Π, R.
    This sets the framework (Sec. II).
  • domain assumption The four auxiliary constraints are assumed to be independent and second-class on the generic branch (rank ∆=4 in the original construction; analogous nondegeneracy for the modified theory).
    Eqs. (11)-(12) verify this for the original construction; for the Laplacian-modified theory, the constraint algebra is not fully re-derived but assumed via the perturbative count.
  • domain assumption Boundary conditions exclude non-trivial harmonic functions, so D²f=0 implies f=0 for k≠0 and the homogeneous multiplier sector drops out.
    Sec. III, text after Eq. (38). Zero modes are said to be treated separately, but all later claims assume they vanish.
  • standard math The scalar and tensor mode quantization uses the standard Bunch–Davies vacuum and slow-roll approximations.
    Sec. V C; standard practice but an assumption for the power spectra.
  • domain assumption Minimal coupling of a canonical scalar field does not alter the second-class nature of the auxiliary constraints.
    Sec. IV; based on Ref. [20]'s general argument, not re-derived here.

pith-pipeline@v1.3.0-alltime-deepseek · 18937 in / 18156 out tokens · 175955 ms · 2026-08-01T00:47:57.662130+00:00 · methodology

0 comments
read the original abstract

We construct a minimally modified gravity theory that propagates only two tensorial gravitational degrees of freedom around a spatially flat Friedmann--Lema\^itre--Robertson--Walker (FLRW) background and admits a predictive cosmological perturbation theory. We first show that, in the original four-constraint construction, the homogeneous values of the Lagrange multipliers are not fully determined and nevertheless enter the quadratic tensor action, thereby obstructing cosmological predictivity. We remove this ambiguity by coupling the auxiliary constraints to spatial Laplacians of the multipliers. The multiplier sector then drops out of the homogeneous background equations while continuing to constrain the inhomogeneous scalar sector. The tensor dispersion relation generically contains both $k^{2}$ and $k^{4}$ contributions, whereas no propagating gravitational vector or scalar mode is present on the generic branch for which the constraint reduction is nondegenerate. We subsequently study a subclass whose gravitational Hamiltonian density is proportional to the lapse and contains cubic momentum invariants. In this subclass the $k^{4}$ tensor term vanishes and the lapse can be absorbed into a time redefinition at the background level. After coupling a canonical inflaton, the only propagating scalar mode is the inflaton fluctuation, with unit sound speed. For a quadratic potential, departures from general relativity shift the scalar spectral index, while a tensor speed $c_{T}>1$ suppresses the tensor-to-scalar ratio according to $r=16\epsilon_{s}/c_{T}$. The parameter regions compatible with the observational bounds adopted in this work require a superluminal tensor speed; however, very large $c_{T}$ simultaneously reduces the tensor kinetic coefficient and may lower the perturbative cutoff, although determining the actual strong-coupling scale requires a nonlinear analysis.

Figures

Figures reproduced from arXiv: 2607.26031 by Antonio De Felice, Jakkrit Sangtawee, Khamphee Karwan.

Figure 1
Figure 1. Figure 1: The scalar spectral index ns and the combination rcT = 16ϵs as functions of c for several values of A. Each curve terminates when ns leaves the adopted interval 0.965 ± 0.004. the observable predictions provided that it is sufficiently small. The observables are evaluated for modes that cross the Hubble radius 60 e-folds before the end of inflation. Plots of ns and the combination rcT = 16ϵs for various va… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 20 linked inside Pith

  1. [1]

    In this sub- section, primes on perturbations denote derivatives with respect to τ, whereas ϵs, ηs, and ηs2 retain their e-fold definitions

    Curvature perturbation Since the metric curvature perturbation is constrained to vanish in the spatial gauge used above, the gauge- invariant comoving curvature perturbation is carried en- tirely by the scalar-field fluctuation, δϕ=− ˙ϕ N HR=−M P ˜ϕ′ R.(153) Introducing conformal time through dτ= N a dt(154) and defining Qs ≡M 2 Pϵs,(155) the quadratic ac...

  2. [2]

    loss of perturbative control

    Tensor perturbations For either tensor polarization, the quadratic action of the toy model is ST = Z dtd 3x N a3QT " ˙h2 N 2 −c 2 T (∂ih)(∂ih) a2 # ,(163) where QT = M 2 P 8c2 T , c 2 T = 1 + 3 2 ˜f( ˜ξ2 + ˜ξ3).(164) The tensor no-ghost and gradient conditions are therefore QT >0, c 2 T >0.(165) Defining the canonical variable vT,k = a√2QT hk, the leading...

    2018

  3. [3]

    Kosowsky, M

    A. Kosowsky, M. Milosavljevic, and R. Jimenez, Phys. Rev. D66, 063007 (2002), arXiv:astro-ph/0206014

    Pith/arXiv arXiv 2002

  4. [4]

    Durrer, Class

    R. Durrer, Class. Quant. Grav.32, 124007 (2015), arXiv:1506.01907 [astro-ph.CO]

    Pith/arXiv arXiv 2015

  5. [5]

    Sahni, Class

    V. Sahni, Class. Quant. Grav.19, 3435 (2002), arXiv:astro-ph/0202076

    Pith/arXiv arXiv 2002

  6. [6]

    Weinberg, Rev

    S. Weinberg, Rev. Mod. Phys.61, 1 (1989)

    1989

  7. [7]

    A. A. Starobinsky, Phys. Lett. B91, 99 (1980)

    1980

  8. [8]

    A. H. Guth, Phys. Rev. D23, 347 (1981)

    1981

  9. [9]

    B. P. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. Lett.116, 221101 (2016), [Erratum: Phys.Rev.Lett. 121, 129902 (2018)], arXiv:1602.03841 [gr-qc]

    Pith/arXiv arXiv 2016

  10. [10]

    B. P. Abbottet al.(LIGO Scientific, Virgo, Fermi- GBM, INTEGRAL), Astrophys. J. Lett.848, L13 (2017), arXiv:1710.05834 [astro-ph.HE]

    Pith/arXiv arXiv 2017

  11. [11]

    Lovelock, Aequat

    D. Lovelock, Aequat. Math.4, 127 (1970)

    1970

  12. [12]

    Lovelock, J

    D. Lovelock, J. Math. Phys.13, 874 (1972)

    1972

  13. [13]

    Afshordi, D

    N. Afshordi, D. J. H. Chung, and G. Geshnizjani, Phys. Rev. D75, 083513 (2007), arXiv:hep-th/0609150

    Pith/arXiv arXiv 2007

  14. [14]

    Lin and S

    C. Lin and S. Mukohyama, JCAP10, 033, arXiv:1708.03757 [gr-qc]

    Pith/arXiv arXiv

  15. [15]

    Mukohyama and K

    S. Mukohyama and K. Noui, JCAP07, 049, arXiv:1905.02000 [gr-qc]

    Pith/arXiv arXiv 1905

  16. [16]

    Carballo-Rubio, F

    R. Carballo-Rubio, F. Di Filippo, and S. Liberati, JCAP06, 026, [Erratum: JCAP 11, E02 (2018)], arXiv:1802.02537 [gr-qc]

    Pith/arXiv arXiv 2018

  17. [17]

    K. Aoki, C. Lin, and S. Mukohyama, Phys. Rev. D98, 044022 (2018), arXiv:1804.03902 [gr-qc]

    Pith/arXiv arXiv 2018

  18. [18]

    K. Aoki, A. De Felice, S. Mukohyama, K. Noui, M. Oliosi, and M. C. Pookkillath, Eur. Phys. J. C80, 708 (2020), arXiv:2005.13972 [astro-ph.CO]

    Pith/arXiv arXiv 2020

  19. [19]

    Sangtawee and K

    J. Sangtawee and K. Karwan, Phys. Rev. D104, 023511 (2021), arXiv:2103.11463 [astro-ph.CO]

    Pith/arXiv arXiv 2021

  20. [20]

    De Felice, S

    A. De Felice, S. Mukohyama, and M. C. Pookkillath, Phys. Rev. D106, 084050 (2022), arXiv:2206.03338 [gr-qc]

    Pith/arXiv arXiv 2022

  21. [21]

    De Felice, S

    A. De Felice, S. Kumar, S. Mukohyama, and R. C. Nunes, JCAP04, 013, arXiv:2311.10530 [astro-ph.CO]

    Pith/arXiv arXiv

  22. [22]

    Z.-B. Yao, M. Oliosi, X. Gao, and S. Mukohyama, Phys. Rev. D107, 104052 (2023), arXiv:2302.02090 [gr-qc]

    Pith/arXiv arXiv 2023

  23. [23]

    Gao, Phys

    X. Gao, Phys. Rev. D90, 104033 (2014), arXiv:1409.6708 [gr-qc]

    Pith/arXiv arXiv 2014

  24. [24]

    Mukohyama, R

    S. Mukohyama, R. Namba, R. Saitou, and Y. Watanabe, Phys. Rev. D92, 024005 (2015), arXiv:1504.07357 [hep- th]

    Pith/arXiv arXiv 2015

  25. [25]

    Saitou, Phys

    R. Saitou, Phys. Rev. D94, 104054 (2016), arXiv:1604.03847 [hep-th]

    Pith/arXiv arXiv 2016

  26. [26]

    Z.-B. Yao, M. Oliosi, X. Gao, and S. Mukohyama, Phys. Rev. D103, 024032 (2021), arXiv:2011.00805 [gr-qc]

    Pith/arXiv arXiv 2021

  27. [27]

    Mertzios and M

    B. Mertzios and M. Christodoulou, IEEE Transactions on Automatic Control31, 156 (1986)

    1986

  28. [28]

    Aghanimet al.(Planck), Astron

    N. Aghanimet al.(Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

    Pith/arXiv arXiv 2020