Thermodynamics of fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) Theories

S.E. JorƔs Instituto de Fƭsica, Universidade Federal do Rio de Janeiro,
CEP 21941-972 Rio de Janeiro, RJ, Brazil
Abstract

This series of three lectures was presented at ā€œEscola de Cosmologia e GravitaĆ§Ć£oā€ https://cosmosecontexto.org.br/ecg-inscricoes and webcast at https://shorturl.at/2ZSI7 ā€” in Portuguese, but slides in English.

We will go through a brief review on fā¢(R)š‘“š‘…f(R)italic_f ( italic_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)š‘“š‘…f(R)italic_f ( italic_R ) theories.

This a very short review on fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) theory in the metric approach along with a recently developed Thermodynamics interpretation. The main goal here is to detail such new approach, making an analogy to Thermodynamics ā€” which is written in terms of Catastrophe Theory. The latter deals with the coalescence and birth of extrema in an ordinary (polynomial) function as a given set of free parameters are varied. We will see that is precisely what happens at a first-order phase transition and, therefore, it describes exactly the exit from the inflationary regime in the early universe.

Long reviews on fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) can be found in Refs.Ā SotiriouĀ andĀ Faraoni (2010); Sotiriou (2007); De FeliceĀ andĀ Tsujikawa (2010) as well as Chap.Ā 9 of Ref.Ā AmendolaĀ andĀ Tsujikawa (2015). A modern approach to phase transitions can be found in Ref.Ā Callen (1985). Reviews on Catastrophes are available in Refs.Ā Gilmore (1981); PostonĀ andĀ Stewart (1996); Saunders (1980)

I fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) Theories

I.1 Motivation

The motivations for modifying such a successful theory as General Relativity (GR, from now on) can be divided in two (non-exclusive) categories: theoretical and experimental evidences. The former usually relies on arguments from either the absence of renormalizability, or the first corrections from a possible theory of quantum gravity (whatever that might be), or just because that is the way to better appreciate the power and beauty of GR111The latter is my personal choice, I admit.. The latter case ā€” experimental evidence ā€” is also compelling. We believe the universe experienced, in its early stages, a short but tremendously strong phase of almost-exponential expansion ā€” inflation ā€” and/or a bounce (a contraction phase followed by expansion), besides the current mild accelerated expansion. The latter may be ā€œexplainedā€ by an ad-hoc Cosmological Constant Ī›Ī›\Lambdaroman_Ī›, which must be seen either as a new constant of nature (such as Newtonā€™s Gravitational Constant GNsubscriptšŗš‘G_{N}italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT itself, for instance) or as a consequence from a more fundamental principle (such as a suitably renormalized ā€œvacuum energyā€). Even so, that is not the case for the primordial accelerated phase, since a simple Cosmological Constant Ī›Ī›\Lambdaroman_Ī› wouldnā€™t have allowed an exit (graceful or not) from inflation. The universe would not have been through decelerated expansion (radiation- and matter-dominated) phases. In either case, the standard replacement for Ī›Ī›\Lambdaroman_Ī› is a scalar field, with, truth to be told, not-so-standard potential end/or kinetic terms (k-essence) and couplings.

I should say that some authors do investigate if modifications of GR can explain the rotation curves of galaxies (see, for instance, works on MOND Milgrom (1983); MarraĀ etĀ al. (2020)), but we will not pursue this goal here.

Either way, we are currently at a crossroads in GR, just like Newtonā€™s Gravitation (NG) was in the XIX century, at two similar situations, with opposite results. In one of them, we relied on NG to explain small disturbances on the orbit of Uranus, which ultimately led to the discovery of Neptune. On the other hand, the explanation of Mercuryā€™s orbit could not be explained by an inner planet called Vulcan (simply because it was never found at the ā€œrightā€ position) and, ultimately, GR extended NG. Now, like then, we have to decide if we should change the theory or add an unknown component of the universe (a new ā€œplanetā€, particle or scalar field, perhaps?) but keep the theory (GR).

I.2 Modified Einstein Equations

Here, we will focus on the so-called fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) theory in the metric frame, where the metric gĪ¼ā¢Ī½subscriptš‘”šœ‡šœˆg_{\mu\nu}italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT is the one independent variable. We just briefly mention the Palatini Olmo (2011) and Metric-Affine SotiriouĀ andĀ Liberati (2007); Hernandez-ArboledaĀ etĀ al. (2023) approaches, according to which the Affine Connection Ī“Ī¼ā¢Ī½Ī±subscriptsuperscriptĪ“š›¼šœ‡šœˆ\Gamma^{\alpha}_{\mu\nu}roman_Ī“ start_POSTSUPERSCRIPT italic_Ī± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT is an independent object as well.

Of course, if fā¢(R)=Rāˆ’2ā¢Ī›š‘“š‘…š‘…2Ī›f(R)=R-2\Lambdaitalic_f ( italic_R ) = italic_R - 2 roman_Ī›, then we have GR with a standard Cosmological Constant. Therefore, we shall look for non-linear functions only. To be more accurate, we replace the standard Einstein-Hilbert Lagrangian (plus an extra term that describes the matter/radiation/field content of the universe)

ā„’Gā¢R=āˆ’gā¢12ā¢Īŗ2ā¢(Rāˆ’2ā¢Ī›)+ā„’msubscriptā„’šŗš‘…š‘”12superscriptšœ…2š‘…2Ī›subscriptā„’š‘š{\cal L}_{GR}=\sqrt{-g}\frac{1}{2\kappa^{2}}(R-2\Lambda)+{\cal L}_{m}caligraphic_L start_POSTSUBSCRIPT italic_G italic_R end_POSTSUBSCRIPT = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( italic_R - 2 roman_Ī› ) + caligraphic_L start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT (1)

by

ā„’Gā¢R=āˆ’gā¢12ā¢Īŗ2ā¢[fā¢(R)āˆ’2ā¢Ī›]+ā„’m,subscriptā„’šŗš‘…š‘”12superscriptšœ…2delimited-[]š‘“š‘…2Ī›subscriptā„’š‘š{\cal L}_{GR}=\sqrt{-g}\frac{1}{2\kappa^{2}}[f(R)-2\Lambda]+{\cal L}_{m},caligraphic_L start_POSTSUBSCRIPT italic_G italic_R end_POSTSUBSCRIPT = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ italic_f ( italic_R ) - 2 roman_Ī› ] + caligraphic_L start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT , (2)

where Īŗ2ā‰”1/Mpā¢l2ā‰”8ā¢Ļ€ā¢Gsuperscriptšœ…21superscriptsubscriptš‘€š‘š‘™28šœ‹šŗ\kappa^{2}\equiv 1/M_{pl}^{2}\equiv 8\pi Gitalic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ā‰” 1 / italic_M start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ā‰” 8 italic_Ļ€ italic_G. The former yields the Einstein Equations

RĪ¼ā¢Ī½āˆ’12ā¢gĪ¼ā¢Ī½ā¢R+Ī›ā¢gĪ¼ā¢Ī½=Īŗ2ā¢TĪ¼ā¢Ī½,subscriptš‘…šœ‡šœˆ12subscriptš‘”šœ‡šœˆš‘…Ī›subscriptš‘”šœ‡šœˆsuperscriptšœ…2subscriptš‘‡šœ‡šœˆR_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R+\Lambda g_{\mu\nu}=\kappa^{2}T_{\mu\nu},italic_R start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_R + roman_Ī› italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT = italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT , (3)

while the latter yields

RĪ¼ā¢Ī½ā¢fā€²āˆ’12ā¢gĪ¼ā¢Ī½ā¢f+gĪ¼ā¢Ī½ā¢ā–”ā¢fā€²āˆ’āˆ‡Ī¼āˆ‡Ī½ā”fā€²=Īŗ2ā¢TĪ¼ā¢Ī½,subscriptš‘…šœ‡šœˆsuperscriptš‘“ā€²12subscriptš‘”šœ‡šœˆš‘“subscriptš‘”šœ‡šœˆā–”superscriptš‘“ā€²subscriptāˆ‡šœ‡subscriptāˆ‡šœˆsuperscriptš‘“ā€²superscriptšœ…2subscriptš‘‡šœ‡šœˆR_{\mu\nu}f^{\prime}-\frac{1}{2}g_{\mu\nu}f+g_{\mu\nu}\,\Box f^{\prime}-\nabla% _{\mu}\nabla_{\nu}f^{\prime}=\kappa^{2}T_{\mu\nu},italic_R start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_f + italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT ā–” italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT - āˆ‡ start_POSTSUBSCRIPT italic_Ī¼ end_POSTSUBSCRIPT āˆ‡ start_POSTSUBSCRIPT italic_Ī½ end_POSTSUBSCRIPT italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT = italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT , (4)

which looks (and actually is) much more complicated! For starters, one notices that it is a fourth-order differential equation in the metric (we recall the reader that the Ricci scalar Rš‘…Ritalic_R itself depends on second-order derivatives on gĪ¼ā¢Ī½subscriptš‘”šœ‡šœˆg_{\mu\nu}italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT). As such, it bears an extra degree of freedom, that will become even more explicit shortly.

For now, we can make Eq.Ā (4) more palatable if it is rewritten as

