This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Counterexamples to the maximum force conjecture

Aden Jowsey and Matt Visser ID
Abstract

Dimensional analysis shows that the speed of light and Newton’s constant of gravitation can be combined to define a quantity F=c4GNF_{*}={c^{4}\over G_{N}} with the dimensions of force (equivalently, tension). Then in any physical situation we must have Fphysical=fFF_{\mathrm{physical}}=f\;F_{*}, where the quantity ff is some dimensionless function of dimensionless parameters. In many physical situations explicit calculation yields f=𝒪(1)f={\mathcal{O}}(1), and quite often f14f\leq{1\over 4}. This has lead multiple authors to suggest a (weak or strong) maximum force/maximum tension conjecture. Working within the framework of standard general relativity, we will instead focus on counter-examples to this conjecture, paying particular attention to the extent to which the counter-examples are physically reasonable. The various counter-examples we shall explore strongly suggest that one should not put too much credence into any universal maximum force/maximum tension conjecture. Specifically, fluid spheres on the verge of gravitational collapse will generically violate the weak (and strong) maximum force conjectures. If one wishes to retain any general notion of “maximum force” then one will have to very carefully specify precisely which forces are to be allowed within the domain of discourse.


Date: Wednesday 3 February 2021; Tuesday 2 March 2021; -ed


Keywords: maximum force; maximum tension; general relativity.

1 Introduction

The maximum force/maximum tension conjecture was independently mooted some 20 years ago by Gary Gibbons [1] and Christoph Schiller [2]. At its heart one starts by noting that in (3+1) dimensions the quantity

F=c4GN1.2×1044 N\displaystyle F_{*}=\frac{c^{4}}{G_{N}}\approx 1.2\times 10^{44}\hbox{ N} (1.1)

has the dimensions of force (equivalently, tension). Here cc is the speed of light in vacuum, and GNG_{N} is Newton’s gravitational constant. Thereby any physical force can always be written in the form

Fphysical=fF,F_{\mathrm{physical}}=f\;F_{*}, (1.2)

where the quantity ff is some dimensionless function of dimensionless parameters. In very many situations [1, 2, 3, 4] explicit calculations yield f14f\leq{1\over 4}, though sometimes numbers such as f12f\leq{1\over 2} also arise [5]. Specifically, Yen Chin Ong [5] formulated strong and weak versions of the conjecture:

  1. 1.

    Strong form:  f14f\leq{1\over 4}.

  2. 2.

    Weak form:  f=𝒪(1)f={\mathcal{O}}(1).

Note that F=EPlanck/LPlanckF_{*}=E_{\mathrm{Planck}}/L_{\mathrm{Planck}} can also be interpreted as the Planck force, though it is not intrinsically quantum as the various factors of \hbar cancel, at least in (3+1) dimensions. Furthermore it is sometimes interesting [6] to note that the Einstein equations

Gab=8πGNc4Tab,G_{ab}=8\pi\;{G_{N}\over c^{4}}\;T_{ab}, (1.3)

can be written in terms of FF_{*} as

Tab=F8πGab.T_{ab}={F_{*}\over 8\pi}\;G_{ab}. (1.4)

When recast in this manner, maximum forces conjectures have tentatively been related to Jacobson’s entropic derivation of the Einstein equations [7].

Considerable work has also gone into attempts at pushing various modifications of the maximum force conjecture beyond the framework of standard general relativity [8, 9]. Overall, while there is little doubt that the quantity FF_{*} is physically important, we feel that the precise status of the maximum force conjecture is much less certain, and is less than universal.

We shall investigate these conjectures within the context of standard general relativity, focussing on illustrative counter-examples based on simple physical systems, analyzing the internal forces, and checking the extent to which the counter-examples are physically reasonable. Specifically, we shall consider static spherically symmetric fluid spheres [10, 11, 12, 13, 14, 15, 16, 17, 18], and investigate both radial and equatorial forces. We shall also include an analysis of the speed of sound, and the relevant classical energy conditions, specifically the dominant energy condition (DEC), see [19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. We shall see that even the most elementary static spherically symmetric fluid sphere, Schwarzschild’s constant density star, raises significant issues for the maximum force conjecture. Other models, such as the Tolman IV solution and its variants are even worse. Generically, we shall see that any prefect fluid sphere on the verge of gravitational collapse will violate the weak (and strong) maximum force conjectures. Consequently, if one wishes to retain any truly universal notion of “maximum force” then one will at the very least have to very carefully delineate precisely which forces are to be allowed within the domain of discourse.

2 Spherical symmetry

Consider spherically symmetric spacetime, with metric given in Schwarzschild curvature coordinates:

ds2=gttdt2+grrdt2+r2(dθ2+sin2θdϕ2).ds^{2}=g_{tt}\;dt^{2}+g_{rr}\;dt^{2}+r^{2}(d\theta^{2}+\sin^{2}\theta\;d\phi^{2}). (2.1)

We do not yet demand pressure isotropy, and for the time being allow radial and transverse pressures to differ, that is prptp_{r}\neq p_{t}.

Pick a spherical surface at some specified value of the radial coordinate rr. Define

Fr(r)=pr(r)𝑑A=4πpr(r)r2.F_{r}(r)=\int p_{r}(r)\;dA=4\pi\;p_{r}(r)\;r^{2}. (2.2)

This quantity simultaneously represents the compressive force exerted by outer layers of the system on the core, and the supporting force exerted by the core on the outer layers of the system.

Consider any equatorial slice through the system and define the equatorial force by

Feq=pt(r)𝑑A=2π0Rsgrrpt(r)r𝑑r.F_{eq}=\int p_{t}(r)\;dA=2\pi\int^{R_{s}}_{0}\sqrt{g_{rr}}\;p_{t}(r)\;rdr. (2.3)

This quantity simultaneously represents the force exerted by the lower hemisphere of the system on the upper hemisphere, and the force exerted by the upper hemisphere of the system on the lower hemisphere. Here RsR_{s} is the location of the surface of the object (potentially taken as infinite). As we are investigating with spherically symmetric systems, the specific choice of hemisphere is irrelevant.

3 Perfect fluid spheres

3.1 Generalities

The perfect fluid condition excludes pressure anisotropy so that radial and transverse pressures are set equal: p(r)=pr(r)=pt(r)p(r)=p_{r}(r)=p_{t}(r). Once this is done, the radial and equatorial forces simplify

Fr(r)=p(r)𝑑A=4πp(r)r2;F_{r}(r)=\int p(r)\;dA=4\pi\;p(r)\;r^{2}; (3.1)
Feq=p(r)𝑑A=2π0Rsgrrp(r)r𝑑r.F_{eq}=\int p(r)\;dA=2\pi\int^{R_{s}}_{0}\sqrt{g_{rr}}\;p(r)\;rdr. (3.2)

Additionally, we shall impose the conditions that pressure is positive and decreases as one moves outwards with zero pressure defining the surface of the object [10, 11, 12, 13, 14, 15, 16, 17, 18].111There is a minor technical change in the presence of a cosmological constant, the surface is then defined by p(Rs)=pΛp(R_{s})=p_{\Lambda}. Similarly density is positive and does not increase as one moves outwards, though density need not be and typically is not zero at the surface [10, 11, 12, 13, 14, 15, 16, 17, 18].

We note that for the radial force we have by construction

Fr(0)=0;Fr(Rs)=0;and forr(0,Rs):Fr(r)>0.F_{r}(0)=0;\qquad F_{r}(R_{s})=0;\qquad\hbox{and for}\quad r\in(0,R_{s}):\;F_{r}(r)>0. (3.3)

In particular in terms of the central pressure p0p_{0} we have the particularly simple bound

Fr(r)<4πp0Rs2.F_{r}(r)<4\pi\;p_{0}\;R_{s}^{2}. (3.4)

This suggests that in general an (extremely) weak version of the maximum force conjecture might hold for the radial force, at least within the framework outlined above, and as long as the central pressure is finite. Unfortunately without some general relationship between central pressure p0p_{0} and radius RsR_{s} this bound is less useful than one might hope. For the strong version of the maximum force conjecture no such simple argument holds for FrF_{r}, and one must perform a case-by-case analysis. For the equatorial force FeqF_{eq} there is no similar argument of comparable generality, and one must again perform a case-by-case analysis.

Turning now to the classical energy conditions [19, 20, 21, 22, 23, 24, 25, 26, 27, 28], they add extra restrictions to ensure various physical properties remain well-behaved. For our perfect fluid solutions, these act as statements relating the pressure pp and the density ρ\rho given by the stress-energy tensor Tμ^ν^T_{\hat{\mu}\hat{\nu}}. Since, (in view of our fundamental assumptions that pressure and density are both positive), the null, weak, and strong energy conditions, (NEC, WEC, SEC) are always automatically satisfied, we will only be interested in the dominant energy condition (DEC). In the current context the dominant energy condition only adds the condition |p|ρ|p|\leq\rho. But since in the context of perfect fluid spheres, the pressure is always positive, it is more convenient to simply write this as

pρ1;that ispρ.\frac{p}{\rho}\leq 1;\qquad\hbox{that is}\qquad p\leq\rho. (3.5)

The best physical interpretation of the DEC is that it guarantees that any timelike observer with 4-velocity VaV^{a} will observe a flux Fa=TabVbF^{a}=T^{ab}\,V_{b} that is non-spacelike (either timelike or null) [25]. However, it should be pointed out that the DEC, being the strongest of the classical energy conditions, is also the easiest to violate — indeed there are several known situations in which the classical DEC is violated by quantum effects [20, 21, 22, 23, 24, 25, 26, 27, 28].

The DEC is sometimes [somewhat misleadingly] interpreted in terms of the speed of sound not being superluminal: naively vs2=p/ρ1v_{s}^{2}=\partial p/\partial\rho\leq 1; whence pρρsurface<ρp\leq\rho-\rho_{\mathrm{surface}}<\rho. But the implication is only one-way, and in addition the argument depends on extra technical assumptions to the effect that the fluid sphere is well-mixed with a unique barotropic equation of state p(ρ)p(\rho) holding throughout the interior. To clarify this point, suppose the equation of state is not barotropic, so that p=p(ρ,zi)p=p(\rho,z_{i}), with the ziz_{i} being some collection of intensive variables, (possibly chemical concentrations, entropy density, or temperature). Then we have

dpdr=pρdρdr+ipzidzidr=vs2(ρ,zi)dρdr+ipzidzidr.{dp\over dr}={\partial p\over\partial\rho}\;{d\rho\over dr}+\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}=v_{s}^{2}(\rho,z^{i})\;{d\rho\over dr}+\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}. (3.6)

