Quantum Haplodynamics, Dark Matter and Dark Energy

In quantum haplodynamics (QHD) the weak bosons, quarks and leptons are bound states of fundamental constituents, denoted as haplons. The confinement scale of the associated gauge group SU(2)_h is of the order of $\Lambda_h\simeq 0.3$ TeV. One scalar state has zero haplon number and is the resonance observed at the LHC. In addition, there exist new bound states of haplons with no counterpart in the SM, having a mass of the order of 0.5 TeV up to a few TeV. In particular, a neutral scalar state with haplon number 4 is stable and can provide the dark matter in the universe. The QHD, QCD and QED couplings can unify at the Planck scale. If this scale changes slowly with cosmic time, all of the fundamental couplings, the masses of the nucleons and of the DM particles, including the cosmological term (or vacuum energy density), will evolve with time. This could explain the dark energy of the universe.


Introduction
The standard model (SM) of strong and electroweak interactions may not be the final theory of the universe. The dark matter (DM) and dark energy (DE) are fundamental problems awaiting for an explanation [1]. In this Letter we consider the possible impact of the chiral gauge theory QHD (quantum haplodynamics [2]) on these problems. QHD is a theory of bound states for all the SM particles and the dark matter particles.
If the QCD and QHD couplings take definite values around the Planck scale, and this scale is permitted to slowly change with time, the QCD and QHD couplings should change also with time. Hints that the electromagnetic fine structure constant α em might change with the cosmic evolution are reported in the literature [3]. If α em changes in time, we expect that all fundamental coupling constants change in time, including the gravity constant [4,5]. Since the gravity constant G N determines the Planck mass M P = G −1/2 N , one expects that M P depends also on time and slowly evolves with the cosmic expansion.
The present framework obviously implies a link between gravity and particle physics. We shall come back to the cosmological implications after describing the essentials of QHD.