RĪ¼ā¢Ī½āˆ’12ā¢gĪ¼ā¢Ī½ā¢R=Īŗ2ā¢TĪ¼ā¢Ī½+[RĪ¼ā¢Ī½āˆ’12ā¢gĪ¼ā¢Ī½ā¢Rāˆ’(RĪ¼ā¢Ī½ā¢fā€²āˆ’12ā¢gĪ¼ā¢Ī½ā¢f+gĪ¼ā¢Ī½ā¢ā–”ā¢fā€²āˆ’āˆ‡Ī¼āˆ‡Ī½ā”fā€²)].subscriptš‘…šœ‡šœˆ12subscriptš‘”šœ‡šœˆš‘…superscriptšœ…2subscriptš‘‡šœ‡šœˆdelimited-[]subscriptš‘…šœ‡šœˆ12subscriptš‘”šœ‡šœˆš‘…subscriptš‘…šœ‡šœˆsuperscriptš‘“ā€²12subscriptš‘”šœ‡šœˆš‘“subscriptš‘”šœ‡šœˆā–”superscriptš‘“ā€²subscriptāˆ‡šœ‡subscriptāˆ‡šœˆsuperscriptš‘“ā€²R_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R=\kappa^{2}T_{\mu\nu}+\Bigg{[}R_{\mu\nu}-% \frac{1}{2}g_{\mu\nu}R-\bigg{(}R_{\mu\nu}f^{\prime}-\frac{1}{2}g_{\mu\nu}f+g_{% \mu\nu}\,\Box f^{\prime}-\nabla_{\mu}\nabla_{\nu}f^{\prime}\bigg{)}\Bigg{]}.italic_R start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_R = italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT + [ italic_R start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_R - ( italic_R start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT italic_f + italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT ā–” italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT - āˆ‡ start_POSTSUBSCRIPT italic_Ī¼ end_POSTSUBSCRIPT āˆ‡ start_POSTSUBSCRIPT italic_Ī½ end_POSTSUBSCRIPT italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ) ] . (5)

The term inside the square brackets on the right-hand side can be read as an effective energy-momentum tensor for the so-called ā€œcurvature fluidā€ Īŗ2ā¢TĪ¼ā¢Ī½(c)superscriptšœ…2superscriptsubscriptš‘‡šœ‡šœˆš‘\kappa^{2}T_{\mu\nu}^{(c)}italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_c ) end_POSTSUPERSCRIPT ā€” since it is pure geometry. Such form is more than adequate (as opposed to an alternate form that divides the previous equation by fā€²superscriptš‘“ā€²f^{\prime}italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT, that will not be adopted here) because the energy-momentum tensor of the standard components and its conservation does not need to change. Consequently, the curvature fluid is conserved by itself as well.

One can further define the energy density ĻcsubscriptšœŒš‘\rho_{c}italic_Ļ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT and pressure pcsubscriptš‘š‘p_{c}italic_p start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT for such fluid, as

Īŗ2ā¢Ļcsuperscriptšœ…2subscriptšœŒš‘\displaystyle\kappa^{2}\rho_{c}italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_Ļ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰”12ā¢(fā€²ā¢Rāˆ’f)āˆ’3ā¢Hā¢fā€²Ė™+3ā¢H2ā¢(1āˆ’fā€²)absent12superscriptš‘“ā€²š‘…š‘“3š»Ė™superscriptš‘“ā€²3superscriptš»21superscriptš‘“ā€²\displaystyle\equiv\frac{1}{2}\left(f^{\prime}R-f\right)-3H\dot{f^{\prime}}+3H% ^{2}(1-f^{\prime})ā‰” divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT italic_R - italic_f ) - 3 italic_H overĖ™ start_ARG italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT end_ARG + 3 italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ) (6)
Īŗ2ā¢pcsuperscriptšœ…2subscriptš‘š‘\displaystyle\kappa^{2}p_{c}italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰”fā€²ĀØ+2ā¢Hā¢fā€²Ė™āˆ’(2ā¢HĖ™+3ā¢H2)ā¢(1āˆ’fā€²)+12ā¢(fāˆ’fā€²ā¢R).absentĀØsuperscriptš‘“ā€²2š»Ė™superscriptš‘“ā€²2Ė™š»3superscriptš»21superscriptš‘“ā€²12š‘“superscriptš‘“ā€²š‘…\displaystyle\equiv\ddot{f^{\prime}}+2H\dot{f^{\prime}}-(2\dot{H}+3H^{2})(1-f^% {\prime})+\frac{1}{2}(f-f^{\prime}R).ā‰” overĀØ start_ARG italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT end_ARG + 2 italic_H overĖ™ start_ARG italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT end_ARG - ( 2 overĖ™ start_ARG italic_H end_ARG + 3 italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ( 1 - italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ) + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_f - italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT italic_R ) . (7)

One can then proceed to investigate the system as if it were described by GR in the presence of two fluids: the curvature fluid defined above and the standard matter/radiation fluid. We point out that the curvature fluid does not need to satisfy the standard (strong, weak, null) energy conditions, since it is not an actual physical fluid.

I.3 The Frames

In this section, we will make it explicit the extra degree of freedom inherent to fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) theories. We start over from the Lagrangian (2) and write it as

ā„’g=āˆ’gā¢12ā¢Īŗ2ā¢[Ļ•ā¢Rā¢(Ļ•)āˆ’Wā¢(Ļ•)]+ā„’msubscriptā„’š‘”š‘”12superscriptšœ…2delimited-[]italic-Ļ•š‘…italic-Ļ•š‘Šitalic-Ļ•subscriptā„’š‘š{\cal L}_{g}=\sqrt{-g}\frac{1}{2\kappa^{2}}\bigg{[}\phi R(\phi)-W(\phi)\bigg{]% }+{\cal L}_{m}caligraphic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ italic_Ļ• italic_R ( italic_Ļ• ) - italic_W ( italic_Ļ• ) ] + caligraphic_L start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT (8)

where

Ļ•italic-Ļ•\displaystyle\phiitalic_Ļ• ā‰”fā€²ā¢(R)absentsuperscriptš‘“ā€²š‘…\displaystyle\equiv f^{\prime}(R)ā‰” italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ( italic_R ) (9)
Wā¢(Ļ•)š‘Šitalic-Ļ•\displaystyle W(\phi)italic_W ( italic_Ļ• ) ā‰”Ļ•ā¢Rā¢(Ļ•)āˆ’fā¢[Rā¢(Ļ•)]absentitalic-Ļ•š‘…italic-Ļ•š‘“delimited-[]š‘…italic-Ļ•\displaystyle\equiv\phi R(\phi)-f[R(\phi)]ā‰” italic_Ļ• italic_R ( italic_Ļ• ) - italic_f [ italic_R ( italic_Ļ• ) ] (10)

which is simply a Legendre Transform, that replaces Rš‘…Ritalic_R by Ļ•italic-Ļ•\phiitalic_Ļ• as the independent variable. For that, we require a definite sign for fā€²ā€²ā‰”d2ā¢f/dā¢R2superscriptš‘“ā€²ā€²superscriptš‘‘2š‘“š‘‘superscriptš‘…2f^{\prime\prime}\equiv d^{2}f/dR^{2}italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT ā‰” italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_f / italic_d italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (usually taken positive), except, perhaps, at particular values of Rš‘…Ritalic_R.

We now perform a conformal transformation where the conformal factor is Ļ•italic-Ļ•\phiitalic_Ļ• itself:

gĪ¼ā¢Ī½ā†’g~Ī¼ā¢Ī½ā‰”Ļ•ā‹…gĪ¼ā¢Ī½,ā†’subscriptš‘”šœ‡šœˆsubscript~š‘”šœ‡šœˆā‹…italic-Ļ•subscriptš‘”šœ‡šœˆg_{\mu\nu}\rightarrow\tilde{g}_{\mu\nu}\equiv\phi\cdot g_{\mu\nu},italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT ā†’ over~ start_ARG italic_g end_ARG start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT ā‰” italic_Ļ• ā‹… italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT , (11)

and define a new field

Ļ•~ā‰”+32ā¢Īŗ2ā¢lnā”Ļ•,~italic-Ļ•32superscriptšœ…2italic-Ļ•\tilde{\phi}\equiv+\sqrt{\frac{3}{2\kappa^{2}}}\ln\phi,over~ start_ARG italic_Ļ• end_ARG ā‰” + square-root start_ARG divide start_ARG 3 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG roman_ln italic_Ļ• , (12)

which requires that fā€²ā‰”Ļ•>0superscriptš‘“ā€²italic-Ļ•0f^{\prime}\equiv\phi>0italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ā‰” italic_Ļ• > 0, already assumed when we performed the conformal transformation (11). Note that, since the field Ļ•italic-Ļ•\phiitalic_Ļ• is dimensionless and Īŗā‰”1/Mpā¢lšœ…1subscriptš‘€š‘š‘™\kappa\equiv 1/M_{pl}italic_Īŗ ā‰” 1 / italic_M start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT, then Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG has the dimensions of Mpā¢lsubscriptš‘€š‘š‘™M_{pl}italic_M start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT.

Upon the definitions above, Eq.Ā (2) is written as

ā„’~g=āˆ’g~ā¢[12ā¢Īŗ2ā¢R~āˆ’12ā¢g~Ī±ā¢Ī²ā¢āˆ‚Ī±Ļ•~ā¢āˆ‚Ī²Ļ•~āˆ’Vā¢(Ļ•~)]+ā„’mā¢[Ļˆ,1Ļ•ā¢g~Ī¼ā¢Ī½]subscript~ā„’š‘”~š‘”delimited-[]12superscriptšœ…2~š‘…12superscript~š‘”š›¼š›½subscriptš›¼~italic-Ļ•subscriptš›½~italic-Ļ•š‘‰~italic-Ļ•subscriptā„’š‘ššœ“1italic-Ļ•subscript~š‘”šœ‡šœˆ\tilde{\cal L}_{g}=\sqrt{-\tilde{g}}\bigg{[}\frac{1}{2\kappa^{2}}\tilde{R}-% \frac{1}{2}\tilde{g}^{\alpha\beta}\partial_{\alpha}\tilde{\phi}\partial_{\beta% }\tilde{\phi}-V(\tilde{\phi})\bigg{]}+{\cal L}_{m}\bigg{[}\psi,\frac{1}{\phi}% \tilde{g}_{\mu\nu}\bigg{]}over~ start_ARG caligraphic_L end_ARG start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT = square-root start_ARG - over~ start_ARG italic_g end_ARG end_ARG [ divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG over~ start_ARG italic_R end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG over~ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_Ī± italic_Ī² end_POSTSUPERSCRIPT āˆ‚ start_POSTSUBSCRIPT italic_Ī± end_POSTSUBSCRIPT over~ start_ARG italic_Ļ• end_ARG āˆ‚ start_POSTSUBSCRIPT italic_Ī² end_POSTSUBSCRIPT over~ start_ARG italic_Ļ• end_ARG - italic_V ( over~ start_ARG italic_Ļ• end_ARG ) ] + caligraphic_L start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_Ļˆ , divide start_ARG 1 end_ARG start_ARG italic_Ļ• end_ARG over~ start_ARG italic_g end_ARG start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT ] (13)

