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

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) color charge. It may therefore be used to explore the stability and confinement of gauge singlet hadrons, though not to describe scattering processes that require color 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 1GeV-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1~{\mathrm{GeV}}^{-1}$$\end{document}. We revisit the strong gravity theory and investigate the strong gravity field equations in the presence of a mixing term which induces an effective strong cosmological constant, Λf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _{f}$$\end{document}. This introduces a strong de Sitter radius for strongly interacting fermions, producing a confining bubble, which allows us to identify Λf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _{f}$$\end{document} with the ‘bag constant’ of the MIT bag model, B≃2×1014gcm-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B \simeq 2 \times 10^{14}~{\mathrm{g}}~{\mathrm{cm}}^{-3}$$\end{document}. 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 generalized 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 Λf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _f$$\end{document}, 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.


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 linearized gravity, and the quantum theory of closed bosonic strings naturally includes a candidate graviton [1]. 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 G E 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 / √ G 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 [2,3].
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) [4]. Strong support for this idea is provided by the fact that the quantization of linearized gravity, obtained by substituting g μν = η μν + h μν , where η μν is the Minkowski metric and h μν 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, h μν = 0, h μν,ν = 0, h μ μ = 0 [4]. 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 [5][6][7].
Classically, the (non-vacuum) Einstein field equations, G μν = κ T μν , where κ = 8π G/c 4 , relate a covariant geometrical quantity, the Einstein tensor G μν = R μν − (1/2)Rg μν , to a covariant physical quantity, the conserved energymomentum tensor T μν , via the proportionality (coupling) constant κ [8]. However, the derivation of the Einstein field equations does not place any restriction whatsoever on the numerical value of the constant κ. For the canonical gravita-tional interaction, this must be recovered from the Newtonian limit of the theory [8]. Noting the existence of strongly interacting, massive, spin-2 meson states (such as the f-meson), and arguing by analogy with the quantization of linearized Einstein gravity, it was proposed in [9][10][11][12]] 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 μν , which determines the properties of the strong gravity field, as well as a new strong coupling constant, determined to be of order 1 GeV −1 to ensure consistency with the known physics of strong interactions, were introduced [9][10][11][12]. 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 respectively. The dimensionless strong gravity coupling is taken to be of the same order of magnitude as the strong interaction coupling, giving This is equivalent to G f 10 38 G = 6.67×10 30 cm 3 g −1 s −2 for a strong interaction scale of m 10 −19 m Pl 10 −24 g, where m Pl = √h c/G 10 −5 g is the Planck mass. It is interesting to note that the ratio G f /G yields 'Dirac's number', i.e.
the same dimensionless quantity that formed the basis of his Large Number Hypothesis (LNH) [13][14][15][16]. (See also [17][18][19] 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 Sect. 4 Hence, in order to describe strong interactions, the following Lagrangian density was proposed [9][10][11][12]: (4) where L f = R f ( f ) and R f ( f ) is the scalar curvature constructed from f μν and its derivatives. The term L f g describes the interaction between f-mesons and gravitons and L m is the Lagrangian for both strongly interacting and non-strongly interacting matter. If one drops the interaction term L f g and considers that the dominant term in the field equation obtained from (4) is given by L f , one obtains an Einsteintype equation for the f μν field, with the strong coupling constant k 2 f in place of k 2 g . This explains the name strong gravity given to the theory. The physical and mathematical properties of the strong gravity model were investigated in [20][21][22][23][24]. For a comprehensive review of the early results in strong gravity, and the corresponding references, see [25].
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 [26,27]. In this model, QCD in the IR region is approximated by the exchange of a dressed two-gluon phenomenological field G μν (x) = B a μ B b ν η ab , where η ab is a color-SU(3) metric, and B a μ is the dressed gluon field. This model produces color 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 it 'predicts' scaling. The effective Ricci tensor R μν constructed from the field G μν induces Einsteinian dynamics. Alternative approaches were proposed in [28][29][30] and [31], respectively. The relations between Yang-Mills fields and Riemannian geometry were also investigated in [32,33], 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 [34] states that, for stable, compact, charge-neutral objects of mass M and radius R, the condition must be satisfied. If the bound (5) is violated, collapse to a black hole becomes inevitable. In [35] the general-relativistic maximum mass of a stable compact astrophysical stellartype object was found to be of the order of 3.2M . 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 ( = 0) [36], of the electric charge of the sphere [37], and of an anisotropic interior pressure distribution [38]. Sharp bounds on the mass/radius ratio for neutral and charged compact objects, with both isotropic and anisotropic pressure distributions, in the presence of = 0, were also obtained in [39][40][41][42][43].
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 [44] 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 general-ized Buchdahl inequality in the presence of > 0, given by 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.
Though the derivation of this condition is somewhat involved [44], its physical meaning is intuitively obvious. The dark energy density is given by 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 ρ ρ min = ρ /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 ρ = 5.971 × 10 −30 g cm −3 [45][46][47][48][49]. This corresponds to a value of = 1.114 × 10 −56 cm −2 for the cosmological constant.
The minimum-mass/radius bound, in the absence of dark energy (i.e., for = 0), was extended for the case of charged objects in [50], where it was shown that a stable object with charge Q must obey the relation This result was extended to include the presence of a cosmological constant ( = 0), and the mass/radius ratio of a stable charged object was found to obey the relation to leading order in Q 2 /R. For = 0, this relation recovers an earlier result due to Bekenstein [51], 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 [52][53][54], 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 where and q Pl =hc is the Planck charge. Evaluated for Q = ±e, this gives where α e = e 2 /q 2 Pl 1/137 is the usual fine structure constant and m e is the electron mass, which is equivalent to 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 destabilize its Compton wavelength [52][53][54]. We also note the close similarity between the physical picture of the electron, modeled as an extended charged fluid sphere [51], used to derive Eq. (14) in [52], and Dirac's 'extensible' model of the electron, proposed in [55]. Equation (14) may also be rewritten as yielding an upper bound for which is consistent with the current best fit value inferred from cosmological observations [45][46][47][48][49]. Interestingly, Eq. (15) was previously derived via three different methods (see [56][57][58]) and, for the mass scale m = m e /α e , is equivalent to Zel'dovich's estimate of , based on his reformulation of Dirac's LNH for an asymptotically de Sitter Universe [59,60]. 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 they are discussed further in the context of the strong gravity model, with a 'strong cosmological constant' f , in Sect. 4.
The particle physics and cosmological implications of the mass scale M T = (h 2 √ /G) 1/3 , which corresponds to taking Q 2 → q 2 Pl in Eq. (11), were considered in [61], 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 = m 2 Pl /M T = c(h/G 2 √ ) 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 [62][63][64]. 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 [65][66][67]. 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 [68], 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 selfinteracting 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 [9][10][11][12], 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 energymomentum 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, 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 behavior 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 L f g , as well as of the geometric properties of the hadrons. The physical implications of our results are also discussed.
This paper is organized as follows. The TOV equation for strong gravity is derived in Sect. 2. The generalized Buchdahl inequality, and the resulting upper and lower bounds on the mass/radius ratio of strongly interacting particles, are derived in Sect. 3. The strong gravity mass gap, and its implications for holography, are discussed in Sect. 4. In Sect. 5, we combine the mass bounds obtained in the strong gravity model with MLURs motivated by quantum gravity research, exchanging G → G f and → 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 2 × 10 14 g cm −3 , then gives rise to new mass bounds for both neutral and charged strongly interacting particles. Section 6 contains a brief summary and discussion of our main conclusions, of outstanding problems, and of prospects for future work.

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 Einsteintype 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 1 GeV −1 , which replaces the Newtonian coupling k g 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 generalized 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 [8].

Strong gravity
As stated in the Introduction, the Lagrangian for the interacting strong field metric f μν and gravitational metric g μν can be constructed as (16) where the corresponding 'volume element' Here the first term represents the standard general-relativistic Lagrangian for the gravitational field, while the second is its strong interaction analog, obtained by replacing k g by k f and g μν by f μν . To give the 'elementary' particles mass, as well as their weak gravitational interaction, a mixing term between the strong and weak gravitational fields, L f g , is needed. 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 μν and its derivative only, whereas the former may depend (generically) on both g μν and f μν . Hence, although the strong gravity metric minimally couples to all forms of matter (see below), strongly interacting particles and nonstrongly interacting particles 'feel' the curvature of the strong metric differently. A simple covariant mixing term was proposed in [21] and is given by where M is a constant with the dimension of mass. For later use, the full dimension of the mass mixing parameter is given by M 2 → M 2 c 2 /h 2 , the inverse Compton wavelength squared.
In the limit in which the gravitational field may be ignored, g μν → η μν , the field equations of the strong gravity theory can be written as where 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 {t, r, θ, φ} the line element with respect to the strong gravity metric f μν is assumed to be of the form where d 2 = dθ 2 +sin 2 θ dφ 2 is the line element for the unit 2-sphere and ν(r ), λ(r ) and (r ) are arbitrary functions of the radial coordinate. Furthermore, we assume that the fluid can be described by the standard energy-momentum tensor where ρc 2 is the fluid energy density, P r is the radial pressure and P ⊥ denotes the tangential pressure. Ignoring the weak gravitational interaction, the field equations for the strong gravity field coupled to the matter fluid are, therefore, where the raising and lowering of space-time-like indices is performed via the tensors f νλ and f μν , respectively, satisfying the condition Note the universality of coupling between strong gravity and any form of matter energy-momentum tensor.

Strong gravity field equations
The stress tensor of the mixing term between massive and massless gravitons can be computed straightforwardly, giving For simplicity, and also to make T t (s) t interpretable as a proper density, we fix the gauge so that [21] = 2r 2 3 , e (ν+λ) = = constant > 0.
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 [52]. With this choice of gauge, there is an anisotropic Poincaré stress associated with the θ -component of the energy-momentum tensor, given by Moreover, from the second of Eq. (27) it follows that, since (ν + λ ) = 0 (where a prime denotes differentiation with respect to r ), the functions ν and λ satisfy the condition at all points inside the strongly interacting fluid sphere. The field equations of the strong gravity model may then be written in the form where the last equation is simply the conservation law, 0 = f κμ ∇ κ T μν , with respect to the strong gravity metric f μν .

The Tolman-Oppenheimer-Volkoff equation in strong gravity
Equation (31) can be directly integrated to give where the accumulated mass inside radius r is defined by Substituting ν 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 Note that, by taking into account that ν + λ = 0, Eq. (37) can also be written as which, combined with Eq. (35), gives an alternative definition of the mass of the fluid,

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 generalized 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 easily be obtained.

The Buchdahl inequality
In order to obtain the generalized Buchdahl inequality for strong gravity, we define the following Buchdahl variables: where From Eqs. (35), (37) and (38), we then have where we have defined 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 . Further manipulation then leads to To separate positive and negative terms on the right-hand side of Eq. (45), we introduce two new quantities, γ and γ − , defined as From Eq. (35), it follows that, since e λ ≥ 2/3 and e ν ≥ 0, γ − is always negative. On the other hand, γ is positive definite for P ⊥ > P r −27M 2 /16 3/2 k 2 f . For a static charged sphere, P ⊥ = −P r = Q 2 (r )/2r 4 [52], where Q(r ) is the accumulated charge, this condition is valid as long as the pressure from the mass mixing contribution satisfies Note that the mixing term could be negative. Equation (45) can be further simplified by defining to obtain Assuming monotonically decreasing profiles for both the density ρ(r ) and γ (r ), it follows that for all r > r the conditions must hold at all points inside the compact object. We then immediately obtain w (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 ψ with respect to z: For Using condition (51), we find where and is the effective cosmological constant of the strong gravity model.
Finally, we substitute Eqs. (54), (57) and (37) into the inequality (53) and divide by ζ to obtain the Buchdahl inequality in strong gravity as where D ≡ P ⊥ − P r + 27 16 Equation (58) is valid for all r inside the strong gravity particle. Moreover, its validity does not depend on the sign of D.
3.2 The minimum-and maximum-mass/radius ratios of hadrons

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 D = 0, corresponding to the condition Moreover, we assume that the effective pressure also vanishes at the surface of the massive particle, so that By evaluating Eq. (58) at the surface of the hadron r = R, using the conditions (61), we obtain where we have denoted M eff = m eff (R), leading to the well-known result 2M eff /R ≤ 8/9 [34]. 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 eff and radius R can be written as For G f = 6.67×10 30 cm 3 g −1 s −2 , this relation is obviously satisfied in the case of proton, with mass m p = 1.672 × 10 −24 g and classical radius r p = 0.875 × 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 eff = 0. In the quasi-isotropic limit D = 0, Eq. (58) gives the upper mass-radius bound where we have defined the mean density of the compact object as ρ eff = M eff /R 3 . We assume that the matter radial pressureP r vanishes at the surface of the strong gravity par-ticle, and thus we obtain for the surface effective pressure the expression Taking into account that where we have defined the mean fluid density as , we obtain the following upper limit for the ordinary matter mass/radius ratio of a stable compact object in strong gravity: Next, we consider the anisotropic case with D = 0. We define the function f (M eff , R, f , D) as Assuming again that the effective pressure vanishes at the surface of the compact object, P eff ≡ 0, Eq. (58) leads to the following general restriction on the mass/radius ratio for a spherical hadronic fluid: By taking into account the definition of the total effective mass as given by Eq. (66), we immediately find With the use of the Taylor series expansion of the function arcsin x/x − 1, for small values of the argument, we can approximate the Then the maximum-mass bound for compact objects in strong gravity with vanishing surface effective pressure can be reformulated as Finally, we consider the case of strong gravity compact objects with vanishing surface radial pressure, i.e. with P r (R) = 0, and P eff = 0, given by Eq. (65). In this case, we obtain for the maximum-mass/radius bounds the expressions and respectively. The maximum-mass/radius bound for compact objects is generally obtained in the constant density regime, with ρ eff ρ 0 = constant. Hence the total mass of the compact object can be approximated as M 0 = 4πρ 0 R 3 /3. Therefore, in the mass/radius ratio bounds obtained above, we can approximate the effective mean density as ρ eff M 0 /R 3 ρ 0 = constant, a relation that is satisfied by the maximum-mass objects with a very good approximation. Therefore, in all the above results, the ratio f / ρ eff 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.

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 where we have assumed that k 2 f P eff /3 + M 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 eff (R) = −3 f /k 2 f , and therefore P 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 eff → −P eff , Eq. (78) can be written as By introducing a new variable v, defined as Eq. (79) takes the form where we have denoted and respectively. Then, by squaring Eq. (81), we can reformulate the corresponding inequality as where we have denoted and In the following analysis, we keep only the first order terms in both p (depending on f ) and q (depending on D) in the expressions involving square roots. Since v ≥ 0, Eq. (85) is However, one can easily check that the condition v ≥ 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 ≥ 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 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: By taking into account the explicit expression for the effective mass of the hadron, Eq. (89) may be rewritten as 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 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, M, and the arbitrary constant , which fixes the value of the metric tensor coefficients inside the hadrons. For f ≡ 0, the minimum-mass bound simply reduces to the positivity requirement for the bare mass, M 0 ≥ 0.
By taking into account the explicit expression for the strong gravity cosmological constant and by assuming that (k 2 f /9)D(R)R 2 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 This equation imposes strong (pun intended) constraints on the mass parameter M, which controls the strength of the mixing between the f μν and g μν fields, i.e.
where ρ min = M (min) 0 /R 3 is the density corresponding to the minimum-mass hadron.

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 μν , by analogy with general relativity, as [69,70] where the vector ξ i is a time translation Killing field, ξ s denotes its value at S, and [K ] is the jump across the shell of the trace of the extrinsic curvature of S, assumed to be embedded in the 2-space t = constant. E M = S T k i ξ i √ −gdS k and E F are the energies of the ordinary matter and of the strong gravitational field in the f μν metric, respectively. This definition of the total energy is manifestly coordinate invariant. In the case of static spherical symmetry, for both the g μν and the f μν fields, we obtain the total energy for the hadron, from Eq. (93) [70]: 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 f → . Next, for the sake of convenience, we rescale the metric tensor component so that (2/3)e −λ → e −λ . Equation (35), which may be expressed in terms of the effective mass of the strong gravity object M eff , can then be written as Such a rescaling is always possible and is performed explicitly via coordinate transformation r 2 → r 2 = 2r 2 /3.
Note that the mass per unit radial distance transforms as M(r )/r = 2M(r )/3r . For convenience, we will also redefine f → f /3. Then, by taking into account the relation ν+λ = 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 or, equivalently, With the use of Eq. (58), we find the following upper limit for the total energy of the compact hadronic object in strong gravity: In the case of a vanishing strong cosmological constant, f → 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 [69,70]: For a quasi-isotropic matter distribution with D = 0, with the assumption of vanishing radial pressure P r (R) = 0, we obtain Equation (97) also allows one 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 ∂ E/∂ R = 0. This gives the following mass/radius relation, as a function of the strong cosmological constant, which is valid for 2 f R 2 > 1: The two values of the mass differ by a quantity

The mass gap in strong gravity
In Sect. 3.2.1, we have considered the upper bound on the mass/radius ratio in the situation where P eff (R) = 0 and D = 0, P eff = 0. In Sect. 3.2.2, the lower bound is derived when we set P r (R) = 0 in order to highlight the effect of anisotropic parameter D on the bound. In this section, we will consider the most generic case, without making any assumptions as regards the value of pressure at the surface of the object, and we will simply define P eff (R) ≡ P eff . In order to simplify our formalism, we define the additional dimensionless quantities which allows us to express the inequality (58) at r = R as This may be written in the following, explicitly quadratic form: The inequality can be satisfied if for for B < 0. For B > 0, the mass bounds cease to exist. The condition B < 0 leads to the following constraint on F: Another condition for the existence of mass bounds is B 2 > 4C, which is trivially satisfied for a > 0. However, for a < 0, B 2 > 4C requires To simplify these results, we consider the case when 4C/B 2 1 and also a, F 1, to obtain Under these conditions, the contribution from the anisotropic stress, F = k 2 f 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 The nontrivial minimum-mass/radius exists when P eff < −2 f /k 2 f . This is valid even for 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 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, λ C , equal to R in Eq. (114) with P eff , f , D = 0. We then have in the standard units. We define M max to be the corresponding 'Planck mass' of the strong gravity. Whereas, for M > m Pl , fundamental particles inevitably collapse to form black holes, strongly interacting particles with M > M max inevitably undergo a deconfinement phase transition. The corresponding length scale, given by a Compton relation, R min =h/M 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 M (i.e. the mass mixing term is proportional to M 2 ), while the mass gap itself, given by Eq. (113), is proportional to If we set the mass gap equal to the mass of the π -meson, m gap = m π 140 MeV, R R f /2 and M = 2 GeV, the strength of strong gravity becomes M max 3.6 GeV. (117) Any quantum particle with a strong gravity interaction will inevitably collapse to form a strong gravity 'black hole' once the mass exceeds M max . The corresponding Hawking temperature of the strong gravity black hole is 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 μν 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 a spherical object in AdS space and the deconfinement temperature of the dual gauge matter is well known [71][72][73][74]. 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 space-time. It is interesting to note similar features of the strong gravity model.
In [52,68], 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 the 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.) 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 + ρ = 0, which originates from the choice = 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 becomes physically relevant and we need to allow to depend on the radial coordinate r , setting = const. → (r ). We then have resulting in an increasing function (r ) ≥ 1 for a positive (P r + ρ) matter profile. The mass bound analysis above can then be repeated with replacing ρ in Eq. (31) and by setting b = 0. The resulting minimum-and maximum-mass/radius ratio bounds are exactly the same as in Eqs. (113) and (114) under the replace- 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 = r/ √ 3/2 and redefine f → f /3 so that our metric is in the conventional form.

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 = 0, before extending these to the strong gravity case via the substitutions G → G f , → f .

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 λ C , or reduced Compton wavelength, k −1 C , given by 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, where, for future reference, we define the (reduced) Planck mass and length scales m Pl = √h c/G = 2.176 × 10 −5 g, and l Pl = hG/c 3 = 1.616 × 10 −33 cm, respectively, and the (reduced) first/second de Sitter mass and length scales, denoted using unprimed/primed quantities, respectively, as Note that in previous work [52,68,75] the de Sitter scales defined in Eq. (123) were also referred to as the first and second Wesson scales, following the pioneering work [76]. The physical interpretations of these scales are discussed in detail in [75]. For now, we simply note that the first and second de Sitter scales are related via and that the numerical value of the minimum mass, m 10 −3 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 ν ≤ 0.23 eV [45][46][47][48]. To within numerical factors of order unity, m is also the unique mass scale for which the Compton radius of the particle is equal to its gravitational turn-around radius, in the presence of dark energy [53,54,77]. 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 Schwarzschildde Sitter metric, which gives rise to an effective Newtonian potential of the form This, in turn, gives the effective gravitational field strength [78] which changes sign at r = r grav . However, to within numerical factors of order unity, Eq. (125) remains rigorously valid in full general relativity [77]. 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 m dS = (c 2 /G) √ 3/ , which is comparable to the present day mass of the Universe [52,75].

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 μν field gives where m sP and l sP denote the strong gravity Planck mass and strong gravity Planck length, respectively. Note that m sP and l sP are equivalent to M max and R min , defined in Sect. 4. We here relabel these quantities for the sake of easy comparison with results from standard general relativity. Likewise, analogs of the de Sitter scales may be obtained by replacing G → G f and → f in Eq. (123). In addition, based on purely dimensional arguments, we may define two additional mass scales, and their corresponding lengths, by mixing and matching {G, , G f , f }.
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 M 2 in the mixing Lagrangian L f g . This appears in the strong gravity field equations in combination with the 'geometric' parameter , through the definition of the strong cosmological constant, (56). By analogy with Eq. (8), we define the energy density associated with f as 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 where B 2 × 10 14 g cm −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 [21], 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 μν field gives 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 f 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 1.675 × 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 λ n = λ C (m n ) = 2.100 × 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, σ n = πr 2 n 10 −24 cm −2 , yielding r n 5.642×10 −13 cm 26.864×λ n . Identifying the radius R in Eq. (140) with the scattering radius of the particle r scat , yields an alternative estimate of the lightest neutral hadron in the strong gravity theory, which we denoteM f . This is given by Finally, we note that, due to the approximate numerical coincidence of m sP and m sd , we have respectively, from now on.

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 Refs. [52][53][54]. (Note that Eq. (15) was also derived, using different methods, in [56,58] and that its cosmological implications were investigated in [79], while similar expressions, equivalent to replacing m e /α e ↔ m, were obtained in [59,60,80].) As we shall see, the approach taken in [52][53][54], 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 → G f and → f . Having obtained these, we again identify the energy density associated with f , ρ f ≡ f c 2 /(8π G f ), with the 'bag constant' of the MIT bag model, B 2×10 14 g cm −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 = ±2e/3, this is found to be of the same order of magnitude as the mass 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 emphasize that, when estimating 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 ρ n 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 [81,82], an MLUR of the form where α , β = const., was proposed. Here x 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 p = M v and λ C k −1 C =h/(Mc), whereas the second represents a recoil term, due to the emission of the photon [83]. The third is the 'gravitational uncertainty', which is assumed to be of the order of the Schwarzschild radius r S l 2 Pl /λ C [81,82]. The second line is obtained by minimizing the first with respect to v, Minimizing the expression on the second line of Eq. (149) with respect to M then yields and hence The quantity 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 as the minimum mass of a stable, strongly interacting, quantum mechanical particle with charge Q. Here M sT is the strong gravity analog of the critical holographic mass scale M T = (h 2 √ /G) 1/3 investigated in [61].
Hence, ( f /G f ) 1/3 8.334 × ( √ /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 analogs of the hypothetical MLURs for an asymptotically de Sitter Universe derived in [52][53][54] -this should be the minimum possible mass of a stable, strongly interacting, quantum mechanical particle of charge Q = ±(2/3)e.
Reinserting the inequality stemming from Eq. (161) then yields 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), destabilizing the Compton wavelength. (See analogous arguments made in [52][53][54] for the canonical gravitational case.) Equation (166) may also be rewritten as 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 cell(s) , given by Eq. (162). Thus, we obtain the analog 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' f may be expressed in 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 dS /l Pl ) 2 = 1.030 × 10 122 , which may be interpreted as the ratio of the number of 'cells' in the threedimensional bulk space to the number of Planck sized 'bits' on the two-dimensional de Sitter boundary, we can construct the strong gravitational analog, This implies that approximately 10-10 2 'strong gravity cells', each with linear dimension comparable to the classical up quark radius, r u R cell(s) 10 −14 cm, can be packed within the strong de Sitter radius, l sd l f 10 −13 cm. Taking the nucleon radius obtained from scattering cross sections data, r n l sP l f 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 ∼ 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 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) and (173), linking the physics of strongly interacting particles to the electroweak scale, gravity and dark energy.

5.3
The Hagedorn temperature for minimum-mass particles, the expanding Universe, and deconfinement As shown in [75], M 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 λ 3 C (m ) l 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 analog 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 ∼h/(M 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 pairproduction 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 cell new fundamental 'cells' of space are also produced, with total volume V = n cell V cell V . In [52], 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 R cell (l 2 Pl l dS ) 1/3 3.500 × 10 −13 cm r e = e 2 /(m e c 2 ) = 2.818 × 10 −13 cm.
In addition, it is clear that the fundamental field giving 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 behavior 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 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 eff ≤ 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 ρ = −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 colored objects is localized due to the 'curvature' of the strong metric field, which creates a kind of 'bag' [21]. By interpreting the energy density of the strong gravity cosmological constant as the bag constant of QCD [92], it follows that strong gravity imposes the following classical lower bound on the minimum-mass/radius ratio of a hadron: This is related to the bag constant B, via B ρ f f c 2 /(8π G f ) (135). By assuming a hadronic radius of the order of R = 0.8 fm (comparable to the proton radius [93]), and taking the estimate f 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 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 ≥ 3.37 × 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 D, with P r P ⊥ , we cannot exclude (in principle) the situation in which these conditions do not hold. If, for example γ r < γ (r ), ∀r > r , then, for a hadron in strong gravity with monotonically decreasing density, the condition d 2 dξ 2 > 0, ∀r, holds in place of Eq. (52). This situation corresponds to a tangential pressure-dominated hadronic structure, with the tangential pressure 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. (182) with the Compton wavelength λ C . The mass scale thus obtained is roughly comparable to the mass of the neutron, m n 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 → G f and → 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 [52][53][54]. The formal similarity between general relativity and the strong gravity theory again allowed us to replace G → G f and → f , yielding analogous MLURs for strongly interacting particles. Identifying R ( x total ) min , and evaluating the strong gravity MLUR Eq. (160) for Q = ±2e/3, we obtained an estimate of the mass of the up quark, m u 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 [75]. 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 pairproduction 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 [75], we obtained the resulting 'Hagedorn temperature' (T ) of the dark energy field, which was found to be comparable to the present day temperature of the CMB (181). 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 U 1/ √ and M , or, equivalently, T CMB T . 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 color 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 color charge interactions.