Baryogenesis in Lorentz-violating gravity theories

Lorentz-violating theories of gravity typically contain constrained vector fields. We show that the lowest-order coupling of such vectors to $\mathrm{U}(1)$-symmetric scalars can naturally give rise to baryogenesis in a manner akin to the Affleck-Dine mechanism. We calculate the cosmology of this new mechanism, demonstrating that a net $B-L$ can be generated in the early Universe, and that the resulting baryon-to-photon ratio matches that which is presently observed. We discuss constraints on the model using solar system and astrophysical tests of Lorentz violation in the gravity sector. Generic Lorentz-violating theories can give rise to the observed matter-antimatter asymmetry without violating any current bounds.


I. INTRODUCTION
Why is there so much more matter than antimatter? One most likely cannot appeal to initial conditions, as these would be washed away by inflation. The standard model can provide such an asymmetry during the electroweak phase transition, but cannot produce enough to accommodate observations [1]. It seems probable, then, that a dynamical generation mechanism, or baryogenesis, arises from new physics beyond the standard model. 1 In this paper we point out that Lorentz violation might play a key role in this new physics. While Lorentz invariance is extraordinarily well tested in the matter sector, the possibility of gravitational Lorentz violation remains relatively unconstrained. If boosts are broken but rotational invariance is maintained-i.e., if gravity picks out a preferred rest frame-then the low-energy physics is described by Einstein-aether theory, a vector-tensor theory in which the vector field is constrained to have a fixed, timelike norm. This is the general effective field theory when boosts are broken [8]; for example, a special case of Einstein-aether arises in the low-energy limit of Hořava-Lifschitz gravity [9,10], a putative UV completion of general relativity which relies on the existence of a preferred foliation.
We demonstrate that if a U(1) B−L scalar couples to the vector of Einstein-aether theory, then the lowest-order interactions between the two can lead to baryogenesis. This operates in a manner qualitatively similar to Affleck-Dine baryogenesis [11] or the recent model of Ref. [12], in which the U(1) B−L symmetry is broken at early times due to a tachyonic mass proportional the Hubble parameter H appearing in the effective scalar potential. 2 In purely metric theories this is difficult to achieve, as H is not a spacetime scalar. Breaking boosts cures this difficulty, and indeed * sakstein@physics.upenn.edu † adamsol@physics.upenn.edu 1 For reviews of baryogenesis, see, e.g., Refs. [1][2][3][4][5][6][7]. 2 In the former model, this arises due to a coupling of the scalar to the inflaton, while in the latter, it arises from a Weyl coupling to dark matter.
in Einstein-aether H is simply proportional to the divergence of the timelike vector field. As Einstein-aether is the most general low-energy effective theory for broken boosts, our conclusion can be stated as follows: if the Universe contains a U(1) B−L scalar with softly broken symmetry and spontaneous Lorentz violation, a working baryogenesis mechanism comes for free.
This paper is organized as follows. In section II we introduce the model of scalar-aether baryogenesis and discuss known constraints on the theory. In section III we derive the cosmology and verify that this model can yield the observed baryon-to-photon ratio with sensible parameters, and we conclude in section IV.