QHD
In QHD all of the SM particles (except the photon and the gluons) are bound states of the fundamental constituents called haplons, h, and their antiparticles. A first model of this type was introduced in 1981 (see [6,7,8]). Here we extend it and assume that the QHD chiral gauge group is the unitary left-right group SU (2) L × SU (2) R , which we will denote SU (2) h for short. All species of haplons h are SU (2) doublets, hence each one has two internal states h i represented by the SU (2) h quantum number i = 1, 2 . Rotations among these states are performed by the exchange of two sets of massless SU (2) h gauge bosons X r L,R µ (r = 1, 2, 3) for each chirality.
There are six haplon flavors, two of them are electrically charged chiral spinors (χ = α, β) and four are charged scalars S. One scalar (ℓ) has electric charge (+1/2) and carries leptonic flavor. The other three scalars have charge (-1/6) and carry color: c k = R, G, B ("red, green, blue"). In Table I we indicate the relevant quantum numbers.
The complete QHD gauge group is SU (3) c × SU (2) h × U (1) em , with the coupling constants (g s , g h , e). There is a gauge coupling g h = g L h , g R h for each chiral factor in SU (2) h . The QHD part of the interaction Lagrangian involving a generic chiral haplon χ (χ L or χ R ) and a scalar haplon S reads where D µ ij = δ ij ∂ µ + ig L,R h X r L,R µ (σ r /2) ij is the SU (2) h covariant derivative and σ r are the Pauli matrices. The other covariant derivative is coincident with the previous one if S is any scalar haplon.
For colored scalar haplons c k , however, we have additional terms where the covariant derivative is D µ ij = δ ij ∂ µ + ig s A µ a (λ a /2) ij and involves the gluons A µ a and Gell-Mann matrices λ a . Only the scalar haplons c k interact with the gluons, since these are the only colored constituents. Besides the fine structure constant α em = e 2 /4π there is the strong coupling α s = g 2 s /4π and its SU (2) h counterpart α h = g 2 h /4π, which is also strong at energies µ Λ h , but asymptotically free well above it. Here Λ h denotes generically any of the two confining scales Λ L h and Λ R h associated to the chiral factors of SU (2) h .
From the various haplon flavors the bound states of QHD can be constructed. Only for energies µ well above Λ h these states break down into the fundamental haplons. The weak gauge bosons are s-wave bound states of left-handed haplons α or β and their antiparticles: W + =βα, W − =ᾱβ and W 3 = ᾱα −ββ / √ 2. The neutral weak boson mixes with the photon (similar to the mixing between the photon and the neutral ρ-meson. One obtains the physical Z-boson with a mass slightly heavier than the W -boson: where the mixing parameter m = e F W /M W is related to the W decay constant F W [2]. The confinement scale Λ L h for SU (2) L defines the Fermi scale G Owing to the SU (2) L ×SU (2) R chiral structure of QHD, besides the three observed weak bosons there are also vector bosons that couple to the righthanded leptons and quarks. We assume that the righthanded confining scale Λ R h is higher so that the masses of the new vector bosons lie well above 1 TeV. They might be observed in the new experiments at the LHC. The observed scalar resonance is a p-wave excitation of the Z. Here we will not discuss these aspects of QHD in detail [2].
The leptons and quarks are themselves bound states. They are composed of a chiral haplon (α or β) and a scalar haplon: ℓ for leptons and c k for quarks. The electron and its neutrino have the structure ν = (αl) and e − = (βl), which is consistent with the quantum numbers of Table I. Similarly, the up and down quarks (with c k color) are given by: u = (αc k ) and d = (βc k ).
In QHD the first generation of leptons and quarks describe the ground states of the fermionscalar bound states; the second and third generation must be dynamical excitations. Likewise the c-quark and t-quark families are excitations of the first quark generation. Compared to the QHD mass scale the masses of the observed leptons and quarks are essentially zero.
The outcome is an effective theory equivalent to the electroweak SM in good approximation.
However, new matter content is predicted. In particular, the simplest neutral bound state of the four scalars with haplon number H = 4 is a stable color singlet spinless boson: D = (lRGB). It is stable due to haplon number conservation, which is similar to the conservation of baryon number.

A New Dark Matter Candidate
The mass of the D boson is expected to be in the region of a few TeV. It can be produced together with its antiparticle by the LHC-accelerator, and it can be observed by the large missing energy. We interpret it as the particle providing the DM in the universe. The properties of this DM particle are similar to a "Weakly Interacting Massive Particle" (WIMP), but it can be much more elusive concerning the interactions with nuclei. After the Big Bang the universe is not only filled with a gas of quarks and antiquarks, but also with a relativistic gas of D-bosons and the corresponding antiparticles, which annihilate into other particles. Due to the CP-violation there is an asymmetry in the number of D-bosons and anti-D-bosons. After freezout, at temperatures roughly ∼ 1/20 of the D-boson mass, the D particles abandon the equilibrium. A relic density remains today, a gas of D bosons, forming the dark matter. If we use the average density of matter in our galaxy, we find that there should be O(100) D-particles per cubic meter.
In contrast to more conventional WIMP's, the QHD candidate for DM can escape more easily the recent, highly restrictive, bounds obtained from scattering of DM particles off nuclei [9]. It is not difficult to estimate the cross section for the D-boson off a nucleon N (of mass m N ). It should be of order where G In QCD the proton mass is given by a few times the value of the confining scale Λ QCD , with only a tiny contribution from the quark masses. Similarly, in QHD the masses of the bound states are proportional to Λ h , although here the spectroscopy is richer and the proportionality factor is bigger for the more deeply bound states. Whereas for a weak boson such factor is of order one, in the case of the DM candidate D = (lRGB) it can be much larger. This way we can naturally explain the small cross-section of the DM particles with ordinary matter. For B D ≃ 5− 10 TeV the scattering cross-section of D-bosons off nucleons, Eq. (3), can be reduced to the level of 10 −45 cm 2 . It is compatible with the current bounds [9].

Unification at the Planck Scale
The QHD, QCD and QED couplings might unify at the Planck scale. We can check it at one-loop level, starting from their low-energy values and using the renormalization group equations (RGE's) to compute the running of these parameters. For SU (N ) groups (N > 1) one has [10]: Here we have α i = α h , α s (n f and n s are the number of fermion flavors and scalars). For the U (1) coupling α em we have a similar formula as (4), but in this case Here N h = 2 for SU (2) h . The electric charges Q are defined in Table I. For energies below Λ h we have to replace N h in b 1 with N c = 3 (or 1) and use the electric charges of the quarks (leptons) rather than those of the haplons. For the fine structure constant α em we extrapolate its value from low energies to the Planck scale M P ≃ 1.22 × 10 19 GeV. At the mass of the Z-boson we have α −1 em (M Z ) = 127.94 ± 0.014. From the mass scale of the Z-boson, µ = M Z , until a scale well above Λ h , say µ ∼ 2 TeV, we use the RGE, taking into account the charges of the three charged leptons and of the five quarks, not including the top-quark: α −1 em (2 TeV) = 123.57 . From ∼ 2 TeV up to the Planck mass we take into account the electric charges of the two spin-1/2 haplons and of the four scalar haplons (cf. Table  I). The result is: equivalently: α em (M P ) = 0.008727 (cf . Table II). A similar procedure can be followed to compute the QCD coupling constant at various energies. The accurate measurement of this constant at the Z pole yields: α s (M Z ) = 0.1184 ± 0.0007. At the Fermi scale Λ h ∼ 0.3 GeV we find α s (Λ h ) = 0.1010. Well above 1 TeV up to the Planck scale the renormalization proceeds via haplon-pairs (as indicated in Table II).
For the SU (2) h group, we focus here on the lefthanded sector and assume Λ h ≃ 0.3 TeV. We find e.g. α h (2TeV) = 0.62, and eventually at the Planck energy: α h (M P ) ≃ 0.030. From Table  II we see that the three couplings approach each other at the Planck scale. The details of the unification will depend on the particular GUT group and can be affected by Clebsch -Gordan coefficients of O(1).
We note that SU (3) × SU (2) L × SU (2) R × U (1) is a natural breakdown step for GUT groups such as e.g. SO (10). In our case we do not have spontaneous symmetry breaking (SSB), the breaking is always meant to be dynamical. The complete QHD group can thus be naturally linked to the GUT framework without generating unconfined vacuum energy.
If the three couplings come close at the Planck scale, interesting consequences can be derived in connection to the time variation of the fundamental constants, of which hints in the literature appear quite often [11]. An exact unification is not essential -we only need that the three couplings take fixed values at or around M P .

Time evolution of Fundamental "Constants"
A cosmic time change of Newton's constant G N (and hence of M P ) is conceivable in the same way as one admits a possible time change of α em [3,11]. If the QED, QCD as well as the QHD coupling constants emerge at the Planck epoch, their primeval values should be very close and not be time-dependent. Since the Planck energy changes in time, there must be a time evolution of the gauge couplings at lower energies, say around the confining scale of the weak bosons, Λ h ∼ 300 GeV. At the same time the masses of all the particles (including of course the baryons and the D-bosons) are forced to slowly evolve with the cosmic expansion since their binding energies are functions of the coupling strengths.
We can estimate the time change of G N in QHD. We use the approximate time variation of α em suggested in a typical measurement where the current value of the QED coupling is compared with that of a quasar some 12 billion years ago [3]: ∆α em /α em = (−0.54 ± 0.12) × 10 −5 .
From the RGE's and setting µ = M P we can obtain the time variation (indicated by a dot) of the Planck scale. Since b 1 = −14/9 in this case, we finḋ It follows: ∆M P /M P ≃ 0.0027 or ∆G/G ≃ −0.0054. The time variation of the non-Abelian gauge couplings α i (i.e. α s and α h ) at an arbitrary scale µ below M P is also determined: where b N was defined in (4). SinceṀ P /M P is fixed from (7), the above equation enables us to compute the cosmic time variation of the QCD and QHD couplings within the last 12 billion years at any desired energy well above Λ h , e.g. at µ 1 = 2 TeV (cf. Table II): Using the definition of the corresponding confining scales Λ i (viz. Λ QCD , Λ h ) we can check from the above formulas that their cosmic time evolution [4] is renormalization group invariant and is directly tied to the cosmic evolution of M P itself: Numerically, ∆Λ i /Λ i ≃ 3 × 10 −3 for the indicated period.

QHD and Dark Energy
The unification of the QCD and QHD couplings can have a nontrivial significance for the combined framework of particle physics and General Relativity (GR). It suggests a cosmic evolution of all the masses in the universe, both of the nuclei and of the DM particles. This can be perfectly compatible with GR. In order to preserve the Bianchi identity that is satisfied by the Einstein tensor of the gravitational field equations (∇ µ G µν = 0), the time evolution of the masses can be compensated for by the time variation of one or more fundamental gravitational parameters, typically the gravitational constant G N , or the cosmological constant Λ, or both [4,5].
With the help of the Friedmann-Lemaître-Robertson-Walker (FLRW) metric one finds that the most general local conservation law preserving the Bianchi identity, corresponding to an isotropic and homogeneous dust matter fluid with density ρ m , reads: where a is the scale factor and the primes denote derivatives of the various quantities with respect to it. ρ Λ is the vacuum energy density. From Eq. (11) we can understand how a theory of bound states at low energies can lead to a general evolution of the vacuum energy density ρ Λ and Newton's coupling G N in combination with the particle masses.
The predicted time evolution of the particle masses in QHD can be parameterized as ρ m ∼ a −3(1−ν) [4], where the presence of |ν| ≪ 1 denotes a very small departure from the standard conservation law ∼ a −3 . Such a departure is not viewed here as a loss or an excess in the number of particles in a comoving volume (beyond the normal dilution law), but rather as a change in the value of their masses.
As a result of the anomalous matter conservation law ρ m ∼ a −3(1−ν) , Eq. (11) implies that a dynamical response will be generated from the parameters of the gravitational sector G N and ρ Λ . This is how the dark energy emerges: it is the change of the vacuum energy density triggered by the time evolution of all the masses in the universe. Details will be provided elsewhere.

Conclusions
The framework we have outlined here proposes a new formulation of GUT's involving gravity ab initio (in contrast to more conventional formulations). It also proposes a new candidate for the dark matter that is not in conflict with the recent, highly restrictive, bounds for the scattering of DM particles off nuclei. The neutral and stable D-boson is the particle of the dark matter.
QHD is not based on the conventional SSB mechanism and it does not lead to a large contribution to the cosmological term. The DE appears here as the tiny (but observable) dynamical change of the vacuum energy density of the expanding background, and hence is a part of the generic response of GR to the cosmic time variation of the masses of all the stable baryons and dark matter particles in the universe.
These ideas can be tested by future astrophysical and laboratory tests in quantum optics, which are expected to detect potential proton mass variations of order 10 −14 yr −1 [12], hence at the level of the expected running of the cosmological parameters.