where

Vā¢(Ļ•~)ā‰”Rā¢(Ļ•)ā¢Ļ•āˆ’fā¢[Rā¢(Ļ•)]2ā¢Īŗ2ā¢Ļ•2š‘‰~italic-Ļ•š‘…italic-Ļ•italic-Ļ•š‘“delimited-[]š‘…italic-Ļ•2superscriptšœ…2superscriptitalic-Ļ•2V(\tilde{\phi})\equiv\frac{R(\phi)\phi-f[R(\phi)]}{2\kappa^{2}\phi^{2}}italic_V ( over~ start_ARG italic_Ļ• end_ARG ) ā‰” divide start_ARG italic_R ( italic_Ļ• ) italic_Ļ• - italic_f [ italic_R ( italic_Ļ• ) ] end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_Ļ• start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (14)

and the determinant g~~š‘”\tilde{g}over~ start_ARG italic_g end_ARG and the Ricci scalar R~~š‘…\tilde{R}over~ start_ARG italic_R end_ARG are calculated from the new metric g~Ī¼ā¢Ī½subscript~š‘”šœ‡šœˆ\tilde{g}_{\mu\nu}over~ start_ARG italic_g end_ARG start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT.

The Lagrangian (13) is the one of GR with an extra scalar field Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG, which carries the extra degree of freedom, as previously announced. One should notice that the field Ļˆšœ“\psiitalic_Ļˆ (that represents matter/radiation) does not follow geodesics in the new frame, although its equation of motion should be expressed in terms of the new metric. Accordingly, it is not conserved if the covariant derivative is calculated in terms of g~Ī¼ā¢Ī½subscript~š‘”šœ‡šœˆ\tilde{g}_{\mu\nu}over~ start_ARG italic_g end_ARG start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT, since this is supposed to hold in the previous metric only. One can also notice that Ļˆšœ“\psiitalic_Ļˆ is not minimally coupled to the metric g~Ī¼ā¢Ī½subscript~š‘”šœ‡šœˆ\tilde{g}_{\mu\nu}over~ start_ARG italic_g end_ARG start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT since Ļ•italic-Ļ•\phiitalic_Ļ• shows up

This is the so-called Einstein frame. Its ā€œequivalenceā€ to the Jordan frame (before the conformal transformation) is a topic of debate in the current literature, along with the question about which one is the ā€œphysicalā€ one (whatever that means). Here, we will rely on the conservation of the energy-momentum tensor of matter/radiation in the Jordan frame in order to define that all the observations are supposed to be done in that frame. One can, obviously, use the Einstein frame to make calculations, but, at the end of the day, one should go to the Jordan frame in order to make predictions on observable quantities. With that in mind, we point out that eventual divergences in the EF are not always bad, if they are mapped into finite quantities in the Jordan frame.

One final piece of information: the aforementioned ā€œequivalenceā€ between the frames is robust when it comes to inflation. Indeed, it is possible to show Kaiser (1995) that if the slow-roll parameters hold in one of the frames, then they also hold in the other one.

I.4 Constraints

One can always write Eq.Ā (2) as

ā„’f=āˆ’gā¢12ā¢Īŗ2ā¢fā¢(R)=āˆ’gā¢12ā¢Īŗ2ā¢(R+Ļµā¢Ī”ā¢(R))=āˆ’gā¢12ā¢Īŗ2ā¢Rā¢(1+Ļµā¢Ī”ā¢(R)R),subscriptā„’š‘“š‘”12superscriptšœ…2š‘“š‘…š‘”12superscriptšœ…2š‘…italic-ĻµĪ”š‘…š‘”12superscriptšœ…2š‘…1italic-ĻµĪ”š‘…š‘…{\cal L}_{f}=\sqrt{-g}\frac{1}{2\kappa^{2}}f(R)=\sqrt{-g}\frac{1}{2\kappa^{2}}% \bigg{(}R+\epsilon\Delta(R)\bigg{)}=\sqrt{-g}\frac{1}{2\kappa^{2}}R\bigg{(}1+% \epsilon\frac{\Delta(R)}{R}\bigg{)},caligraphic_L start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_f ( italic_R ) = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( italic_R + italic_Ļµ roman_Ī” ( italic_R ) ) = square-root start_ARG - italic_g end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_R ( 1 + italic_Ļµ divide start_ARG roman_Ī” ( italic_R ) end_ARG start_ARG italic_R end_ARG ) , (15)

which defines Ī”ā¢(R)Ī”š‘…\Delta(R)roman_Ī” ( italic_R ) and Ļµitalic-Ļµ\epsilonitalic_Ļµ, where the latter indicates the amplitude of the deviation Ī”ā¢(R)Ī”š‘…\Delta(R)roman_Ī” ( italic_R ) from the Einstein-Hilbert Lagrangian ā„’Eā¢Hā‰”Rsubscriptā„’šøš»š‘…{\cal L}_{EH}\equiv Rcaligraphic_L start_POSTSUBSCRIPT italic_E italic_H end_POSTSUBSCRIPT ā‰” italic_R. The standard constraints on fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) theories are almost always listed then as PogosianĀ andĀ Silvestri (2008)

limRā†’āˆžĪ”R=0subscriptā†’š‘…Ī”š‘…0\displaystyle\lim_{R\to\infty}\frac{\Delta}{R}=0roman_lim start_POSTSUBSCRIPT italic_R ā†’ āˆž end_POSTSUBSCRIPT divide start_ARG roman_Ī” end_ARG start_ARG italic_R end_ARG = 0 (16)
Ļ•ā‰”fā€²ā‰”1+Ī”ā€²>0italic-Ļ•superscriptš‘“ā€²1superscriptĪ”ā€²0\displaystyle\phi\equiv f^{\prime}\equiv 1+\Delta^{\prime}>0italic_Ļ• ā‰” italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ā‰” 1 + roman_Ī” start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT > 0 (17)
fā€²ā€²ā‰”Ī”ā€²ā€²>0.superscriptš‘“ā€²ā€²superscriptĪ”ā€²ā€²0\displaystyle f^{\prime\prime}\equiv\Delta^{\prime\prime}>0.italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT ā‰” roman_Ī” start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT > 0 . (18)

Nevertheless, they should be taken with a grain of salt. Let us look at them, one by one, and see why.

The first one states that the deviations from GR should be small when the Ricci scalar is large, i.e, in the primordial universe. That relies on the fact that we indeed understand the evolution of the universe from nucleosynthesis on as one of the many success of GR. Nevertheless, this claim is poorly stated. The reason is threefold. First, one aims exactly at non-negligible deviations from GR when modified theories are applied to explain the primordial inflationary era. Secondly, even when the extra terms Ī”ā¢(R)Ī”š‘…\Delta(R)roman_Ī” ( italic_R ) are negligible, the equations of motion (EoM) are still higher-order differential equations. As such, there is no guarantee that their solutions are going to be closer to the ones from GR. Actually, it can be easily seen that this is not the case if one adds higher-order terms at the EoM of linear harmonic oscillator. (except with strong fine tuning of the initial conditions). An even simpler case (a homework) is to compare solutions of the equation Ļµā¢xĀØ+xĖ™=0italic-ĻµĀØš‘„Ė™š‘„0\epsilon\ddot{x}+\dot{x}=0italic_Ļµ overĀØ start_ARG italic_x end_ARG + overĖ™ start_ARG italic_x end_ARG = 0 in the limit Ļµā†’0ā†’italic-Ļµ0\epsilon\to 0italic_Ļµ ā†’ 0 to the ones from the ā€œoriginalā€ equation xĖ™=0Ė™š‘„0\dot{x}=0overĖ™ start_ARG italic_x end_ARG = 0. Obviously, in modified theories of gravity, the nonlinearities only make it worse. Finally, the Ricci scalar is not necessarily large in the early universe ā€” we recall the reader that the trace of Eq.Ā (4) is now

3ā¢ā–”ā¢fRā¢(R)+Rā¢fRā¢(R)āˆ’2ā¢fā¢(R)=Īŗ2ā¢T3ā–”subscriptš‘“š‘…š‘…š‘…subscriptš‘“š‘…š‘…2š‘“š‘…superscriptšœ…2š‘‡3\square f_{R}(R)+Rf_{R}(R)-2f(R)=\kappa^{2}T3 ā–” italic_f start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ( italic_R ) + italic_R italic_f start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ( italic_R ) - 2 italic_f ( italic_R ) = italic_Īŗ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T (19)

(where Tš‘‡Titalic_T is the trace of TĪ¼ā¢Ī½subscriptš‘‡šœ‡šœˆT_{\mu\nu}italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT) which is not an algebraic relation anymore (as in GR). Therefore, one could still have a large Tš‘‡Titalic_T (as it is expected in the early universe), but a small Rš‘…Ritalic_R.

The second constraint usually implies that it is required for an attractive gravity. It relies on dividing Eq.Ā (4) by fā€²superscriptš‘“ā€²f^{\prime}italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT and noticing that it would change the sign of the Gravitational coupling GšŗGitalic_G. Nevertheless, that mapping would disrupt the very conservation of TĪ¼ā¢Ī½subscriptš‘‡šœ‡šœˆT_{\mu\nu}italic_T start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT and we will not follow it here. However, it is indeed necessary to perform the conformal transformation to the Einstein frame (which is desired, but it is far from being an actual constraint).

Finally, the third constraint is well-posed. It can be shown that perturbations of the flat space (Minkowski) are stable if and only if fā€²ā€²>0superscriptš‘“ā€²ā€²0f^{\prime\prime}>0italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT > 0. Here, we follow Ref.Ā SotiriouĀ andĀ Faraoni (2010). Let us expand Eq.Ā (19) around Minkwoski metric Ī·Ī¼ā¢Ī½subscriptšœ‚šœ‡šœˆ\eta_{\mu\nu}italic_Ī· start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT, i.e,