Then, (noting that dρ/drd\rho/dr is non-positive as one moves outwards), enforcing the speed of sound to not be superluminal implies

dpdrdρdr+ipzidzidr.{dp\over dr}\geq{d\rho\over dr}+\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}. (3.7)

Integrating this from the surface inwards we have

p(r)ρ(r)ρ(Rs)+irRsp(ρ,zi)zidzidr𝑑r.p(r)\leq\rho(r)-\rho(R_{s})+\sum_{i}\int_{r}^{R_{s}}{\partial p(\rho,z^{i})\over\partial z^{i}}\;{dz^{i}\over dr}\;dr. (3.8)

Consequently, unless one either makes an explicit barotropic assumption p/zi=0\partial p/\partial z^{i}=0, or otherwise at the very least has some very tight control over the partial derivatives p/zi\partial p/\partial z^{i}, one simply cannot use an assumed non-superluminal speed of sound to deduce the DEC. Neither can the DEC be used to derive a non-superluminal speed of sound, at least not without many extra and powerful technical assumptions. We have been rather explicit with this discussion since we have seen considerable confusion on this point. Finally we note that there is some disagreement as to whether or not the DEC is truly fundamental [21, 22, 23, 24].

3.2 Schwarzschild’s constant density star

We shall now consider a classic example of perfect fluid star, Schwarzschild’s constant density star [29], (often called the Schwarzschild interior solution), which was discovered very shortly after Schwarzschild’s original vacuum solution [30], (often called the Schwarzschild exterior solution).

It is commonly argued that Schwarzschild’s constant density star is “unphysical” on the grounds that it allegedly leads to an infinite speed of sound. But this is a naive result predicated on the physically unreasonable hypothesis that the star is well-mixed with a barotropic equation of state p=p(ρ)p=p(\rho). To be very explicit about this, all realistic stars are physically stratified with non-barotropic equations of state p=p(ρ,zi)p=p(\rho,z_{i}), with the ziz_{i} being some collection of intensive variables, (possibly chemical concentrations, entropy density, or temperature). We have already seen that

dpdr=pρdρdr+ipzidzidr=vs2(ρ,zi)dρdr+ipzidzidr.{dp\over dr}={\partial p\over\partial\rho}\;{d\rho\over dr}+\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}=v_{s}^{2}(\rho,z^{i})\;{d\rho\over dr}+\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}. (3.9)

Thence for a constant density star, dρ/dr=0d\rho/dr=0, we simply deduce

dpdr=ipzidzidr.{dp\over dr}=\sum_{i}{\partial p\over\partial z^{i}}\;{dz^{i}\over dr}. (3.10)

This tells us nothing about the speed of sound, one way or the other — it does tell us that there is a fine-tuning between the pressure pp and the intensive variables ziz^{i}, but that is implied by the definition of being a “constant density star”. We have been rather explicit with this discussion since we have seen considerable confusion on this point. Schwarzschild’s constant density star is not “unphysical”; it may be “fine-tuned” but it is not a priori “unphysical”.

Specifically, the Schwarzschild interior solution describes the geometry inside a static spherically symmetric perfect fluid constant density star with radius RsR_{s} and mass MM by the metric:

ds2=14(312MRs12Mr2Rs3)2dt2+(12Mr2Rs3)1dr2+r2dΩ2.ds^{2}=-\frac{1}{4}\left(3\sqrt{1-\frac{2M}{R_{s}}}-\sqrt{1-\frac{2Mr^{2}}{R_{s}^{3}}}\right)^{2}dt^{2}+\left(1-\frac{2Mr^{2}}{R_{s}^{3}}\right)^{-1}dr^{2}+r^{2}d\Omega^{2}. (3.11)

Here we have adopted geometrodynamic units (c1c\to 1, GN1G_{N}\to 1). Calculating the non-zero orthonormal stress-energy components from the Einstein equations applied to this metric yields:

Tt^t^\displaystyle T_{\hat{t}\hat{t}} =\displaystyle= ρ=3M4πRs3;\displaystyle\rho=\frac{3M}{4\pi R_{s}^{3}}; (3.12)
Tr^r^\displaystyle T_{\hat{r}\hat{r}} =\displaystyle= Tθ^θ^=Tϕ^ϕ^=p=ρ12Mr2Rs312MRs312MRs12Mr2Rs3.\displaystyle T_{\hat{\theta}\hat{\theta}}=T_{\hat{\phi}\hat{\phi}}=p=\rho\;\;\frac{\sqrt{1-\frac{2Mr^{2}}{R_{s}^{3}}}-\sqrt{1-\frac{2M}{R_{s}}}}{3\sqrt{1-\frac{2M}{R_{s}}}-\sqrt{1-\frac{2Mr^{2}}{R_{s}^{3}}}}. (3.13)

This gives us the relation between density and pressure, as well as demonstrating the perfect fluid condition (p=pr=ptp=p_{r}=p_{t}), and also verifying that the density is (inside the star) a position independent constant. In these geometrodynamic units both density and pressure have units 1/(length)2, while forces are dimensionless. Note that the pressure does in fact go to zero at rRsr\to R_{s}, so RsR_{s} really is the surface of the “star”. Rewriting the relation between pressure and density in terms of the simplified dimensionless quantities χ=2MRs\chi=\frac{2M}{R_{s}} and y=r2Rs2y=\frac{r^{2}}{R_{s}^{2}} we see

p=ρ1χy1χ31χ1χy.\displaystyle p=\rho\;\;\frac{\sqrt{1-\chi y}-\sqrt{1-\chi}}{3\sqrt{1-\chi}-\sqrt{1-\chi y}}. (3.14)

Here 0y10\leq y\leq 1, and 0χ<890\leq\chi<\frac{8}{9}. The first of these ranges is obvious from the definition of yy, while the second comes from considering the central pressure at y=0y=0:

p0=ρ11χ31χ1.p_{0}=\rho\;\;\frac{1-\sqrt{1-\chi}}{3\sqrt{1-\chi}-1}. (3.15)

Demanding that the central pressure be finite requires χ<89\chi<\frac{8}{9}. (This is actually a rather more general result of general relativistic stellar dynamics, not restricted to constant density, see various discussions of the Buchdahl–Bondi bound [31, 32].)

3.2.1 Radial Force

The radial force FrF_{r} as defined by equation (3.1) can be combined with the pressure-density relation given by equation (3.14), giving:

Fr=4πpr2=4πρRs2y1χy1χ31χ1χy=32χy1χy1χ31χ1χy.F_{r}=4\pi pr^{2}=4\pi\rho R_{s}^{2}\;y\;\frac{\sqrt{1-\chi y}-\sqrt{1-\chi}}{3\sqrt{1-\chi}-\sqrt{1-\chi y}}=\frac{3}{2}\chi y\;\frac{\sqrt{1-\chi y}-\sqrt{1-\chi}}{3\sqrt{1-\chi}-\sqrt{1-\chi y}}. (3.16)

As advertised in both abstract and introduction, this quantity is indeed a dimensionless function of dimensionless variables. Furthermore this quantity is defined on the bounded range 0y10\leq y\leq 1, 0χ<890\leq\chi<\frac{8}{9}. To find if FrF_{r} itself is bounded we analyse the multi-variable derivative for critical points.

For χFr\partial_{\chi}F_{r} we find:

χFr=3y2({4χ(3+y)}1χ1χy{4χ(3+5y4χy)}1χ1χy(31χ1χy)2).\partial_{\chi}F_{r}=-{3y\over 2}\left(\frac{\{4-\chi(3+y)\}\sqrt{1-\chi}\sqrt{1-\chi y}-\{4-\chi(3+5y-4\chi y)\}}{\sqrt{1-\chi}\sqrt{1-\chi y}\;(3\sqrt{1-\chi}-\sqrt{1-\chi y})^{2}}\right). (3.17)

