Thermodynamics of Theories
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 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 theories.
This a very short review on 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 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 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 , which must be seen either as a new constant of nature (such as Newtonās Gravitational Constant 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 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 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 theory in the metric frame, where the metric 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 is an independent object as well.
Of course, if , 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)
(1) |
by
(2) |
where . The former yields the Einstein Equations
(3) |
while the latter yields
(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 itself depends on second-order derivatives on ). 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
(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ā ā since it is pure geometry. Such form is more than adequate (as opposed to an alternate form that divides the previous equation by , 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 and pressure for such fluid, as
(6) | ||||
(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 theories. We start over from the Lagrangian (2) and write it as
(8) |
where
(9) | ||||
(10) |
which is simply a Legendre Transform, that replaces by as the independent variable. For that, we require a definite sign for (usually taken positive), except, perhaps, at particular values of .
We now perform a conformal transformation where the conformal factor is itself:
(11) |
and define a new field
(12) |
which requires that , already assumed when we performed the conformal transformation (11). Note that, since the field is dimensionless and , then has the dimensions of .
Upon the definitions above, Eq.Ā (2) is written as
(13) |
where
(14) |
and the determinant and the Ricci scalar are calculated from the new metric .
The Lagrangian (13) is the one of GR with an extra scalar field , which carries the extra degree of freedom, as previously announced. One should notice that the field (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 , since this is supposed to hold in the previous metric only. One can also notice that is not minimally coupled to the metric since 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
(15) |
which defines and , where the latter indicates the amplitude of the deviation from the Einstein-Hilbert Lagrangian . The standard constraints on theories are almost always listed then as PogosianĀ andĀ Silvestri (2008)
(16) | ||||
(17) | ||||
(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 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 in the limit to the ones from the āoriginalā equation . 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
(19) |
(where is the trace of ) which is not an algebraic relation anymore (as in GR). Therefore, one could still have a large (as it is expected in the early universe), but a small .
The second constraint usually implies that it is required for an attractive gravity. It relies on dividing Eq.Ā (4) by and noticing that it would change the sign of the Gravitational coupling . Nevertheless, that mapping would disrupt the very conservation of 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 . Here, we follow Ref.Ā SotiriouĀ andĀ Faraoni (2010). Let us expand Eq.Ā (19) around Minkwoski metric , i.e,
(20) | ||||
(21) |
and assume Eq.Ā (15), which yields
(22) |
Since can be as small as needed, the term inside the parenthesis is always positive. The sign of the effective mass squared for is the same of , 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 so that the picture described above is fulfilled. The authors define the quantities
(23) | ||||
(24) |
which measure the deviation from GR ā for which and .
The existence of a MD phase with the required constraints is translated to
(25) |
while a de Sitter attractor requires either
(29) | |||
(30) | |||
(31) |
That paper has shown that simple expressions for , 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 (its standard value) or to HerreraĀ etĀ al. (2017).
A second approach is to model the value of as a function of wavenumber 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 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 ; in the latter, the time . 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
(32) |
or the plain integration of the energy density (with or without the proper-volume factor ) bounded by the starās surface at
(33) | ||||
(34) |
and the āsurface-redshift massā, estimated from the light redshift (emitted from the surface to infinity) as
(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 , 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 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 ā 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):
(36) |
where the constants and 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 , 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.
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 , there is a volume range where
(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 , whose variation is given by
(38) |
Therefore,
(39) | ||||
(40) |
This means that in the unstable region.
One can always write Eq.Ā (36) as
(41) |
As a cubic equation with real coefficients, there is always (i.e, for any ) at least one real root, i.e, one pressure and one temperature for a given volume . On the other hand, for , Eq.Ā (41) yields three real roots: three values for (one of them in the aforementioned unstable region and, therefore, non physical) for a a given pressure and temperature . At the critical temperature , the system has critical pressure and critical Volume . Above , there is no difference between the states and , 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:
(42) | ||||
(43) | ||||
(44) |
upon which Eq.Ā (41) is written as
(45) | ||||
(46) | ||||
(47) |
The plot of , i.e, the equilibrium configurations as a function of the parameters 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 of the vdW gas. . A simple and general expression444One can get rid of a possible cubic term by a suitable shift on . that does the job is
(48) |
where we have made explicit the dependence on the control parameters . 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:
(49) |
where is one of the extrema, given, of course, by
(50) |
which is the so-called Equilibrium Set. The solution of Eqs.Ā (49) and (50) is the so-called Bifurcation Set
(51) |
In the parameter space , 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 .
In Fig.Ā 4, one can see the equilibrium surface (i.e, the values of the equilibrium points) as a function of the control parameters . 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 , i.e., toward the closest fold branch, then the largest (smallest) values of coalesce.
Accordingly, we can plot the so-called swallowtail curve also as a function of ā see Fig.Ā 4. one can see that the highest 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 that coalesces with the one corresponding to the intermediate value , i.e, the local maximum. Outside the fold, only the global minimum survives.
III 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 (in the Jordan frame) and, via a conformal transformation given by , arrives at the Einstein frame, where the physics is described by GR and an extra scalar field subject to a potential ), completely defined by the function 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 in the Jordan frame. We assume
(52) |
For now, the extra (constant) parameters and are simply a matter of generalizing the potential , while keeping it as simple as possible. One might argue against the introduction of a cosmological constant , 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 :
(53) | ||||
(54) |
We have plotted the above expression for the potential (52) in Fig.Ā 5, for different values of .
We first notice a few important characteristics: First and more obvious: the function is multi valued. Its actual value depends on while it evolves towards the bottom of its potential. Secondly, for all and only in the middle, lower branch, indicating instability of the theory. Thirdly, that very branch is absent for , 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 , the final step is almost automatic: to make the correspondence between itself and the Gibbs Energy . This is clearly seen if we plot the 3D behavior of and the Gibbs Energy ) in Fig.Ā 6.
The effective volume is the variable ācanonically conjugatedā to the effective pressure , i.e, since
(55) |
We have, therefore, the full expressions for the Gibbs Energy , the pressure and the Volume , in terms of the field :
(56) | ||||
(57) | ||||
(58) |
which allows us to plot the 3D behavior of in Fig.Ā 6. The unstable region, where (notice that the axis points down) corresponds to higher values of the Gibbs Energy for a fixed .
One can also calculate the Helmholtz energy
(59) | ||||
(60) |
and the Entropy
(61) | ||||
(62) |
There are two important notes from the result above. First, we recall the reader that, as the field sets itself at the bottom of its potential (), the volume assumes its minimum value ), 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
(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:
(64) |
The behavior of for five different values of is shown in Fig.Ā 7, which bears strong resemblance to a vdW gas555Nevertheless, here one obtains in the high-temperature limit, instead of the standard ideal-gas behavior .. Even though the equations of state are not exactly the same, they do describe the same phenomena, as we will now see.
The interesting Physics happens around the region where , 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 ) into the right-most region in Fig.Ā 7. If we follow the -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 () 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 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 is , 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 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
- SotiriouĀ andĀ Faraoni (2010) T.Ā P.Ā SotiriouĀ andĀ V.Ā Faraoni,Ā Reviews of Modern PhysicsĀ 82,Ā 451 (2010),Ā arXiv:0805.1726 [gr-qc] .
- Sotiriou (2007) T.Ā P.Ā Sotiriou,Ā āModified actions for gravity: Theory and phenomenology,āĀ (2007),Ā arXiv:0710.4438 [gr-qc] .
- De FeliceĀ andĀ Tsujikawa (2010) A.Ā De FeliceĀ andĀ S.Ā Tsujikawa,Ā Living Reviews in RelativityĀ 13,Ā 3 (2010),Ā arXiv:1002.4928 [gr-qc] .
- AmendolaĀ andĀ Tsujikawa (2015) L.Ā AmendolaĀ andĀ S.Ā Tsujikawa,Ā Dark EnergyĀ (2015).
- Callen (1985) H.Ā B.Ā Callen,Ā Thermodynamics and an Introduction to Thermostatistics, 2nd EditionĀ (1985).
- Gilmore (1981) R.Ā Gilmore,Ā Catastrophe Theory for Scientists and EngineersĀ (Dover,Ā New York,Ā 1981).
- PostonĀ andĀ Stewart (1996) T.Ā PostonĀ andĀ I.Ā Stewart,Ā Catastrophe Theory and Its Applications,Ā Dover books on mathematicsĀ (Dover Publications,Ā 1996).
- Saunders (1980) P.Ā T.Ā Saunders,Ā An Introduction to Catastrophe TheoryĀ (Cambridge University Press,Ā 1980).
- Milgrom (1983) M.Ā Milgrom,Ā Astrophys.Ā J.Ā Ā 270,Ā 371 (1983).
- MarraĀ etĀ al. (2020) V.Ā Marra, D.Ā C.Ā Rodrigues, Ā andĀ A.Ā O.Ā F.Ā deĀ Almeida,Ā Monthly Notices of the Royal Astronomical SocietyĀ 494,Ā 2875ā2885 (2020).
- Olmo (2011) G.Ā J.Ā Olmo,Ā International Journal of Modern Physics DĀ 20,Ā 413 (2011),Ā arXiv:1101.3864 [gr-qc] .
- SotiriouĀ andĀ Liberati (2007) T.Ā P.Ā SotiriouĀ andĀ S.Ā Liberati,Ā Annals of PhysicsĀ 322,Ā 935 (2007),Ā arXiv:gr-qc/0604006 [gr-qc] .
- Hernandez-ArboledaĀ etĀ al. (2023) A.Ā Hernandez-Arboleda, D.Ā C.Ā Rodrigues, J.Ā D.Ā Toniato, Ā andĀ A.Ā Wojnar,Ā āPalatini gravity tests in the weak field limit: Solar system, seismology and galaxies,āĀ (2023),Ā arXiv:2306.04475 [gr-qc] .
- Kaiser (1995) D.Ā I.Ā Kaiser,Ā Phys.Ā Rev.Ā DĀ 52,Ā 4295 (1995),Ā arXiv:astro-ph/9408044 [astro-ph] .
- PogosianĀ andĀ Silvestri (2008) L.Ā PogosianĀ andĀ A.Ā Silvestri,Ā Phys. Rev. DĀ 77,Ā 023503 (2008).
- AmendolaĀ etĀ al. (2007) L.Ā Amendola, R.Ā Gannouji, D.Ā Polarski, Ā andĀ S.Ā Tsujikawa,Ā Physical Review DĀ 75 (2007),Ā 10.1103/physrevd.75.083504.
- HerreraĀ etĀ al. (2017) D.Ā Herrera, I.Ā Waga, Ā andĀ S.Ā JorĆ”s,Ā Physical Review DĀ 95 (2017),Ā 10.1103/physrevd.95.064029.
- BertschingerĀ andĀ Zukin (2008) E.Ā BertschingerĀ andĀ P.Ā Zukin,Ā Phys.Ā Rev.Ā DĀ 78,Ā 024015 (2008),Ā arXiv:0801.2431 [astro-ph] .
- BorisovĀ etĀ al. (2012) A.Ā Borisov, B.Ā Jain, Ā andĀ P.Ā Zhang,Ā Phys.Ā Rev.Ā DĀ 85,Ā 063518 (2012),Ā arXiv:1102.4839 [astro-ph.CO] .
- Batista (2021) R.Ā C.Ā Batista,Ā UniverseĀ 8,Ā 22 (2021).
- PretelĀ etĀ al. (2023) J.Ā M.Ā Z.Ā Pretel, S.Ā E.Ā JorĆ”s, R.Ā R.Ā R.Ā Reis, S.Ā B.Ā Duarte, Ā andĀ J.Ā D.Ā V.Ā ArbaƱil,Ā arXiv e-printsĀ ,Ā arXiv:2308.00203 (2023),Ā arXiv:2308.00203 [gr-qc] .
- SbisĆ Ā etĀ al. (2020) F.Ā SbisĆ , P.Ā O.Ā Baqui, T.Ā Miranda, S.Ā E.Ā JorĆ”s, Ā andĀ O.Ā F.Ā Piattella,Ā Physics of the Dark UniverseĀ 27,Ā 100411 (2020),Ā arXiv:1907.08714 [gr-qc] .
- KhouryĀ andĀ Weltman (2004) J.Ā KhouryĀ andĀ A.Ā Weltman,Ā Phys.Ā Rev.Ā DĀ 69,Ā 044026 (2004),Ā arXiv:astro-ph/0309411 [astro-ph] .
- vanĀ der Waals (ates) 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)).
- PeraltaĀ andĀ JorĆ” s (2020) C.Ā PeraltaĀ andĀ S.Ā JorĆ” s,Ā Journal of Cosmology and Astroparticle PhysicsĀ 2020,Ā 053 (2020).
- MagnanoĀ andĀ SokoÅowski (1994) G.Ā MagnanoĀ andĀ L.Ā M.Ā SokoÅowski,Ā Phys.Ā Rev.Ā DĀ 50,Ā 5039 (1994),Ā arXiv:gr-qc/9312008 [gr-qc] .