gĪ¼ā¢Ī½subscriptš‘”šœ‡šœˆ\displaystyle g_{\mu\nu}italic_g start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT =Ī·Ī¼ā¢Ī½+hĪ¼ā¢Ī½absentsubscriptšœ‚šœ‡šœˆsubscriptā„Žšœ‡šœˆ\displaystyle=\eta_{\mu\nu}+h_{\mu\nu}= italic_Ī· start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT + italic_h start_POSTSUBSCRIPT italic_Ī¼ italic_Ī½ end_POSTSUBSCRIPT (20)
Rš‘…\displaystyle Ritalic_R =T+Ī“ā¢R,absentš‘‡š›æš‘…\displaystyle=T+\delta R,= italic_T + italic_Ī“ italic_R , (21)

and assume Eq.Ā (15), which yields

Ī“ā¢RĀØ+12ā¢Ī”ā€²ā€²ā¢(1Ļµāˆ’Ī”ā€²)ā¢Ī“ā¢Rā‰ˆ0.ĀØš›æš‘…12superscriptĪ”ā€²ā€²1italic-ĻµsuperscriptĪ”ā€²š›æš‘…0\ddot{\delta R}+\frac{1}{2\Delta^{\prime\prime}}\bigg{(}\frac{1}{\epsilon}-% \Delta^{\prime}\bigg{)}\delta R\approx 0.overĀØ start_ARG italic_Ī“ italic_R end_ARG + divide start_ARG 1 end_ARG start_ARG 2 roman_Ī” start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT end_ARG ( divide start_ARG 1 end_ARG start_ARG italic_Ļµ end_ARG - roman_Ī” start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ) italic_Ī“ italic_R ā‰ˆ 0 . (22)

Since Ļµitalic-Ļµ\epsilonitalic_Ļµ can be as small as needed, the term inside the parenthesis is always positive. The sign of the effective mass squared for Ī“ā¢Rš›æš‘…\delta Ritalic_Ī“ italic_R is the same of Ī”ā¢"Ī”"\Delta"roman_Ī” ", which should be positive for stable perturbations. We will see further below that the same constraint can be obtained from different and independent calculations.

I.5 Practical Applications

I.5.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 ā€” followed by the current accelerated Dark-Energy Dominated (DED) phase. One could also add the inflationary era at the very beginning, but let us postpone that discussion for now, because it involves the (p)reheating process.

Using the language of dynamical systems, we should then require that the RD phase is not an attractor (otherwise, the universe would never have left it). Analogously, the MD phase has to be a saddle point (i.e, with at least one stable direction and another unstable one), so that the universe will be drawn towards it (after RD, along the stable direction) but move on to the DED phase (along the unstable direction). Usually, one requires that the last phase is an attractor (i.e, only negative eigenvalues) so that the initial conditions do not have to be fine tuned.

Using that approach, Ref.Ā AmendolaĀ etĀ al. (2007) has determined the necessary conditions for any fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) so that the picture described above is fulfilled. The authors define the quantities

mš‘š\displaystyle mitalic_m ā‰”Rā¢fā€²ā€²fā€²absentš‘…superscriptš‘“ā€²ā€²superscriptš‘“ā€²\displaystyle\equiv\frac{Rf^{\prime\prime}}{f^{\prime}}ā‰” divide start_ARG italic_R italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT end_ARG start_ARG italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT end_ARG (23)
rš‘Ÿ\displaystyle ritalic_r ā‰”āˆ’Rā¢fā€²f,absentš‘…superscriptš‘“ā€²š‘“\displaystyle\equiv-\frac{Rf^{\prime}}{f},ā‰” - divide start_ARG italic_R italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT end_ARG start_ARG italic_f end_ARG , (24)

which measure the deviation from GR ā€” for which m=0š‘š0m=0italic_m = 0 and r=āˆ’1š‘Ÿ1r=-1italic_r = - 1.

The existence of a MD phase with the required constraints is translated to