For yFr\partial_{y}F_{r} we find:

yFr=3χ2({4χ(3+y)}1χy{45χy}1χ1χy(31χ1χy)2).\partial_{y}F_{r}=-{3\chi\over 2}\left(\frac{\{4-\chi(3+y)\}\sqrt{1-\chi y}-\{4-5\chi y\}\sqrt{1-\chi}}{\sqrt{1-\chi y}\;(3\sqrt{1-\chi}-\sqrt{1-\chi y})^{2}}\right). (3.18)

In particular we see that

χχFryyFr=3χ2y2(1χy)1χ(31χ1χy)2).\chi\partial_{\chi}F_{r}-y\partial_{y}F_{r}={3\chi^{2}y\over 2}\left(\frac{\sqrt{1-\chi y})}{\sqrt{1-\chi}\;(3\sqrt{1-\chi}-\sqrt{1-\chi y})^{2}}\right). (3.19)

To have a critical point, χFr=yFr=0\partial_{\chi}F_{r}=\partial_{y}F_{r}=0, we certainly require χy=0\chi y=0. So either χ=0\chi=0 or y=0y=0. But for y=0y=0, and χ(0,89)\chi\in(0,{8\over 9}) we have

yFr3χ(3χ+41χ4)2(31χ1)2>0.\partial_{y}F_{r}\quad\longrightarrow\quad{3\chi(3\chi+4\sqrt{1-\chi}-4)\over 2(3\sqrt{1-\chi}-1)^{2}}>0. (3.20)

In contrast, for χ=0\chi=0, and y(0,1)y\in(0,1), we have χFr0\partial_{\chi}F_{r}\to 0. So the only critical points lie on one of the boundary segments:

χFr=yFr=0χ=0.\partial_{\chi}F_{r}=\partial_{y}F_{r}=0\quad\iff\quad\chi=0. (3.21)

Therefore to find the maxima of Fr(χ,r)F_{r}(\chi,r) we must inspect all four of the boundary segments of the viable region. Along three of the boundary segments we can see that the three lines corresponding to χ=0\chi=0, y=0y=0, and y=1y=1 all give Fr(χ,r)=0F_{r}(\chi,r)=0, leaving only χ89\chi\rightarrow\frac{8}{9} to be investigated.

We note

limχ89Fr(χ,y)=4y3(98y1)398y.\lim_{\chi\rightarrow\frac{8}{9}}F_{r}(\chi,y)={4y\over 3}\;{(\sqrt{9-8y}-1)\over 3-\sqrt{9-8y}}. (3.22)

Inserting this into the partial derivative yFr\partial_{y}F_{r} reveals:

limχ89yFr=43(1+198y).\lim_{\chi\rightarrow\frac{8}{9}}\partial_{y}F_{r}=-\frac{4}{3}\left(1+\frac{1}{\sqrt{9-8y}}\right). (3.23)

This is a strictly negative function in the range 0y10\leq y\leq 1.

Thus the maximum of Fr(χ,y)F_{r}(\chi,y) can be found by taking the limit limy0\lim_{y\rightarrow 0} giving:

(Fr)max=limy0limχ89Fr=2.(F_{r})_{\mathrm{max}}=\lim_{y\rightarrow 0}\lim_{\chi\rightarrow\frac{8}{9}}F_{r}=2. (3.24)

This is therefore bounded, with the radial force approaching its maximum at the centre of a fluid star which is on the verge of collapse. This force violates the strong maximum force conjecture, though it satisfies the weak maximum force conjecture. This limit can easily be seen graphically in Figure 1.

Refer to caption
Figure 1: Radial force Fr(χ,y)F_{r}(\chi,y) for the interior Schwarzschild solution. Note Fr(χ,y)F_{r}(\chi,y) is bounded above by 2 in the region of interest y[0,1]y\in[0,1], χ[0,8/9)\chi\in[0,8/9).

3.2.2 Equatorial force

Using equation (3.2) and the metric defined in equation (3.11), with the relabelling of the previous subsection in terms of χ\chi and yy gives:

Feq(χ)=38χ0111χy(1χy1χ31χ1χy)𝑑y.\displaystyle F_{eq}(\chi)=\frac{3}{8}\;\chi\;\hbox{\Large$\displaystyle\int$}^{1}_{\!\!\!\!\!0}\frac{1}{\sqrt{1-\chi y}}\;\left(\frac{\sqrt{1-\chi y}-\sqrt{1-\chi}}{3\sqrt{1-\chi}-\sqrt{1-\chi y}}\right)dy. (3.25)

The integral evaluates to:

Feq(χ)\displaystyle F_{eq}(\chi) =\displaystyle= 34[1χy+21χln(31χ1χy)]|y=0y=1.\displaystyle\frac{3}{4}\left.\left[\sqrt{1-\chi y}+2\sqrt{1-\chi}\;\ln(3\sqrt{1-\chi}-\sqrt{1-\chi y})\right]\right|_{y=0}^{y=1}. (3.26)

Ultimately

Feq(χ)\displaystyle F_{eq}(\chi) =\displaystyle= 34[1χ{1+ln(44χ)2ln(31χ1)}1].\displaystyle\frac{3}{4}\left[\sqrt{1-\chi}\left\{1+\ln(4-4\chi)-2\ln(3\sqrt{1-\chi}-1)\right\}-1\right]. (3.27)

However, due to the presence of the ln(31χ1)-\ln(3\sqrt{1-\chi}-1) term in this equation, it can be seen that as χ89\chi\rightarrow\frac{8}{9}, Feq(χ)+F_{eq}(\chi)\rightarrow+\infty. Indeed

Feq(χ)=ln(89χ)2+𝒪(1),F_{eq}(\chi)={\ln({8\over 9}-\chi)\over 2}+{\mathcal{O}}(1), (3.28)

implying that the equatorial force in this space-time will grow without bound as the star approaches the critical size, (just prior to gravitational collapse), in violation of both the strong and weak maximum force conjectures.

So while the interior Schwarzschild solution has provided a nice example of a bounded radial force, Fr(y,χ)F_{r}(y,\chi), it also clearly provides an explicit counter-example, where the equatorial force Feq(χ)F_{eq}(\chi) between two hemispheres of the fluid star grows without bound.

3.2.3 DEC

Imposing the DEC (equation 3.5) within the fluid sphere we would require:

pρ=1χy1χ31χ1χy1\frac{p}{\rho}=\frac{\sqrt{1-\chi y}-\sqrt{1-\chi}}{3\sqrt{1-\chi}-\sqrt{1-\chi y}}\leq 1 (3.29)

That is

1χy41χ,\sqrt{1-\chi y}\leq 4\sqrt{1-\chi}, (3.30)

whence

1χy16(1χ).1-\chi y\leq 16(1-\chi). (3.31)

Applying the boundary conditions of 0χ890\leq\chi\leq\frac{8}{9}, 0y10\leq y\leq 1, we have a solution range:

(0χ34,0y1)(34<χ89,43χy1).\left(0\leq\chi\leq\frac{3}{4},\quad 0\leq y\leq 1\right)\quad\quad\bigcup\quad\quad\left(\frac{3}{4}<\chi\leq\frac{8}{9},\quad 4-\frac{3}{\chi}\leq y\leq 1\right). (3.32)

See figures 3 and 3.

Refer to caption
Figure 2: pρ\frac{p}{\rho} in first range
Refer to caption
Figure 3: pρ\frac{p}{\rho} in second range

Within the first region 0χ34,0y10\leq\chi\leq\frac{3}{4},\quad 0\leq y\leq 1, the radial force is maximised at:

χ=34,y=16(55)0.46;Fr=316(51)0.23<14.\chi=\frac{3}{4},\quad y=\frac{1}{6}(5-\sqrt{5})\approx 0.46;\quad\rightarrow\quad F_{r}=\frac{3}{16}(\sqrt{5}-1)\approx 0.23<{1\over 4}. (3.33)

Under these conditions the strong maximum force conjecture is satisfied. This can be seen visually in figure 4.

Refer to caption
Figure 4: Radial force Fr(χ,y)F_{r}(\chi,y) for the interior Schwarzschild solution in region 1 (0χ34,0y1)\left(0\leq\chi\leq\frac{3}{4},\quad 0\leq y\leq 1\right) where the DEC is satisfied.

Within the second region 34<χ89,34χy1\frac{3}{4}<\chi\leq\frac{8}{9},\quad 3-\frac{4}{\chi}\leq y\leq 1, the radial force is maximised at:

χ=89,y=58;Fr=56>14.\chi=\frac{8}{9},\quad y=\frac{5}{8};\quad\rightarrow\quad F_{r}=\frac{5}{6}>{1\over 4}. (3.34)

Under these conditions the strong maximum force conjecture is violated, though the weak maximum force conjecture is satisfied. This can be seen visually in figure 5.