II. MODEL
The model we will consider is the constrained vector (or "aether") u µ of Einstein-aether theory coupled to a new U(1) B−L scalar φ. At leading order, the most general action we can write down is where is the most general kinetic term for the Lorentz-violating vector u µ , and κ is a Lagrange multiplier that ensures that the vector is timelike and of fixed norm, u µ u µ = −m 2 . We have included some U(1) B−L -violating terms proportional to ε in order to generate a net B − L. The last line is the leading-order interaction one can write down between φ and u µ given the symmetries. 3 This is the general low-energy effective theory with φ, broken boosts, and softly-broken U(1) B−L . The coupling between u µ and φ can give rise to baryogenesis using a mechanism akin to (but distinct from) that of Affleck and Dine [11] or similar generalizations [12]. In particular, the effective potential for φ is so that ∇ µ u µ acts as a tachyonic mass term for φ. We can see this explicitly by considering a homogeneous and isotropic cosmological setting, so that the Universe is described by a Freidmann-Lemaître-Robertson-Walker metric (in cosmic time) The on-shell condition u µ u µ = −m 2 and symmetry imply that the vector must be of the form [19][20][21] u µ = (m, 0, 0, 0) .
In this case, the divergence of u µ is ∇ µ u µ = m∂ t ln √ −g = 3mH, so that the effective potential becomes (6) This potential leads to baryogenesis similarly to the well-known Affleck-Dine mechanism. In the early Universe, the U(1) B−L symmetry is broken as the tachyonic mass term is more important than the bare mass, while at late times the symmetry is restored. During the broken-symmetry phase, the motion of the angular component of the scalar generates a net B − L due to the symmetry-breaking terms εφ 4 + ε † φ † 4 . When the symmetry is restored, the B − L is stored in the field, and can be transferred to the standard model through sphaleron processes [22], although one must first transfer the asymmetry to left-handed standard model particles. The details of the transfer were discussed for models such as ours in Ref. [12], where the neutrino portal [23,24] was identified as one promising mechanism. We note that our model differs quantitatively from Affleck-Dine: the tachyonic mass scales like mH rather than H 2 , which can lead to novel and interesting new features. In what follows, we will calculate the cosmology of this model, paying special attention to the generation of a net B − L.
In this subsection, we briefly summarize observational and theoretical constraints on the parameters in our model. Most of these will apply to Einstein-aether theory or to its coupling to a real scalar. We will use the notation Experimental constraints on Einstein-aether theory tend to place upper bounds on the aether vacuum expectation value (VEV) m, with the result c i m 2 ≪ M 2 Pl for generic values of the c i parameters. We note that any of these constraints can be weakened or removed entirely by tuning the c i parameters, although these tunings cannot all be done simultaneously. We refer the reader to Sec. V.D of Ref. [16] for a more comprehensive summary of constraints on the aether.
The strongest constraints come from gravitational Cerenkov radiation: high-energy cosmic rays could lose energy to subluminal aether-graviton modes, leading to a degradation in cosmic ray propagation which has not been observed, constraining m/M Pl < 3 × 10 −8 [25]. These constraints can be avoided by tuning the c i or allowing for superluminal propagation in aether-graviton modes; since this is an explicitly Lorentz-breaking theory, superluminality may not be as deadly as one normally expects. The preferred-frame parameters α 1,2 in the parametrized post-Newtonian formalism are modified by the aether, constraining m/M Pl < 6 × 10 −4 in the absence of tuning c i [26,27]. Note that the tuning which eliminates gravitationalCerenkov radiation (c 3 = −c 1 , c 2 = c 1 /(1 − 2c 1 )) also sets α 1 = 0, so that the dominant constraint comes from α 2 , in which case the strongest constraint on m is rather mild, m/M Pl < ∼ 10 −2 .
Under the assumption that c i m 2 ≪ M 2 Pl , we are justified in ignoring the mixing with gravity [20], in which case there are a few constraints on the c i from flat space perturbation theory. In the vector sector, the absence of ghosts requires c 1 > 0 [20]; coupling a scalar to ∇ µ u µ , as in this paper, does not modify the vector perturbations around flat space [16]. The no-ghost condition for the spin-0 piece of u µ is the same, while gradient stability requires c 123 > 0. Some authors require the spin-0 aether mode to propagate subluminally, which would imply c 123 < c 1 [20], although the scalar coupling relaxes this bound to c 123 < c 1 + x for some x > 0 [16]. If we require the sound speed of tensors to be subluminal then we would require c 13 > 0 [20].
The aether-scalar coupling can lead to a gradient instability, placing an upper bound on α [16]. Writing φ(x) = 1 √ 2 ρ(x)e iθ(x) , the real scalar ρ interacts with the aether through a potential V (θ, ρ) = 1 2 m 2 φ ρ 2 + 1 4 λρ 4 − 1 6 αΘρ 2 , where Θ ≡ ∇ µ u µ . Gradient stability around flat space where V Θρ = ∂ Θ ∂ ρ V and V ρρ = ∂ 2 ρ V are evaluated at the background values of Θ and ρ. Applying this to our potential, we find that the constraint is trivially satisfied in the unbroken symmetry phase (ρ = 0), while in the broken symmetry phase (ρ =ρ) we have, in the k → 0 limit, 5 where we have assumed m φ ≪ρ, as we will throughout this paper. This constraint places a mild upper bound on the coupling, α < ∼ √ λ. We expect λ < ∼ O(1), otherwise the theory is strongly coupled. Only the combination αm is relevant for baryogenesis, and this constraint then implies that we cannot simultaneously take m < M Pl to satisfy the constraints above whilst having αm parametrically larger. 4 While the analysis of Ref. [16], in contrast to our model, assumed a single real scalar, the field θ decouples from ρ and u µ around the background ∂µθ = 0. In principle a non-zero ∂µθ = cµ could modify the constraint by shifting the effective mass for ρ fluctuations, m 2 φ → m 2 φ + c 2 . 5 Strictly speaking the k → 0 limit is not physical, as we will be dealing with cosmological spacetimes, which only resemble flat space for k ≫ H. A more exact condition is obtained by setting k =H, whereH is the value of the Hubble parameter below which the symmetry is restored, as discussed in the next section. This only significantly relaxes the constraint (8) if m φ ≫ αm, in which case the constraint becomes α 2 ≤ √ 18c 123 λm φ /m.