{m(rā‰ˆāˆ’1)ā‰ˆ0+Ā anddā¢mdā¢r|rā‰ˆāˆ’1>āˆ’1Ā ,\left\{\begin{tabular}[]{l}$m(r\approx-1)\approx 0^{+}$ ~{} and\\ $\frac{dm}{dr}\big{|}_{r\approx-1}>-1$ ,\end{tabular}\right.{ start_ROW start_CELL italic_m ( italic_r ā‰ˆ - 1 ) ā‰ˆ 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and end_CELL end_ROW start_ROW start_CELL divide start_ARG italic_d italic_m end_ARG start_ARG italic_d italic_r end_ARG | start_POSTSUBSCRIPT italic_r ā‰ˆ - 1 end_POSTSUBSCRIPT > - 1 , end_CELL end_ROW (25)

while a de Sitter attractor requires either

{m=āˆ’rāˆ’1,3āˆ’12<m<1,dā¢mdā¢r|rā‰ˆāˆ’2<āˆ’1\displaystyle\left\{\begin{tabular}[]{l}$m=-r-1$,\\ $\frac{\sqrt{3}-1}{2}<m<1$,\\ $\frac{dm}{dr}\big{|}_{r\approx-2}<-1$\end{tabular}\right.{ start_ROW start_CELL italic_m = - italic_r - 1 , end_CELL end_ROW start_ROW start_CELL divide start_ARG square-root start_ARG 3 end_ARG - 1 end_ARG start_ARG 2 end_ARG < italic_m < 1 , end_CELL end_ROW start_ROW start_CELL divide start_ARG italic_d italic_m end_ARG start_ARG italic_d italic_r end_ARG | start_POSTSUBSCRIPT italic_r ā‰ˆ - 2 end_POSTSUBSCRIPT < - 1 end_CELL end_ROW (29)
oror\displaystyle{\rm or}roman_or (30)
0<mā¢(r=āˆ’1)ā‰¤1.0š‘šš‘Ÿ11\displaystyle 0<m(r=-1)\leq 1.0 < italic_m ( italic_r = - 1 ) ā‰¤ 1 . (31)

That paper has shown that simple expressions for fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ), such as simple power laws just will not work.

The growth rate of perturbations can be affected, of course, by modifications on gravity. There is a couple of mainstream approached for this problem. One may, for instance, realize that there are two opposite regimes, known as weak- and strong-field limits. They correspond to situations where the chameleon effect (see next section) is in action or not, which can be modeled by an effective Newtonā€™s constant equal to GNsubscriptšŗš‘G_{N}italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT (its standard value) or to 43ā¢GN43subscriptšŗš‘\frac{4}{3}G_{N}divide start_ARG 4 end_ARG start_ARG 3 end_ARG italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT HerreraĀ etĀ al. (2017).

A second approach is to model the value of GNsubscriptšŗš‘G_{N}italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT as a function of wavenumber kš‘˜kitalic_k BertschingerĀ andĀ Zukin (2008). Being dependent on scale, gravity will produce different effects at different distances from a central mass and, then, the so-called shell crossing BorisovĀ etĀ al. (2012) ā€” a top-hat overdensity will then not retain its shape.

The final goal in both approaches is to obtain the number of observed objects in a given mass range, as a function of its redshift. That number can then be used to constraint the free parameters of your chosen model.

A short review can be found in Ref.Ā Batista (2021).

I.5.2 Relativistic Stars

I briefly mention that, as a serious candidate for a full theory of gravity, a given fā¢(R)š‘“š‘…f(R)italic_f ( italic_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š‘Ÿritalic_r; in the latter, the time tš‘”titalic_t. Both cover a large range of values of densities (from the core to the radial infinity, in stars) and have plenty of data available.

The catch is the uncertainty with respect to the equation of state of the standard (sometimes not much so) matter that makes up the star. Using the aforementioned ā€œcurvature fluidā€ analogy, there is clearly some degeneracy between the contribution from such a term and from a different matter composition. Nevertheless, one may still provide constraints on particular models, testing both the maximum mass (and minimal radius) and the stability of such configurations (usually against radial perturbations). Indeed, the latter are able to yield the strongest constraints on the free parameters of your theory 2020JCAP...11..048P; 2022JCAP...09..058P; PretelĀ etĀ al. (2023).

One caveat is the many non-equivalent definitions of mass SbisĆ Ā etĀ al. (2020): the astrophysical ā€” the one measured by an observer at radial infinity, that can be calculated from an asymptotic expansion of the metric components, given by

Mastroā‰”limrā†’āˆžr2ā¢GNā¢(1āˆ’1grā¢r),subscriptš‘€astrosubscriptā†’š‘Ÿš‘Ÿ2subscriptšŗš‘11subscriptš‘”š‘Ÿš‘ŸM_{\rm astro}\equiv\lim_{r\to\infty}\frac{r}{2G_{N}}\bigg{(}1-\frac{1}{g_{rr}}% \bigg{)},italic_M start_POSTSUBSCRIPT roman_astro end_POSTSUBSCRIPT ā‰” roman_lim start_POSTSUBSCRIPT italic_r ā†’ āˆž end_POSTSUBSCRIPT divide start_ARG italic_r end_ARG start_ARG 2 italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT end_ARG ( 1 - divide start_ARG 1 end_ARG start_ARG italic_g start_POSTSUBSCRIPT italic_r italic_r end_POSTSUBSCRIPT end_ARG ) , (32)

or the plain integration of the energy density (with or without the proper-volume factor āˆ’gš‘”\sqrt{-g}square-root start_ARG - italic_g end_ARG) bounded by the starā€™s surface at r=rāˆ—š‘Ÿsubscriptš‘Ÿr=r_{*}italic_r = italic_r start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT

MĻsubscriptš‘€šœŒ\displaystyle M_{\rho}italic_M start_POSTSUBSCRIPT italic_Ļ end_POSTSUBSCRIPT ā‰”4ā¢Ļ€ā¢āˆ«0rāˆ—Ļā¢(r)ā¢r2ā¢š‘‘rabsent4šœ‹superscriptsubscript0subscriptš‘ŸšœŒš‘Ÿsuperscriptš‘Ÿ2differential-dš‘Ÿ\displaystyle\equiv 4\pi\int_{0}^{r_{*}}\rho(r)r^{2}\,drā‰” 4 italic_Ļ€ āˆ« start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_Ļ ( italic_r ) italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_r (33)
Mpropsubscriptš‘€prop\displaystyle M_{\rm prop}italic_M start_POSTSUBSCRIPT roman_prop end_POSTSUBSCRIPT ā‰”4ā¢Ļ€ā¢āˆ«0rāˆ—Ļā¢(r)ā¢r2ā¢āˆ’grā¢rā¢š‘‘rabsent4šœ‹superscriptsubscript0subscriptš‘ŸšœŒš‘Ÿsuperscriptš‘Ÿ2subscriptš‘”š‘Ÿš‘Ÿdifferential-dš‘Ÿ\displaystyle\equiv 4\pi\int_{0}^{r_{*}}\rho(r)r^{2}\,\sqrt{-g_{rr}}\,drā‰” 4 italic_Ļ€ āˆ« start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_Ļ ( italic_r ) italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG - italic_g start_POSTSUBSCRIPT italic_r italic_r end_POSTSUBSCRIPT end_ARG italic_d italic_r (34)

and the ā€œsurface-redshift massā€, estimated from the light redshift zāˆ—subscriptš‘§z_{*}italic_z start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT (emitted from the surface to infinity) as

Msrā‰”rāˆ—2ā¢GNā¢zāˆ—ā¢(2+zāˆ—)(1+zāˆ—)2.subscriptš‘€srsubscriptš‘Ÿ2subscriptšŗš‘subscriptš‘§2subscriptš‘§superscript1subscriptš‘§2M_{\rm sr}\equiv\frac{r_{*}}{2G_{N}}\frac{z_{*}(2+z_{*})}{(1+z_{*})^{2}}.italic_M start_POSTSUBSCRIPT roman_sr end_POSTSUBSCRIPT ā‰” divide start_ARG italic_r start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_G start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT end_ARG divide start_ARG italic_z start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT ( 2 + italic_z start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT ) end_ARG start_ARG ( 1 + italic_z start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (35)

Although all of them coincide in GR (except, obviously for the proper mass, whose difference to the others yields the binding energy), they are not the same in fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ), due to the simple fact that the scalar degree of freedom does leaks out from the star. One should then be careful to properly define the star surface: does it happen when the baryonic pressure vanishes or when the total fluid pressure (including then the ā€œcurvature fluidā€) does?

I.5.3 Local gravity constraints

No review (however short) on fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) would be complete without mentioning (however briefly) the chameleon effect KhouryĀ andĀ Weltman (2004). 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 the scalar field depend on the energy density of the standard matter, due to the aforementioned coupling between them in the Einstein Frame. The larger the average local density, the heavier the field is and the shorter is its range. Since we are surrounded by a dark-matter halo, all the classical solar-system tests will still hold.

I.6 Summary

We finish this section by stating how long we have gone so far: one must come up with an ingenious function fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) ā€” that could either be completely ad hoc (aiming then to pick important characteristics and their outcome to observable quantities) or from first principles (e.g, some low-energy limit of a ā€œquantum gravityā€ theory). Such function must follow theoretical (self-consistency, non-divergences) and observational constraints (background and perturbed cosmology, local gravity and solar system). A very useful tool is the conformal transformation to the Einstein frame, where the Lagrangian is similar to GR and, therefore, more palatable.

In the next Section we will open a large parenthesis to review some basic Thermodynamics, but in a not-so-ordinary approach, given by Catastrophe Theory. The final Section will then make the link between the first two.

II Thermodynamics

In the following subsections we will use the van der Waals gas as a typical description of a system that goes through a first-order phase transition. Iin spite of its analytical simplicity, it bears the fundamental characteristics of such transformation and, at the same time, allows for full analytical calculations.

II.1 Phase transitions

The equation of state for a non-ideal gas was introduced by van der Waals in 1873 vanĀ der Waals (ates):

(P+aV2)ā¢(Vāˆ’b)=Rā¢Tš‘ƒš‘Žsuperscriptš‘‰2š‘‰š‘š‘…š‘‡\left(P+\frac{a}{V^{2}}\right)\left(V-b\right)=R\,T( italic_P + divide start_ARG italic_a end_ARG start_ARG italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ( italic_V - italic_b ) = italic_R italic_T (36)

where the constants aš‘Žaitalic_a and bš‘bitalic_b are introduced to take into account, respectively, a residual interaction between the gas molecules (therefore, decreasing the effective pressure measured) and the volume taken by them. Three of such curves, for different temperatures, are plotted in Fig.Ā 1. One can see that it resembles the ideal gas behavior for large temperatures Tš‘‡Titalic_T, as expected. For low temperatures, though, there are two distinguishable regimes: for larger volumes, the pressure is low, while the opposite happens for small volumes. The former region describes a gas-like behavior, while the latter, where the pressure increases fast with a slight decrease in volume, resembles a liquid state, where the interactions between the particles can no longer be neglected.

Refer to caption
Figure 1: Plot of the pressure Ɨ\timesƗ volume for 3 temperatures; the lower (upper) curve corresponds to temperatures below (above) the critical temperature Tcsubscriptš‘‡š‘T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. In the middle curve (T=Tcš‘‡subscriptš‘‡š‘T=T_{c}italic_T = italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT), the two extrema (where dā¢P/dā¢V=0š‘‘š‘ƒš‘‘š‘‰0dP/dV=0italic_d italic_P / italic_d italic_V = 0) coalesce at the inflection point (where d2ā¢P/dā¢V2=0superscriptš‘‘2š‘ƒš‘‘superscriptš‘‰20d^{2}P/dV^{2}=0italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_P / italic_d italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0).

We recall the reader that, as any equation of state in Thermodynamics, it describes equilibrium configurations, but not necessarily stable ones. For instance, for temperatures below a critical value Tcsubscriptš‘‡š‘T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, there is a volume range where

dā¢Pdā¢V>0.š‘‘š‘ƒš‘‘š‘‰0\frac{dP}{dV}>0.divide start_ARG italic_d italic_P end_ARG start_ARG italic_d italic_V end_ARG > 0 . (37)

That indicates an unstable equilibrium, since, is the gas was compressed, i.e, its volume would be decreased by an external force, then its pressure would also decrease. In other words, the gas would not try to restore the previous configuration and ultimately it would collapse222Of course, the opposite argument applies for an expanding gas in this volume range..

A very important piece of information is obtained when one writes the Gibbs Energy Gā¢(P,T)šŗš‘ƒš‘‡G(P,T)italic_G ( italic_P , italic_T ), whose variation is given by

dā¢Gā¢(P,T)=Vā‹…dā¢Pāˆ’Sā‹…dā¢T.š‘‘šŗš‘ƒš‘‡ā‹…š‘‰š‘‘š‘ƒā‹…š‘†š‘‘š‘‡dG(P,T)=V\cdot dP-S\cdot dT.italic_d italic_G ( italic_P , italic_T ) = italic_V ā‹… italic_d italic_P - italic_S ā‹… italic_d italic_T . (38)

Therefore,

dā¢Gdā¢Pš‘‘šŗš‘‘š‘ƒ\displaystyle\frac{dG}{dP}divide start_ARG italic_d italic_G end_ARG start_ARG italic_d italic_P end_ARG =Vandabsentš‘‰and\displaystyle=V\qquad{\rm and}= italic_V roman_and (39)
d2ā¢Gdā¢P2superscriptš‘‘2šŗš‘‘superscriptš‘ƒ2\displaystyle\frac{d^{2}G}{dP^{2}}divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_G end_ARG start_ARG italic_d italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG =dā¢Vdā¢P.absentš‘‘š‘‰š‘‘š‘ƒ\displaystyle=\frac{dV}{dP}.= divide start_ARG italic_d italic_V end_ARG start_ARG italic_d italic_P end_ARG . (40)

This means that d2ā¢G/dā¢P2>0superscriptš‘‘2šŗš‘‘superscriptš‘ƒ20d^{2}G/dP^{2}>0italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_G / italic_d italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT > 0 in the unstable region.

One can always write Eq.Ā (36) as

V3āˆ’(b+Rā¢TP)ā¢V2āˆ’(aP)ā¢Vāˆ’aā¢bP=0.superscriptš‘‰3š‘š‘…š‘‡š‘ƒsuperscriptš‘‰2š‘Žš‘ƒš‘‰š‘Žš‘š‘ƒ0V^{3}-\left(b+\frac{RT}{P}\right)V^{2}-\left(\frac{a}{P}\right)V-\frac{ab}{P}=0.italic_V start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - ( italic_b + divide start_ARG italic_R italic_T end_ARG start_ARG italic_P end_ARG ) italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( divide start_ARG italic_a end_ARG start_ARG italic_P end_ARG ) italic_V - divide start_ARG italic_a italic_b end_ARG start_ARG italic_P end_ARG = 0 . (41)

As a cubic equation with real coefficients, there is always (i.e, for any Tš‘‡Titalic_T) at least one real root, i.e, one pressure Pš‘ƒPitalic_P and one temperature Tš‘‡Titalic_T for a given volume Vš‘‰Vitalic_V. On the other hand, for T<Tcā‰”8ā¢a/(27ā¢bā¢R)š‘‡subscriptš‘‡š‘8š‘Ž27š‘š‘…T<T_{c}\equiv 8a/(27bR)italic_T < italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰” 8 italic_a / ( 27 italic_b italic_R ), Eq.Ā (41) yields three real roots: three values for Vš‘‰Vitalic_V (one of them in the aforementioned unstable region and, therefore, non physical) for a a given pressure Pš‘ƒPitalic_P and temperature T<Tcš‘‡subscriptš‘‡š‘T<T_{c}italic_T < italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. At the critical temperature Tā¢cš‘‡š‘Tcitalic_T italic_c, the system has critical pressure Pcā‰”a/(27ā¢b2)subscriptš‘ƒš‘š‘Ž27superscriptš‘2P_{c}\equiv a/(27b^{2})italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰” italic_a / ( 27 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and critical Volume Vcā‰”3ā¢bsubscriptš‘‰š‘3š‘V_{c}\equiv 3bitalic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰” 3 italic_b. Above Tcsubscriptš‘‡š‘T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, there is no difference between the states Aš“Aitalic_A and BšµBitalic_B, i.e, one can go from the former to the latter without going through first-order phase transition.

A simple variable change is usually made so that we can deal only with dimensionless quantities:

pš‘\displaystyle pitalic_p ā‰”P/Pcāˆ’1absentš‘ƒsubscriptš‘ƒš‘1\displaystyle\equiv P/P_{c}-1ā‰” italic_P / italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - 1 (42)
tš‘”\displaystyle titalic_t ā‰”T/Tcāˆ’1absentš‘‡subscriptš‘‡š‘1\displaystyle\equiv T/T_{c}-1ā‰” italic_T / italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - 1 (43)
xš‘„\displaystyle xitalic_x ā‰”Vc/Vāˆ’1,absentsubscriptš‘‰š‘š‘‰1\displaystyle\equiv V_{c}/V-1,ā‰” italic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT / italic_V - 1 , (44)

upon which Eq.Ā (41) is written as

x3superscriptš‘„3\displaystyle x^{3}italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT +Ī±ā¢x+Ī²=0whereš›¼š‘„š›½0where\displaystyle+\alpha x+\beta=0\qquad{\rm where}+ italic_Ī± italic_x + italic_Ī² = 0 roman_where (45)
Ī±š›¼\displaystyle\alphaitalic_Ī± ā‰”8ā¢t+p3absent8š‘”š‘3\displaystyle\equiv\frac{8t+p}{3}ā‰” divide start_ARG 8 italic_t + italic_p end_ARG start_ARG 3 end_ARG (46)
Ī²š›½\displaystyle\betaitalic_Ī² ā‰”8ā¢tāˆ’2ā¢p3.absent8š‘”2š‘3\displaystyle\equiv\frac{8t-2p}{3}.ā‰” divide start_ARG 8 italic_t - 2 italic_p end_ARG start_ARG 3 end_ARG . (47)

The plot of xš‘„xitalic_x, i.e, the equilibrium configurations as a function of the parameters {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² } can be seen in Fig.Ā 51. The region with multiple (to be more precise, three) solutions can be easily seen and its boundary will be determined later on. For now, it suffices to remember that one of the solutions is the unstable one. Outside that region, there is one and only one equilibrium solution.

II.2 Catastrophe theory

The equilibrium solutions can be seen as the ones that extremize a suitable potential energy333Not to be confused with the volume Vš‘‰Vitalic_V of the vdW gas. Vā¢(x)š‘‰š‘„V(x)italic_V ( italic_x ). A simple and general expression444One can get rid of a possible cubic term by a suitable shift on xš‘„xitalic_x. that does the job is

Vā¢(x;Ī±,Ī²)ā‰”14ā¢x4+Ī±2ā¢x2+Ī²ā¢x,š‘‰š‘„š›¼š›½14superscriptš‘„4š›¼2superscriptš‘„2š›½š‘„V(x;\alpha,\beta)\equiv\frac{1}{4}x^{4}+\frac{\alpha}{2}x^{2}+\beta x,italic_V ( italic_x ; italic_Ī± , italic_Ī² ) ā‰” divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_x start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT + divide start_ARG italic_Ī± end_ARG start_ARG 2 end_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_Ī² italic_x , (48)

where we have made explicit the dependence on the control parameters {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² }. In the region with three extrema, there are three Real solutions to Eq.Ā (41). It is easy to see that one of them is necessarily a local maximum, another one is a local minimum and the last one is the global minimum ā€” of course, the two minima could be degenerated, but we will come back to this point in a second. Outside that region, there is only one minimum. The boundary is defined, then, as the location in the parameter space where the local maximum coalesces to the local minimum. Mathematically, it is defined by the so-called Singularity Set:

d2ā¢Vdā¢x2|xāˆ—=0evaluated-atsuperscriptš‘‘2š‘‰š‘‘superscriptš‘„2subscriptš‘„0\frac{d^{2}V}{dx^{2}}\bigg{|}_{x_{*}}=0divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_V end_ARG start_ARG italic_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG | start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 0 (49)

where xāˆ—subscriptš‘„x_{*}italic_x start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT is one of the extrema, given, of course, by

dā¢Vdā¢x|xāˆ—=0,evaluated-atš‘‘š‘‰š‘‘š‘„subscriptš‘„0\frac{dV}{dx}\bigg{|}_{x_{*}}=0,divide start_ARG italic_d italic_V end_ARG start_ARG italic_d italic_x end_ARG | start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 0 , (50)

which is the so-called Equilibrium Set. The solution of Eqs.Ā (49) and (50) is the so-called Bifurcation Set

4ā¢Ī±3+27ā¢Ī²2=0.4superscriptš›¼327superscriptš›½204\alpha^{3}+27\beta^{2}=0.4 italic_Ī± start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 27 italic_Ī² start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 . (51)

In the parameter space {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² }, the function (48) has only one extremum (a minimum) outside the curve (51). Inside, there are 3 extrema, 2 of which are minima and 1, a maximum. On the curve itself, two of them coalesce ā€” which ones is set by the sign of Ī²š›½\betaitalic_Ī².

Refer to caption
Refer to captionRefer to captionRefer to captionRefer to captionRefer to captionRefer to caption
Figure 2: Plot of the fold given by Eq.Ā (51) and the corresponding form of the potential energy (48) in each region. The middle column of insets indicates the form of the potential on the bifurcation set, as indicated by the arrows. Inside the set, the potential has 3 extrema; outside, only one.
Refer to caption
Figure 3: Plot of the equilibrium positions, given by Eq.Ā (50) as a function of the control parameters {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² }. The black solid curve is the projection of the fold (51) onto that surface. The vertical red and green straight lines are arbitrary, except for both being inside the fold. Each one indicate the number (3) of equilibrium solutions for each pair {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² } and which ones coalesce and which one survives when the fold is crossed.

In Fig.Ā 4, one can see the equilibrium surface (i.e, the values of the equilibrium points) as a function of the control parameters {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² }. The fold is projected onto this surface. One can see that, inside the fold, vertical lines (fixed values of the control parameters) cross the surface three times, each for a given equilibrium solution. If we shift the red (green) line towards larger (smaller) values of Ī²š›½\betaitalic_Ī², i.e., toward the closest fold branch, then the largest (smallest) values of xeā¢qsubscriptš‘„š‘’š‘žx_{eq}italic_x start_POSTSUBSCRIPT italic_e italic_q end_POSTSUBSCRIPT coalesce.

Accordingly, we can plot the so-called swallowtail curve also as a function of {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² } ā€” see Fig.Ā 4. one can see that the highest Vā¢(xeā¢q)š‘‰subscriptš‘„š‘’š‘žV(x_{eq})italic_V ( italic_x start_POSTSUBSCRIPT italic_e italic_q end_POSTSUBSCRIPT ) is always unstable, as expected from a general principle of minimizing the energy. Here, as in the previous plot, when we shift the red (green) vertical line towards the closet fold branch, it is always the equilibrium point corresponding to the highest Vā¢(xeā¢q)š‘‰subscriptš‘„š‘’š‘žV(x_{eq})italic_V ( italic_x start_POSTSUBSCRIPT italic_e italic_q end_POSTSUBSCRIPT ) that coalesces with the one corresponding to the intermediate value Vā¢(xeā¢q)š‘‰subscriptš‘„š‘’š‘žV(x_{eq})italic_V ( italic_x start_POSTSUBSCRIPT italic_e italic_q end_POSTSUBSCRIPT ), i.e, the local maximum. Outside the fold, only the global minimum survives.

Refer to caption
Figure 4: Plot of the potential (48), calculated at the equilibrium positions xeā¢qsubscriptš‘„š‘’š‘žx_{eq}italic_x start_POSTSUBSCRIPT italic_e italic_q end_POSTSUBSCRIPT given by Eq.Ā (50), as a function of the control parameters {Ī±,Ī²}š›¼š›½\{\alpha,\beta\}{ italic_Ī± , italic_Ī² }. The black solid curve is the projection of the fold (51) onto that surface. Vertical red and green lines as in the previous figure.

III fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) Theories and Thermodynamics

Here we follow Ref.Ā PeraltaĀ andĀ JorĆ” s (2020). We have described the traditional discussion: one comes up with an ingenious non-linear function fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) (in the Jordan frame) and, via a conformal transformation given by Ļ•ā‰”fā€²ā¢(R)italic-Ļ•superscriptš‘“ā€²š‘…\phi\equiv f^{\prime}(R)italic_Ļ• ā‰” italic_f start_POSTSUPERSCRIPT ā€² end_POSTSUPERSCRIPT ( italic_R ), arrives at the Einstein frame, where the physics is described by GR and an extra scalar field Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG subject to a potential V(Ļ•~V(\tilde{\phi}italic_V ( over~ start_ARG italic_Ļ• end_ARG), completely defined by the function fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) chosen at the beginning.

Here, we take the opposite direction: we start at the Einstein frame with the simplest potential and look for the corresponding fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) in the Jordan frame. We assume

Vā¢(Ļ•~)=12ā¢m~Ļ•2ā¢(Ļ•~āˆ’a)2+Ī›.š‘‰~italic-Ļ•12superscriptsubscript~š‘šitalic-Ļ•2superscript~italic-Ļ•š‘Ž2Ī›V(\tilde{\phi})=\frac{1}{2}\tilde{m}_{\phi}^{2}(\tilde{\phi}-a)^{2}+\Lambda.italic_V ( over~ start_ARG italic_Ļ• end_ARG ) = divide start_ARG 1 end_ARG start_ARG 2 end_ARG over~ start_ARG italic_m end_ARG start_POSTSUBSCRIPT italic_Ļ• end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Ī› . (52)

For now, the extra (constant) parameters aš‘Žaitalic_a and Ī›Ī›\Lambdaroman_Ī› are simply a matter of generalizing the potential Vš‘‰Vitalic_V, while keeping it as simple as possible. One might argue against the introduction of a cosmological constant Ī›Ī›\Lambdaroman_Ī›, since it seems to be incompatible with one of the core motivations for modifying GR. It will be, however, essential to what follows.

Following a previously established procedure MagnanoĀ andĀ Sokołowski (1994), one arrives at the following parametric expressions, where Ī²ā‰”2/3š›½23\beta\equiv\sqrt{2/3}italic_Ī² ā‰” square-root start_ARG 2 / 3 end_ARG:

fā¢(Ļ•~)š‘“~italic-Ļ•\displaystyle f(\tilde{\phi})italic_f ( over~ start_ARG italic_Ļ• end_ARG ) =e2ā¢Ī²ā¢Ļ•~ā¢[2ā¢Vā¢(Ļ•~)+2ā¢Ī²āˆ’1ā¢dā¢Vā¢(Ļ•~)dā¢Ļ•~]andabsentsuperscripte2š›½~italic-Ļ•delimited-[]2š‘‰~italic-Ļ•2superscriptš›½1dš‘‰~italic-Ļ•š‘‘~italic-Ļ•and\displaystyle={\rm e}^{2\beta\tilde{\phi}}\left[2V(\tilde{\phi})+2\beta^{-1}% \frac{{\rm d}V(\tilde{\phi})}{d\tilde{\phi}}\right]\quad{\rm and}= roman_e start_POSTSUPERSCRIPT 2 italic_Ī² over~ start_ARG italic_Ļ• end_ARG end_POSTSUPERSCRIPT [ 2 italic_V ( over~ start_ARG italic_Ļ• end_ARG ) + 2 italic_Ī² start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT divide start_ARG roman_d italic_V ( over~ start_ARG italic_Ļ• end_ARG ) end_ARG start_ARG italic_d over~ start_ARG italic_Ļ• end_ARG end_ARG ] roman_and (53)
Rā¢(Ļ•~)š‘…~italic-Ļ•\displaystyle R(\tilde{\phi})italic_R ( over~ start_ARG italic_Ļ• end_ARG ) =eĪ²ā¢Ļ•~ā¢[4ā¢Vā¢(Ļ•~)+2ā¢Ī²āˆ’1ā¢dā¢Vā¢(Ļ•~)dā¢Ļ•~].absentsuperscripteš›½~italic-Ļ•delimited-[]4š‘‰~italic-Ļ•2superscriptš›½1dš‘‰~italic-Ļ•š‘‘~italic-Ļ•\displaystyle=\,{\rm e}^{\beta\tilde{\phi}}\,\left[4V(\tilde{\phi})+2\beta^{-1% }\frac{{\rm d}V(\tilde{\phi})}{d\tilde{\phi}}\right].= roman_e start_POSTSUPERSCRIPT italic_Ī² over~ start_ARG italic_Ļ• end_ARG end_POSTSUPERSCRIPT [ 4 italic_V ( over~ start_ARG italic_Ļ• end_ARG ) + 2 italic_Ī² start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT divide start_ARG roman_d italic_V ( over~ start_ARG italic_Ļ• end_ARG ) end_ARG start_ARG italic_d over~ start_ARG italic_Ļ• end_ARG end_ARG ] . (54)

We have plotted the above expression for the potential (52) in Fig.Ā 5, for different values of {Ī›,a}Ī›š‘Ž\{\Lambda,a\}{ roman_Ī› , italic_a }.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: Plots of fā¢(R)ƗRš‘“š‘…š‘…f(R)\times Ritalic_f ( italic_R ) Ɨ italic_R for different values of the free parameters Ī›Ī›\Lambdaroman_Ī› and aš‘Žaitalic_a, as indicated in each panel. Note the absence of the unstable region (where fā€²ā€²<0superscriptš‘“ā€²ā€²0f^{\prime\prime}<0italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT < 0) for large Ī›Ī›\Lambdaroman_Ī›. Different values of aš‘Žaitalic_a simply change the axis scales.

We first notice a few important characteristics: First and more obvious: the function fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) is multi valued. Its actual value depends on Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG while it evolves towards the bottom of its potential. Secondly, fā¢(R)>0š‘“š‘…0f(R)>0italic_f ( italic_R ) > 0 for all Rš‘…Ritalic_R and fā€²ā€²<0superscriptš‘“ā€²ā€²0f^{\prime\prime}<0italic_f start_POSTSUPERSCRIPT ā€² ā€² end_POSTSUPERSCRIPT < 0 only in the middle, lower branch, indicating instability of the theory. Thirdly, that very branch is absent for Ī›ā©¾Ī›c=15/16Ī›subscriptĪ›š‘1516\Lambda\geqslant\Lambda_{c}=15/16roman_Ī› ā©¾ roman_Ī› start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 15 / 16, which indicates it is a control parameter, such as temperature.

Besides, if we look at Eq.Ā (2) as the Lagrangian of a relativistic fluid, described by its pressure Pš‘ƒPitalic_P, the final step is almost automatic: to make the correspondence between Rš‘…Ritalic_R itself and the Gibbs Energy GšŗGitalic_G. This is clearly seen if we plot the 3D behavior of fā¢(R,Ī›)š‘“š‘…Ī›f(R,\Lambda)italic_f ( italic_R , roman_Ī› ) and the Gibbs Energy Gā¢(P,T)šŗš‘ƒš‘‡G(P,T)italic_G ( italic_P , italic_T )) in Fig.Ā 6.

Refer to caption
Refer to caption
Figure 6: Plots of (left panel) fā¢(R,Ī›)š‘“š‘…Ī›f(R,\Lambda)italic_f ( italic_R , roman_Ī› ), given by Eqs.Ā (53) and (54) with a=Ī›=0š‘ŽĪ›0a=\Lambda=0italic_a = roman_Ī› = 0, and (right panel) Gā¢(P,T)šŗš‘ƒš‘‡G(P,T)italic_G ( italic_P , italic_T ) for de vdW gas. Two of the axes in the latter are inverted to allow a clear comparison to the former.

The effective volume Vš‘‰Vitalic_V is the variable ā€œcanonically conjugatedā€ to the effective pressure Pš‘ƒPitalic_P, i.e, since

dā¢Gā¢(P,T)=Vā‹…dā¢Pāˆ’Sā‹…dā¢T.š‘‘šŗš‘ƒš‘‡ā‹…š‘‰š‘‘š‘ƒā‹…š‘†š‘‘š‘‡dG(P,T)=V\cdot dP-S\cdot dT.italic_d italic_G ( italic_P , italic_T ) = italic_V ā‹… italic_d italic_P - italic_S ā‹… italic_d italic_T . (55)

We have, therefore, the full expressions for the Gibbs Energy GšŗGitalic_G, the pressure Pš‘ƒPitalic_P and the Volume Vš‘‰Vitalic_V, in terms of the field Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG:

Gšŗ\displaystyle Gitalic_G =eĪ²ā¢Ļ•~ā¢(2ā¢(Ļ•~āˆ’a)ā¢[Ī²ā¢(Ļ•~āˆ’a)+1]Ī²+4ā¢T)absentsuperscriptš‘’š›½~italic-Ļ•2~italic-Ļ•š‘Ždelimited-[]š›½~italic-Ļ•š‘Ž1š›½4š‘‡\displaystyle=e^{\beta\tilde{\phi}}\left(\frac{2(\tilde{\phi}-a)[\beta(\tilde{% \phi}-a)+1]}{\beta}+4T\right)= italic_e start_POSTSUPERSCRIPT italic_Ī² over~ start_ARG italic_Ļ• end_ARG end_POSTSUPERSCRIPT ( divide start_ARG 2 ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) [ italic_Ī² ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) + 1 ] end_ARG start_ARG italic_Ī² end_ARG + 4 italic_T ) (56)
Pš‘ƒ\displaystyle Pitalic_P =e2ā¢Ī²ā¢Ļ•~ā¢(2ā¢(Ļ•~āˆ’a)Ī²+(Ļ•~āˆ’a)2+2ā¢T)absentsuperscriptš‘’2š›½~italic-Ļ•2~italic-Ļ•š‘Žš›½superscript~italic-Ļ•š‘Ž22š‘‡\displaystyle=e^{2\beta\tilde{\phi}}\left(\frac{2(\tilde{\phi}-a)}{\beta}+(% \tilde{\phi}-a)^{2}+2T\right)= italic_e start_POSTSUPERSCRIPT 2 italic_Ī² over~ start_ARG italic_Ļ• end_ARG end_POSTSUPERSCRIPT ( divide start_ARG 2 ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) end_ARG start_ARG italic_Ī² end_ARG + ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_T ) (57)
Vš‘‰\displaystyle Vitalic_V =expā”(āˆ’Ī²ā¢Ļ•~)ā‡”Ļ•~=āˆ’1Ī²ā¢logā”(V),formulae-sequenceabsentš›½~italic-Ļ•ā‡”~italic-Ļ•1š›½š‘‰\displaystyle=\exp(-\beta\tilde{\phi})\quad\Leftrightarrow\quad\tilde{\phi}=-% \frac{1}{\beta}\log(V),= roman_exp ( - italic_Ī² over~ start_ARG italic_Ļ• end_ARG ) ā‡” over~ start_ARG italic_Ļ• end_ARG = - divide start_ARG 1 end_ARG start_ARG italic_Ī² end_ARG roman_log ( italic_V ) , (58)