Refer to caption
Figure 5: Radial force Fr(χ,y)F_{r}(\chi,y) for the interior Schwarzschild solution in region 2 (34<χ89,43χy1)\left(\frac{3}{4}<\chi\leq\frac{8}{9},\quad 4-\frac{3}{\chi}\leq y\leq 1\right) where the DEC is satisfied.

Turning to the equatorial force, we see that the integrand used to define integral for Feq(χ)F_{eq}(\chi) satisfies the DEC only within the range 0χ340\leq\chi\leq\frac{3}{4}. Using the result for Feq(χ)F_{eq}(\chi) given above, equation (3.27), we have:

(Feq)max,DEC=Feq(χ=3/4)=38(2ln21)0.1448603854<14.(F_{eq})_{\mathrm{max,DEC}}=F_{eq}(\chi=3/4)={3\over 8}(2\ln 2-1)\approx 0.1448603854<{1\over 4}. (3.35)

This now satisfies the strong maximum force conjecture.

3.2.4 Summary

Only if we enforce the DEC can we then make Schwarzschild’s constant density star satisfy the weak and strong maximum force conjectures. Without adding the DEC Schwarzschild’s constant density star will violate both the weak and strong maximum force conjectures. Since it is not entirely clear that the DEC represents fundamental physics [21, 22, 23, 24], it is perhaps a little sobering to see that one of the very simplest idealized stellar models already raises issues for the maximum force conjecture. We shall soon see that the situation is even worse for the Tolman IV model (and its variants).

3.3 Tolman IV solution

The Tolman IV solution is another perfect fluid star space-time, however it does not have the convenient (albeit fine-tuned) property of constant density like the interior Schwarzschild solution. The metric can be written in the traditional form [10]:

ds2=(1+r2A2)dt2+1+2r2A2(1r2R2)(1+r2A2)dr2+r2dΩ2.ds^{2}=-\left(1+\frac{r^{2}}{A^{2}}\right)dt^{2}+\frac{1+\frac{2r^{2}}{A^{2}}}{\left(1-\frac{r^{2}}{R^{2}}\right)\left(1+\frac{r^{2}}{A^{2}}\right)}dr^{2}+r^{2}d\Omega^{2}. (3.36)

Here AA and RR represent some arbitrary scale factors with units of length. This metric yields the orthonormal stress-energy tensor:

Tt^t^\displaystyle T_{\hat{t}\hat{t}} =\displaystyle= ρ=18π6r4+(7A2+2R2)r2+3A2(A2+R2)R2(A2+2r2)2;\displaystyle\rho=\frac{1}{8\pi}\frac{6r^{4}+(7A^{2}+2R^{2})r^{2}+3A^{2}(A^{2}+R^{2})}{R^{2}(A^{2}+2r^{2})^{2}}; (3.37)
Tr^r^\displaystyle T_{\hat{r}\hat{r}} =\displaystyle= Tθ^θ^=Tϕ^ϕ^=p=18πR2A23r2R2(A2+2r2).\displaystyle T_{\hat{\theta}\hat{\theta}}=T_{\hat{\phi}\hat{\phi}}=p=\frac{1}{8\pi}\frac{R^{2}-A^{2}-3r^{2}}{R^{2}(A^{2}+2r^{2})}. (3.38)

This demonstrates the non-constancy of the energy-density ρ\rho as well as the perfect fluid conditions. Physically, the surface of a fluid star is defined as the zero pressure point, which now is at:

Rs=R2A23.R_{s}=\sqrt{\frac{R^{2}-A^{2}}{3}}. (3.39)

And likewise we can find the surface density (ρ\rho at R0R_{0}):

ρs=34π2A2+R2R2(A2+2R2).\rho_{s}=\frac{3}{4\pi}\frac{2A^{2}+R^{2}}{R^{2}(A^{2}+2R^{2})}. (3.40)

The central pressure and central density are

p0=18πR2A2R2A2;ρ0=18π3(R2+A2)R2A2.p_{0}={1\over 8\pi}{R^{2}-A^{2}\over R^{2}A^{2}};\qquad\rho_{0}={1\over 8\pi}{3(R^{2}+A^{2})\over R^{2}A^{2}}. (3.41)

Moving forwards, we will likewise calculate the radial and equatorial forces in this space-time.

3.3.1 Radial force

Using the previously defined radial force equation, (2.2), we can write the radial force for the Tolman IV space-time as:

Fr(r,R,A)=4πprr2=r22R2(R2A23r2A2+2r2).F_{r}(r,R,A)=4\pi p_{r}r^{2}=\frac{r^{2}}{2R^{2}}\left(\frac{R^{2}-A^{2}-3r^{2}}{A^{2}+2r^{2}}\right). (3.42)

Defining y=r2/Rs2y=r^{2}/R_{s}^{2} and a=A2/R2a=A^{2}/R^{2} we have y(0,1)y\in(0,1) and a(0,1)a\in(0,1). The radial force then simplifies to

Fr(a,y)=(1a)2y(1y)6a+4(1a)y.F_{r}(a,y)={(1-a)^{2}y(1-y)\over 6a+4(1-a)y}. (3.43)

Note this is, as expected, a dimensionless function of dimensionless variables.

The multivariable derivatives are:

aFr(a,y)=y(1y)(1a)(2ay3a2y3)2(2ay3a2y)2.\partial_{a}F_{r}(a,y)={y(1-y)(1-a)(2ay-3a-2y-3)\over 2(2ay-3a-2y)^{2}}. (3.44)
yFr(a,y)=(1a)2(2ay26ay2y23a)2(2ay3a2y)2.\partial_{y}F_{r}(a,y)={(1-a)^{2}(2ay^{2}-6ay-2y^{2}-3a)\over 2(2ay-3a-2y)^{2}}. (3.45)

For both derivatives to vanish, (within the physical region), we require a=1a=1. However a=1a=1 actually minimises the function with Fr(a,y)=0F_{r}(a,y)=0. So we need to look at the boundaries of the physical region. Both y=0y=0 and y=1y=1 also minimise the function with Fr(a,y)=0F_{r}(a,y)=0. We thus consider a=0a=0:

Fr(0,y)=14(1y),F_{r}(0,y)=\frac{1}{4}(1-y), (3.46)

where then it is clear that the function is maximised at a=0a=0, y0y\to 0, which corresponds (Fr)max=14(F_{r})_{\mathrm{max}}={1\over 4}. This can be seen visually in figure 6. Thus Fr(a,y)F_{r}(a,y) for the Tolman IV solution is compatible with the strong maximum force conjecture, but as for the Schwarzschild constant density star, we shall soon see that the equatorial force does not behave as nicely.

Refer to caption
Figure 6: Fr(a,y)F_{r}(a,y) for the Tolman IV solution.

3.3.2 Equatorial force

Using equation (3.2) for this space-time, and combining it with radial surface result of equation (3.39), we obtain:

Feq=2π0Rsgrrp(r)r𝑑r=140R2A231+2r2A2(1r2R2)(1+r2A2)R2A23r2R2(A2+2r2)rdr.F_{eq}=2\pi\int^{R_{s}}_{0}\sqrt{g_{rr}}\;p(r)\;r\,dr=\frac{1}{4}\hbox{\Large$\displaystyle\int$}^{\sqrt{\frac{R^{2}-A^{2}}{3}}}_{\!\!\!\!\!0}\sqrt{\frac{1+\frac{2r^{2}}{A^{2}}}{\left(1-\frac{r^{2}}{R^{2}}\right)\left(1+\frac{r^{2}}{A^{2}}\right)}}\;\frac{R^{2}-A^{2}-3r^{2}}{R^{2}(A^{2}+2r^{2})}\;r\,dr. (3.47)

As an integral this converges, however the resultant function is intractable. Instead, we will opt for a simpler approach by finding a simple bound. Since the radial coordinate is physically bound by 0rRs=R2A23<R0\leq r\leq R_{s}=\sqrt{\frac{R^{2}-A^{2}}{3}}<R, we find that in that range:

grr=11r2R21+2r2A21+r2A21+2r2A21+r2A21.\displaystyle g_{rr}=\frac{1}{1-\frac{r^{2}}{R^{2}}}\frac{1+\frac{2r^{2}}{A^{2}}}{1+\frac{r^{2}}{A^{2}}}\geq\frac{1+\frac{2r^{2}}{A^{2}}}{1+\frac{r^{2}}{A^{2}}}\geq 1. (3.48)

This is actually a much more general result; for any perfect fluid sphere we have

grr=112m(r)/r,g_{rr}={1\over 1-2m(r)/r}, (3.49)

where m(r)m(r) is the Misner–Sharp quasi-local mass.

So as long as m(r)m(r) is positive, which is guaranteed by positivity of the density ρ(r)\rho(r), we have grr>1g_{rr}>1, and so in all generality we have

Feq>2π0Rsp(r)r𝑑r.F_{eq}>2\pi\int^{R_{s}}_{0}\;p(r)\;r\,dr. (3.50)

For the specific case of Tolman IV we can write

Feq\displaystyle F_{eq} >\displaystyle> 140R2A23R2A23r2R2(A2+2r2)r𝑑r.\displaystyle\frac{1}{4}\hbox{\Large$\displaystyle\int$}^{\sqrt{\frac{R^{2}-A^{2}}{3}}}_{0}\frac{R^{2}-A^{2}-3r^{2}}{R^{2}(A^{2}+2r^{2})}rdr. (3.51)

