The QCD mass gap and quark deconfinement scales as mass bounds in strong gravity

TL;DR

This paper explores a strong gravity model inspired by Salam's theory to understand confinement, mass gaps, and deconfinement scales in QCD, deriving bounds on hadron properties and linking them to a strong cosmological constant.

Contribution

It introduces a modified strong gravity framework with a mixing term that models confinement and mass gaps, providing bounds on hadron stability and phase transitions.

Findings

Identifies a strong cosmological constant as the bag constant in QCD.

Derives mass/radius bounds for stable hadronic objects.

Links deconfinement to an upper mass/radius bound in strong gravity.

Abstract

Though not a part of mainstream physics, Salam's theory of strong gravity remains a viable effective model for the description of strong interactions in the gauge singlet sector of QCD, capable of producing particle confinement and asymptotic freedom, but not of reproducing interactions involving $SU(3)$ colour charge. It may therefore be used to explore the stability and confinement of gauge singlet hadrons, though not to describe scattering processes that require colour interactions. It is a two-tensor theory of both strong interactions and gravity, in which the strong tensor field is governed by equations formally identical to the Einstein equations, apart from the coupling parameter, which is of order $1 \ \rm GeV^{-1}$. We revisit the strong gravity theory and investigate the strong gravity field equations in the presence of a mixing term which induces an effective {\it strong cosmological constant}, $\Lambda_{f}$. This introduces a {\it strong de Sitter radius} for strongly interacting fermions, producing a confining bubble, which allows us to identify $\Lambda_{f}$ with the `bag constant' of the MIT bag model, $B \simeq 2 \times 10^{14} \rm gcm^{-3}$. Assuming a static, spherically symmetric geometry, we derive the strong gravity TOV equation, which describes the equilibrium properties of compact hadronic objects. From this, we determine the generalised Buchdahl inequalities for a strong gravity `particle', giving rise to upper and lower bounds on the mass/radius ratio of stable, compact, strongly interacting objects. We show, explicitly, that the existence of the lower mass bound is induced by the presence of $\Lambda_f$, producing a mass gap, and that the upper bound corresponds to a deconfinement phase transition. The physical implications of our results for holographic duality in the context of the AdS/QCD and dS/QCD correspondences are also discussed.

Equations

$$k_g^2=8\pi G/c^4,k_f^2=8\pi G_f/c^4,$$

$${\alpha }_s=\frac{G_f{m}^2}{\hslash c}\simeq 1.$$

$$G_f/G\simeq 10^{38},$$

$$\mathcal{L}=-\frac{1}{2k_g^2}\sqrt{-g}R(g)-\frac{1}{2k_f^2}\sqrt{-f}{R}_f(f)+{\mathcal{L}}_{fg}+{\mathcal{L}}_m$$

$$\frac{2GM}{c^{2}R}\le \frac{8}{9},$$

$$\sqrt{1-\frac{2GM}{R{c}^2}-\frac{\Lambda R^2}{3}}\ge \frac{1}{3}-\frac{\Lambda c^2}{12\pi G\rho },$$

$$\frac{2GM}{R{c}^2}\ge \frac{\Lambda }{6}{R}^2⟺\rho =\frac{3M}{4\pi R^3}\ge {\rho }_{\min }:=\frac{\Lambda c^2}{16\pi G}.$$

$${\rho }_{\Lambda }=-p_{\Lambda }/c=\frac{\Lambda c^2}{8\pi G}.$$

$$\frac{2M{c}^2}{R}\ge \frac{3}{2}\frac{Q^2}{R^2}\left[1-\frac{Q^2}{36R^2}+\mathcal{O}{\left(\frac{Q^2}{R^2}\right)}^4\right].$$

$$M\ge \frac{3}{4}\frac{Q^2}{R{c}^2}+\frac{\Lambda R^3{c}^2}{12G}.$$

$$M\gtrsim M_{\mathrm{Q}}=2^{-1/3}{\alpha }_{\mathrm{Q}}(m_{\mathrm{Pl}}^2{m}_{\mathrm{dS}}{)}^{1/3},$$

$${\alpha }_{\mathrm{Q}}=Q^2/q_{\mathrm{Pl}}^2,$$

$$M$$

$$\frac{Q^2}{M}$$

$$\Lambda \lesssim \frac{l_{\mathrm{Pl}}^4}{r_e^6}=\frac{3{\hslash }^2{G}^2{m}_e^6{c}^6}{e^{12}}\simeq 1.366\times 10^{-56}{\mathrm{cm}}^{-2},$$

$$\mathcal{L}=-\frac{1}{2k_g^2}\sqrt{-g}R(g)-\frac{1}{2k_f^2}\sqrt{-f}{R}_f(f)+{\mathcal{L}}_{fg}+{\mathcal{L}}_m$$

$${\mathcal{L}}_{fg}$$

$$R_{\mu \nu }(f)-\frac{1}{2}{f}_{\mu \nu }R(f)=k_f^2{T}_{\mu \nu }^{(s)},$$

$$k_f^2{T}_{\mu \nu }^{(s)}=-\frac{1}{2}{\mathcal{M}}^2\left(f^{\kappa \lambda }-{\eta }^{\kappa \lambda }\right)\left({\eta }_{\kappa \nu }{\eta }_{\lambda \nu }-{\eta }_{\mu \nu }{\eta }_{\kappa \lambda }\right)\frac{\sqrt{-\eta }}{\sqrt{-f}}.$$

$$d{q}^2=f_{\mu \nu }d{x}^{\mu }d{x}^{\nu }=e^{\nu (r)}{c}^{2}d{t}^2-e^{\lambda (r)}d{r}^2-\Sigma (r)d{\Omega }_2,$$

$$T_{\nu }^{\mu }=\text{diag}(\rho c^2,-P_r,-P_{\perp },-P_{\perp }),$$

$$R_{\nu }^{\mu }(f)-\frac{1}{2}{f}_{\nu }^{\mu }R(f)=k_f^2(T_{\nu }^{\mu (s)}+T_{\nu }^{\mu }),$$

$$f_{\mu \nu }{f}^{\nu \lambda }={\delta }_{\mu }^{\lambda }.$$

$$k_f^2{T}_t^{t(s)}$$

$$k_f^2{T}_r^{r(s)}$$

$$k_f^2{T}_{\theta }^{\theta (s)}$$

$$\Sigma =\frac{2r^2}{3},e^{(\nu +\lambda )}=\Delta =\mathrm{constant}>0.$$

$$P_r+\rho c^2=0,$$

$$k_f^2{T}_{\theta }^{\theta (s)}=\frac{9}{8}\frac{{\mathcal{M}}^2}{\sqrt{\Delta }}\left(\frac{3}{2}-(e^{-\nu }+e^{-\lambda })\right).$$

$${\nu }^{\prime }+{\lambda }^{\prime }=0,$$

$$\frac{e^{-\lambda }{\lambda }^{\prime }}{r}+\frac{(\frac{3}{2}-e^{-\lambda })}{r^2}$$

Full text

Keywords: strong interaction, strong gravity, mass gap, confinement/deconfinement, mass bounds, holography, bag model, AdS/QCD, dS/QCD, Dirac Large Number Hypothesis

1. I Introduction

2. II Tolman-Oppenheimer-Volkoff Equation in the Strong Gravity Model

3. II.1 Strong gravity

4. II.2 Strong gravity field equations

5. II.3 The Tolman-Oppenheimer-Volkoff equation in strong gravity

6. III The Buchdahl inequality and the minimum and maximum mass/radius ratios of stable compact objects in strong gravity

7. III.1 The Buchdahl inequality

8. III.2 The minimum and maximum mass/radius ratios of hadrons

9. III.2.1 The upper bound on the mass/radius ratio of hadrons

10. III.2.2 The lower bound on the mass/radius ratio of hadrons

11. III.3 The energy localization problem in strong gravity

12. IV The mass gap in strong gravity

13. V Quantum mechanical implications of the classical mass/radius ratio bounds in strong gravity

14. V.1 Quantum mass bounds in standard general relativity

15. V.2 Quantum mass bounds in strong gravity

16. V.2.1 Mass bounds for neutral particles

17. V.2.2 Mass bounds for charged particles

18. V.3 The Hagedorn temperature for minimum-mass particles, the expanding Universe, and deconfinement

19. VI Discussions and final remarks

I Introduction

One of the most intriguing aspects of short-distance physics is that the strong interactions of hadrons in the infrared (IR) regime exhibit certain features bearing a close resemblance to gravity. For example, string theory – originally proposed as a theory of strongly interacting hadrons – can be reinterpreted as a theory of linearised gravity, and the quantum theory of closed bosonic strings naturally includes a candidate graviton Green:2012pqa . On the other hand, it is generally assumed that, in low energy physics, the gravitational interaction plays a negligible role. However, it is important to note that the strength of the gravitational interaction increases with energy, the coupling being proportional to $GE^{2}$, where $E$ is the total energy of the particle (including rest mass) and $G$ is Newton’s constant. Hence, gravitational interactions become more and more important at higher energies. In fact, if the particle energy exceeds $E=ec^{2}/\sqrt{G}\simeq 10^{18}$ GeV, the gravitational interaction is stronger than the electromagnetic interaction and, at energies of the order of $10^{19}$ GeV, it is as strong as the strong nuclear interaction. The important role of gravitation in fundamental particle physics was also pointed out in 1 ; 2 .

At the quantum level, the gravitational interaction is expected to be mediated by massless spin-2 bosons (gravitons) in a way that is analogous to the mediation of the electromagnetic interaction by massless spin-1 bosons (photons) 3 . Strong support for this idea is provided by the fact that the quantisation of linearised gravity, obtained by substituting $g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}$, where $\eta_{\mu\nu}$ is the Minkowski metric and $h_{\mu\nu}$ is an arbitrary perturbation, into the vacuum Einstein equations (neglecting second order quantities), leads to the well-known Pauli-Fierz equations for massless spin-2 particles, $\Box h_{\mu\nu}=0$, $h_{\mu\nu,\nu}=0$, $h_{\mu}^{\mu}=0$ 3 . Thus, from a quantum theoretical point of view, the long range of the gravitational force is a consequence of the masslessness of the mediating particles. This result can be extended to the massive graviton case, and hence it follows that a non-linear self-interacting spin-2 field can also be described by Einstein’s field equations 3a ; 3b ; 3c .

Classically, the (non-vacuum) Einstein field equations, $G_{\mu\nu}=\kappa T_{\mu\nu}$, where $\kappa=8\pi G/c^{4}$, relate a covariant geometrical quantity, the Einstein tensor $G_{\mu\nu}=R_{\mu\nu}-(1/2)Rg_{\mu\nu}$, to a covariant physical quantity, the conserved energy-momentum tensor $T_{\mu\nu}$, via the proportionality (coupling) constant $\kappa$ 4 . However, the derivation of the Einstein field equations does not place any restriction whatsoever on the numerical value of the constant $\kappa$. For the canonical gravitational interaction, this must be recovered from the Newtonian limit of the theory 4 . Noting the existence of strongly interacting, massive, spin-2 meson states (such as the f-meson), and arguing by analogy with the quantisation of linearised Einstein gravity, it was proposed in sg1 ; sg1a ; sg1b ; sg1c that a short-range ‘strong gravity’ interaction may be responsible for the properties of such elementary particles at a microscopic level.

Thus, a new ‘metric’ tensor $f_{\mu\nu}$, which determines the properties of the strong gravity field, as well as a new strong coupling constant, determined to be of order $1\ \rm GeV^{-1}$ to ensure consistency with the known physics of strong interactions, were introduced sg1 ; sg1a ; sg1b ; sg1c . Throughout the rest of this paper, we follow Salam’s original (though slightly unconventional) notation, denoting the canonical gravitational and strong gravity coupling constants as

[TABLE]

respectively. The dimensionless strong gravity coupling is taken to be of the same order of magnitude as the strong interaction coupling, giving

[TABLE]

This is equivalent to $G_{f}\simeq 10^{38}G=6.67\times 10^{30}\rm cm^{3}g^{-1}s^{-2}$ for a strong interaction scale of $m\simeq 10^{-19}m_{\rm Pl}\simeq 10^{-24}$ g, where $m_{\rm Pl}=\sqrt{\hbar c/G}\simeq 10^{-5}$ g is the Planck mass.

It is interesting to note that the ratio $G_{f}/G$ yields ‘Dirac’s number’, i.e.

[TABLE]

the same dimensionless quantity that formed the basis of his Large Number Hypothesis (LNH) Dirac1937 ; Dirac1938 ; Dirac1974 ; Dirac1979 . (See also MenaMarugan:2001qn ; Ray:2007cc ; Nassif:2015sdu for contemporary viewpoints of the LNH.) The possible implications of this, and its relation to more recent results that imply a relation between the cosmological constant and physics at the electroweak scale, à la Dirac, are discussed below, and at length in Sec. IV

Hence, in order to describe strong interactions, the following Lagrangian density was proposed sg1 ; sg1a ; sg1b ; sg1c

[TABLE]

where $\mathcal{L}_{f}=R_{f}(f)$ and $R_{f}(f)$ is the scalar curvature constructed from $f_{\mu\nu}$ and its derivatives. The term $\mathcal{L}_{fg}$ describes the interaction between f-mesons and gravitons and $\mathcal{L}_{m}$ is the Lagrangian for both strongly interacting and non-strongly interacting matter. If one drops the interaction term $\mathcal{L}_{fg}$ and considers that the dominant term in the field equation obtained from (4) is given by $\mathcal{L}_{f}$, one obtains an Einstein-type equation for the $f_{\mu\nu}$ field, with the strong coupling constant $k_{f}^{2}$ in place of $k_{g}^{2}$. This explains the name strong gravity given to the theory. The physical and mathematical properties of the strong gravity model were investigated in sg2 ; sg3 ; sg4 ; sg5 ; sg6 . For a comprehensive review of the early results in strong gravity, and the corresponding references, see sg7 .

An alternative attempt to describe the physics of strong interactions using a geometric, general relativity-inspired picture is the reformulation of Yang-Mills theory proposed in the so-called ‘chromogravity’ model cg1 ; cg2 . In this model, QCD in the IR region is approximated by the exchange of a dressed two-gluon phenomenological field $\mathcal{G}_{\mu\nu}(x)=B_{\mu}^{a}B_{\nu}^{b}\eta_{ab}$, where $\eta_{ab}$ is a color-SU(3) metric, and $B_{\mu}^{a}$ is the dressed gluon field. This model produces colour confinement, explains the successful features of the hadronic string, predicts the spectrum of baryons and mesons with their Regge trajectories, justifies the interacting boson model, and also ‘predicts’ scaling. The effective Ricci tensor $\mathcal{R}_{\mu\nu}$ constructed from the field $\mathcal{G}_{\mu\nu}$ induces Einsteinian dynamics. Alternative approaches were proposed in Lun1 ; Lun2 ; John and Pav , respectively. The relations between Yang-Mills fields and Riemannian geometry were also investigated in Haag1 ; Haag2 , where it was shown that it is possible to define gauge invariant variables in the Hilbert space of Yang-Mills theories that manifestly implement Gauss’ law on physical states.

From a ‘true’ gravitational perspective, the conditions under which upper and/or lower mass bounds exist for different physical systems is of fundamental importance in theoretical general relativity and relativistic astrophysics. A classic result by Buchdahl Buch states that, for stable, compact, charge-neutral objects of mass $M$ and radius $R$, the condition

[TABLE]

must be satisfied. If the bound (5) is violated, collapse to a black hole becomes inevitable. In Rho1 the general relativistic maximum mass of a stable compact astrophysical stellar type object was found to be of the order of $3.2\rm M_{\odot}$. This result was obtained with the use of the principle of causality, requiring that in the dense matter nuclear matter the speed of sound cannot exceed the speed of light, and of Le Chatelier’s principle. The Buchdahl limit (5) has been extended to include the effects of a nonzero cosmological constant ($\Lambda\neq 0$) g1 , of the electric charge of the sphere g2 , and of an anisotropic interior pressure distribution g3 . Sharp bounds on the mass/radius ratio for neutral and charged compact objects, with both isotropic and anisotropic pressure distributions, in the presence of $\Lambda\neq 0$, were also obtained in g4 ; g5 ; g6 ; g7 ; g8 .

If the existence of an upper mass bound for stellar type structures seems to be a reasonable physical requirement of general relativity, the possible existence of a minimum mass is less obvious. In m1 it was shown that the presence of a positive cosmological constant implies the existence of a minimum classical mass and of a minimum density in nature. These results rigorously follow from the generalized Buchdahl inequality in the presence of $\Lambda>0$, given by

[TABLE]

which implies the existence of a lower bound for the mass/radius ratio or, equivalently, the density of a stable, charge-neutral, gravitating compact object, i.e.

[TABLE]

Though the derivation of this condition is somewhat involved m1 , its physical meaning is intuitively obvious. The dark energy density is given by

[TABLE]

Hence, Eq. (7) simply states that spherical objects with densities significantly lower than the dark energy density have insufficient self-gravity to overcome the effects of dark energy repulsion. For $\rho\lesssim\rho_{\rm min}=\rho_{\Lambda}/2$, the classical radius $R$ becomes unstable and dark energy repulsion blows the object apart. For future convenience, we note that the current experimental value of the vacuum energy density, inferred from observations of high-redshift type 1A supernovae (SN1A), Large Scale Structure (LSS) data from the Sloan Digital Sky Survey (SDSS) and Cosmic Microwave Background (CMB) data from the Planck satellite, is is $\rho_{\Lambda}=5.971\times 10^{-30}\;{\rm g\;cm^{-3}}$ Betoule:2014frx ; Planckresults . This corresponds to a value of $\Lambda=1.114\times 10^{-56}\;\;\rm cm^{-2}$ for the cosmological constant.

The minimum mass/radius bound, in the absence of dark energy (i.e., for $\Lambda=0$), was extended for the case of charged objects in m2 , where it was shown that a stable object with charge $Q$ must obey the relation

[TABLE]

This result was extended to to include the presence of a cosmological constant ($\Lambda\neq 0$), and the mass/radius ratio of a stable charged object was found to obey the relation

[TABLE]

to leading order in $Q^{2}/R$. For $\Lambda=0$, this relation recovers an earlier result due to Bekenstein Bekenstein:1971ej , which demonstrated that the expression for the classical radius of a charged particle – obtained from equating its rest mass with its electrostatic potential energy – remains rigorously valid in general relativity.

In B2 and LakePaterek1 ; Lake1 , Eq. (10) was combined with minimum length uncertainty relations (MLURs), obtained from gravitational extensions of canonical quantum theory, leading to a minimum mass bound for stable, charged, quantum mechanical and gravitating compact objects of the form

[TABLE]

where

[TABLE]

and $q_{\rm Pl}=\hbar c$ is the Planck charge. Evaluated for $Q=\pm e$, this gives

[TABLE]

where $\alpha_{e}=e^{2}/q_{\rm Pl}^{2}\simeq 1/137$ is the usual fine structure constant and $m_{e}$ is the electron mass, which is equivalent to

[TABLE]

According to this relation, if the electron were any less massive (with the same charge $e$) or more highly charged (with the same mass $m_{e}$) a combination of electrostatic and dark energy repulsion would destabilise its Compton wavelength B2 ; LakePaterek1 ; Lake1 . We also note the close similarity between the physical picture of the electron, modelled as an extended charged fluid sphere Bekenstein:1971ej , used to derive Eq. (14) in B2 , and Dirac’s ‘extensible’ model of the electron, proposed in Dirac:1962iy . Equation (14) may also be rewritten as

[TABLE]

yielding an upper bound for $\Lambda$ which is consistent with the current best fit value inferred from cosmological observations Planckresults ; Betoule:2014frx . Interestingly, the relation (15) was previously derived via three different methods (see Nottale1993 ; m3 ; Beck:2008rd ) and, for the mass scale $m=m_{e}/\alpha_{e}$, is equivalent to Zel’dovich’s estimate of $\Lambda$, based on his reformulation of Dirac’s LNH for an asymptotically de Sitter Universe Zel'dovich:1968zz . These results suggest the existence of a deep connection between gravity, the presence of a positive cosmological constant, and the stability of fundamental particles, and are discussed further in the context of the strong gravity model, with a ‘strong cosmological constant’ $\Lambda_{f}$, in Sec. IV.

The particle physics and cosmological implications of the mass scale $M_{T}=(\hbar^{2}\sqrt{\Lambda}/G)^{1/3}$, which corresponds to taking $Q^{2}\rightarrow q^{2}_{\rm Pl}$ in Eq. (11), were considered in B4 , where, based on an MLUR, it was shown that a black hole with age comparable to the age of the Universe would stop radiating and form a relic state when its mass reaches the dual mass scale $M_{T}^{{}^{\prime}}=m_{\rm Pl}^{2}/M_{T}=c(\hbar/G^{2}\sqrt{\Lambda})^{1/3}$. Moreover, it was shown that a holographic relation exists between the entropy and horizon area of the remnant black hole in generic dimensions.

Though, in the present work, we derive mass bounds in the context of the original strong gravity theory, based on the analogy between general relativity and strong interactions, we note that, in recent years, many theories of modified gravity have been proposed in the literature mg1 ; mg2 ; mg3 . In general, these aim to solve the problems posed by modern observational cosmology without the need to posit the existence of exotic states of matter and energy, i.e. dark matter and dark energy d1 ; d2 ; d3 . Theoretically, such approaches may also be extended to the physics of strong interactions: if modified gravity theories possess desirable properties from a cosmological perspective, could modified strong gravity theories possess desirable properties from a particle physics point of view?

Though beyond the scope of this paper, we note that upper and lower bounds on the mass/radius ratio of stable compact objects in modified gravity theories were obtained in B3 , in which modifications of the canonical gravitational dynamics were described by an effective contribution to the matter energy-momentum tensor. As an application of the general formalism developed therein, compact bosonic objects, described by scalar-tensor gravitational theories with self-interacting scalar potentials, and charged compact objects were considered. For Higgs type potentials, it was found that the mass bounds can be expressed in terms of the value of the potential at the surface of the compact object. The general implications of minimum mass bounds for the gravitational stability of fundamental particles and for the existence of holographic duality between bulk and boundary degrees of freedom were also investigated.

It is the goal of this work to investigate the existence of mass bounds in the strong gravity model proposed in sg1 ; sg1a ; sg1b ; sg1c , and to discuss the relevance of these bounds for hadronic physics and cosmology via the holographic principle. To prove the existence of both minimum and maximum mass bounds, we consider a static, spherically symmetric ‘geometry’ for the strong gravity metric, together with the Einstein gravitational field equations, in which the matter energy-momentum tensor consists of two components: ordinary matter, described thermodynamically by its energy density and anisotropic pressure distribution, and a mixing term. With a specific choice of metric tensor, the coefficients of the contribution from the mixing term take the form of an effective strong cosmological constant, $\Lambda_{f}$, whose repulsive (or attractive) force is ‘felt’ only inside the strongly interacting matter.

After determining the effective Einstein field equations, the Tolman-Oppenheimer-Volkoff (TOV) equation describing the equilibrium properties of the strong gravity system is obtained. With the use of this equation, and adopting some physically reasonable assumptions about the behaviour of the physical and geometrical quantities, we derive the generalized Buchdahl inequality, which is valid at all points inside the compact objects. By evaluating this bound on the surface of the hadronic ‘particle’, we therefore obtain both upper and lower bounds of the mass/radius ratio of the hadrons in the strong gravity model. These bounds depend on the mass parameter (i.e. coupling) in the mixing Lagrangian $\mathcal{L}_{fg}$, as well as of the geometric properties of the hadrons. The physical implications of our results are also discussed.

This paper is organised as follows. The TOV equation for strong gravity is derived in Sec. II. The generalised Buchdahl inequality, and the resulting upper and lower bounds on the mass/radius ratio of strongly interacting particles, are derived in Sec. III. The strong gravity mass gap, and its implications for holography, are discussed in Sec. IV. In Sec. V, we combine the mass bounds obtained in the strong gravity model with MLURs motivated by quantum gravity research, exchanging $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$ where necessary. Identifying the ‘strong dark energy density’ with the bag constant of the MIT bag model, which is of the order of the nuclear density, $B\simeq 2\times 10^{14}\ {\rm gcm^{-3}}$, then gives rise to new mass bounds for both neutral and charged strongly interacting particles. Section VI contains a brief summary and discussion of our main conclusions, of outstanding problems, and of prospects for future work.

II Tolman-Oppenheimer-Volkoff Equation in the Strong Gravity Model

In the present Section we briefly review the physical basis and mathematical formalism of the strong gravity model, in which it is assumed that a tensor field, obeying an Einstein-type equation, plays a fundamental role in strong interaction physics. In this approach, strong interactions are governed by a set of field equations formally identical to the Einstein equations, apart from the coupling parameter $k_{f}\simeq 1$ GeV*-1*, which replaces the Newtonian coupling $k_{g}\simeq 10^{-19}$ GeV*-1*. Under the assumption of static spherical symmetry we write down the field equations of the model in the presence of an energy-momentum tensor containing an anisotropic fluid term, and derive the generalised Tolman-Oppenheimer-Volkoff (TOV) equation describing the equilibrium properties of stable, compact, hadronic objects. Throughout, we use the sign conventions and the definitions of the geometric tensors given in 4 .

II.1 Strong gravity

As stated in the Introduction, the Lagrangian for the interacting strong field metric $f_{\mu\nu}$ and gravitational metric $g_{\mu\nu}$ can be constructed as

[TABLE]

where the corresponding ‘volume element’ $\sqrt{-f}$ is defined via $f={\rm det}f_{\mu\nu}$. Here, the first term represents the standard general relativistic Lagrangian for the gravitational field, while the second is its strong interaction analogue, obtained by replacing $k_{g}$ by $k_{f}$ and $g_{\mu\nu}$ by $f_{\mu\nu}$. To give the ‘elementary’ particles mass, as well as their weak gravitational interaction, a mixing term between the strong and weak gravitational fields, $\mathcal{L}_{fg}$, is needed. $\mathcal{L}_{m}$ represents the matter Lagrangian for both strongly interacting matter and non-strongly interacting matter, where it is assumed that the latter contains terms in $g_{\mu\nu}$ and its derivative only, whereas the former may depend (generically) on both $g_{\mu\nu}$ and $f_{\mu\nu}$. Hence, although the strong gravity metric minimally couples to all forms of matter (see below), strongly interacting particles and non-strongly interacting particles ‘feel’ the curvature of the strong metric differently. A simple covariant mixing term was proposed in sg3 and is given by

[TABLE]

where $\mathcal{M}$ is a constant with the dimension of mass. For later use, the full dimension of the mass mixing parameter is given by $\mathcal{M}^{2}\to\mathcal{M}^{2}c^{2}/\hbar^{2}$, the inverse Compton wavelength squared.

In the limit in which the gravitational field may be ignored, $g_{\mu\nu}\rightarrow\eta_{\mu\nu}$, the field equations of the strong gravity theory can be written as

[TABLE]

where

[TABLE]

In the following, we will consider the effect of the strong gravity interaction for a static sphere filled with strongly interacting matter fluid. In spherical polar coordinates $\left\{t,r,\theta,\phi\right\}$ the line element with respect to the strong gravity metric $f_{\mu\nu}$ is assumed to be of the form

[TABLE]

where $d\Omega_{2}=d\theta^{2}+\sin^{2}\theta d\phi^{2}$ is the line-element for the unit $2$-sphere and $\nu(r)$, $\lambda(r)$ and $\Sigma(r)$ are arbitrary functions of the radial coordinate. Furthermore, we assume that the fluid can be described by the standard energy-momentum tensor

[TABLE]

where $\rho c^{2}$ is the fluid energy density, $P_{r}$ is the radial pressure and $P_{\perp}$ denotes the tangential pressure.

Ignoring the weak gravitational interaction, the field equations for the strong gravity field coupled to the matter fluid are, therefore,

[TABLE]

where the raising and lowering of spacetime-like indices is performed via the tensors $f^{\nu\lambda}$ and $f_{\mu\nu}$, respectively, satisfying the condition

[TABLE]

Note the universality of coupling between strong gravity and any form of matter energy-momentum tensor.

II.2 Strong gravity field equations

The stress tensor of the mixing term between massive and massless gravitons can be computed straightforwardly, giving

[TABLE]

For simplicity, and also to make $T^{t(s)}_{t}$ interpretable as proper density, we fix the gauge so that sg3

[TABLE]

As a result, $T^{r(s)}_{r}=T^{t(s)}_{t}c^{2}$, giving rise to an equation of state, $P^{(s)}_{r}=-\rho^{(s)}c^{2}$, which is characteristic for the cosmological constant. The particular gauge choice, $\Delta={\rm const.}$, forces the equation of state of matter to satisfy

[TABLE]

since the sum of Eqs. (31) and (32) is identically zero. Such an equation of state is uncommon for ordinary matter, but it is satisfied by a $U(1)$ gauge field in the Coulomb gauge, for example, in the description of a static charged sphere B2 . With this choice of gauge, there is an anisotropic Poincaré stress associated with the $\theta$-component of the energy-momentum tensor, given by

[TABLE]

Moreover, from the second of Eqs. (27) it follows that, since $\left(\nu^{\prime}+\lambda^{\prime}\right)\Delta=0$ (where a prime denotes differentiation with respect to $r$), the functions $\nu$ and $\lambda$ satisfy the condition

[TABLE]

at all points inside the strongly interacting fluid sphere. The field equations of the strong gravity model may then be written in the form

[TABLE]

where the last equation is simply the conservation law, $0=f^{\kappa\mu}\nabla_{\kappa}T_{\mu\nu}$, with respect to the strong gravity metric $f_{\mu\nu}$.

II.3 The Tolman-Oppenheimer-Volkoff equation in strong gravity

Eq. (31) can be directly integrated to give

[TABLE]

where the accumulated mass inside radius $r$ is defined by

[TABLE]

and $\Omega_{2}=\int d\Omega_{2}=4\pi$. After substituting $M_{0}$ into Eq. (32), we have

[TABLE]

Substituting $\nu^{\prime}$ into the conservation law, we obtain the TOV equation for a fluid sphere in the strong gravity model, in the presence of anisotropic stresses, as

[TABLE]

Note that, by taking into account that $\nu+\lambda=0$, Eq. (37) can also be written as

[TABLE]

which, combined with Eq. (35), gives an alternative definition of the mass of the fluid,

[TABLE]

III The Buchdahl inequality and the minimum and maximum mass/radius ratios of stable compact objects in strong gravity

In the present section we derive the generalised Buchdahl inequality that constrains the values of the mass and pressure, as well as the geometric quantities of the strong gravity field (i.e. the strong cosmological constant) at an arbitrary point $r$ inside a compact object. From this inequality, both lower and upper bounds on the mass/radius ratio of a static, spherically symmetric object, interacting according to the strong gravity law, can be easily obtained.

III.1 The Buchdahl inequality

In order to obtain the generalised Buchdahl inequality for strong gravity, we define the following Buchdahl variables

[TABLE]

where

[TABLE]

From Eqs. (35), (37) and (38), we then have

[TABLE]

where we have defined

[TABLE]

which denotes the effective radial pressure. From here on, a prime indicates differentiation with respect to $x$, though this convention does not apply to the notation for the radial coordinate $r^{\prime}$. Further manipulation then leads to

[TABLE]

To separate positive and negative terms on the right-hand side of Eq. (45), we introduce two new quantities, $\gamma$ and $\gamma_{-}$, defined as

[TABLE]

From Eq. (35), it follows that, since $e^{\lambda}\geq 2/3$ and $e^{\nu}\geq 0$, $\gamma_{-}$ is always negative. On the other hand, $\gamma$ is positive definite for $P_{\perp}>P_{r}-27\mathcal{M}^{2}/16\Delta^{3/2}k_{f}^{2}$. For a static charged sphere, $P_{\perp}=-P_{r}=Q^{2}(r)/2r^{4}$ B2 , where $Q(r)$ is the accumulated charge, this condition is valid as long as the pressure from the mass mixing contribution satisfies

[TABLE]

Note that the mixing term could be negative.

Equation (45) can be further simplified by defining

[TABLE]

to obtain

[TABLE]

Assuming monotonically decreasing profiles for both the density $\rho(r)$ and $\gamma(r)$, it follows that for all $r>r^{\prime}$ the conditions

[TABLE]

must hold at all points inside the compact object. We then immediately obtain $w^{\prime}(x)<0$ and, thus, it follows that the right-hand side of Eq. (50) is always negative. Using the mean value theorem, we therefore obtain the following inequalities for the first and second derivatives of $\psi$ with respect to $z$,

[TABLE]

For $\psi(0)=\zeta(0)-\eta(0)>0$, this leads to

[TABLE]

Using condition (51), we find

[TABLE]

where

[TABLE]

and

[TABLE]

is the effective cosmological constant of the strong gravity model.

Similarly,

[TABLE]

Finally, we substitute Eqs. (LABEL:app1), (LABEL:app2) and (37) into the inequality (53) and divide by $\zeta$ to obtain the Buchdahl inequality in strong gravity as

[TABLE]

where

[TABLE]

Equation (58) is valid for all $r$ inside the strong gravity particle. Moreover, its validity does not depend on the sign of $\mathcal{D}$.

III.2 The minimum and maximum mass/radius ratios of hadrons

III.2.1 The upper bound on the mass/radius ratio of hadrons

As a simple application of the Buchdahl inequality in strong gravity (58), we consider the quasi-isotropic limit $\mathcal{D}=0$, corresponding to the condition

[TABLE]

Moreover, we assume that the effective pressure also vanishes at the surface of the massive particle, so that

[TABLE]

By evaluating Eq. (58) at the surface of the hadron $r=R$, using the conditions (III.2.1), we obtain

[TABLE]

where we have denoted $M_{{\rm eff}}=m_{{\rm eff}}(R)$, leading to the well-known result $2M_{{\rm eff}}/R\leq 8/9$ Buch . This shows that, in strong gravity, the maximum possible mass/radius ratio for hadrons should be constrained (at least approximately) by a Buchdahl-type relation. Written in a dimensional form, the Buchdahl inequality for strongly interacting particles of mass $M_{\rm eff}$ and radius $R$ can be written as

[TABLE]

For $G_{f}=6.67\times 10^{30}\ \rm cm^{3}g^{-1}s^{-2}$, this relation is obviously satisfied in the case of proton, with mass $m_{p}=1.672\times 10^{-24}$ g and classical radius $r_{p}=0.875\times 10^{-13}$ cm, such that $2G_{f}m_{p}/c^{2}r_{p}=0.288$. Interestingly, a particle radius around 3.2 times smaller than $r_{p}$ would make the proton unstable from the point of view of strong interactions.

Next we consider the case $P_{\rm eff}\neq 0$. In the quasi-isotropic limit $\mathcal{D}=0$, Eq. (58) gives the upper mass-radius bound

[TABLE]

where we have defined the mean density of the compact object as $\left<\rho_{\rm eff}\right>=M_{\rm eff}/R^{3}$. We assume that the matter radial pressure $P_{r}$ vanishes at the surface of the strong gravity particle, and thus we obtain for the surface effective pressure the expression

[TABLE]

Taking into account that

[TABLE]

[TABLE]

[TABLE]

where we have defined the mean fluid density as $\left<\rho_{0}\right>=M_{0}/R^{3}$, with $M_{0}\equiv M_{0}(R)$, we obtain the following upper limit for the ordinary matter mass/radius ratio of a stable compact object in strong gravity,

[TABLE]

Next, we consider the anisotropic case with $\mathcal{D}\neq 0$. We define the function $f(M_{{\rm eff}},R,\Lambda_{f},\mathcal{D})$ as

[TABLE]

Assuming again that the effective pressure vanishes at the surface of the compact object, $P_{\rm eff}\equiv 0$, Eq. (58) leads to the following general restriction on the mass/radius ratio for a spherical hadronic fluid,

[TABLE]

By taking into account the definition of the total effective mass as given by Eq. (66), we immediately find

[TABLE]

With the use of the Taylor series expansion of the function $\arcsin x/x-1$,

[TABLE]

for small values of the argument, we can approximate the function $f\left(M_{{\rm eff}},R,\Lambda_{f},\mathcal{D}\right)$ as

[TABLE]

Then the maximum mass bound for compact objects in strong gravity with vanishing surface effective pressure can be reformulated as

[TABLE]

Finally, we consider the case of strong gravity compact objects with vanishing surface radial pressure, i.e. with $P_{r}(R)=0$, and $P_{\rm eff}\neq 0$, given by Eq. (65). In this case, we obtain for the maximum mass/radius bounds the expressions

[TABLE]

and

[TABLE]

respectively.

The maximum mass/radius bound for compact objects is generally obtained in the constant density regime, with $\rho_{\rm eff}\simeq\rho_{0}={\rm constant}$. Hence the total mass of the compact object can be approximated as $M_{0}=4\pi\rho_{0}R^{3}/3$. Therefore, in the mass/radius ratio bounds obtained above, we can approximate the effective mean density as $\left<\rho_{\rm eff}\right>\simeq M_{0}/R^{3}\simeq\rho_{0}={\rm constant}$, a relation that is satisfied by the maximum mass objects with a very good approximation. Therefore, in all the above results, the ratio $\Lambda_{f}/\left<\rho_{\rm eff}\right>$ can then be approximated as a constant. Hence, it follows that, generally, the right-hand sides of the upper bounds on the mass/radius ratio can be regarded as independent of the masses of the compact objects.

III.2.2 The lower bound on the mass/radius ratio of hadrons

On the vacuum boundary of the anisotropic fluid distribution, $r=R$, Eq. (58) takes the general form

[TABLE]

where we have assumed that $k_{f}^{2}P_{\mathrm{eff}}/3+M_{\mathrm{eff}}/R^{3}>0$. In addition, we assume that the surface radial pressure $P_{r}$ is negligibly small on the hadron’s surface, so that $P_{\mathrm{eff}}(R)=-3\Lambda_{f}/k^{2}_{f}$, and therefore $P_{\mathrm{eff}}(R)<0$. Again using the fact that, for small values of the argument, the function $\arcsin x/x-1$ can be approximated using Eq. (73), and performing the replacement $P_{\text{eff}}\to-P_{\text{eff}}$, Eq. (78) can be written as

[TABLE]

By introducing a new variable $v$, defined as

[TABLE]

Eq. (79) takes the form

[TABLE]

where we have denoted

[TABLE]

and

[TABLE]

respectively. Then, by squaring Eq. (81), we can reformulate the corresponding inequality as

[TABLE]

or, equivalently,

[TABLE]

where we have denoted

[TABLE]

and

[TABLE]

In the following analysis, we keep only the first order terms in both $p$ (depending on $\Lambda_{f}$) and $q$ (depending on $\mathcal{D}$) in the expressions involving square roots. Since $v\geq 0$, Eq. (85) is satisfied if $v\leq v_{1}$ and $v\geq v_{2}$, or $v\geq v_{1}$ and $v\leq v_{2}$. However, one can easily check that the condition $v\geq v_{1}$ contradicts the upper bound on the mass/radius given by Eq. (72). Therefore, it follows that Eq. (85) is identically satisfied if and only if, for all values of the physical parameters determining the total mass of the hadronic particle, the condition $v\geq v_{2}$ holds. This result is equivalent to the existence of a minimum bound for the mass/radius ratio of particles in strong gravity, which is given by

[TABLE]

By explicitly substituting the expressions for $p$, $q$ and $v$, as defined above, we then obtain the following (alternative) form of lower bound for the mass/radius ratio of hadronic particles in strong gravity,

[TABLE]

By taking into account the explicit expression for the effective mass of the hadron, Eq. (89) may be rewritten as

[TABLE]

Hence, we see that the presence of the anisotropic matter distribution weakens the lower bound on the hadron mass. Nonetheless, in the strong gravity theory, there still exists an absolute minimum mass for hadrons in nature. If the surface anisotropy, described by the coefficient $\mathcal{D}$, can be neglected, the existence of a minimum mass is determined by the presence of the effective strong gravity cosmological constant only. This is constructed from the mass parameter of the model, $\mathcal{M}$, and the arbitrary constant $\Delta$, which fixes the value of the metric tensor coefficients inside the hadrons. For $\Lambda_{f}\equiv 0$, the minimum mass bound simply reduces to the positivity requirement for the bare mass, $M_{0}\geq 0$.

By taking into account the explicit expression for the strong gravity cosmological constant and by assuming that $(k^{2}_{f}/9)\mathcal{D}(R)R^{2}\ll 1$, i.e. assuming that the pressure anisotropy vanishes at the vacuum boundary of the hadron, the minimum mass/radius ratio of hadronic particles, given by Eq. (90), can be written as

[TABLE]

This equation imposes strong (pun intended) constraints on the mass parameter $\mathcal{M}$ which controls the strength of the mixing between the $f_{\mu\nu}$ and $g_{\mu\nu}$ fields, i.e.

[TABLE]

where $\rho_{{\rm min}}=M_{0}^{({\rm min})}/R^{3}$ is the density corresponding to the minimum mass hadron.

III.3 The energy localization problem in strong gravity

An important issue in general relativity is the problem of energy localization. Tentatively, we assume that the total effective energy in strong gravity can be described in a similar way as in canonical Einstein gravity. This assumption allows us to derive explicit limits on the total energy of compact hadronic objects. Hence, we define the total energy inside an equipotential surface $S$, which includes the contribution from the strong tensor field $f_{\mu\nu}$, by analogy with general relativity, as 32 ; 33

[TABLE]

where the vector $\xi^{i}$ is a Killing field of time translation, $\xi_{s}$ denotes its value at $S$, and $\left[K\right]$ is the jump across the shell of the trace of the extrinsic curvature of $S$, assumed to be embedded in the 2-space $t=\mathrm{constant}$. $E_{M}=\int_{S}T_{i}^{k}\xi^{i}\sqrt{-g}dS_{k}$ and $E_{F}$ are the energies of the ordinary matter and of the strong gravitational field in the $f_{\mu\nu}$ metric, respectively. This definition of the total energy is manifestly coordinate invariant. In the case of static spherical symmetry, for both the $g_{\mu\nu}$ and $f_{\mu\nu}$ fields, we obtain the total energy for the hadron, from Eq. (93), as 33

[TABLE]

where, as usual, $[\;\;]_{S}$ denotes the jump across the surface $S$. We also make the fundamental assumption that the metric outside the strong gravity system is of de Sitter type, under the replacement $\Lambda_{f}\rightarrow\Lambda$.

Next, for the sake of convenience, we rescale the metric tensor component so that $(2/3)e^{-\lambda}\rightarrow e^{-\lambda}$. Eq. (35), which may be expressed in terms of the effective mass of the strong gravity object $M_{\rm eff}$, can then be written as

[TABLE]

Such a rescaling is always possible, and is performed explicitly via coordinate transformation $r^{2}\rightarrow r^{\prime 2}=2r^{2}/3$. Note that the mass per unit radial distance transforms as $M(r^{\prime})/r^{\prime}=2M(r)/3r$. For convenience, we will also redefine $\Lambda_{f}\to\Lambda_{f}/3$.

Then, by taking into account the relation $\nu+\lambda=0$, as well as the definition of the strong gravity cosmological constant, the total energy of a compact self-gravitating object may be written as

[TABLE]

or, equivalently,

[TABLE]

With the use of Eq. (58), we find the following upper limit for the total energy of the compact hadronic object in strong gravity,

[TABLE]

In the case of a vanishing strong cosmological constant, $\Lambda_{f}\rightarrow 0$, and also assuming that the matter pressure is zero at the vacuum boundary of the object, from Eq. (98) we obtain the strong gravity equivalent of the standard upper energy bound 32 ; 33

[TABLE]

For a quasi-isotropic matter distribution with $\mathcal{D}=0$, with the assumption of vanishing radial pressure $P_{r}(R)=0$, we obtain

[TABLE]

Eq. (97) also allows to obtain a mass/radius relation for hadronic objects in strong gravity by requiring that the particle is in its minimum energy state, corresponding to $\partial E/\partial R=0$. This gives the following mass/radius relation, as a function of the strong cosmological constant, which is valid for $2\Lambda_{f}R^{2}>1$,

[TABLE]

The two values of the mass differ by a quantity

[TABLE]

IV The mass gap in strong gravity

In Section III.2.1, we have considered the upper bound on the mass/radius ratio in the situation where $P_{\rm eff}(R)=0$ and $\mathcal{D}=0,P_{\rm eff}\neq 0$. In Section III.2.2, the lower bound is derived when we set $P_{r}(R)=0$ in order to highlight the effect of anisotropic parameter $\mathcal{D}$ on the bound. In this section, we will consider the most generic case, without making any assumptions about the value of pressure at the surface of the object, and will simply define $P_{\rm eff}(R)\equiv P_{\rm eff}$.

In order to simplify our formalism, we define the additional dimensionless quantities

[TABLE]

which allows us to express the inequality (58) at $r=R$ as

[TABLE]

This may be written in the following, explicitly quadratic form,

[TABLE]

where

[TABLE]

The inequality can be satisfied if

[TABLE]

for

[TABLE]

for $\mathcal{B}<0$. For $\mathcal{B}>0$, the mass bounds cease to exist. The condition $\mathcal{B}<0$ leads to the following constraint on $F$,

[TABLE]

Another condition for the existence of mass bounds is $\mathcal{B}^{2}>4\mathcal{C}$, which is trivially satisfied for $a>0$. However, for $a<0$, $\mathcal{B}^{2}>4\mathcal{C}$ requires

[TABLE]

To simplify these results, we consider the case when $4\mathcal{C}/\mathcal{B}^{2}\ll 1$ and also $a,F\ll 1$, to obtain

[TABLE]

Under these conditions, the contribution from the anisotropic stress, $F=k^{2}_{f}\mathcal{D}R^{2}/4$, only appears at the next-to-leading order for the lower bound. Finally, the mass bounds for small $a$, $b$ and $F$ are

[TABLE]

The nontrivial minimum mass/radius exists when $P_{\rm eff}<-2\Lambda_{f}/k^{2}_{f}$. This is valid even for $\Lambda_{f}<0$.

Mixing conventional massless gravity with massive ‘gravity’ from the covariant interaction term results in the minimum and maximum mass bounds for any static spherical configuration within the framework of the theory. If we take the QCD glueball as the massive spin-2 state to be mixed with the massless graviton, the model will predict the mass gap of any object composed of particles that couple to the glueball. This does not solve the mass gap problem explicitly, since we assume the glueball mass a priori. However, the mass gap of the glueball is then transmitted to other particles via the universal interaction of the strong gravity field. The mass gap is given by the minimum mass and it is proportional to $\mathcal{M}^{2}/k^{2}_{f}$.

That the mass is generated by mixing with strong gravity induced by the glueball and not by chiral symmetry breaking is a remarkable aspect of this mechanism. It is a universal way to generate mass for a stable static object in the strong interaction. On the other hand, the upper limit represents the maximum mass of the QCD sphere at a given radius (i.e. maximum density) before it undergoes ‘gravitational collapse’ to form a quark-gluon plasma, which cannot be contained within a static sphere. This strong gravitational collapse is nothing but the deconfinement phase transition in the strong interaction. The critical density predicted by this strong gravity model is proportional to $1/k^{2}_{f}$ at the leading order.

However, we may ask, what determines the strong gravity coupling, $k^{2}_{f}$, of the gauge singlet massive ‘graviton’ interaction? As in conventional general relativity, where the Planck mass defines the energy scale at which quantum effects become comparable to those of classical gravity, the analogous mass scale for the strong gravity theory can be used to determine the coupling $k_{f}$ by setting the Compton wavelength of the particle, $\lambda_{\rm C}$, equal to $R$ in Eqn. (LABEL:maxb) with $P_{\text{eff}},\Lambda_{f},\mathcal{D}=0$. We then have

[TABLE]

in the standard units. We define $M_{\rm max}$ to be the corresponding ‘Planck mass’ of the strong gravity. Whereas, for $M>m_{\rm Pl}$, fundamental particles inevitably collapse to form black holes, strongly interacting particles with $M>M_{\rm max}$ inevitably undergo a deconfinement phase transition. The corresponding length scale, given by a Compton relation, $R_{\rm min}=\hbar/M_{\rm max}c$, may be referred to as the strong gravity ‘Planck length’.

The mass gap generation mechanism discussed in this paper assumes the glueball mass to be proportional to $\mathcal{M}$ (i.e. the mass mixing term is proportional to $\mathcal{M}^{2}$), while the mass gap itself, given by Eq. (113), is proportional to $\mathcal{M}^{2}R^{3}/k^{2}_{f}$, yielding

[TABLE]

If we set the mass gap equal to the mass of the $\pi$-meson, $m_{\text{gap}}=m_{\pi}\simeq 140$ MeV, $R\simeq R_{f}/2$ and $\mathcal{M}=2$ GeV, the strength of strong gravity becomes

[TABLE]

Any quantum particle with a strong gravity interaction will inevitably collapse to form a strong gravity ‘black hole’ once the mass exceeds $M_{\rm{max}}$. The corresponding Hawking temperature of the strong gravity black hole is

[TABLE]

This is the Hagedorn temperature, i.e. the maximal possible temperature of the nuclear matter, before the phase transition occurs. Hence, we can identify this temperature with the deconfinement temperature of the hadron. In this picture, strong gravitational collapse is the deconfinement of the strong nuclear interaction. The strong gravity field $f_{\mu\nu}$ becomes zero/infinity at the ‘horizon’ and reverses sign inside the ‘black hole’. We interpret this as the non-existence of the glueball in the deconfined phase.

From the point of view of holographic duality, this looks very familiar. The correspondence between the maximum mass of spherical object in AdS space and the deconfinement temperature of the dual gauge matter is well known Witten:1998zw ; Danielsson:1999zt ; Arsiwalla:2010bt ; B5 . Furthermore, the thermal phase of an AdS black hole is argued to be dual to the thermal phase of deconfined gauge matter living on the boundary of the AdS space. Hence, gravitational collapse in the bulk AdS space corresponds to the deconfinement phase transition of the gauge matter on the AdS boundary. This relationship is holographic in nature since it relates two theories living in different dimensions of spacetime. It is interesting to note similar features of the strong gravity model.

In B2 ; B3 , it was argued that the minimum mass bound should be interpreted as the minimum density required for the nuclear matter to maintain its static configuration without evaporating into a hadron gas. However, in light of the strong gravity model presented here, the minimum mass bound in the bulk AdS could also very well be interpreted as corresponding to the minimum mass of any nuclear particle which is stable under the strong interaction. (Note, however, that such particles may still be unstable with respect to other interactions, such as weak decay, etc.) In this picture, the minimum mass bound in the bulk simply is the mass gap in the strong interaction on the boundary.

Finally, before concluding this Section, we comment on the requirement that $P_{r}+\rho=0$, which originates from the choice $\Delta={\rm const.}$ in the ‘gauge fixing’ of the interaction term of the strong gravity Lagrangian. Even though a $U(1)$ charged sphere satisfies this equation of state, generic matter does not obey this condition. In this case, the quantity $\Delta$ becomes physically relevant and we need to allow $\Delta$ to depend on the radial coordinate $r$, setting $\Delta={\rm const.}\rightarrow\Delta(r)$. We then have

[TABLE]

resulting in an increasing function $\Delta(r)\geq 1$ for a positive $(P_{r}+\rho)$ matter profile. The mass bound analysis above can then be repeated with

[TABLE]

replacing $\rho$ in Eq. (31) and by setting $b=0$. The resulting minimum and maximum mass/radius ratio bounds are exactly the same as in Eqns. (113) and (LABEL:maxb) under the replacement $\Lambda_{f}\to\Lambda_{f}^{\rm(min)},\Lambda_{f}^{\rm(max)}$ where

[TABLE]

respectively. The existence of mass gap is generic in this kind of model.

For convenience, here and henceforth we rewrite the metric in a rescaled coordinate $r^{\prime}=r/\sqrt{3/2}$ and redefine $\Lambda_{f}\to\Lambda_{f}/3$ so that our metric is in the conventional form.

V Quantum mechanical implications of the classical mass/radius ratio bounds in strong gravity

In the present Section, we investigate the quantum mechanical implications of the mass/radius bounds in the strong gravity model. We begin with a brief discussion of the quantum implications of mass bounds in conventional general relativity, in the presence of a cosmological constant $\Lambda\neq 0$, before extending these to the strong gravity case via the substitutions $G\rightarrow G_{f}$, $\Lambda\rightarrow\Lambda_{f}$.

V.1 Quantum mass bounds in standard general relativity

For fundamental particles, viewed as stable compact objects, the radius of an uncharged particle may be identified with the Compton wavelength $\lambda_{\rm C}$, or reduced Compton wavelength, $k_{\rm C}^{-1}$, given by

[TABLE]

For order of magnitude relations, we use these two expressions interchangeably from here on. The combination of Eqs. (8) and (121) then implies the existence of minimum mass for a stable, charge-neutral, quantum mechanical and gravitating compact object,

[TABLE]

where, for future reference, we define the (reduced) Planck mass and length scales $m_{\rm Pl}=\sqrt{\hbar c/G}=2.176\times 10^{-5}\ {\rm g}$, and $l_{\rm Pl}=\sqrt{\hbar G/c^{3}}=1.616\times 10^{-33}\ {\rm cm}$, respectively, and the (reduced) first/second de Sitter mass and length scales, denoted using unprimed/primed quantities, respectively, as

[TABLE]

Note that, in previous work B1 ; B2 ; B3 , the de Sitter scales defined in Eq. (V.1) were also referred to as the first and second Wesson scales, following the pioneering work Wesson:2003qn . The physical interpretations of these scales are discussed in detail in B1 . For now, we simply note that the first and second de Sitter scales are related via

[TABLE]

and that the numerical value of the minimum mass, $m_{\Lambda}\simeq 10^{-3}\rm eV$ is consistent with current experimental bounds on the mass of the electron neutrino, the lightest known neutral particle, obtained from the Planck satellite data, $m_{\nu}\leq 0.23\ \rm eV$ Planckresults . To within numerical factors of order unity, $m_{\Lambda}$ is also the unique mass scale for which the Compton radius of the particle is equal to its gravitational turn-around radius,

[TABLE]

in the presence of dark energy Bhattacharya:2016vur ; LakePaterek1 ; Lake1 . This represents the radius beyond which the repulsive effects of the background dark energy dominate over the attractive force of canonical gravity. Eq. (125) may be obtained by considering the Newtonian limit of general relativity for the Schwarzschild-de Sitter metric, which gives rise to an effective Newtonian potential of the form

[TABLE]

This, in turn, gives the effective gravitational field strength Hobson:2006se

[TABLE]

which changes sign at $r=r_{\rm grav}$. However, to within numerical factors of order unity, the expression (125) remains rigorously valid in full general relativity Bhattacharya:2016vur . This gives a neat way of reinterpreting the stability bound (7). If the quantum mechanical (i.e. Compton) radius of the particle lay outside its gravitational turn-around radius, it would clearly be unstable due to dark energy repulsion. Interestingly, we may also ask, for what mass is the turn-around radius equal to the Schwarzschild radius? The answer is $M\simeq m_{\rm dS}^{\prime}=(c^{2}/G)\sqrt{3/\Lambda}$, which is comparable to the present day mass of the Universe B1 ; B2 .

V.2 Quantum mass bounds in strong gravity

V.2.1 Mass bounds for neutral particles

The Planck mass and length scales are obtained by equating the Compton wavelength of a quantum mechanical particle with the Schwarzschild radius induced by its classical gravitational field (ignoring numerical factors of order unity). In the strong gravity model, an analogous construction using the Schwarzschild radius of the $f_{\mu\nu}$ field gives

[TABLE]

where $m_{\rm sP}$ and $l_{\rm sP}$ denote the strong gravity Planck mass and strong gravity Planck length, respectively. Note that $m_{\rm sP}$ and $l_{\rm sP}$ are equivalent to $M_{\rm max}$ and $R_{\rm min}$, defined in Sec. IV. We here relabel these quantities for the sake of easy comparison with results from standard general relativity. Likewise, analogues of the de Sitter scales may be obtained by replacing $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$ in Eq. (V.1). In addition, based on purely dimensional arguments, we may define two additional mass scales, and their corresponding lengths, by mixing and matching $\left\{G,\Lambda,G_{f},\Lambda_{f}\right\}$.

To investigate the physical meaning (if any) of these scales, and of other mass/length scales constructed using the strong gravity model parameters, we must first consider the physical interpretation of the mass term $\mathcal{M}^{2}$ in the mixing Lagrangian $\mathcal{L}_{fg}$. This appears in the strong gravity field equations in combination with the ‘geometric’ parameter $\Delta$, through the definition of the strong cosmological constant, $\Lambda_{f}=\mathcal{M}^{2}/4\Delta^{3/2}$ (56). By analogy with Eq. (8), we define the energy density associated with $\Lambda_{f}$ as

[TABLE]

This may be related to the equation of state for deconfined quark matter, obtained from perturbation theory in QCD, as follows. Neglecting quark masses in the first order perturbation, the relation between the pressure and energy density of nuclear matter is given by

[TABLE]

where $B\simeq 2\times 10^{14}\rm gcm^{-3}$ is the difference in energy density between the perturbative and the non-perturbative QCD vacuums, and is of the order of the nuclear density. This model is known as the MIT ‘bag’ model and the constant $B$ is called the bag constant. When nuclear matter is compressed to sufficiently high density, a phase transition is thought to occur which converts confined hadronic matter into free, three-flavor (strange) quark matter. The collapse of the quark fluid is described by the bag model equation of state (130).

In sg3 , it was already pointed out in that the QCD bag constant effectively resembles a cosmological constant for strongly interacting matter. Qualitatively at least, it is not difficult to understand how the effective potential of the strong gravity model mimics the effective (‘bag-type’) potential obtained from QCD. The Newtonian limit of the Schwarzschild-de Sitter type metric for the $f_{\mu\nu}$ field gives

[TABLE]

[TABLE]

by analogy with Eqs. (126)-(127). However, here, the $f_{\mu\nu}$ ‘vacuum’ corresponds to the presence of strongly interacting matter and the ‘cosmological constant’ $\Lambda_{f}$ is generated by the interaction term in the strong gravity Lagrangian, $\mathcal{L}_{fg}$. Hence, the strong force is attractive on small scales, $r\leq r_{\rm grav(s)}$, where $r_{\rm grav(s)}$ is the strong gravity turn-around radius,

[TABLE]

repulsive on intermediate scales, $r_{\rm grav(s)}\leq r\leq R$, and quickly tends to zero in the true vacuum ($r\gtrsim R$), where the density of the strongly interacting matter also falls quickly to zero. Here, $r_{\rm sS}$ denotes the string gravity Schwarzschild radius.

However, if we may identify the outer radius of the strongly interacting ‘particle’ with the strong gravity turn-around radius $R\simeq r_{\rm grav(s)}$, the repulsive phase is never realised. The effective potential is strongly attractive over short distances and is (in principle) capable of countering the effects of electrostatic repulsion if the matter is also charged. It then vanishes at $R\simeq r_{\rm grav(s)}$ and remains zero outside the particle. Furthermore, if the analysis presented above is modified to include $r$-dependence in the geometric parameter $\Delta$ in Eqn. (35), i.e. such that

[TABLE]

we may generate a confining potential, $\Phi_{\rm s}\propto r$, for $r\gtrsim r_{\rm grav(s)}\simeq R$. This possibility is discussed further in Sec. VI. In the present section, we assume an effective potential of the form (132), which holds up to $r\simeq r_{\rm grav}$, given by (133). This allows us to treat the ‘bag’ of the MIT bag model within the context of the strong gravity theory.

Technically, since $B$ has dimensions of density, it plays the role of the ‘vacuum’ energy density associated with $\Lambda_{f}$ (129), which is nonzero inside the strongly interacting particles. We note also that, for $\rho=B\simeq{\rm const.}$, we have $\rho=-p/c\simeq{\rm const.}$, a dark energy-type equation of state. Hence, we may identify

[TABLE]

where $\rho_{\Lambda}=\Lambda c^{2}/(8\pi G_{f})=5.971\times 10^{-30}{\rm gcm^{-3}}$, $\Lambda=1.114\times 10^{-52}\ {\rm cm^{-3}}$ are the current best fit values obtained from cosmological observations Planckresults ; Betoule:2014frx . For $G_{f}\simeq 10^{38}G$ (3), this implies

[TABLE]

For this value of $\Lambda_{f}$ we obtain, by analogy with Eq. (V.1),

[TABLE]

where the subscript ‘sd’ stands for ‘strong-de Sitter’. Remarkably, we note that, for $G_{f}\simeq 10^{38}G$ (3), the strong gravity Planck and de Sitter scales are approximately equal,

[TABLE]

As we shall see, this has important implications for the model as an effective theory, able to mimic confinement.

Additional mixing and matching of $\left\{G,\Lambda,G_{f},\Lambda_{f}\right\}$ also yields

[TABLE]

based on purely dimensional arguments. We christen these the first and second Salam mass/length scales, though their physical interpretations are not investigated here as such an analysis lies beyond the scope of the present work.

Instead, in the analysis that follows, we focus on the strong gravity analogues of mass/length scales which are well defined and understood in standard general relativity, and on additional mass/length scales that may be derived from them using quantum gravity arguments. These include the strong gravity Planck scales (V.2.1), strong gravity de Sitter scales (V.2.1), the strong gravity Schwarzschild and turn-around radii (133) and the analogues of the charge-neutral an charged particle mass bounds, $M\gtrsim M_{\Lambda}$ (122) and $M\gtrsim M_{Q}$ (11), obtained by replacing $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$. In addition, we consider a modified charge-neutral bound, $M\gtrsim\tilde{M}_{\Lambda}$, which corresponds to replacing $G\rightarrow G_{f}$, $\Lambda\rightarrow\Lambda_{f}$ and identifying the compact radius $R$ with the scattering radius of the particle, rather than its Compton wavelength. As we shall see, this allows us to recover the mass of the neutron as the mass of the lightest, stable, compact, charge-neutral and strongly interacting quantum mechanical particle in nature, according the strong gravity theory.

In strong gravity, the analogue of the stability condition for neutral particles in general relativity with $\Lambda>0$, Eq. (7), is

[TABLE]

where $M_{0}$ is the bare mass of the hadron. (However, from here on, we simply use $M_{0}\rightarrow M$ to denote the bare mass.) This follows directly from Eq. (90) for the quasi-anisotropic case, $\mathcal{D}=0$.

Even though the bounds obtained in previous sections are derived from the classical strong gravity field equations, we would like to extend them to the quantum mechanical regime. Thus, we consider the situation where the minimum bound is saturated by a ‘classical’ particle that is equally quantum mechanical, in the sense that its classical size is equal to its Compton quantum wavelength. The particle will become ‘purely’ quantum if its classical radius becomes even smaller. Setting $R=k_{\rm C}^{-1}$, where $k_{\rm C}^{-1}$ denotes the reduced Compton wavelength (121), Eqn. (140) then gives

[TABLE]

This may also be obtained, to within numerical factors of order unity, by equating the Compton scale with the strong gravity turn-around radius, defined by Eq. (133).

Hence, according to the strong gravity theory, combined with elementary quantum mechanics, $M_{\Lambda_{f}}\simeq 10^{-25}$g should correspond to the mass of the lightest possible stable, compact, charge neutral, strongly interacting and quantum mechanical particle found in nature. With this in mind, we note that this is almost equal to the neutron rest mass $m_{n}\simeq 1.675\times 10^{-24}$ g, though a discrepancy of around one order of magnitude remains. That said, also we note that the neutron is not a fundamental particle. Therefore it is unclear whether $R$ in Eq. (140) should be identified with $\lambda_{n}=\lambda_{\rm C}(m_{n})=2.100\times 10^{-14}$ cm – here we estimate the reduced Compton wavelength but use the standard (ambiguous) notation – or with some other measure of the neutron radius.

In particular, we may instead consider the neutron radius obtained from scattering cross section data, $\sigma_{n}=\pi r_{n}^{2}\simeq 10^{-24}\ {\rm cm^{-2}}$, yielding $r_{n}\simeq 5.642\times 10^{-13}\ {\rm cm}\simeq 26.864\times\lambda_{n}$. Identifying the radius $R$ in Eq. (140) with the scattering radius of the particle $r_{\rm scat}$,

[TABLE]

yields an alternative estimate of the lightest neutral hadron in the strong gravity theory, which we denote $\tilde{M}_{\Lambda_{f}}$. This is given by

[TABLE]

For $\chi_{\rm scat}\simeq 26.864$, the relevant value for the neutron, this yields

[TABLE]

Finally, we note that, due to the approximate numerical coincidence of $m_{\rm sP}$ and $m_{\rm sd}$, we have that

[TABLE]

and similar relations hold for $M=\tilde{M}_{\Lambda_{f}}$. This justifies our earlier assumption that the effective potential (132) holds up to $R\simeq r_{\rm grav(s)}$ (133).

For convenience, we denote the reduced Compton scales associated with $M_{\Lambda}$, $M_{\Lambda_{f}}$ and $\tilde{M}_{\Lambda_{f}}$ as

[TABLE]

[TABLE]

and

[TABLE]

respectively, from now on.

V.2.2 Mass bounds for charged particles

Having considered neutral particles, we now try to combine classical stability bounds for charged strongly interacting fluid spheres with quantum mechanics. To this end, we now (briefly) review the derivations of Eqs. (11)-(15) presented in references B2 ; LakePaterek1 ; Lake1 . (Note that Eq. (15) was also derived, using different methods, in Nottale1993 ; Beck:2008rd and that its cosmological implications were investigated in Wei:2016moy , while similar expressions, equivalent to replacing $m_{e}/\alpha_{e}\leftrightarrow m$, were obtained in Zel'dovich:1968zz and Funkhouser:2005hp .)

As we shall see, the approach taken in B2 ; LakePaterek1 ; Lake1 , based on hypothetical minimum length uncertainty relations (MLURs), may be readily extended to the strong gravity theory, leading to expressions analogous to Eqs. (11)-(15), but with $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$. Having obtained these, we again identify the energy density associated with $\Lambda_{f}$, $\rho_{\Lambda_{f}}\equiv\Lambda_{f}c^{2}/(8\pi G_{f})$, with the ‘bag constant’ of the MIT bag model, $B\simeq 2\times 10^{14}\ {\rm gcm^{-3}}$. This, in turn, allows us to obtain a numerical estimate of the minimum mass of a stable, compact, charged, strongly interacting and quantum mechanical particle. For $Q=\pm 2e/3$, this is found to be of the same order of magnitude as the masses of the lightest known particle of this form, i.e. the mass of the up quark, which is believed to lie in the range 1.7-3.3 MeV.

We emphasise that, when estimating the the mass of the lightest stable, charge neutral and strongly interacting quantum mechanical particle, we expected to obtain an estimate of the neutron mass $m_{n}$, whose density is $\rho_{n}\simeq B$. By contrast, when considering the lightest possible charged and strongly interacting quantum mechanical particle, we expect to obtain an estimate of the lightest known quark mass. This is because there are no known fundamental, charge-neutral, and strongly interacting particles in nature, whereas fundamental charged and strongly interacting particles (i.e. quarks) do exist. However, in considering the mass of a free quark, we must consider the point at which it becomes unconfined, and identify this with the ‘strong dark energy’ density after the phase transition to the quark-gluon plasma.

In Ng:1993jb ; Ng:1994zk , an MLUR of the form

[TABLE]

where $\alpha^{\prime}$, $\beta={\rm const.}$, was proposed. Here, $\Delta x_{\rm total}$ represents the minimum possible uncertainty in the position of a ‘probe’ particle, which is used to measure a distance $r=ct$ by means of the emission and reabsorption of a photon. Thus, it is equal to the minimum possible uncertainty in the measurement of the probe distance $r$. The first term on the top line of Eq. (154) is the standard Heisenberg term, rewritten using the relations $\Delta p=M\Delta v$ and $\lambda_{\rm C}\simeq k_{\rm C}^{-1}=\hbar/(Mc)$, whereas the second represents a recoil term, due to the emission of the photon Salecker:1957be . The third is the ‘gravitational uncertainty’, which is assumed to be of the order of the Schwarzschild radius $r_{\rm S}\simeq l_{\rm Pl}^{2}/\lambda_{\rm C}$ Ng:1993jb ; Ng:1994zk . The second line is obtained by minimising the first with respect to $\Delta v$, giving

[TABLE]

Minimising the expression on the second line of Eq. (149) with respect to $M$ then yields

[TABLE]

and hence

[TABLE]

The $M$ in Eq. (151) represents the optimum mass for the probe particle. This yields the minimum possible uncertainty in the measurement of the probe distance $r$, given by Eq. (152). The probe distance which may be measured with minimum uncertainty, denoted $r_{\rm min}$, is defined via Eq. (151).

The canonical quantum part of the MLUR (149), $(\Delta x_{\rm canon.})_{\rm min}\gtrsim\sqrt{\lambda_{\rm C}r}$, was originally derived by Salecker and Wigner using the gedanken experiment considered above (neglecting the particle’s self-gravity) Salecker:1957be but may also be derived more rigorously by directly solving the Schödinger equation in the Heisenberg picture, before setting $t=r/c$ Calmet:2004mp ; Calmet:2005mh . Though derived via different means, an MLUR of the form (152) was originally obtained by Károlyházy, under the assumption of asymptotically flat space, in Karolyhazy:1966zz ; KFL . In most of the existing quantum gravity literature, the constants $\alpha^{\prime}$ and $\beta$ are assumed to be of order unity, $\alpha^{\prime}$, $\beta\sim\mathcal{O}(1)$ Hossenfelder:2012jw ; Garay:1994en .

However, in B2 ; LakePaterek1 ; Lake1 , it is argued that the introduction of a constant dark energy density, i.e. $\Lambda>0$, and, hence, the existence of a de Sitter horizon, $l_{\rm dS}=\sqrt{3/\Lambda}$, implies a fundamental modification of the MLUR (149), equivalent to the substitution

[TABLE]

where $\beta^{\prime}={\rm const.}\sim\mathcal{O}(1)$. Equation (149) then becomes

[TABLE]

and the analogues of Eqs. (151)-(152) are

[TABLE]

and

[TABLE]

respectively. Here, $R_{\rm cell}$ represents the linear dimension associated with a fundamental ‘cell’ within the de Sitter horizon, yielding a holographic relation between the number of degrees of freedom in the bulk and on the boundary B2 .

In LakePaterek1 it is also shown that the dark energy-modified MLUR, dubbed the ‘dark energy uncertainty principle’ or DE-UP for short, is consistent with the minimum-mass bound $M\gtrsim M_{\Lambda}$, obtained independently in B1 . Requiring every (potentially) observable length scale in the DE-UP, i.e. $r$, $(\Delta x_{\rm canon.})_{\rm min}\gtrsim\sqrt{\lambda_{\rm C}r}$ and $\Delta x_{\rm grav}\simeq\beta^{\prime}l_{\rm Pl}^{2}l_{\rm dS}/(\lambda_{\rm C}r)$, to be super-Planckian leads naturally to Eq. (122). In addition, since Eq. (154) is invariant under simultaneous rescalings of the form

[TABLE]

where $\alpha_{\rm Q}>0$ is a positive real parameter, the minimum uncertainty (152) may also be obtained for rescaled values of $\Delta v_{\rm max}$ and $M$, obtained by applying (V.2.2) to Eqs. (150) and (155).

These results may be combined with Bekenstein’s relation for the for the stability of charged, self-gravitating fluid spheres (Eq. (10) in the limit $R\lesssim r_{\rm grav}(M)$), by identifying $(\Delta x_{\rm total})_{\rm min}\simeq R\gtrsim Q^{2}/(Mc^{2})$. This is equivalent to assuming that the particle simultaneously saturates both the classical and quantum stability bounds, and allows us to solve the resulting equations explicitly, yielding

[TABLE]

together with

[TABLE]

where $\alpha_{\rm Q}=Q^{2}/q_{\rm Pl}^{2}$ (12). Evaluating Eq. (158) for $Q^{2}=e^{2}$ and reinserting the inequality arising from Eq. (10) then yields Eqs. (13)-(15), given in the Introduction.

For strongly interacting particles, the relevant horizon is the strong de Sitter radius $l_{\rm sd}$, defined in Eq. (V.2.1), and the relevant ‘Planck scale’ is $l_{\rm sP}$, given in Eq. (V.2.1). Modifying Eq. (154) to incorporate the ‘gravitational uncertainty’ of the strong gravity metric $f_{\mu\nu}$, by making the substitutions $l_{\rm Pl}\rightarrow l_{\rm sP}$, $l_{\rm dS}\rightarrow l_{\rm sd}$, or equivalently $G\rightarrow G_{f}$, $\Lambda\rightarrow\Lambda_{f}$, yields

[TABLE]

Combining this with the analogue of the generalised Buchdahl-Bekenstein bound for the strongly interacting fluid sphere,

[TABLE]

gives the radius of a holographic strong gravity metric ‘cell’ as

[TABLE]

We interpret the associated mass scale,

[TABLE]

as the minimum mass of a stable, strongly interacting, quantum mechanical particle with charge $Q$. Here $M_{\rm sT}$ is the strong gravity analogue of the critical holographic mass scale $M_{T}=(\hbar^{2}\sqrt{\Lambda/G})^{1/3}$ investigated in B4 .

Evaluating this for $Q=\pm(2/3)e$ gives

[TABLE]

where we have taken the mass of the up quark as $m_{\rm u}\simeq(1.7+3.3)/2=2.5$ MeV $=4.457\times 10^{-27}$ g. The numerical estimates in Eqs. (163) and (164) follow directly from the fact that, assuming $G_{f}\simeq 10^{38}G$ (3) and identifying $B\simeq 2\times 10^{14}\ {\rm gcm^{-3}}\equiv\rho_{\Lambda_{f}}$ (135), we have $\Lambda_{f}\simeq 3.351\times 10^{81}\Lambda$ (136), so that

[TABLE]

Hence, $(\sqrt{\Lambda_{f}}/G_{f})^{1/3}\simeq 8.334\times(\sqrt{\Lambda}/G)^{1/3}$, so that the estimate for the minimum mass given by Eq. (164) is only around one order of magnitude higher than the electron mass. According to the strong gravity model – combined with analogues of the hypothetical MLURs for an asymptotically de Sitter Universe derived in B2 ; LakePaterek1 ; Lake1 – this should be the minimum possible mass of a stable, strongly interacting, quantum mechanical particle of charge $Q=\pm(2/3)e$.

Reinserting the inequality stemming from Eq. (161) then yields

[TABLE]

According to this relation, if the up quark were any less massive (with the same charge $+2e/3$) or more highly charged (with the same mass $m_{u}$) a combination of canonical quantum pressure and electrostatic repulsion would overcome the strong force attraction (at the outer regions of its mass distribution), destabilising the Compton wavelength. (See analogous arguments made in B2 ; LakePaterek1 ; Lake1 for the canonical gravitational case). Equation (166) may also be rewritten as

[TABLE]

where

[TABLE]

is the classical radius of the up the quark. Numerically, this is of the order of the size of a fundamental strong gravity ‘cell’ $R_{\rm cell(s)}$, given by Eq. (162). Thus, we obtain the analogue of the maximum charge-squared to mass bound in a dark energy Universe, Eq. (14), for strongly interacting matter, and have demonstrated that the ‘strong cosmological constant’ $\Lambda_{f}$ may be expressed in a a form analogous to Eq. (15), i.e. in terms of the relevant ‘Planck length’ and the (classical) radius of a particle that saturates the upper charge-squared to mass ratio bound.

Finally, we note that, in addition to the dimensionless constant $N=(l_{\rm dS}/l_{\rm Pl})^{2}=1.030\times 10^{122}$, which may be interpreted as the ratio of the the number of ‘cells’ in the three-dimensional bulk space to the number of Planck sized ‘bits’ on the two-dimensional de Sitter boundary, we can construct the strong gravitational analogue,

[TABLE]

This implies that approximately $10-10^{2}$ ‘strong gravity cells’, each with linear dimension comparable to the classical up quark radius, $r_{u}\simeq R_{\rm cell(s)}\simeq 10^{-14}$ cm, can be packed within the strong de Sitter radius, $l_{\rm sd}\simeq l_{\Lambda_{f}}\simeq 10^{-13}$ cm. Taking the nucleon radius obtained from scattering cross sections data, $r_{n}\simeq l_{\rm sP}\simeq\tilde{l}_{\Lambda_{f}}\simeq 10^{-14}$ cm, this implies that $1-10$ such quarks can exist in a nucleon bound state. This is (obviously) consistent with known physics, but it is interesting to note that such a requirement may also be viewed as a holographic relation for the ‘de Sitter’ horizon of the strong gravity metric. In the holographic picture, the density yielding $\sim 308$ cells per cubic femtometer may be interpreted as the critical density, above which deconfinement will occur.

In addition, using purely dimensional arguments, we may also define the dimensionless constants

[TABLE]

By analogy with Eq. (V.2.1), we christen these the first and second ‘Salam numbers’, respectively, though their physical meaning (if any) remains unclear and further investigations lie beyond the scope of the present work.

Returning again to the Dirac-type relation $G_{f}\simeq 10^{38}G$ (3), we now speculate that

[TABLE]

Combining this with Eqs. (167) and (135), we then have

[TABLE]

where we have used $m_{u}\simeq 4.457\times 10^{-27}$g, as before. Substituting for $\Lambda$ from Eq. (15) then gives

[TABLE]

where we have used the previous value of $m_{u}$ together with the standard value of $m_{e}$. Thus, we have used the strong gravity theory to ‘derive’ two new large number coincidences, Eqs. (172)-(173), linking the physics of strongly interacting particles to the electroweak scale, gravity and dark energy.

V.3 The Hagedorn temperature for minimum-mass particles, the expanding Universe, and deconfinement

As shown in B1 , $M_{\Lambda}$ may also be interpreted as the effective mass of a dark energy particle. In this picture, the dark energy field is composed of a ‘sea’ of quantum particles, each occupying a volume $V_{\Lambda}\simeq\lambda_{\rm C}^{3}(m_{\Lambda})\simeq l_{\Lambda}^{3}$. Based on this, we now consider an alternative interpretation of the dark energy density and resultant late-time accelerated expansion of the Universe. Though speculative, this interpretation gives rise to a number of interesting results, and it is clear that its analogue in the strong gravity model may be relevant for hadronic physics.

If the dark energy particles are charge neutral and are their own antiparticles, then, under these conditions (in which the average inter-particle distance is comparable to the Compton wavelength $\sim\hbar/(M_{\Lambda}c)$), standard quantum theory implies that they will readily pair-produce. However, this is impossible without a concomitant expansion of space itself. In this picture, otherwise ‘empty’ space is full of dark energy particles, which give rise to an effective constant energy density on large scales. Borrowing a term from basic chemistry to describe this state, we may say that the space is ‘saturated’, and remains so as more space is produced to ‘absorb’ the newly created particles. Thus, the creation of dark energy particles and of space-like quanta go hand in hand.

It is straightforward to see that, if the probability of pair-production remains constant, the scale factor of the Universe $a(t)$ will grow exponentially, since the number of particles produced by a given volume, per unit time, is proportional to the volume itself. Let us assume that, together with the pair-production of a single dark energy particle, $n_{\rm cell}$ new fundamental ‘cells’ of space are also produced, with total volume $V=n_{\rm cell}V_{\rm cell}\simeq V_{\Lambda}$. In B2 , it was already shown that, if there exists a holographic relation between the bulk and the boundary of our (asymptotically) de Sitter Universe, such fundamental cells must have linear dimension of order

[TABLE]

This is also the Compton radius associated with the critical mass scale $M_{T}=(\hbar^{2}\sqrt{\Lambda}/G)^{1/3}$, investigated in B4 , and the Schwarzschild radius of the dual mass, $M_{T}^{{}^{\prime}}=m_{\rm Pl}^{2}/M_{T}=c(\hbar/G^{2}\sqrt{\Lambda})^{1/3}$. Equation (174) ensures that the number of cells in the bulk is equal to the number of Planck sized bits on the de Sitter boundary, i.e.

[TABLE]

Next, let us suppose that the probability of a single cell of space ‘pair-producing’ within a time $t_{\rm Pl}=l_{\rm Pl}/c$, due to the presence of the dark energy density, is given by

[TABLE]

This leads naturally to a de Sitter-type expansion, modelled by the differential equation

[TABLE]

or, equivalently,

[TABLE]

In this picture, the macroscopic dark energy energy density $\rho_{\Lambda}$ remains approximately constant, in spite of spatial expansion, the additional (positive) mass-energy of a newly created dark energy particle being exactly counterbalanced by the additional (negative) energy contained in its gravitational field. This may be shown explicitly by considering the Komar energy (see B1 ; B2 ; B3 ).

However, if this picture is correct, we may expect ‘empty’ three-dimensional space to exhibit granularity on scales $\sim l_{\Lambda}$. It is therefore particularly intriguing that recent experiments provide tentative hints of fluctuations in the gravitational field strength on scales comparable to $l_{\Lambda}=\hbar/(M_{\Lambda}c)$, which is of order $0.1$ mm Perivolaropoulos:2016ucs ; Antoniou:2017mhs . Though many theoretical models may account for this, including those exhibiting spatial variation of Newton’s constant $G$, the results presented above imply that the ‘granularity’ of the dark energy density, due to the presence of effective dark energy particles on sub-millimetre scales, cannot be discounted a priori.

In this model, the number of holographic spatial cells created when one dark energy particle is pair-produced is $n_{\rm cell}\simeq(\hbar G\Lambda/3c^{3})^{1/4}=N^{1/4}\simeq 3.186\times 10^{30}$. As already noted in B1 , this number is also the multiplying factor that naturally generates a sequence of mass scales between $m_{\rm dS}$ and $m^{\prime}_{\rm dS}$, i.e. $m^{\prime}_{\rm dS}=N^{1/4}M^{\prime}_{\Lambda}=N^{1/2}m_{\rm Pl}=N^{3/4}M_{\Lambda}=Nm_{\rm dS}$, where $M^{\prime}_{\Lambda}\equiv m_{\rm Pl}^{2}/M_{\Lambda}$.

In addition, it is clear that the fundamental field giving rise to the dark energy density (whatever its precise nature may be) remains ‘trapped’ in a Hagedorn phase. Any attempt to further compress (i.e. heat) the ‘sea’ of dark energy particles – even if such compression results simply from random quantum fluctuations – results in pair-production rather than increased kinetic energy. The saturation condition implies the existence of not-so-UV cut-off for the vacuum field modes, given by $\lambda_{\rm DE}\simeq\lambda_{\rm C}(M_{\Lambda})=l_{\Lambda}$, yielding

[TABLE]

Thus, the temperature associated with the field remains constant, on large scales, and is comparable to the present day temperature of the CMB B1 ,

[TABLE]

Here the factor of $(8\pi)^{-1}$ is included in the definition of $T_{\Lambda}$ by analogy with the standard expression for the Hawking temperature, yielding $T_{\Lambda}(M_{\Lambda})=T_{\rm H}(M^{\prime}_{\Lambda})$, where $M^{\prime}_{\Lambda}=m_{\rm Pl}^{2}/M_{\Lambda}$ again denotes the dual mass.

Though this too may seem like another incredible coincidence, we note that, in the dark energy model implied by the DE-UP (154), it is simply a restatement of the standard coincidence problem of cosmology, i.e. the Universe begins a phase of accelerated expansion when $r_{\rm U}\simeq l_{\rm dS}$, at which point $\Omega_{\rm M}\simeq\Omega_{\Lambda}$ and, hence, $T_{\rm CMB}\simeq T_{\Lambda}$. The question remains, why do we live at precisely this epoch? However, no new coincidences are required, in order to explain Eq. (180).

The implications of this picture for the dual strong gravity model are self-evident. In this case, the mass associated with the dark energy ‘cell’ $M_{\Lambda}$ is replaced by the nucleon mass $m_{n}\simeq M_{\Lambda_{f}}\simeq\tilde{M}_{\Lambda_{f}}$. When compressed beyond the nuclear density, a free quark fluid is formed, and further attempts at compression (i.e. heating) simply result in the production of more strongly interacting matter (quark-gluon plasma). The free strange quark matter remains ‘locked’ in a Hagedorn phase and the temperature of the plasma remains constant, given by

[TABLE]

This is consistent with previous estimates for the temperature of the deconfinement transition obtained in Sec. IV.

VI Discussions and final remarks

In the present paper, we have considered the mass/radius ratio bounds for spherical compact objects with anisotropic pressure distributions in the strong gravity model, which represents an attempt to describe the gauge singlet sector of the strong interaction using a ‘geometric’ theory, based on analogy with general relativity. Though the theory cannot describe the $SU(3)$ colour charge sector of QCD, it remains a viable candidate for an effective theory which may be used to model stability conditions and confinement in strongly interacting particles.

Strong gravity is a two tensor theory, in which the canonical (weak) gravitational field is described by the usual spacetime metric $g_{\mu\nu}$ and the strong interaction is described by an additional metric-type tensor $f_{\mu\nu}$. The strong gravity action contains the usual Einstein-Hilbert term $(1/2k_{f}^{2})R(g)\sqrt{-g}$, where $k_{g}^{2}=8\pi G/c^{4}$ plus an additional ‘copy’ constructed from the strong tensor field, $(1/2k_{f}^{2})R(f)\sqrt{-f}$, where $k_{f}^{2}=8\pi G_{f}/c^{4}$ and $G_{f}\simeq 10^{38}G$ is the strong gravity coupling constant. The strong gravity Lagrangian also includes an interaction term $\mathcal{L}_{fg}$ and the standard matter term $\mathcal{L}_{m}$. The interaction term is proportional to $\mathcal{M}^{2}$, where $\mathcal{M}$ is the mass of the spin-2 ‘strong graviton’.

An appropriate choice of ‘gauge’ for the mixing term generates an effective strong gravitational constant, $\Lambda_{f}\propto\mathcal{M}^{2}$. Hence, in the strong gravity theory, there exist analogues of many results that can be derived from canonical general relativity with a ‘true’ cosmological constant term, $\Lambda$. These include mass bounds obtained using the appropriate Buchdahl-type inequalities for the physical system under consideration. Therefore, in strong gravity, it is possible to obtain explicit inequalities giving upper and lower bounds on the ratio $M_{{\rm eff}}/R$, where $M_{\rm eff}$ is the effective mass of a compact object, including the contribution from $\Lambda_{f}$. Alternatively, these may be reformulated as bounds on $M_{0}/R$, where $M_{0}$ is the bare mass, and the inequality is written as an explicit function of the strong cosmological constant and the (pressure) anisotropy parameter $\mathcal{D}$, which also depends on $\Lambda_{f}$.

As is the case for compact objects in general relativity, in the presence of a nonzero cosmological constant ($\Lambda\neq 0$), we found two universal limits (upper and lower) for the mass/radius ratio of strongly interacting particles. However, due to the presence of the strong cosmological constant ($\Lambda_{f}\neq 0$) and of the anisotropies in the pressure distribution, the physical and geometric properties of such hadronic-type compact objects are significantly modified within the particle interior, as compared to their counterparts in standard general relativity. Both the upper and lower mass/radius ratios depend sensitively on the values of $\Lambda_{f}$ and $\mathcal{D}$ at the surface of the hadron.

In addition, different physical models for the mixing term that generates the effective strong cosmological constant can lead to very different mass/radius relations. It is a general feature of the behaviour of the physical and geometrical parameters of anisotropic objects in strong gravity that the increase in mass is proportional to the deviations from isotropy, described in our approach by the function $\mathcal{D}$. Since these deviations from isotropy are arbitrary, there are no mathematical or theoretical restrictions that restrict the radii of hadronic-type strong gravity structures, which may therefore extend up to the ‘apparent horizon’ of the strong tensor field. Under the assumption that the exterior of the hadronic objects is described by the strong Schwarzschild metric, this corresponds (approximately) to the mass $M_{{\rm eff}}\leq c^{2}R/(2G_{f})$.

Usually, the upper mass/radius bound can be obtained when we assume the ultra-stiff equation of state, i.e. constant density inside the sphere. For the minimum mass/radius bound, a straightforward interpretation (in the de Sitter case) is the situation where we have constant pressure equal to the vacuum pressure throughout the sphere so that the pressure is always balanced inside and outside. Constant pressure implies $\rho=-P$ inside the sphere and this is simply the vacuum density. The minimum mass/radius ratio actually implies the minimum density of the object, which is equal to the vacuum density due to the presence of dark energy. It actually implies that no particle should have smaller density than the vacuum density.

The Schwarzschild-de Sitter-type solution of the strong gravity field equations describes a microscopic system embedded in an ordinary, flat space-time, in which the mass of compact coloured objects is localized due to the ‘curvature’ of the strong metric field, which creates a kind of ‘bag’ sg3 . By interpreting the energy density of the strong gravity cosmological constant as the bag constant of QCD Shur , it follows that strong gravity imposes the following classical lower bound on the minimum mass/radius ratio of a hadron,

[TABLE]

This is related to the bag constant $B$, via $B\simeq\rho_{\Lambda_{f}}\simeq\Lambda_{f}c^{2}/(8\pi G_{f})$ (135). By assuming a hadronic radius of the order of $R=0.8$ fm (comparable to the proton radius Pohl ), and taking the estimate $\Lambda_{f}\simeq 10^{25}$ cm*-2*, obtained from the identification (135) together with Eq. (3), we obtain a lower bound on the bare mass of a hadron as $M_{0}\gtrsim 10^{-24}$ g, a value which is of the same order of magnitude as the mass of a nucleon. Smaller particle radii, of the order of $0.1$ fm, will give considerable lower hadron masses, of the order of $M_{0}\geq 3.37\times 10^{-28}$ g.

All results regarding the mass/radius ratios for anisotropic hadronic objects have been obtained by assuming the basic conditions (51) and (52). However, for an arbitrarily large anisotropy parameter $\mathcal{D}$, with $P_{r}\gg P_{\perp}$, we can not exclude (in principle) the situation in which these conditions do not hold. If, for example

[TABLE]

then, for a hadron in strong gravity with monotonically decreasing density, the condition

[TABLE]

holds in place of Eq. (52). This situation corresponds to a tangential pressure-dominated hadronic structure, with the tangential increasing inside the compact object. In this case, we obtain a restriction on the minimum mass/radius ratio so that, for this hypothetical, ultra-compact hadronic particle, $4/9$ is an absolute lower bound for the value of $G_{f}M/(c^{2}R)$.

In addition, we have investigated possible quantum mechanical implications of the strong gravity model, for both neutral and charged particles. For neutral particles, the quantum minimum mass bound follows by identifying the classical radius $R$ in Eq. (181) with the Compton wavelength $\lambda_{\rm C}$. The mass scale thus obtained is roughly comparable to the mass of the neutron, $m_{n}\simeq 10^{-24}$g, the lightest known stable, compact, charge-neutral and strongly interacting particle.

To treat charged hadronic objects, we combined classical stability bounds for charged compact objects in strong gravity, obtained by substituting $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$ into their general-relativistic counterparts, with hypothetical minimum length uncertainty relations (MLURs). These, in turn, were based on MLURs obtained by considering canonical gravitational ‘corrections’ to the standard Heisenberg uncertainty principle, including the effects of dark energy and the existence of a de Sitter horizon B2 ; LakePaterek1 ; Lake1 . The formal similarity between general relativity and the strong gravity theory again allowed us to replace $G\rightarrow G_{f}$ and $\Lambda\rightarrow\Lambda_{f}$, yielding analogous MLURs for strongly interacting particles. Identifying $R\simeq(\Delta x_{\rm total})_{\rm min}$, and evaluating the strong gravity MLUR Eq. (160) for $Q=\pm 2e/3$, we obtained an estimate of the mass of the up quark, $m_{u}\simeq 10^{-27}$g, the lightest known stable, compact, charged and strongly interacting particle. This estimate is equivalent to a new large number coincidence, ‘derived’ from the quantum mechanical MLUR/strong gravity model, which relates the nuclear density, dark energy density, and physics at the electroweak scale, Eq. (172).

The formal equivalence between the mathematical structure of the strong gravity theory and canonical general relativity also permits us to draw parallels between the strong gravity model of quark deconfinement and the expansion of the Universe in the particle ‘sea’ model of dark energy, proposed in B1 . In the former, an expanding deconfined quark matter remains ‘trapped’ in a Hagedorn phase, in which further compression of the quark-gluon plasma, even if this arises as a result of random quantum fluctuations, leads to pair-production rather than increased temperature. By interpreting the minimum mass of a stable, compact, charge-neutral, and quantum mechanical object as the mass of an effective dark energy particle B1 , we obtained the resulting ‘Hagedorn temperature’ ($T_{\Lambda}$) of the dark energy field, which was found to be comparable to the present day temperature of the CMB (180). In this model, such a ‘coincidence’ is not really a new coincidence at all, but simply a restatement of the standard coincidence problem in cosmology, whereby a phase of accelerated expansion begins when $r_{\rm U}\simeq 1/{\sqrt{\Lambda}}$ and $\Omega_{\rm M}\simeq\Omega_{\Lambda}$, or, equivalently, $T_{\rm CMB}\simeq T_{\Lambda}$.

Finally, we note again that, since there is no explicit SU(3) gauge symmetry in the strong gravity field equations, these may describe only the gauge singlet sector of the strong interaction, mediated by massless and massive spin-2 particles, coupled to the energy-momentum tensor of the strongly interacting matter. Hence, strong gravity is not expected to replace QCD, but may be used to describe interactions involving only gauge singlet states, using a gravitational type formalism, though not the sector including colour charges. It is therefore justified to use strong gravity to explore the stability and confinement of gauge singlet mesons and baryons, but not scattering processes that require colour charge interactions.

Acknowledgements.

We would like to thank to the anonymous referee for comments and suggestions that helped us to improve our manuscript. P.B. is supported in part by the Thailand Research Fund (TRF) and Chulalongkorn University under grant RSA5780002. M. L. was supported by a Naresuan University Research Fund Individual Research Grant during the preparation of this work and by the gracious hospitality of Sun Yat-Sen University and Nanyang Technological University.