III. COSMOLOGY
In this section we will assume a homogeneous and isotopic cosmological background and derive a simple estimate (19) for the baryon-to-photon ratio generated by our model. In order to verify that the approximations we make are valid, we also solve the equations of motion numerically, with the results plotted in figures 1 and 2.
We can see from the potential (6) that the U(1) symmetry is broken at early times and restored at late times, when When H >H there is a time-dependent symmetrybreaking minimum at where we have assumed that the small symmetryviolating terms (εφ 4 + ε † φ † 4 ) are negligible, or, equivalently, have chosen the coefficients ε so that this is the case at early times. The field tracks this minimum nearly adiabatically until H ≈H, at which point the symmetry is restored and the field begins to oscillate around the symmetry-restored minimum at φ = 0. When this occurs, the U(1) B−Lviolating terms play an important role. We would like these to become important around the time that the symmetry is restored. Expanding ε =εe iψ , the correction to the potential is where φ = |φ|e iθ . The B − L charge density is so we see that the motion of the angular field is responsible for generating a net B − L. This means that θ should not sit at its minimum in the early Universe, and, indeed, one expects it to be frozen at some initial value due to Hubble damping. We would like it to begin rolling around the time of symmetry restoration in order to generate a net B −L before the field settles into the new minimum at φ = 0, which will be the case if the canonicallynormalized field's mass 6 m 2 θ = ∆V eff θθ /|φ| 2 is of order H 2 . If m θ <H the angular field will not roll after symmetry restoration and no B − L will be generated. Similarly, if m θ >H the field starts rolling long before symmetry restoration, and the value of B − L is set by tuning the initial conditions. Setting m θ ∼H implies that where we have used |φ| = |φ min | at the time of symmetry restoration.
where we have omitted factors of order unity. Using equations (9), (10), and (13) we can estimate the B − L conserved charge density as We do not directly observe n B−L , but rather the baryon-to-photon ratio n b = n B−L /s, where s is the entropy density. This introduces some model dependence; for concreteness, and to minimize the number of free parameters, we will focus on a minimal model in which the Universe reheats instantaneously after inflation, and the U(1) B−L symmetry is restored shortly thereafter. The assumption of instantaneous reheating yields where T R is the reheat temperature. Combining this with equation (17) we find We see that this new mechanism can produce the observed baryon-to-photon ratio, n b ∼ 10 −10 , with sensible choices for the reheat temperature and model parameters. Note that the parameters α and m only appear in the combination αm, while it is m 2 (in combination with the c i ) which is constrained by experimental tests of Lorentz violation, as discussed in section II A. We expect λ < ∼ O(1), otherwise the theory is strongly coupled, and as discussed above, we should have α 2 < ∼ O(10)λ to ensure gradient stability around flat space.
In order to verify the approximations we have made above, we have numerically integrated the scalar field equations assuming a radiation-dominated Universe. In figure 1 we plot the motion of the complex scalar for two different models. One can see the behavior we predicted qualitatively above: the field tracks its time-dependent minimum at early times before the angular field begins to roll when the symmetry is restored, giving rise to a spiral trajectory. In figure 2 we plot the baryon-to-photon ratio n B−L /s for the same models. One can see that our numerical results agree well with our prediction (19).

IV. CONCLUSIONS
In this paper we have studied baryogenesis in Lorentzviolating theories of gravity, which, at low energies, are naturally described by a constrained vector so that there is a preferred frame. Baryogenesis requires a field charged under U(1) B−L , the simplest choice being a complex scalar. We have demonstrated here that the lowest-order interaction between the scalar and vector can give rise to a tachyonic mass term for the scalar proportional to the Hubble parameter so that the U(1) B−L symmetry is broken at early times. Inverse phase transitions such as these can generate a net B − L through the coherent motion of the scalar when the symmetry is restored at late times, and we have shown here that Lorentz-violating theories can successfully generate the observed baryon-to-photon ratio using this phenomenon. Furthermore, this can be achieved for parameter choices that are not ruled out by current constraints. Our theory differs from the quintessential paradigm-the Affleck-Dine mechanismin that the tachyonic mass is proportional to H rather than H 2 , which gives rise to new features and a qualitatively different cosmology, which we have calculated in detail. Lorentz-violating gravity theories continue to be important in the study of dark energy and quantum gravity. Here, we have shown that they may also shed light on the origin of the mater-antimatter asymmetry.