Now make the substitutions y=r2/Rs2y=r^{2}/R_{s}^{2} and a=A2/R2a=A^{2}/R^{2}. We find

Feq>πrs201p𝑑y=1801(1a)2(1y)2(1a)y+3a𝑑y.F_{eq}>\pi r_{s}^{2}\int^{1}_{0}p\;dy={1\over 8}\int^{1}_{0}\frac{\left(1-a\right)^{2}\left(1-y\right)}{2\left(1-a\right)y+3a}\;dy. (3.52)

This integral yields

Feq>{(a+2)ln[(2y3)a2y]32(1a)y16}|01.F_{eq}>\left.\left\{{(a+2)\ln[(2y-3)a-2y]\over 32}-{(1-a)y\over 16}\right\}\right|_{0}^{1}. (3.53)

Thence

Feq>(a+2)[ln(2+a)ln(3a)32(1a)16.F_{eq}>{(a+2)[\ln(2+a)-\ln(3a)\over 32}-{(1-a)\over 16}. (3.54)

Under the limit a0a\rightarrow 0 we find that the term log(3a)+-\log(3a)\rightarrow+\infty. So the inequality (3.54) diverges to infinity, demonstrating that the equatorial force in the Tolman-IV space-time can be made to violate the weak maximum force conjecture.

Thus, as in the case of the interior Schwarzschild solution, we have shown that the radial force is bounded (and in this case obeys both the weak and strong maximum force conjectures). However, the equatorial force can be made to diverge to infinity and act as a counter example to both weak and strong conjectures.

3.3.3 DEC

To see if the DEC is satisfied over the range of integration for the equatorial force, we inquire as to whether or not

pρ=(A2+2r2)(R2A23r2)3A4+7A2r2+6r4+(3A2+2r2)R21?\frac{p}{\rho}=\frac{(A^{2}+2r^{2})(R^{2}-A^{2}-3r^{2})}{3A^{4}+7A^{2}r^{2}+6r^{4}+(3A^{2}+2r^{2})R^{2}}\leq 1\;? (3.55)

It is straightforward to check that this inequality will always hold in the physical region. Using the definitions a=A2/R2a=A^{2}/R^{2} and z=r2/R2z=r^{2}/R^{2}, so that a(0,1)a\in(0,1), and z(0,1a3)z\in(0,{1-a\over 3}), we can write this as

pρ1=2(2a2+6az+a+6z2)(7a+2)z+3a(a+1)+6z2<0,{p\over\rho}-1=-{2(2a^{2}+6az+a+6z^{2})\over(7a+2)z+3a(a+1)+6z^{2}}<0, (3.56)

which is manifestly negative. So adding the DEC does not affect or change our conclusions. Indeed, we have already seen that the equatorial force diverges in the limit of a0a\to 0 implying A0A\rightarrow 0. Applying this limit to the ratio p/ρp/\rho gives:

limA0pρ=16r23r2+R2=16y1+3y1.\lim_{A\rightarrow 0}\frac{p}{\rho}=1-\frac{6r^{2}}{3r^{2}+R^{2}}=1-{6y\over 1+3y}\leq 1. (3.57)

Again, this is always true within any physical region, so we verify that adding the DEC does not change our conclusions.

3.3.4 Summary

For the Tolman IV solution, while the radial force is bounded (and obeys both the weak and strong maximum force conjectures), the equatorial force can be made to diverge to infinity in certain parts of parameter space (A0A\to 0) and acts as a counter-example to both weak and strong maximum force conjectures. For the Tolman IV solution, adding the DEC does not save the situation, the violation of both weak and strong maximum force conjectures is robust.

3.4 Buchdahl–Land spacetime: ρ=ρs+p\rho=\rho_{s}+p

The Buchdahl–Land spacetime is a special case of the Tolman IV spacetime, corresponding to the limit A0A\rightarrow 0 (equivalently a0a\rightarrow 0). It is sufficiently simple that it is worth some discussion in its own right. The Tolman IV metric (with a re-scaled time coordinate tAtt\rightarrow At) can be written:

ds2=(A2+r2)dt2+1+2r2A2(1r2R2)(1+r2A2)dr2+r2dΩ2.ds^{2}=-(A^{2}+r^{2})dt^{2}+\frac{1+\frac{2r^{2}}{A^{2}}}{\left(1-\frac{r^{2}}{R^{2}}\right)\left(1+\frac{r^{2}}{A^{2}}\right)}dr^{2}+r^{2}d\Omega^{2}. (3.58)

Under the limit A0A\rightarrow 0, this becomes:

ds2=r2dt2+2R2R2r2dr2+r2dΩ2.ds^{2}=-r^{2}dt^{2}+\frac{2R^{2}}{R^{2}-r^{2}}dr^{2}+r^{2}d\Omega^{2}. (3.59)

Then the orthonormal stress-energy components are:

Tt^t^\displaystyle T_{\hat{t}\hat{t}} =\displaystyle= ρ=116π(1r2+3R2);\displaystyle\rho=\frac{1}{16\pi}\left(\frac{1}{r^{2}}+\frac{3}{R^{2}}\right);
Tr^r^\displaystyle T_{\hat{r}\hat{r}} =\displaystyle= Tθ^θ^=Tϕ^ϕ^=p=116π(1r23R2).\displaystyle T_{\hat{\theta}\hat{\theta}}=T_{\hat{\phi}\hat{\phi}}=p=\frac{1}{16\pi}\left(\frac{1}{r^{2}}-\frac{3}{R^{2}}\right). (3.60)

The surface is located at

Rs=R3;withρs=38πR2=18πRs2.R_{s}=\frac{R}{\sqrt{3}};\qquad\hbox{with}\qquad\rho_{s}={3\over 8\pi R^{2}}={1\over 8\pi R_{s}^{2}}. (3.61)

At the centre the pressure and density both diverge — more on this point later.

We recast the metric as

ds2=r2dt2+2113r2Rs2dr2+r2dΩ2.ds^{2}=-r^{2}dt^{2}+{2\over 1-{1\over 3}{r^{2}\over R_{s}^{2}}}dr^{2}+r^{2}d\Omega^{2}. (3.62)

This is simply a relabelling of equation (3.59). The orthonormal stress-energy tensor is now relabelled as:

Tt^t^\displaystyle T_{\hat{t}\hat{t}} =\displaystyle= ρ=116π(1r2+1Rs2);\displaystyle\rho=\frac{1}{16\pi}\left(\frac{1}{r^{2}}+{1\over R_{s}^{2}}\right);
Tr^r^\displaystyle T_{\hat{r}\hat{r}} =\displaystyle= Tθ^θ^=Tϕ^ϕ^=p=116π(1r21Rs2).\displaystyle T_{\hat{\theta}\hat{\theta}}=T_{\hat{\phi}\hat{\phi}}=p=\frac{1}{16\pi}\left(\frac{1}{r^{2}}-{1\over R_{s}^{2}}\right). (3.63)

Note that

p=ρρs;that isρ=ρs+p.p=\rho-\rho_{s};\qquad\hbox{that is}\qquad\rho=\rho_{s}+p. (3.64)

That is, the Buchdahl–Land spacetime represents a “stiff fluid”. This perfect fluid solution has a naked singularity at r=0r=0 and a well behaved surface at finite radius. The singularity at r=0r=0 is not really a problem as one can always excise a small core region near r=0r=0 to regularize the model.

3.4.1 Radial force

Due to the simplicity of the pressure, the radial force can be easily calculated as:

Fr=14(1r2Rs2).F_{r}=\frac{1}{4}\left(1-{r^{2}\over R_{s}^{2}}\right). (3.65)

The radial force is trivially bounded with a maximum of 14\frac{1}{4} at the centre of the star. This obeys the strong (and so also the weak) maximum force conjecture.

3.4.2 Equatorial force

The equatorial force is:

Feq\displaystyle F_{eq} =\displaystyle= 2π0Rs2113r2Rs2116π(1r21Rs2)r𝑑r.\displaystyle 2\pi\hbox{\Large$\displaystyle\int$}^{{R_{s}}}_{\!\!\!\!\!0}\sqrt{2\over 1-{1\over 3}{r^{2}\over R_{s}^{2}}}\;\;\frac{1}{16\pi}\left(\frac{1}{r^{2}}-{1\over R_{s}^{2}}\right)\;rdr. (3.66)

This is now simple enough to handle analytically. Using the dimensionless variable y=r2/Rs2y=r^{2}/R_{s}^{2}, with range y(0,1)y\in(0,1), we see:

Feq\displaystyle F_{eq} =\displaystyle= 116012113y(1y)dyy.\displaystyle{1\over 16}\hbox{\Large$\displaystyle\int$}^{1}_{\!\!\!\!\!0}\sqrt{2\over 1-{1\over 3}{y}}\;\;\left(1-y\right)\;{dy\over y}. (3.67)

This is manifestly dimensionless, and manifestly diverges to ++\infty. If we excise a small region r<rcorer<r_{\mathrm{core}}, (corresponding to y<ycorey<y_{\mathrm{core}}) to regularize the model, replacing r<rcorer<r_{\mathrm{core}} with some well-behaved fluid ball, then we have the explicit logarithmic divergence

Feq=116lnycore+𝒪(1).F_{eq}=-{1\over 16}\ln y_{\mathrm{core}}+{\mathcal{O}}(1). (3.68)

This violates the weak (and so also the strong) maximum force conjecture.

3.4.3 DEC

The DEC for this space-time is given by:

pρ=ρρsρ=1ρsρ1.\displaystyle\frac{p}{\rho}={\rho-\rho_{s}\over\rho}=1-{\rho_{s}\over\rho}\leq 1. (3.69)

which is always true for positive values of rr, ρs\rho_{s}.

3.4.4 Summary

The Buchdahl–Land spacetime is another weak maximum force conjecture counter-example, one which again obeys the classical energy conditions.

3.5 Scaling solution

The scaling solution is

ds2=r4w1+wdt2+(w2+6w+1(1+w)2)dr2+r2dΩ2.ds^{2}=-r^{\frac{4w}{1+w}}dt^{2}+\left(\frac{w^{2}+6w+1}{(1+w)^{2}}\right)dr^{2}+r^{2}d\Omega^{2}. (3.70)

This produces the following stress energy tensor:

Tt^t^\displaystyle T_{\hat{t}\hat{t}} =\displaystyle= ρ=w2π(w2+6w+1)r2;\displaystyle\rho=\frac{w}{2\pi(w^{2}+6w+1)r^{2}};
Tr^r^\displaystyle T_{\hat{r}\hat{r}} =\displaystyle= Tθ^θ^=Tϕ^ϕ^=p=w22π(w2+6w+1)r2.\displaystyle T_{\hat{\theta}\hat{\theta}}=T_{\hat{\phi}\hat{\phi}}=p=\frac{w^{2}}{2\pi(w^{2}+6w+1)r^{2}}. (3.71)

This perfect fluid solution has a naked singularity at r=0r=0 and does not have a finite surface — it requires rr\rightarrow\infty for the pressure to vanish. Nevertheless, apart from a small region near r=0r=0 and small fringe region near the surface r=Rsr=R_{s}, this is a good approximation to the bulk geometry of a star that is on the verge of collapse [34, 35]. To regularize the model excise two small regions, a core region at r(0,rcore)r\in(0,r_{\mathrm{core}}), and an outer shell at r(rfringe,Rs)r\in(r_{\mathrm{fringe}},R_{s}), replacing them by segments of well-behaved fluid spheres. Note that for r(rcore,rfringe)r\in(r_{\mathrm{core}},r_{\mathrm{fringe}}) we have p/ρ=wp/\rho=w, (and since ρ>0\rho>0 we must have w>0w>0), so the DEC implies w(0,1]w\in(0,1].

3.5.1 Radial force

Using equation (2.2), we find that the radial force is very simply given by:

Fr=2w2w2+6w+1.F_{r}=\frac{2w^{2}}{w^{2}+6w+1}. (3.72)

This is independent of rr and attains a maximum value of 14\frac{1}{4} when w=1w=1, giving a bounded force obeying the strong maximum force conjecture.

3.5.2 Equatorial force

Now, using equation (3.2), the equatorial force can be calculated as:

Feq=rcorerfringew2+6w+1(1+w)2(w2(w2+6w+1)r)dr+𝒪(1).F_{eq}=\hbox{\Large$\displaystyle\int$}^{r_{\mathrm{fringe}}}_{\!\!\!\!\!r_{\mathrm{core}}}\sqrt{\frac{w^{2}+6w+1}{(1+w)^{2}}}\left(\frac{w^{2}}{(w^{2}+6w+1)r}\right)dr+{\mathcal{O}}(1). (3.73)

That is

Feq=w2(1+w)w2+6w+1ln(rfringe/rcore)+𝒪(1),F_{eq}=\frac{w^{2}}{(1+w)\sqrt{w^{2}+6w+1}}\;\ln(r_{\mathrm{fringe}}/r_{\mathrm{core}})+{\mathcal{O}}(1), (3.74)

which trivially diverges logarithmically as either rcore0r_{\mathrm{core}}\to 0 or rfringer_{\mathrm{fringe}}\to\infty, providing a counter-example to weak maximum force conjecture.

3.5.3 Summary

Again we have an explicit model where the radial force FrF_{r} is well-behaved, but the equatorial force FeqF_{eq} provides an explicit counter-example to weak maximum force conjecture. This counter-example is compatible with the DEC.

3.6 TOV equation

Let us now see how far we can push this sort of argument using only the TOV equation for the pressure profile in perfect fluid spheres — we will (as far as possible) try to avoid making specific assumptions on the metric components and stress-energy. The TOV equation is

dp(r)dr={ρ(r)+p(r)}{m(r)+4πp(r)r3}r2{12m(r)/r}.{dp(r)\over dr}=-{\{\rho(r)+p(r)\}\{m(r)+4\pi p(r)r^{3}\}\over r^{2}\{1-2m(r)/r\}}. (3.75)

3.6.1 Radial force

From the definition of radial force Fr=4πpr2F_{r}=4\pi pr^{2}, we see that at the maximum of FrF_{r} we must have

(2pr+r2p)|rmax=0.\left.(2pr+r^{2}p^{\prime})\right|_{r_{\mathrm{max}}}=0. (3.76)

Thence, at the maximum

(Fr)max=(4πpr2)max=2π(r3p)|rmax.(F_{r})_{\mathrm{max}}=(4\pi pr^{2})_{\mathrm{max}}=-2\pi\left.(r^{3}p^{\prime})\right|_{r_{\mathrm{max}}}. (3.77)

In particular, now using the TOV at the location rmaxr_{\mathrm{max}} of the maximum of FrF_{r}:

(Fr)max=2π[(ρ+p)r(m+4πpr3)(12m/r)]|rmax.(F_{r})_{\mathrm{max}}=2\pi\left.\left[{(\rho+p)r(m+4\pi pr^{3})\over(1-2m/r)}\right]\right|_{r_{\mathrm{max}}}. (3.78)

Let us define the two parameters

χ=[2m(r)r]rmax=2m(rmax)rmax,andw=[p(r)ρ(r)]rmax=p(rmax)ρ(rmax).\chi=\left[2m(r)\over r\right]_{r_{\mathrm{max}}}={2m(r_{\mathrm{max}})\over r_{\mathrm{max}}},\quad\hbox{and}\quad w=\left[p(r)\over\rho(r)\right]_{r_{\mathrm{max}}}={p(r_{max})\over\rho(r_{max})}. (3.79)

Then

(Fr)max=12(4πp(1+1w)r2(χ/2+4πpr2)1χ)|rmax.(F_{r})_{\mathrm{max}}={1\over 2}\left.\left({4\pi p(1+{1\over w})r^{2}(\chi/2+4\pi pr^{2})\over 1-\chi}\right)\right|_{r_{\mathrm{max}}}. (3.80)

Simplifying, we see:

(Fr)max=12(Fr)max[1+1/w][(Fr)max+χ/2)1χ(F_{r})_{\mathrm{max}}={1\over 2}{(F_{r})_{\mathrm{max}}\;[1+1/w]\;[(F_{r})_{\mathrm{max}}+\chi/2)\over 1-\chi} (3.81)

Discarding the unphysical solution (Fr)max=0(F_{r})_{\mathrm{max}}=0, we find

(Fr)max=4wχ5wχ2(1+w)=2w1+wχ1+5w1+w.(F_{r})_{\mathrm{max}}={4w-\chi-5w\chi\over 2(1+w)}={2w\over 1+w}-\chi\;{1+5w\over 1+w}. (3.82)

The physical region corresponds to 0χ<10\leq\chi<1, while w>0w>0. Furthermore we have (Fr)max>0(F_{r})_{\mathrm{max}}>0, whence 4wχ5wχ>04w-\chi-5w\chi>0, implying χ<4w/(1+5w)<4/5\chi<4w/(1+5w)<4/5. That is, at the location rmaxr_{\mathrm{max}} of the maximum of FrF_{r} we have

[2m(r)r]rmax=2m(rmax)rmax<45.\left[2m(r)\over r\right]_{r_{\mathrm{max}}}={2m(r_{\mathrm{max}})\over r_{\mathrm{max}}}<{4\over 5}. (3.83)

This is not the Buchdahl–Bondi bound, it is instead a bound on the compactness of the fluid sphere at the internal location rmaxr_{\mathrm{max}} where FrF_{r} is maximized.

Observe that (Fr)max(F_{r})_{\mathrm{max}} is maximized when χ=0\chi=0 and w=w=\infty, when (Fr)max2(F_{r})_{\mathrm{max}}\to 2. This violates the strong conjecture maximum force but not the weak maximum force conjecture. If we impose the DEC then w1w\leq 1, and (Fr)max(F_{r})_{\mathrm{max}} is maximized when χ=0\chi=0 and w=1w=1, when (Fr)max1(F_{r})_{\mathrm{max}}\to 1. This still violates the strong maximum force conjecture but not the weak maximum force conjecture. Consequently the weak conjecture for FrF_{r} generically holds for any prefect fluid sphere satisfying the TOV.

3.6.2 Equatorial force

As we have by now come to expect, dealing with the equatorial force will be considerably trickier. In view of the non-negativity of the Misner–Sharp quasi-local mass we have:

Feq=2π0Rsgrrpr𝑑r=2π0Rs112m(r)/rpr𝑑r>2π0Rspr𝑑r.F_{eq}=2\pi\int_{0}^{R_{s}}\sqrt{g_{rr}}\;p\;rdr=2\pi\int_{0}^{R_{s}}{1\over\sqrt{1-2m(r)/r}}\;p\;rdr>2\pi\int_{0}^{R_{s}}p\,r\;dr. (3.84)

To make the integral 0Rspr𝑑r\int_{0}^{R_{s}}p\,r\;dr converge it is sufficient to demand p(r)=o(1/r2)p(r)=o(1/r^{2}). However, for stars on the verge of gravitational collapse it is known that p(r)K/r2p(r)\sim K/r^{2}, see for instance [34, 35]. More specifically, there is some core region r(0,rcore)r\in(0,r_{\mathrm{core}}) designed to keep the central pressure finite but arbitrarily large, a large scaling region r(rcore,rfringe)r\in(r_{\mathrm{core}},r_{\mathrm{fringe}}) where pK/r2p\sim K/r^{2}, and an outer fringe r(rfringe,Rs)r\in(r_{\mathrm{fringe}},R_{s}) where one has p(r)p(Rs)=0p(r)\to p(R_{s})=0. Then we have the identity

0Rspr𝑑r=0rcorepr𝑑r+rcorerfringepr𝑑r+rfringeRspr𝑑r.\int_{0}^{R_{s}}p\,r\;dr=\int_{0}^{r_{\mathrm{core}}}p\,r\;dr+\int_{r_{\mathrm{core}}}^{r_{\mathrm{fringe}}}p\,r\;dr+\int_{r_{\mathrm{fringe}}}^{R_{s}}p\,r\;dr. (3.85)

But under the assumed conditions this implies

0Rspr𝑑r=𝒪(1)+[rcorerfringeKr𝑑r+𝒪(1)]+𝒪(1).\int_{0}^{R_{s}}p\,r\;dr={\mathcal{O}}(1)+\left[\int_{r_{\mathrm{core}}}^{r_{\mathrm{fringe}}}{K\over r}\;dr+{\mathcal{O}}(1)\right]+{\mathcal{O}}(1). (3.86)

Thence

0Rspr𝑑r=Kln(rfringe/rcore)+𝒪(1).\int_{0}^{R_{s}}p\,r\;dr=K\ln\left({r_{\mathrm{fringe}}}/r_{\mathrm{core}}\right)+{\mathcal{O}}(1). (3.87)

Finally

Feq>2πKln(rfringe/rcore)+𝒪(1).F_{eq}>2\pi K\ln\left({r_{\mathrm{fringe}}}/r_{\mathrm{core}}\right)+{\mathcal{O}}(1). (3.88)

This can be made arbitrarily large for a star on the verge of gravitational collapse, so the weak and strong maximum force conjectures are both violated.

Note that technical aspects of the argument are very similar to what we saw for the exact scaling solution to the Einstein equations, but the physical context is now much more general.

3.6.3 Summary

We see that the weak maximum force conjecture generically holds for the radial force FrF_{r} when considering perfect fluid spheres satisfying the TOV. In contrast we see that the weak maximum force conjecture fails for the equatorial force FeqF_{eq} when considering perfect fluid spheres satisfying the TOV that are close to gravitational collapse.

4 Discussion

With the notion a natural unit of force F=FPlanck=c4/GNF_{*}=F_{\mathrm{Planck}}=c^{4}/G_{N} in hand, one can similarly define a natural unit of power [36, 37, 38, 39, 40]

P=PPlanck=c5GN=1 Dyson3.6×1052 W,P_{*}=P_{\mathrm{Planck}}={c^{5}\over G_{N}}=1\hbox{ Dyson}\approx 3.6\times 10^{52}\hbox{ W}, (4.1)

a natural unit of mass-loss-rate

(m˙)=(m˙)Planck=c3GN4.0×1035 kg/s,(\dot{m})_{*}=(\dot{m})_{\mathrm{Planck}}={c^{3}\over G_{N}}\approx 4.0\times 10^{35}\hbox{ kg/s}, (4.2)

and even a natural unit of mass-per-unit-length

(m)=(m)Planck=c2GN1.36×1027 kg/m.(m^{\prime})_{*}=(m^{\prime})_{\mathrm{Planck}}={c^{2}\over G_{N}}\approx 1.36\times 10^{27}\hbox{ kg/m}. (4.3)

Despite being Planck units, all these concepts are purely classical (the various factors of \hbar cancel, at least in (3+1) dimensions).

Indeed, consider the classical Stoney units which pre-date Planck units by some 20 years [41, 42, 43], and use GNG_{N}, cc, and Coulomb’s constant e24πϵ0e^{2}\over 4\pi\epsilon_{0}, instead of GNG_{N}, cc, and Planck’s constant \hbar. Then we have F=FPlanck=FStoneyF_{*}=F_{\mathrm{Planck}}=F_{\mathrm{Stoney}}. Similarly we have P=PPlanck=PStoneyP_{*}=P_{\mathrm{Planck}}=P_{\mathrm{Stoney}}, (m˙)=(m˙)Planck=(m˙)Stoney(\dot{m})_{*}=(\dot{m})_{\mathrm{Planck}}=(\dot{m})_{\mathrm{Stoney}}, and (m)=(m)Planck=(m)Stoney(m^{\prime})_{*}=(m^{\prime})_{\mathrm{Planck}}=(m^{\prime})_{\mathrm{Stoney}}. Based ultimately on dimensional analysis, any one of these quantities might be used to advocate for a maximality conjecture: maximum luminosity [36, 37, 38, 39, 40], maximum mass-loss-rate, or maximum mass-per-unit-length. The specific counter-examples to the maximum force conjecture that we have discussed above suggest that it might also be worth looking for specific counter-examples to these other conjectures [39].

5 Conclusions

Through the analysis of radial and equatorial forces within perfect fluid spheres in general relativity, we have produced a number of counter-examples to both the strong and weak forms of the maximum force conjecture. These counter-examples highlight significant issues with the current phrasing and understanding of this conjecture, as merely specifying that forces are bounded within the framework of general relativity is manifestly a falsehood. As such, should one wish some version of the maximum force conjecture to be considered viable as a potential physical principle, it must be very clearly specified as to what types of forces they pertain to.

Acknowledgments

MV was supported by the Marsden Fund, via a grant administered by the Royal Society of New Zealand.

References

  • [1] G. W. Gibbons, “The Maximum Tension Principle in General Relativity”,
    Foundations of Physics 32 (2002), 1891
    doi: 10.1023/A:1022370717626. [arXiv:hep-th/0210109 [hep-th]].
  • [2] C. Schiller, “Motion Mountain — A Hike Beyond Space and Time Along the Concepts of Modern Physics” (http://www.motionmountain.net, 1997-2004),
    section 7: Maximum force — a simple principle encompassing general relativity.
  • [3] J. D. Barrow and G. W. Gibbons,
    “Maximum Tension: with and without a cosmological constant”,
    Mon. Not. Roy. Astron. Soc. 446 (2015) no.4, 3874–3877
    doi: 10.1093/mnras/stu2378 [arXiv:1408.1820 [gr-qc]].
  • [4] J. D. Barrow and N. Dadhich, “Maximum Force in Modified Gravity Theories”,
    Phys. Rev. D 102 (2020) no.6, 064018 doi: 10.1103/PhysRevD.102.064018 [arXiv:2006.07338 [gr-qc]].
  • [5] Y.C. Ong,
    “GUP-Corrected Black Hole Thermodynamics and the Maximum Force Conjecture”,
    Physics Letters B 785 (2018), 217
    doi: 10.1016/j.physletb.2018.08.065 [arXiv:1809.00442 [gr-qc]].
  • [6] C. Schiller, “General Relativity and Cosmology Derived From Principle of Maximum Power or Force”, Int J Theor Phys 44 (2005), 1629 doi: 10.1007/s10773-005-4835-2 [arXiv:physics/0607090 [physics.gen-ph]].
  • [7] T. Jacobson, “Thermodynamics of space-time: The Einstein equation of state”,
    Phys. Rev. Lett. 75 (1995), 1260-1263 doi: 10.1103/PhysRevLett.75.1260 [arXiv:gr-qc/9504004 [gr-qc]].
  • [8] M. P. Dabrowski and H. Gohar, “Abolishing the maximum tension principle”,
    Phys. Lett. B 748 (2015), 428-431 doi: 10.1016/j.physletb.2015.07.047 [arXiv:1504.01547 [gr-qc]].
  • [9] Y. L. Bolotin and V. A. Cherkaskiy, “Principle of maximum force and holographic principle: two principles or one?”, [arXiv:1507.02839 [gr-qc]].
  • [10] M. S. R. Delgaty and K. Lake, “Physical acceptability of isolated, static, spherically symmetric, perfect fluid solutions of Einstein’s equations”, Comput. Phys. Commun. 115 (1998) 395 doi: 10.1016/S0010-4655(98)00130-1 [arXiv:gr-qc/9809013 [gr-qc]]
  • [11] S. Rahman and M. Visser, “Space-time geometry of static fluid spheres”,
    Class. Quant. Grav. 19 (2002), 935-952 doi: 10.1088/0264-9381/19/5/307 [arXiv:gr-qc/0103065 [gr-qc]].
  • [12] D. Martin and M. Visser, “Algorithmic construction of static perfect fluid spheres”, Phys. Rev. D 69 (2004), 104028 doi: 10.1103/PhysRevD.69.104028 [arXiv:gr-qc/0306109 [gr-qc]].
  • [13] D. Martin and M. Visser,
    “Bounds on the interior geometry and pressure profile of static fluid spheres”,
    Class. Quant. Grav. 20 (2003), 3699-3716 doi: 10.1088/0264-9381/20/16/311 [arXiv:gr-qc/0306038 [gr-qc]].
  • [14] P. Boonserm, M. Visser and S. Weinfurtner,
    “Generating perfect fluid spheres in general relativity”,
    Phys. Rev. D 71 (2005), 124037 doi: 10.1103/PhysRevD.71.124037 [arXiv:gr-qc/0503007 [gr-qc]].
  • [15] P. Boonserm, M. Visser and S. Weinfurtner,
    “Solution generating theorems for perfect fluid spheres”,
    J. Phys. Conf. Ser. 68 (2007), 012055 doi: 10.1088/1742-6596/68/1/012055 [arXiv:gr-qc/0609088 [gr-qc]].
  • [16] P. Boonserm, M. Visser and S. Weinfurtner,
    “Solution generating theorems: Perfect fluid spheres and the TOV equation”,
    doi: 10.1142/9789812834300_0388 [arXiv:gr-qc/0609099 [gr-qc]].
  • [17] P. Boonserm and M. Visser, “Buchdahl-like transformations in general relativity”, Thai Journal of Mathematics 5 # 2 (2007) 209–223.
  • [18] P. Boonserm and M. Visser, “Buchdahl-like transformations for perfect fluid spheres”, Int. J. Mod. Phys. D 17 (2008), 135-163 doi: 10.1142/S0218271808011912 [arXiv:0707.0146 [gr-qc]].
  • [19] See https://en.wikipedia.org/wiki/Energy_condition.
  • [20] M. Visser, “Lorentzian wormholes: From Einstein to Hawking”,
    (AIP press — now Springer Verlag, Woodbury, USA, 1995).
  • [21] C. Barceló and M. Visser, “Twilight for the energy conditions?”,
    Int. J. Mod. Phys. D 11 (2002), 1553-1560 doi:10.1142/S0218271802002888 [arXiv:gr-qc/0205066 [gr-qc]].
  • [22] M. Visser, “Energy conditions in the epoch of galaxy formation”,
    Science 276 (1997), 88-90 doi:10.1126/science.276.5309.88 [arXiv:1501.01619 [gr-qc]].
  • [23] M. Visser,
    “Gravitational vacuum polarization. 2: Energy conditions in the Boulware vacuum”, Phys. Rev. D 54 (1996), 5116-5122 doi:10.1103/PhysRevD.54.5116 [arXiv:gr-qc/9604008 [gr-qc]].
  • [24] C. Cattöen and M. Visser, “Cosmological milestones and energy conditions”,
    J. Phys. Conf. Ser. 68 (2007), 012011 doi:10.1088/1742-6596/68/1/012011 [arXiv:gr-qc/0609064 [gr-qc]].
  • [25] P. Martín-Moruno and M. Visser, “Classical and semi-classical energy conditions”, Fundam. Theor. Phys. 189 (2017) 193–213, (formerly Lecture Notes in Physics), doi:10.1007/978-3-319-55182-1_9 [arXiv:1702.05915 [gr-qc]].
  • [26] P. Martín-Moruno and M. Visser,
    “Semiclassical energy conditions for quantum vacuum states”,
    JHEP 09 (2013), 050 doi:10.1007/JHEP09(2013)050 [arXiv:1306.2076 [gr-qc]].
  • [27] P. Martín-Moruno and M. Visser,
    “Classical and quantum flux energy conditions for quantum vacuum states”,
    Phys. Rev. D 88 (2013) no.6, 061701 doi:10.1103/PhysRevD.88.061701 [arXiv:1305.1993 [gr-qc]].
  • [28] P. Martín-Moruno and M. Visser, “Semi-classical and nonlinear energy conditions”, Proceedings of the 14th Marcel Grossmann Meeting, Rome, 12–18 July 2015. doi:10.1142/9789813226609_0126 [arXiv:1510.00158 [gr-qc]].
  • [29] K. Schwarzschild, “Über das Gravitationsfeld einer Kugel aus inkompressibler Flüssigkeit nach der Einsteinschen Theorie”, (On the gravitational field of a sphere of incompressible fluid according to Einstein’s theory), Sitzungsberichte der Königlich-Preussischen Akademie der Wissenschaften, Berlin, 7(1916) 424–434.
  • [30] K. Schwarzschild, “Über das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie”, (On the gravitational field of a mass point according to Einstein’s theory), Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften, Berlin, 7 (1916) 189–196.
  • [31] H. A. Buchdahl, “General Relativistic Fluid Spheres”, Phys. Rev. 116 (1959), 1027 doi: 10.1103/PhysRev.116.1027
  • [32] H. Bondi, “Massive spheres in general relativity”,
    Proc. Roy. Soc. Lond. A 282 (1964) 303-317 doi: 10.1098/rspa.1964.0234
  • [33] H. A. Buchdahl and W. Land, “The relativistic incompressible sphere”,
    Journal of the Australian Mathematical Society, 8(1) (1968),6
    doi: 10.1017/S1446788700004559.
  • [34] B. Kent Harrison, Kip S. Thorne, Masami Wakano, and John Archibald Wheeler,
    Gravitation theory and gravitational collapse,
    (University of Chicago Press, Chicago, 1965)
  • [35] M. Visser and N. Yunes, “Power laws, scale invariance, and generalized Frobenius series: Applications to Newtonian and TOV stars near criticality”,
    Int. J. Mod. Phys. A 18 (2003), 3433-3468 doi:10.1142/S0217751X03013892 [arXiv:gr-qc/0211001 [gr-qc]].
  • [36] Freeman Dyson, “Gravitational machines”,
    Chapter 12 in Interstellar Communication,
    A. G. W. Cameron, editor, New York: Benjamin Press, 1963.
    https://www.ifa.hawaii.edu/~barnes/ast242_s14/Dyson_Machines.pdf
    Also awarded 4th prize in the 1962 Gravity Research Foundation essay contest:
    https://www.gravityresearchfoundation.org/year https://static1.squarespace.com/static/5852e579be659442a01f27b8/t/5873d777e4fcb5de501305f7/1483986808648/dyson.pdf
  • [37] J. D. Barrow and G. W. Gibbons,
    “Maximum magnetic moment to angular momentum conjecture”,
    Phys. Rev. D 95 (2017) no.6, 064040 doi:10.1103/PhysRevD.95.064040 [arXiv:1701.06343 [gr-qc]].
  • [38] Footnote 5 in reference [37]: In reply to an enquiry by Christoph Schiller, Dyson replied on 14th Feb 2011: “It is not true that I proposed the formula c5/Gc^{5}/G as a luminosity limit for anything. I make no such claim. Perhaps this notion arose from a paper that I wrote in 1962 with the title, “Gravitational Machines”, published as Chapter 12 in the book, “Interstellar Communication” edited by Alastair Cameron, [New York, Benjamin, 1963]. Equation (11) in that paper is the well-known formula 128V10/5Gc5128V^{10}/5Gc^{5} for the power in gravitational waves emitted by a binary star with two equal masses moving in a circular orbit with velocity VV. As VV approaches its upper limit cc, this gravitational power approaches the upper limit 128c5/5G128c^{5}/5G. The remarkable thing about this upper limit is that it is independent of the masses of the stars. It may be of some relevance to the theory of gamma-ray bursts.”
  • [39] V. Cardoso, T. Ikeda, C. J. Moore and C. M. Yoo,
    “Remarks on the maximum luminosity”,
    Phys. Rev. D 97 (2018) no.8, 084013 doi: 10.1103/PhysRevD.97.084013 [arXiv:1803.03271 [gr-qc]].
  • [40] J. P. Bruneton, “Notes on several phenomenological laws of quantum gravity”, [arXiv:1308.4044 [gr-qc]].
  • [41] See https://en.wikipedia.org/wiki/Stoney_units
  • [42] G. Johnstone Stoney, “On The Physical Units of Nature”,
    Phil.Mag. 11 (1881) 381–391.
  • [43] G. Johnstone Stoney, “On The Physical Units of Nature”,
    Scientific Proceedings of the Royal Dublin Society, 3 (1883) 51–60.