which allows us to plot the 3D behavior of Gā¢(P,T)šŗš‘ƒš‘‡G(P,T)italic_G ( italic_P , italic_T ) in Fig.Ā 6. The unstable region, where d2ā¢G/dā¢P2>0superscriptš‘‘2šŗš‘‘superscriptš‘ƒ20d^{2}G/dP^{2}>0italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_G / italic_d italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT > 0 (notice that the GšŗGitalic_G axis points down) corresponds to higher values of the Gibbs Energy for a fixed Tš‘‡Titalic_T.

One can also calculate the Helmholtz energy

Fā¢(T,V)š¹š‘‡š‘‰\displaystyle F(T,V)italic_F ( italic_T , italic_V ) ā‰”Gāˆ’Pā‹…Vabsentšŗā‹…š‘ƒš‘‰\displaystyle\equiv G-P\cdot Vā‰” italic_G - italic_P ā‹… italic_V (59)
=1Vā¢[(a+1Ī²ā¢logā”V)2+2ā¢T]absent1š‘‰delimited-[]superscriptš‘Ž1š›½š‘‰22š‘‡\displaystyle=\frac{1}{V}\bigg{[}\left(a+\frac{1}{\beta}\log V\right)^{2}+2T% \bigg{]}= divide start_ARG 1 end_ARG start_ARG italic_V end_ARG [ ( italic_a + divide start_ARG 1 end_ARG start_ARG italic_Ī² end_ARG roman_log italic_V ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_T ] (60)

and the Entropy

Sā¢(T,V)š‘†š‘‡š‘‰\displaystyle S(T,V)italic_S ( italic_T , italic_V ) ā‰”āˆ’āˆ‚Fāˆ‚T|Vabsentevaluated-atš¹š‘‡š‘‰\displaystyle\equiv-\left.\frac{\partial F}{\partial T}\right|_{V}ā‰” - divide start_ARG āˆ‚ italic_F end_ARG start_ARG āˆ‚ italic_T end_ARG | start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT (61)
=āˆ’2V.absent2š‘‰\displaystyle=-\frac{2}{V}.= - divide start_ARG 2 end_ARG start_ARG italic_V end_ARG . (62)

There are two important notes from the result above. First, we recall the reader that, as the field Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG sets itself at the bottom of its potential (Ļ•~=a~italic-Ļ•š‘Ž\tilde{\phi}=aover~ start_ARG italic_Ļ• end_ARG = italic_a), the volume assumes its minimum value Vā†’expā”(āˆ’Ī²ā¢a)ā†’š‘‰š›½š‘ŽV\to\exp(-\beta a)italic_V ā†’ roman_exp ( - italic_Ī² italic_a )), which corresponds to a decrease in Entropy. That indicates that the system is not complete, i.e, we are not looking only at the full system. Indeed, we have not included matter nor radiation, which should then absorb the latent heat released in the phase transition from inflation to a radiation-dominated universe (see discussion further down). 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.Ā (56) and redefine the Gibbs Energy it as

