Mirror QCD and Cosmological Constant

An analog of Quantum Chromo Dynamics (QCD) sector known as mirror QCD (mQCD) can affect the cosmological evolution due to a non-trivial contribution to the Cosmological Constant analogous to that induced by the ground state in non-perturbative QCD. In this work, we explore a plausible hypothesis for trace anomalies cancellation between the usual QCD and mQCD. Such an anomaly cancellation between the two gauge theories, if exists in Nature, would lead to a suppression or even elimination of their contributions to the Cosmological Constant. The trace anomaly compensation condition and the form of the non-perturbative mQCD coupling constant in the infrared limit have been proposed by analysing a partial non-perturbative solution of the Einstein--Yang-Mills equations of motion.


I. INTRODUCTION
The ground state of Yang-Mills (YM) theories plays a critical role in both Particle Physics and Cosmology. In particular, the gluon condensate in Quantum Chromo Dynamics (QCD) largely determines non-trivial properties of the topological QCD vacuum and is responsible e.g. for the color confinement effects and hadron mass generation which can be understood beyond the Perturbation Theory (PT) only (for a comprehensive review on the QCD vacuum, see e.g. Refs. [1][2][3][4] and references therein). The gluon condensate directly influences properties of the quark-gluon plasma and its hadronisation, as well as dynamics of the QCD phase transition. On the other hand, YM condensates have various implications in the cosmological evolution ranging from the Cosmic Inflation [5][6][7] to the phenomenon of late-time acceleration and the Dark Energy (DE) [8][9][10] (see also Refs. [11][12][13][14][15][16]).
Currently, the Cosmological Constant (CC) with the vacuum equation of state w ≡ p/ǫ = −1 is a preferred scenario for the late-time acceleration epoch supported by a wealth of recent observations provided that w = −1.006 ± 0.045 (see e.g. Refs. [17,18]). Despite of many DE/CC models existing in the literature, there is not a compelling resolution of the CC problem i.e. why the CC term is small and positive as well as why the CC term is non-zeroth and exists at all. From the Quantum Field Theory (QFT) viewpoint, the ground state energy density of the universe should account for a bulk of various contributions from existing quantum fields at energy scales ranging from the Quantum Gravity (Planck) scale, M PL ≃ 1.2 · 10 19 GeV, down to the QCD confinement scale, Λ QCD ≃ 0.1 GeV. Even such relatively well-known vacuum subsystems of the Standard Model (SM) as the Higgs and quark-gluon condensates exceed by far the observed cosmological constant which is often considered as a severe problem [19,20] (for recent reviews on this topic, see e.g. Refs. [21][22][23] and references therein). Also it is well known, that an every field in the universe forms a divergent perturbative vacuum contribution, which is usually cut off at the Planck scale. The cancellation of these contributions may need the introduction of additional bosonic and fermionic fields putting important constraints for the particle spectrum [24][25][26].
In this work, we discuss the problem connected with formation of big non-perturbative vacuum contributions on the hadronic scale after the QCD phase transition, assuming that the contributions from higher scales are already compensated. In the case of confined QCD with color SU(3) gauge symmetry, there is a rather unique (negative-valued) contribution to the ground state energy of the universe originating from the non-perturbative quantum fluctuations of the quark and gluon fields [1, 2, 27, 28], ǫ QCD < 0. Given the fact that the CC term observed in astrophysical measurements is very small (and positive-valued), one must eliminate the QCD vacuum contribution, ǫ QCD , with an unprecedented accuracy over forty decimal digits. A dynamical mechanism for such a gross cancellation of vacua terms is yet theoretically unknown although several possible scenarios have been discussed so far e.g. in Refs. [9,21,[29][30][31]. This work is devoted to making a further step in exploring this possibility in quantum YM theories with a non-trivial ground state. Clearly, in order to cancel the QCD vacuum contribution, ǫ QCD , formed during the QCD phase transition epoch, a positive contribution to the vacuum energy density should be formed at the same QCD energy scale Λ QCD . Where could such an extra contribution originate from?
Here we suggest a new scenario of compensation realized by means of a hidden (mirror) sector of particles [32] which correspond to the extra non-Abelian gauge group and whose possible interaction with the visible SM sectors is strongly suppressed.
In particular, a class of models known as Neutral Naturalness theories has been proposed in the literature [33][34][35] as a promising solution of the naturalness problem in the SM protecting the weak scale from large radiative corrections. Various phenomenological implications of such a "Mirror World" concept have been discussed e.g. in Ref. [36]. In particular, the mirror color SU(3) gauge group is typically assumed to be a symmetry describing the confined phase in full analogy with ordinary QCD revealing interesting signatures at the Large Hadron Collider due to e.g. a mixing of mirror glueballs with the Higgs boson [37].
Quite naturally, the quantum vacua contributions from the "Mirror World" should contribute to the CC on the same footing as known vacua since "mirror" particles are expected to gravitate in the same way as the usual ones. We argue that mirror QCD (mQCD) sector, if exists, should affect the cosmological expansion, in particular, via an extra non-trivial "mirror gluon" contribution into the ground state energy of the universe. Below, we will demonstrate that under certain conditions the mirror gluon condensate can contribute to the energy density of the universe with positive sign and thus may, in principle, eliminate the negative QCD vacuum effect yielding a vanishingly small CC term. Attributing the positive vacuum energy contribution to the mQCD sector non-interacting with quarks and gluons in ordinary QCD, one may therefore resolve the issue of why such a positive-valued condensate energy density does not emerge in QCD hadron physics and affects the CC-term only. Then the observed CC can, in principle, be formed as a remnant of the gluon condensate cancellation in expanding universe (e.g. due to an uncompensated quantum gravity correction to the QCD ground state energy) [9], which appears to be remarkably consistent with the observed CC and with the Zeldovich scaling relation [38]. The exact compensation of the QCD vacuum effect by means of the mirror gluon condensate is therefore the central point to the observable smallness of the CC.

II. QCD AND MIRROR QCD VACUA COMPENSATION
The condensate in QCD is formed by the contributions of gluon and quark nonperturbative quantum fluctuations Usually, the dominant gluon contribution is given by means of the trace anomaly relation in QCD which to one-loop order reads [1, 2] where b = 9 is the first (one-loop) coefficient of the negative perturbative β-function in SU(3) gluodynamics with three light flavors with Λ QCD being the QCD scale parameter (the normalization of β corresponds to Ref. [39]). The formation of the chromomagnetic gluon condensate F 2 > 0 is typically considered at characteristic momentum scales µ g inverse to the correlation length l g , i.e. µ g ∼ l −1 g ≃ 1.2 GeV [1, 2], where the perturbative QCD still provides a realistic estimate. This validates the use of one-loop approximated expression (2.3).
One would like to explain why such a big negative contribution (2.1), which is responsible for a variety of well-known phenomena in hadron physics and is rather unique for QCD, does not affect the cosmological expansion at late times. Provided that the observed CC-term density is tiny compared to the QCD vacuum density (2. 3), the latter should be almost totally eliminated during the QCD phase transition epoch. Which mechanism could be responsible for that? A mirror copy of QCD may generate a similar gluon contribution to the trace anomaly proportional to the corresponding β-function in mQCD, In mQCD framework mirror quarks can be much heavier than in ordinary QCD [35]. Applying the idea, that mQCD is similar in main features to usual QCD, this means, that in mQCD the vacuum is formed only by mirror gluon contribution with pure gluonic βfunction, as long as the heavy quark condensates [40] are compensated by quark part of β-function, i.e.
ǫ mQCD = ǫ mQCD gluon . (2.7) A possible cancellation of QCD and mQCD vacuum may ensure a required smallness of the observable CC density which means that the corresponding mirror gluon condensate should provide a positive contribution to the vacuum density, i.e. ǫ mQCD > 0. We suppose, that mQCD gluon condensate can compensate both gluon and quark condensates of usual QCD.
Adopting the traditional hypothesis that the mQCD sector of mirror quarks and gluons is confined but is not (or very weakly) interacting with the observed SM sectors [37] and considering only chromomagnetic condensates, the compensation condition (2.8) can be satisfied if and only if the mQCD β-function is positive, i.e.
which is not realized in the perturbative mQCD regime due to Eq. (2.4). It is, however, possible to achieve the positivity of the non-perturbative β-function provided that at the characteristic scale of the QCD gluon condensate formation, µ g , the mQCD sector is in deeply non-perturbative regime. The latter condition can be satisfied if the mQCD scale parameter is large, i.e.
such that the mirror gluon condensate would be in deeply non-perturbative regime by the moment in the cosmological evolution when its density gets precisely cancelled with the QCD contribution. Note that the compensation conditions for the QCD and mQCD contributions (2.8) and (2.9), if indeed realized in nature, may be one of the most important implications of the mirror QCD in Cosmology yielding a vanishing CC-term and thus providing a dynamical way to resolve the CC problem.
Can the sign of the β-function in mQCD become positive in the non-perturbative regime? In order to answer this question, one has to employ a proper formalism which extends the effective action approach in a gauge theory beyond the perturbativity domain. Indeed, the compensation (2.8) emerges as a long-distance phenomenon and thus should hold beyond the PT.

(3.2)
Depending on the sign of the invariant J, one distinguishes the chromoelectric J > 0 and chromomagnetic J < 0 YM fields. The effective YM equation of motion in a non-trivial background metric reads and can be employed beyond the PT as long as the non-perturbative β-function is known.
Remarkably enough, this equation has a simple manifestly non-perturbative ground-state solution with positive β-function which is a complete analog of similar solution (see Ref. [42], Eqs. (13) and (17)) found in the Euclidian case with the negative coupling and the β-function corresponding to the ferromagnetic vacuum. It is also a non-perturbative analog of the perturbative solution [5,29] eliminating the traceless part of the energy-momentum tensor which then takes the following form It is important to point out, that the equations (3.4-3.5) were obtained in the pure YM case when the interaction with other fields can be neglected. In particular, we neglect the mQCD quark current in the right hand side of Eq. (3.3). Let us also stress that we consider the effective YM Lagrangian (3.1) and energy-momentum tensor (3.5) as a classical model [42] which possesses well-known properties of the full quantum theory such as (i) local gauge invariance, (ii) RG evolution and asymptotic freedom, (iii) correct quantum vacuum configurations, and (iv) trace anomaly given by the last term in Eq. (3.5). These provide a sufficient motivation and physics interest in cosmological aspects of the considering effective model.
As the solution (3.4) leads to the energy-momentum tensor of vacuum type (3.5), it immediately follows from the Friedmann equations that the corresponding energy density is constant thus, the contribution of the YM fields has a cosmological constant form. In particular, this can be realised if J = const. Indeed, the solution (3.4) fixes the invariant J to its constant initial value Such solutions were also considered in Refs. [5,29,42] in connection with the spontaneous vacuum magnetisation and in the domain concept of the QCD vacuum [43] (see also the recent paper [44] and references therein). Further, we will apply the solution (3.4) and (3.7) to the the mQCD theory in the non-perturbative regime. The energy-momentum tensor in this case becomes constant as expected whereḡ 2 0 =ḡ 2 (J 0 ). The couplingḡ 2 (J) touches the linear function f (J) =ḡ 2 0 · (J/J 0 ) at the point J 0 (indeed, it has the same value and derivative). And vice versa, ifḡ 2 (J) touches const · J at some point J 0 , then dḡ 2 dJ | J=J 0 =ḡ 2 0 J 0 , which means that Eq. (3.4) is satisfied at J 0 . Indeed, the existence of such a contact point is a necessary and sufficient condition for the solution (3.4) with fixed J = J 0 . This allows us to constrain generic non-perturbative behavior of the correspondingḡ 2 (J). An illustration of the corresponding infrared behavior of the mQCD coupling in consistency with both the non-perturbative asymptotics for the β-function (3.4) and the conventional perturbative regime of asymptotic freedom (2.4) is shown in Fig. 1. A desirable non-monotonic shape of the coupling was earlier discussed in the case of usual QCD [45,46]. From Eq. (3.6) one notices that the mQCD gauge field give a constant vacuum contribution to the energy-momentum tensor in the Einstein equations. Since the QCD and mQCD contributions to the ground-state energy density have opposite signs there is a compelling possibility that they can, in principle, cancel each other at some moment in the cosmological history provided that ǫ mQCD → −ǫ QCD (3.9) in the infrared regime of mQCD. Keepingḡ 2 0 > 0 in both QCD and mQCD, we arrive at the following form of the compensation condition (2.8) which means, that mQCD condensate has to be chromomagnetic. After such a compensation is achieved, only a very small ǫ CC contribution, which could be formed by other vacuum sources ǫ vac and possibly by a non-compensated part of mQCD and QCD vacua contributions ǫ QCD + ǫ mQCD , remains. Of course, a fine-tuning is unavoidable to match the observations, although would not be entirely unreasonable due to the same order magnitude of the QCD and mQCD contributions. Such a vacua alignment, if realised in Nature, can be suggested e.g. by their common Anthropic origins [47]. Then the standard Friedmann equation in the non-stationary FLRW universe 3 κ (a ′ ) 2 a 4 = ǫ mat + ǫ CC , ǫ CC ≡ ǫ QCD + ǫ mQCD + ǫ vac , (3.11) determines the cosmological evolution, a = a(η), driven by the gluon and mirror gluon condensate densities (compensating each other exactly or in part), the matter contribution, ǫ mat , and other possible vacua contributions of a different kind, ǫ vac .
Let us notice for completeness, that Eq. (3.4) allows for a more specific solution apart from constant fields (3.7). This another solution appears, if the non-perturbative β-function becomes constant satisfying Eq. (3.4) in some finite domain which corresponds to the strong coupling regime (saturated behaviour). In this case, Eq. (3.4) can be substituted into the RG equation (ξ is a numerical parameter) which implies that the mQCD coupling is proportional to J in the infrared limit, e.g.ḡ Such behaviour clearly guarantees a constant vacuum energy (see Eq. (3.5)) as well as a possibility for the QCD/mQCD vacua compensation. According to Eq. (3.13) the gauge coupling squared has to approach the linearḡ 2 ∼ J asymptotics in the non-perturbative regime (note, for the constant field solution (3.7)ḡ 2 has to satisfy a much less restrictive constraint and just touches the linear asymptoticsḡ 2 ∼ J at a fixed point). So, within the QCD/mQCD vacua compensation hypothesis, both vacua subsystems should be generated in early universe at close (but different) energy scales and then get compensated during the cosmological QCD phase transition epoch. As was shown above, this can be realised in a deeply non-perturbative regime for the mirror gluon condensate which asymptotically acquires the same absolute value of energy density and opposite sign compared to the QCD gluon one, such that they almost exactly eliminate each other at macroscopically large space-time separations.

IV. SUMMARY
In this Letter, we have considered cosmological consequences of a possible exact cancellation of vacuum energies between QCD and a mirror high-scale copy of QCD in the confined regime. One of the most important implications of such a cancellation is a possibility for the dynamical elimination of the QCD vacuum density contribution to the cosmological constant resulting in its observed smallness (provided that all perturbative vacua are eliminated by some other mechanism, see e.g. Ref. [24][25][26]). By an appropriate fine-tuning of QCD and mirror QCD vacua parameters, the compensation is provided by a partial non-perturbative solution of the Yang-Mills equation of motion corresponding to a positive constant β-function in deeply infrared regime of mirror QCD. Such a dynamical vacua compensation due to different signs of β-functions in QCD and mirror QCD can be yet another implication of the Anthropic Principle [47] realised for QCD and mirror gluon condensates.