REVIEW 4 major objections 5 minor 29 references
Thermodynamics of $f(R)$ Theories
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read If f(R) gravity is read thermodynamically with curvature R as Gibbs energy and the constant Λ as temperature, the reconstructed f(R) from a quadratic Einstein-frame potential is a multivalued function whose unstable branch and van der…
desk verdict Clear and honest lecture notes, but the thermodynamic dictionary is asserted rather than derived, and the homogeneous-mode ambiguity makes the phase-transition language a redescription rather than a prediction. 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 machinery has three parts. The first is the Legendre-transform reconstruction map that sends an Einstein-frame potential to a Jordan-frame f(R): with $\beta=\sqrt{2/3}$, the paper uses $f(\tilde\phi)=e^{2\beta\tilde\phi}[2V(\tilde\phi)+2\beta^{-1}V'(\tilde\phi)]$ and $R(\tilde\phi)=e^{\beta\tilde\phi}[4V(\tilde\phi)+2\beta^{-1}V'(\tilde\phi)]$, which is what makes $f(R)$ multivalued. The second is Catastrophe Theory's fold geometry, summarized by the bifurcation set $4\alpha^3+27\beta^2=0$ for the potential $\frac14 x^4+\frac{\alpha}{2}x^2+\beta x$; it says where two extrema coalesce, which is exactly the swallowtail seen in the reconstructed f(R). The third is the thermodynamic dictionary $R\leftrightarrow G$, $f\leftrightarrow P$, $\Lambda\leftrightarrow T$, from which the paper derives $G(P,T)$, $P(V,T)$, $F(T,V)$ and $S=-2/V$ and plots the spinodal and binodal curves.
What would settle it
Integrate the full Jordan-frame background equations for the reconstructed f(R) of Eqs. (50)-(51) with $V(\phi)=\frac12 m^2(\phi-a)^2+\Lambda$, starting from slow-roll initial conditions, without using any thermodynamic dictionary. If the trajectory exits inflation while always keeping $f''>0$ and never entering the region $dP/dV>0$, the first-order-transition picture is not describing the dynamics; if it crosses into the spinodal region and then settles on the stable small-volume branch, the correspondence is confirmed.
Extended reading notes
Core claim
The author's central claim is that, for the Einstein-frame potential $V(\tilde\phi)=\frac12 m^2(\tilde\phi-a)^2+\Lambda$, the conformal reconstruction of the Jordan-frame Lagrangian $f(R)$ produces a swallowtail-shaped curve in the $(R,f)$ plane: three branches for small $\Lambda$, of which the middle one has $f''<0$ and is therefore unstable. The threshold $\Lambda_c=15/16$ acts like a critical temperature: above it the unstable branch disappears. Once the dictionary $R\leftrightarrow G$, $f\leftrightarrow P$, $\Lambda\leftrightarrow T$ is adopted, the reconstruction formulas give explicit $G(P,T)$, an effective pressure $P(V,T)$, a Helmholtz energy $F(T,V)$, and an entropy $S=-2/V$, and the equation of state reproduces the binodal and spinodal curves of the van der Waals gas. The paper states as its final summary that there is a strong correspondence between f(R) theories in the metric approach and a first-order phase transition as described by Catastrophe Theory, with inflation as the metastable phase.
Load-bearing premise
The load-bearing premise is the dictionary $R\leftrightarrow G$, $f\leftrightarrow P$, and $\Lambda\leftrightarrow T$; the paper asserts this final step as 'almost automatic' rather than deriving it from statistical mechanics, and if it is only a formal analogy the phase-transition language is a redescription of the swallowtail geometry rather than an independent thermodynamic property of f(R) gravity.
Editorial extensions
If this is right
- For $\Lambda<15/16$, the reconstructed theory has a genuine unstable branch $f''<0$; the same condition that destroys stability of Minkowski perturbations also locates the spinodal region of the effective fluid.
- Inflationary initial conditions mapped to the large-volume branch sit below the binodal, so the inflationary phase is metastable by construction and lasts a finite number of e-folds rather than being an attractor.
- The final configuration, with $\tilde\phi$ oscillating at the bottom of the quadratic potential, has $\langle \tilde\omega_\phi\rangle\approx0$, i.e. a matter-dominated phase, so the transition out of inflation is into the standard matter era.
- Because the reconstruction formulas (50)-(51) are general, the swallowtail structure is a feature of the reconstruction map itself, so more complex Einstein-frame potentials should also organize their phases through the same catastrophe-theory geometry.
Reading between the lines
- Editorial extension: if the dictionary is taken literally, the entropy $S=-2/V$ being negative signals that the curvature degree of freedom is an open subsystem; the latent heat of the transition must be carried by the radiation or matter sector, which gives a concrete check: the reheating temperature after inflation should equal the latent heat read off the binodal curve of this equation of state
- Editorial extension: the same formalism could be run in reverse as a classification scheme—any f(R) whose $(R,f)$ curve has a swallowtail would define an effective fluid with two phases, and the control-parameter threshold (the analogue of $\Lambda_c$) would predict whether the theory admits a first-order exit.
- Editorial extension: a testable extension is to compute the curvature perturbation spectrum across the spinodal region; a first-order transition would imprint characteristic non-Gaussianity or bubble signatures in the CMB, whereas their absence would indicate the thermodynamic language is a formal redescription.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript is a set of lecture notes on f(R) gravity and its proposed thermodynamic reinterpretation. After reviewing metric f(R) gravity, the Einstein frame, stability constraints, and a catastrophe-theory description of the van der Waals phase transition, Section III reverses the usual reconstruction: starting from the Einstein-frame potential V(φ) = 1/2 m^2(φ-a)^2 + Λ, the paper uses the reconstruction formulas of Eqs. (50)-(51) to obtain a parametric f(R), identifies R with the Gibbs energy G and f with the pressure P, identifies Λ with the temperature T, and derives an effective equation of state P(V,T) in Eq. (61). The paper claims that the resulting multivalued f(R), the unstable branch with f''<0 for Λ<15/16, and the spinodal/binodal structure of P(V,T) establish a strong correspondence between f(R) theories and first-order phase transitions, with inflation as a metastable phase and its exit as a first-order transition.
Significance. If the thermodynamic dictionary were derived rather than assumed, the paper would offer a striking connection between f(R) gravity and equilibrium thermodynamics, with explicit analytical formulas: the reconstructed f(R), the effective pressure P(V,T), the spinodal and binodal curves, and the identification of the unstable branch with the cusp catastrophe are all obtained in closed form. The paper is also valuable as a pedagogical bridge between two otherwise disjoint literatures. However, the central claim rests on identifications that are asserted in Section III rather than derived from statistical mechanics, and at least one of those identifications is underdetermined by the reconstruction equations. The explicit algebraic construction is internally consistent, but the physical interpretation is not yet established at the level claimed in the final summary.
major comments (4)
- [Section III, Eqs. (50)-(51)] The identification of R with the Gibbs energy G and f(R) with the pressure P is the load-bearing step of the paper, but it is presented only as 'the final step is almost automatic' after Eq. (51). No derivation from an ensemble, a partition function, or a statistical-mechanical principle is given. All subsequent thermodynamic statements, including the phase-transition interpretation, depend on this dictionary. The authors should either derive this mapping from a concrete microscopic model or explicitly reframe the claim as a formal mathematical analogy between the swallowtail geometry of the reconstruction and the cusp catastrophe of a van der Waals fluid. As it stands, the manuscript does not establish that the thermodynamic variables are properties of f(R) gravity rather than a relabeling of the reconstruction parameters.
- [Section III, Eqs. (50)-(51) and Eq. (60)] The reconstruction equations have a homogeneous mode: adding V_h(φ) = C e^{-βφ} to the Einstein-frame potential leaves f(φ) in Eq. (50) unchanged while shifting R(φ) in Eq. (51) by the constant 2C. The parametric curve (f(φ), R(φ)) is therefore defined only up to adding a constant to the would-be Gibbs energy G. Equation (60) explicitly exploits exactly this freedom by adding -2T e^{βa} to G to make the entropy nonnegative, without any principle fixing the constant C(T). Consequently, the entropy S and Helmholtz free energy F are not determined by the f(R) construction; only P(V,T), the spinodal, and the binodal are invariant under this shift. The paper should identify a physical or mathematical selection rule for C(T), or state plainly that the free energy is defined only up to the homogeneous mode and is therefore not a predictive thermodynamic output of the theory.
- [Section III, Eqs. (53)-(54)] The identification of the constant Λ in the quadratic potential (49) with the thermodynamic temperature T is an additional assumption that is never justified. In Eqs. (53)-(54), T appears precisely where Λ appeared in V(φ), but there is no argument connecting a constant shift in a scalar potential to an equilibrium temperature, nor any discussion of why this constant should be the same for all isotherms. If the Λ-T identification is merely a formal bookkeeping device, then the temperature dependence of P(V,T) and the location of the critical point are not physical predictions. The authors should provide a criterion for when a parameter in an Einstein-frame potential can be promoted to a thermodynamic temperature, or alternatively present the result as a mathematical correspondence with the caveat that T is a control parameter, not a thermal temperature.
- [Section III, Eqs. (58)-(60)] The unmodified entropy computed from Eqs. (58)-(59) is S(T,V) = -2/V, which is negative for all positive V. The text acknowledges that this indicates the system is incomplete, and then repairs the Gibbs energy by adding -2T e^{βa} in Eq. (60). This repair is ad hoc: it changes G and S while leaving P, V, and the phase-transition curves invariant, and no physical origin is given for the added term. Since entropy is a central thermodynamic quantity, the need for an arbitrary additive correction undermines the quantitative content of the claimed correspondence. The authors should derive the correction from a microscopic model or from a principle such as extensivity, rather than inserting it by hand to remove a negative entropy.
minor comments (5)
- [Section II A] There is a typo in the opening sentence: 'Iin spite of' should read 'In spite of'.
- [Section I E 2] Reference [21] is written as 'arXiv:2308.00203 [gr-qc]' followed by placeholder question marks '??'; the citation details should be completed.
- [Section II B] The symbol V is used both for the volume in the van der Waals equation and for the potential function in Eq. (45). The footnote disambiguates the two uses, but the notation remains confusing in a few passages, especially when V(x) is plotted against the volume axis.
- [Section I D, Eq. (22)] The linearized perturbation equation around Minkowski space is stated without derivation. The step from Eq. (19) to Eq. (22) is not obvious, and the notation R = T + δR is potentially confusing because R and T are both scalars; a short derivation or a reference to the original calculation would improve clarity.
- [Section III, Eqs. (54) and (61)] The paper notes that P ∝ T V^{-2} in the high-temperature limit, which differs from the ideal-gas behavior P ∝ T V^{-1}. This is a useful caveat, but it deserves a brief discussion of whether the identification of P as a thermodynamic pressure remains physically meaningful when the ideal-gas limit is not recovered.
Circularity Check
Thermodynamic dictionary is underdetermined: identifying P with f and G with R makes the claimed phase-transition correspondence a definitional consequence of the reconstruction, and the full correspondence is deferred to a self-citation.
-
self definitional
[Section III, after Eq. (51) and before Eqs. (53)-(55)]
"Besides, if we look at Eq. (2) as the Lagrangian of a relativistic fluid, described by its pressure P , the final step is almost automatic: to make the correspondence between R itself and the Gibbs Energy G."
The central 'correspondence' to a first-order phase transition is obtained by defining P:=f (Eq. 50) and G:=R (Eq. 51). Under this dictionary, the swallowtail geometry of the reconstructed f(R) for the quadratic potential is, by construction, the swallowtail of G(P,T). The phase-transition language is not a testable consequence of f(R) gravity; it is a restatement of the cusp structure of the Legendre transform. No statistical-mechanical argument fixes these identifications, so the claimed 'strong correspondence' reduces to the chosen dictionary.
-
other
[Section III, after Eq. (59) and Eq. (60)]
"Secondly, we must correct its negative value, since it prevents a physical interpretation in terms of the number of accessible states. For that, it suffices to add an extra term in Eq. (53) and redefine the Gibbs Energy it as [Eq. (60)], which does not spoil the previous results."
The entropy S=-2/V is presented as a prediction, but it comes out negative and is then corrected by adding an arbitrary T-dependent constant -2T e^{βa} to G. This constant leaves P(V,T), the spinodal, and the binodal unchanged, so it is not fixed by the f(R) reconstruction. The thermodynamic entropy and free energy are therefore inputs chosen to make the analogy physical, not outputs of the theory. This exposes the underdetermination of the thermodynamic dictionary.
1 more flagged steps
-
self citation load bearing
[Section III, opening and closing paragraphs]
"Here we follow Ref. [25]. ... We refer the reader to Ref. [25] to a numerical description and the full Thermodynamic correspondence."
Ref. [25] is Peralta & Jorás, a prior paper by the present author. The load-bearing thermodynamic dictionary—the identification of R with G and the 'full Thermodynamic correspondence'—is deferred to this self-citation rather than derived here from independent statistical mechanics or an external benchmark. Since that dictionary is what turns the swallowtail geometry into a phase transition, the central claim leans on an unverified self-citation.
full rationale
The paper's mathematical reconstruction is self-contained: for the quadratic Einstein-frame potential V=1/2(φ-a)^2+Λ, Eqs. (50)-(51) really produce a multivalued f(R) with an unstable branch, and the computed P(V,T) is a genuine cusp-catastrophe equation of state. This part is not fitted. The circularity enters when the paper labels f as pressure and R as Gibbs energy ('the final step is almost automatic') and then presents the resulting swallowtail as a first-order phase transition in f(R) theories. With that dictionary, the phase-transition structure is equivalent by construction to the Legendre reconstruction. The arbitrary T-dependent shift added in Eq. (60) to cure the negative entropy confirms that the thermodynamic quantities are not uniquely determined by f(R). Additionally, the 'full Thermodynamic correspondence' is attributed to the author's own prior work rather than to an independent derivation. I therefore score partial circularity (6), not total (8-10), because the underlying cusp geometry and the explicit equation of state are real, independently checkable mathematics.
Assumptions & free parameters
free parameters (3)
- a =
0 and 0.2 (plotted), arbitrary
- Lambda (identified as temperature T) =
Lambda_c=15/16 critical value; panels use 0, 0.2, 1.2
- m (scalar mass) =
set to 1 implicitly
assumptions (4)
- standard math Reconstruction formulas (50)-(51) from Ref. [26] correctly give the Jordan-frame f(R) corresponding to an Einstein-frame scalar potential.
- domain assumption The scalar field phi and conformal factor f'(R)=e^{beta phi} are positive throughout the relevant branches.
- ad hoc to paper The curvature scalar R can be identified with Gibbs free energy G and the Lagrangian f(R) with pressure P.
- ad hoc to paper The entropy can be corrected by adding the term -2T e^{beta a} to the Gibbs energy.
invented entities (1)
-
Effective thermodynamic variables (P, V, T, G, S) for f(R) gravity
Cite this review
Pith. "Pith review of Thermodynamics of $f(R)$ Theories." pith.science (2026). https://pith.science/paper/HUPPUJ53
@misc{pith2026241118469,
author = {Pith},
title = {Pith review of: Thermodynamics of $f(R)$ Theories},
year = {2026},
howpublished = {\url{https://pith.science/paper/HUPPUJ53}},
note = {Machine review of arXiv:2411.18469}
}
abstract
This series of three lectures was presented at ``Escola de Cosmologia e Gravita\c{c}\~ao" \url{https://cosmosecontexto.org.br/ecg-inscricoes} and webcast at \url{https://shorturl.at/2ZSI7} -- in Portuguese, but slides in English. We will go through a brief review on $f(R)$ theories (in the metric approach) and the usual requirements for successful modifications of General Relativity. Then we will open a large parenthesis to talk about a non-standard approach to Phase Transitions: the Catastrophe Theory. Finally, we will connect the previous lectures to introduce a new Thermodynamic interpretation of $f(R)$ theories.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[25]
Neutron star masses in $R^{2}$-gravity
F. Sbis` a, P. O. Baqui, T. Miranda, S. E. Jor´ as, and O. F. Piattella, Physics of the Dark 23 Universe 27, 100411 (2020), arXiv:1907.08714 [gr-qc]
work page Pith review arXiv 2020
-
[1]
Cosmology When it comes to the background evolution of the universe, one should require that we should “start” in a Radiation-Dominated (RD) universe, pass through a Matter-Dominated (MD) phase — that should last long enough so that matter perturbations can grow — 8 followed by the current accelerated Dark-Energy Dominated (DED) phase. One could also add ...
-
[2]
Relativistic Stars I briefly mention that, as a serious candidate for a full theory of gravity, a given f (R) must be tested in different scenarios, such as, relativistic stars. Spherically-symmetric stars are the perfect equivalence to homogeneous cosmology, since they both have only one independent variable: in the former, the radial coordinate r; in th...
-
[3]
Local gravity constraints No review (however short) on f (R) would be complete without mentioning (however briefly) the chameleon effect [23]. This mechanism explains how the extra degree of freedom does not propagate too far — which would have shown up as an extra force in current experiments. In a nut shell, the authors show that the effective mass of t...
-
[4]
T. P. Sotiriou and V. Faraoni, Reviews of Modern Physics 82, 451 (2010), arXiv:0805.1726 [gr-qc]. 22
arXiv 2010
-
[5]
Modified actions for gravity: Theory and phenomenology,
T. P. Sotiriou, “Modified actions for gravity: Theory and phenomenology,” (2007), arXiv:0710.4438 [gr-qc]
arXiv 2007
-
[6]
A. De Felice and S. Tsujikawa, Living Reviews in Relativity 13, 3 (2010), arXiv:1002.4928 [gr-qc]
arXiv 2010
- [7]
Show all 29 references
-
[8]
H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd Edition (1985)
1985
-
[9]
Gilmore, Catastrophe Theory for Scientists and Engineers (Dover, New York, 1981)
R. Gilmore, Catastrophe Theory for Scientists and Engineers (Dover, New York, 1981)
1981
-
[10]
Poston and I
T. Poston and I. Stewart, Catastrophe Theory and Its Applications , Dover books on mathe- matics (Dover Publications, 1996)
1996
-
[11]
P. T. Saunders, An Introduction to Catastrophe Theory (Cambridge University Press, 1980)
1980
-
[12]
Milgrom, Astrophys
M. Milgrom, Astrophys. J. 270, 371 (1983)
1983
-
[13]
Marra, D
V. Marra, D. C. Rodrigues, and A. O. F. de Almeida, Monthly Notices of the Royal Astro- nomical Society 494, 2875–2885 (2020)
2020
-
[14]
G. J. Olmo, International Journal of Modern Physics D 20, 413 (2011), arXiv:1101.3864 [gr- qc]
2011 arXiv
-
[15]
T. P. Sotiriou and S. Liberati, Annals of Physics 322, 935 (2007), arXiv:gr-qc/0604006 [gr-qc]
2007 arXiv
-
[16]
Palatini f (r) gravity tests in the weak field limit: Solar system, seismology and galaxies,
A. Hernandez-Arboleda, D. C. Rodrigues, J. D. Toniato, and A. Wojnar, “Palatini f (r) gravity tests in the weak field limit: Solar system, seismology and galaxies,” (2023), arXiv:2306.04475 [gr-qc]
2023 arXiv
-
[17]
D. I. Kaiser, Phys. Rev. D 52, 4295 (1995), arXiv:astro-ph/9408044 [astro-ph]
1995 arXiv
-
[18]
Pogosian and A
L. Pogosian and A. Silvestri, Phys. Rev. D 77, 023503 (2008)
2008
-
[19]
Amendola, R
L. Amendola, R. Gannouji, D. Polarski, and S. Tsujikawa, Physical Review D 75 (2007), 10.1103/physrevd.75.083504
2007 doi
-
[20]
Herrera, I
D. Herrera, I. Waga, and S. Jor´ as, Physical Review D95 (2017), 10.1103/physrevd.95.064029
2017 doi
-
[21]
Bertschinger and P
E. Bertschinger and P. Zukin, Phys. Rev. D 78, 024015 (2008), arXiv:0801.2431 [astro-ph]
2008 arXiv
-
[22]
Borisov, B
A. Borisov, B. Jain, and P. Zhang, Phys. Rev. D 85, 063518 (2012), arXiv:1102.4839 [astro- ph.CO]
2012 arXiv
-
[23]
R. C. Batista, Universe 8, 22 (2021)
2021
-
[24]
J. M. Z. Pretel, S. E. Jor´ as, R. R. R. Reis, S. B. Duarte, and J. D. V. Arba˜ nil, arXiv e-prints , arXiv:2308.00203 (2023), arXiv:2308.00203 [gr-qc]
2023 arXiv
-
[26]
Khoury and A
J. Khoury and A. Weltman, Phys. Rev. D 69, 044026 (2004), arXiv:astro-ph/0309411 [astro- ph]
2004 arXiv
-
[27]
J. D. van der Waals, 1873, Ph.D. thesis, Universiteit Leiden (Over de continuiteit van den gas- en vloeistoftoestand (On the Continuity of the Gaseous and Liquid States))
-
[28]
Peralta and S
C. Peralta and S. Jor´ a s, Journal of Cosmology and Astroparticle Physics 2020, 053 (2020)
2020
-
[29]
Magnano and L
G. Magnano and L. M. Soko lowski, Phys. Rev. D50, 5039 (1994), arXiv:gr-qc/9312008 [gr-qc]
1994 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.