G=eĪ²ā¢Ļ•~ā¢(2ā¢(Ļ•~āˆ’a)ā¢[Ī²ā¢(Ļ•~āˆ’a)+1]Ī²+4ā¢T)āˆ’2ā¢Tā¢eĪ²ā¢a,šŗsuperscriptš‘’š›½~italic-Ļ•2~italic-Ļ•š‘Ždelimited-[]š›½~italic-Ļ•š‘Ž1š›½4š‘‡2š‘‡superscriptš‘’š›½š‘ŽG=e^{\beta\tilde{\phi}}\left(\frac{2(\tilde{\phi}-a)[\beta(\tilde{\phi}-a)+1]}% {\beta}+4T\right)-2Te^{\beta a},italic_G = italic_e start_POSTSUPERSCRIPT italic_Ī² over~ start_ARG italic_Ļ• end_ARG end_POSTSUPERSCRIPT ( divide start_ARG 2 ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) [ italic_Ī² ( over~ start_ARG italic_Ļ• end_ARG - italic_a ) + 1 ] end_ARG start_ARG italic_Ī² end_ARG + 4 italic_T ) - 2 italic_T italic_e start_POSTSUPERSCRIPT italic_Ī² italic_a end_POSTSUPERSCRIPT , (63)

which does not spoil the previous results.

