REVIEW 2 major objections 3 minor 15 references
Analytical constraints on gravitational models with a quadratic Weyl tensor
T0 review · 2 major / 3 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read The paper constrains the free coupling of a quadratic Weyl term, showing that viable matter-era cosmology requires $10^{14}H_0^2 \lesssim 1/\alpha \ll M_{\text{cutoff}}^2$.
desk verdict A clean new lower bound on the Weyl-squared coupling from matter-era scalar growth, with a real but addressable sensitivity to initial conditions in the scalar sector. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the quadratic Weyl term $-\alpha C^2$, whose coupling $\alpha$ sets the mass scale of the new degrees of freedom: in the vector and tensor sectors the modes acquire an effective squared mass $1/(2\alpha)$, which converts stability into the condition $\alpha>0$. The scalar sector is reduced to a single fourth-order ordinary differential equation for the gauge-invariant matter overdensity $\delta_m$ (equation A.3), obtained by algebraically eliminating the Bardeen potential $\Phi$ from the coupled second-order system. Expanding $\delta_m = \delta_m^{\mathrm{GR}} + \alpha H_0^2\, \delta_m^{(1)} + \dots$ makes the deviation from general relativity explicit as a source term proportional to $K_0^2(1+z)^2\delta_m^{\mathrm{GR}}$, and the demand that this correction stay small is what yields the lower bound on $1/\alpha$. The same coefficient appears in the vector and tensor stability conditions, so the paper's whole argument hangs on tracking the single dimensionless quantity $\alpha H_0^2$.
What would settle it
Compute the present-day ratio $\delta_m^{\mathrm{QG}}/\delta_m^{\mathrm{GR}}$ by numerically integrating the full fourth-order scalar equation with initial amplitudes and velocities for $\delta_m$ spanning, say, $10^{-7}$ to $10^{-3}$ at the initial epoch used in the paper ($N=-10$), while keeping the background fixed. If the ten-percent crossing moves outside the range $1/\alpha \sim 10^{14}H_0^2$, the paper's bound is not robust; conversely, if the crossing stays put for all these initial conditions, the bound is confirmed as a genuine property of the model.
Extended reading notes
Core claim
On the authors' own terms, the discovery is a pair of inequalities on the coupling $\alpha$ in the action $S = (M_{\mathrm{Pl}}^2/2)\int d^4x\sqrt{-g}(R-2\Lambda-\alpha C^2)+S_{\mathrm{mat}}$. In an FLRW background during matter domination, the vector and tensor perturbations carry a mass term proportional to $1/(2\alpha)$, so $\alpha>0$ is required to avoid classical instabilities; for $\alpha<0$ the modes blow up near the present time. In the scalar sector, the gauge-invariant density contrast $\delta_m$ obeys a fourth-order equation, and expanding around the GR growing mode shows that the correction is controlled by $\alpha H_0^2$ times a redshift-weighted source. Requiring that the Jeans instability be recovered and that the present-day matter power spectrum differ from GR by no more than $10\%$ gives $\bar\alpha = \alpha H_0^2 \lesssim 10^{-14}$, i.e. $10^{14}H_0^2 \lesssim 1/\alpha \ll M_{\text{cutoff}}^2$. The paper verifies this by exact numerical integration of the fourth-order equation and by a first-order perturbative solution.
Load-bearing premise
The load-bearing premise is that setting the density perturbation and all its derivatives to $10^{-5}$ at the initial time is a representative choice for the new scalar mode; if the extra mode starts with a different amplitude or velocity, the value of $1/\alpha$ at which deviations from general relativity reach ten percent could change.
Editorial extensions
If this is right
- If the constraint holds, quadratic Weyl gravity with $0<1/\alpha < 10^{14}H_0^2$ is observationally problematic during matter domination: the extra scalar suppresses the growing density mode and changes the matter power spectrum by more than ten percent.
- Within the allowed window, the vector and tensor sectors are classically stable even though the new modes are ghosts; quantum consistency then requires the fakeon prescription, keeping the theory unitary.
- The two-sided inequality is meaningful only when the effective-field-theory cutoff obeys $M_{\text{cutoff}}^2 \gg 10^{14}H_0^2$, so the result forces the cutoff of the low-energy theory to be far above the Hubble scale.
- Gravitational waves in the subhorizon regime propagate at the speed of light in this model, consistent with current constraints, while scalar perturbations deviate from $\Lambda$CDM by at most ten percent at present.
Reading between the lines
- An extension the authors do not pursue: the same fourth-order scalar equation could be integrated from different initial amplitudes and velocities for the new scalar mode; if the ten-percent power-spectrum criterion is sensitive to those choices, the precise value of the lower bound would shift.
- One testable consequence of the window is that the matter power spectrum in this theory should show a characteristic scale-dependent suppression relative to $\Lambda$CDM controlled by $\alpha H_0^2$; future large-scale-structure surveys that measure $P(k)$ at the percent level could therefore directly measure or bound $\alpha$.
- The matter-era constraint does not by itself settle the behavior of the Weyl-squared term in other epochs; applying the same expansion around the appropriate background solutions during radiation domination or inflation could give independent, possibly tighter windows on $\alpha$.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the effects of a term quadratic in the Weyl tensor, -α C^2, on linear cosmological perturbations during the matter-dominated era. The authors derive the equations of motion for vector, tensor, and scalar perturbations, show that α<0 leads to classical instabilities in the vector and tensor sectors, and require α>0. For the scalar sector, they reduce the two-field system to a fourth-order equation for the matter overdensity δm, use a perturbative expansion in α H_0^2 for small α, and impose a criterion that the matter power spectrum deviate from GR by at most 10% at present. Sampling wavenumbers K0 = k/(a0H0) in [10,100], they conclude that α H_0^2 ≲ 10^-14, i.e., 10^14 H_0^2 ≲ α^-1 ≪ M_cutoff^2. The central claim is that this window of α is required for classical stability and for recovering the growth of structure.
Significance. If the claimed bound holds, it provides a concrete constraint on the parameter space of quadratic Weyl gravity from structure formation, complementing the known stability requirement α>0. The paper's vector and tensor analyses are clear and the reduction of the scalar sector to a single fourth-order equation is a useful technical contribution. The authors are transparent about the heuristic 10% criterion and the perturbative nature of the small-α expansion, and they compare exact and approximate solutions. The central numerical claim, however, is conditional on choices that are not fully justified, so the headline inequality should be interpreted as a conditional result rather than a definitive property of the theory.
major comments (2)
- [Sec. 7.2, Figs. 4-5 captions] The fourth-order scalar equation (A.3) is integrated with all initial values of δ_m and its first three derivatives set to 10^-5 at N=-10, which fixes the amplitude and velocity of the additional scalar degree of freedom in an arbitrary way. Since the theory does not prescribe these initial conditions, the present-day ratio δ_m^QG/δ_m^GR (Fig. 5) is a mixture of the physical particular solution and an initial-condition-dependent homogeneous component, and the resulting bound ᾱ ≲ 10^-14 is conditional on this choice. The paper should either derive these initial conditions from a physical prescription (e.g., inflation or adiabatic vacuum) or explicitly test the sensitivity of the bound to varying the extra-mode amplitudes; without this, the abstract's inequality 10^14 H_0^2 ≲ α^-1 is not robust.
- [Sec. 7.2, paragraph before Fig. 5] The 10% power-spectrum deviation criterion is hand-set ('heuristic criterion') and no observational reference is given. Because the boundary value ᾱ ≲ 10^-14 is read off from where the ratio in Fig. 5 crosses the 10% band, the central numerical claim inherits the arbitrariness of this threshold. Please either connect the 10% level to a specific observational constraint on the matter power spectrum or show how the bound shifts for other thresholds (e.g., 5% and 20%) so that the reader can gauge the robustness of 10^14 H_0^2 ≲ α^-1.
minor comments (3)
- [Fig. 5] The y-axis label reads δ_m^GR(0)/δ_m^QG(0) while the caption says the ratio is between δ_m^QG and δ_m^GR; the text's description of 'suppression' for ᾱ=10^-13 is inconsistent with the plotted ratio being less than unity at large K0. Please clarify the convention.
- [Abstract, Sec. 7] The abstract counts 'two scalar degrees of freedom, δ_m and Φ'; in GR, Φ is not an independent propagating degree of freedom, and the theory adds one new scalar gravitational mode. The phrasing could be adjusted to avoid confusion about the number of new d.o.f.
- [Sec. 7.2, numerical integration] The exact integrations of (A.3) for ᾱ=10^-13 and 10^-14 are performed in a stiff regime (coefficients diverge as ᾱ→0); a brief description of the numerical method and error control would increase confidence in the exact/approximate comparison in Fig. 4.
Circularity Check
No significant circularity: the alpha-bounds are imposed by stability and growth-matching requirements, not derived from assuming the bound.
full rationale
The paper's central claim, 10^14 H0^2 ≲ 1/alpha ≪ Mcutoff^2, is obtained by direct analysis rather than by construction. In the vector and tensor sectors, the sign of alpha is constrained through explicit equations of motion, e.g. Eq. (5.4) has a term Y/(2 alpha_bar) whose sign controls stability, and the tensor equation (6.5) is analyzed numerically for positive and negative alpha_bar. In the scalar sector, the paper expands delta_m = delta_GR + alpha H0^2 delta^(1) in Eqs. (7.6)-(7.13) and then, in Sec. 7.2, imposes a stated heuristic criterion that deviations from GR remain below 10% in the matter power spectrum. The bound alpha_bar ≲ 10^-14 is read off from the numerical scan in Fig. 5 over K0 in [10,100]; it is a constraint imposed on the parameter space, not a fitted parameter renamed as a prediction. The cited prior works with overlapping authorship, [7] and [8], are used as contextual motivation (e.g. Minkowski stability and the Starobinsky background) and are not load-bearing inputs to the FLRW derivation. The arbitrary choice of initial conditions for delta_m and its derivatives (set to 10^-5 at N=-10, as stated in the Fig. 4 and Fig. 5 captions) is a possible robustness limitation that could shift the numerical value of the bound, but it is not a circular step: the paper does not choose those initial conditions to force the stated bound, and the bound is not equivalent to its inputs by definition.
Assumptions & free parameters
free parameters (3)
- Power spectrum deviation threshold =
10 percent
- Initial conditions for scalar perturbations at N=-10 =
delta_m = delta_GR = delta^(1) = 10^-5, derivatives also 10^-5
- Sub-horizon sampling range for K0 =
[10, 100]
assumptions (5)
- standard math The Weyl-squared term vanishes on FLRW backgrounds, so the background dynamics are exactly Lambda CDM.
- domain assumption During matter domination the matter fluid is pressureless dust, with w=0, c_s^2=0, and eta=0.
- domain assumption The effective field theory expansion 1 << 1/|alpha-bar| << (M_cutoff/H0)^2 remains valid.
- ad hoc to paper A deviation of at most 10 percent in the matter power spectrum relative to GR is the correct phenomenological consistency criterion.
- domain assumption The perturbative solution delta_m = delta_GR + alpha H0^2 delta^(1) is valid whenever the inequalities in Eq. (7.13) hold.
Cite this review
Pith. "Pith review of Analytical constraints on gravitational models with a quadratic Weyl tensor." pith.science (2026). https://pith.science/paper/SINDBZI4
@misc{pith2026250415005,
author = {Pith},
title = {Pith review of: Analytical constraints on gravitational models with a quadratic Weyl tensor},
year = {2026},
howpublished = {\url{https://pith.science/paper/SINDBZI4}},
note = {Machine review of arXiv:2504.15005}
}
abstract
We set analytical constraints on the parameter space of models of gravity containing a term quadratic in Weyl curvature $-\alpha C^2$. In this class of models, there are four propagating tensorial degrees of freedom, two vector degrees of freedom, and two scalar degrees of freedom, $\delta_m$ and $\Phi$, corresponding to gauge invariant perturbations in the matter density and the gravitational Bardeen potential, respectively. We consider the era of matter domination, and requiring that growth of perturbations are recovered in the scalar sector and classical instabilities are eliminated in the vector and tensor sectors, we obtain bounds on the free coupling parameter to the quadratic Weyl curvature term, $10^{14}H^2_0 \lesssim \alpha^{-1} \ll M^2_{\text{cutoff}}$, where $M_{\text{cutoff}}$ is the cutoff scale of the low energy effective field theory, and $\alpha^{-1}$ is proportional to the masses of the additional propagating degrees of freedom.
Reference graph
Works this paper leans on
-
[1]
T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, Phys. Rept.513, 1 (2012), arXiv:1106.2476 [astro-ph.CO]
arXiv 2012
-
[2]
K. S. Stelle, Gen. Rel. Grav.9, 353 (1978)
work page 1978
-
[3]
E. S. Fradkin and A. A. Tseytlin, Nucl. Phys. B201, 469 (1982)
work page 1982
-
[4]
A. Hindawi, B. A. Ovrut, and D. Waldram, Phys. Rev. D53, 5583 (1996), arXiv:hep-th/9509142
arXiv 1996
-
[5]
C. Bogdanos, S. Capozziello, M. De Laurentis, and S. Nesseris, Astropart. Phys.34, 236 (2010), arXiv:0911.3094 [gr-qc]
arXiv 2010
-
[6]
K. Hinterbichler and M. Saravani, Phys. Rev. D93, 065006 (2016), arXiv:1508.02401 [hep-th]
arXiv 2016
-
[7]
Gravitational Waves in the Spectral Action of Noncommutative Geometry
W. Nelson, J. Ochoa, and M. Sakellariadou, Phys. Rev. D82, 085021 (2010), arXiv:1005.4276 [hep-th]
work page Pith review arXiv 2010
-
[8]
A. De Felice, R. Kawaguchi, K. Mizui, and S. Tsujikawa, Phys. Rev. D108, 123524 (2023), arXiv:2309.01835 [gr-qc]
arXiv 2023
Show all 15 references
-
[9]
Anselmi, JHEP02, 141, arXiv:1801.00915 [hep-th]
D. Anselmi, JHEP02, 141, arXiv:1801.00915 [hep-th]
-
[10]
V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger, Phys. Rept.215, 203 (1992)
1992
-
[11]
J. M. Bardeen, Phys. Rev. D22, 1882 (1980)
1980
-
[12]
P. J. Peebles,The Large-Scale Structure of the Universe (Princeton University Press, 1980)
1980
-
[13]
Lovelock, J
D. Lovelock, J. Math. Phys.12, 498 (1971)
1971
-
[14]
R. L. Arnowitt, S. Deser, and C. W. Misner, Gen. Rel. Grav.40, 1997 (2008), arXiv:gr-qc/0405109
2008 arXiv
-
[15]
B. F. Schutz and R. Sorkin, Annals Phys.107, 1 (1977). – 15 –
1977
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.