From Eqs.Ā (57) and (58), we can write the equation of state, i.e, the relation between the pressure, volume and temperature:

Pā¢(V,T)=Ī²ā¢(a2ā¢Ī²āˆ’2ā¢a+2ā¢Ī²ā¢T)+(2ā¢aā¢Ī²āˆ’2+logā”V)ā¢logā”VĪ²2ā¢V2.š‘ƒš‘‰š‘‡š›½superscriptš‘Ž2š›½2š‘Ž2š›½š‘‡2š‘Žš›½2š‘‰š‘‰superscriptš›½2superscriptš‘‰2P(V,T)=\frac{\beta\left(a^{2}\beta-2a+2\beta T\right)+(2a\beta-2+\log V)\log V% }{\beta^{2}V^{2}}.italic_P ( italic_V , italic_T ) = divide start_ARG italic_Ī² ( italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_Ī² - 2 italic_a + 2 italic_Ī² italic_T ) + ( 2 italic_a italic_Ī² - 2 + roman_log italic_V ) roman_log italic_V end_ARG start_ARG italic_Ī² start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (64)

The behavior of Pā¢(V,T)š‘ƒš‘‰š‘‡P(V,T)italic_P ( italic_V , italic_T ) for five different values of Tš‘‡Titalic_T is shown in Fig.Ā 7, which bears strong resemblance to a vdW gas555Nevertheless, here one obtains PāˆTā¢Vāˆ’2proportional-toš‘ƒš‘‡superscriptš‘‰2P\propto TV^{-2}italic_P āˆ italic_T italic_V start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT in the high-temperature limit, instead of the standard ideal-gas behavior PāˆTā¢Vāˆ’1proportional-toš‘ƒš‘‡superscriptš‘‰1P\propto TV^{-1}italic_P āˆ italic_T italic_V start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT.. Even though the equations of state are not exactly the same, they do describe the same phenomena, as we will now see.

Refer to caption
Figure 7: Plot of the effective pressure Pš‘ƒPitalic_P as a function of the effective volume Vš‘‰Vitalic_V, for a=aāˆ—ā‰”1/Ī²š‘Žsubscriptš‘Ž1š›½a=a_{*}\equiv 1/\betaitalic_a = italic_a start_POSTSUBSCRIPT āˆ— end_POSTSUBSCRIPT ā‰” 1 / italic_Ī² and different values of temperature: T=Tcā‰”15/16š‘‡subscriptš‘‡š‘1516T=T_{c}\equiv 15/16italic_T = italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ā‰” 15 / 16 (solid thick black); lower (higher) curves, in solid thin gray, correspond to lower (higher) temperatures. The lowest one (bottom not shown) corresponds to Ī›ā‰”T=0Ī›š‘‡0\Lambda\equiv T=0roman_Ī› ā‰” italic_T = 0. The spinodal curve is plotted in dotted red. The binodal curve is plotted in dashed blue. The gray circle indicates the final configuration (Ļ•=aitalic-Ļ•š‘Ž\phi=aitalic_Ļ• = italic_a) for the T=0š‘‡0T=0italic_T = 0 case (higher temperatures correspond to higher final pressures).

The interesting Physics happens around the region where dā¢P/dā¢V>0š‘‘š‘ƒš‘‘š‘‰0dP/dV>0italic_d italic_P / italic_d italic_V > 0, i.e, the unstable one. Fig.Ā 7 shows also the binodal and spinodal lines. The former indicates the limit of existence of unstable configurations (which are exactly the Maxwell Construction, also known as the equal-area principle), while below the latter there are unstable states.

One can map the initial conditions that satisfy the slow-roll conditions (large negative Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG) into the right-most region in Fig.Ā 7. If we follow the (T=0)š‘‡0(T=0)( italic_T = 0 )-curve, the system starts below the binodal line, which is exactly what is needed: the inflationary era should be metastable ā€” i.e, it should last a few (āˆ¼60similar-toabsent60\sim 60āˆ¼ 60) e-folds, not too short (as it would be for an unstable state) and not too long (as in a stable state). The final state, when Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG is oscillating around the bottom of its potential, corresponds to the opposite end of the plot (small volume), in a different phase ā€” for the case at hand, with a quadratic potential (52), the average value of the equation-of-state parameter for the field Ļ•~~italic-Ļ•\tilde{\phi}over~ start_ARG italic_Ļ• end_ARG is āŸØĻ‰~Ļ•āŸ©ā‰”āŸØp~Ļ•/Ļ~Ļ•āŸ©ā‰ˆ0delimited-āŸØāŸ©subscript~šœ”italic-Ļ•delimited-āŸØāŸ©subscript~š‘italic-Ļ•subscript~šœŒitalic-Ļ•0\langle\tilde{\omega}_{\phi}\rangle\equiv\langle\tilde{p}_{\phi}/\tilde{\rho}_% {\phi}\rangle\approx 0āŸØ over~ start_ARG italic_Ļ‰ end_ARG start_POSTSUBSCRIPT italic_Ļ• end_POSTSUBSCRIPT āŸ© ā‰” āŸØ over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_Ļ• end_POSTSUBSCRIPT / over~ start_ARG italic_Ļ end_ARG start_POSTSUBSCRIPT italic_Ļ• end_POSTSUBSCRIPT āŸ© ā‰ˆ 0, i.e, a matter-dominated universe. We refer the reader to Ref.Ā PeraltaĀ andĀ JorĆ” s (2020) to a numerical description and the full Thermodynamic correspondence.

IV Final summary

We have seen a strong correspondence between fā¢(R)š‘“š‘…f(R)italic_f ( italic_R ) theories (in the metric approach) and a first-order phase transition, as described by Catastrophe Theory. That analogy has only been investigated in the inflationary scenario, so far.

Current work in under way for the (p)reheating mechanism that leads to an homogeneous hot dense universe compatible to a radiation-dominated phase followed by a matter-dominated one, that GR describes so nicely. Of course, relativistic stars and perturbation growth are also part of the next steps, as well as more complex potentials in the Einstein frame.

I would like to thank Mario Novello for keeping such a successful School (along with its international counterpart, the Brazilian School of Cosmology and Gravitation) alive and thriving for so many years. I acknowledge financial support from FAPERJ.

References