---
title: Gravitational Baryogenesis Mechanisms
url: https://www.emergentmind.com/topics/gravitational-baryogenesis
type: topic
---

# Gravitational Baryogenesis Mechanisms

Gravitational baryogenesis denotes a class of baryon-asymmetry mechanisms in which the expanding spacetime itself supplies the bias between baryons and antibaryons. In its canonical form, the baryon current couples derivatively to a gravitational scalar—most often the Ricci scalar—so that a time-dependent background curvature induces an effective baryon chemical potential. If baryon-number-violating reactions remain active until a decoupling temperature \(T_D\), a net baryon asymmetry is frozen into the plasma. The mechanism has been embedded in backgrounds ranging from vacuum inflation and cosmological-constant-dominated FRW cosmologies to anisotropic, braneworld, teleparallel, and higher-curvature models, but it also faces a major theoretical challenge: the standard \((\partial_\mu R)J_B^\mu\) coupling generically feeds back into the gravitational equations and can produce higher-order instabilities [1708.09137, 1712.04627].

## 1. Canonical interaction and thermodynamic implementation

The standard operator used in gravitational baryogenesis is
\[
\frac{1}{M_*^2}\int d^4x\,\sqrt{-g}\,(\partial_\mu R)J^\mu,
\]
or equivalently, in some conventions,
\[
{\cal L}_{GBG}=\frac{f}{m_0^2}(\partial_\mu R)J_B^\mu.
\]
Here \(J^\mu\) is the baryon current, \(R\) is the Ricci scalar, and \(M_*\) or \(m_0\) is an effective cutoff scale. In a homogeneous cosmological background only the temporal derivative survives, so the interaction acts as an effective chemical potential,
\[
\mu_B \sim \pm \frac{\dot R}{M_*^2},
\]
with opposite sign for baryons and antibaryons. The same logic is used in extensions where \(R\) is replaced by \(\mathcal G\), \(T\), \(T_G\), \(B\), \(\mathcal T\), \(\mathcal T^2\), or a model-dependent function \(f\) of such quantities [2006.06534, 2103.13328, 1607.00545].

In thermal equilibrium, the induced asymmetry is computed from the usual small-chemical-potential expansion. Several papers write
\[
n_B=\frac{g_b}{6}\mu_B T^2,
\]
while a more detailed expression is
\[
n_b=n_B-n_{\bar B}=\frac{g_b T^3}{6\pi^2}\left(\frac{\pi^2\mu_B}{T}+\left(\frac{\mu_B}{T}\right)^3\right),
\]
whose linear term dominates for \(\mu_B/T\ll 1\). With entropy density
\[
s=\frac{2\pi^2}{45}g_s T^3,
\]
many treatments obtain the standard estimate
\[
\frac{n_B}{s}\simeq -\frac{15g_b}{4\pi^2 g_*}\frac{\dot R}{M_*^2 T_D}\Big|_{T_D},
\]
or closely related expressions with \(g_{\star s}\) and slightly different temperature conventions. The asymmetry is therefore controlled by the geometric time derivative evaluated when baryon-violating processes decouple [1010.4966, 2103.13328, 2012.12965].

A notable feature, emphasized repeatedly in the literature, is that this mechanism can operate in thermal equilibrium. The curvature-induced bias modifies the equilibrium number densities themselves, so the usual intuition that baryogenesis strictly requires conventional out-of-equilibrium dynamics is softened, although freeze-out remains essential because it fixes the final yield when baryon-number violation becomes inefficient [2006.06534, 1712.04627].

## 2. Background dependence, freeze-out, and inflationary realizations

The efficacy of gravitational baryogenesis is highly background dependent. In standard FRW cosmology the familiar obstruction is that for a perfect radiation fluid with \(w=1/3\), one has \(R=0\) and hence \(\dot R=0\); during matter domination \(H\) is typically too small to generate a large asymmetry. This motivates embeddings in nonstandard early-universe backgrounds where curvature remains nonzero and time dependent [2006.06534, 1603.02403].

A prominent example is vacuum inflation, where inflation is driven by a quantum potential rather than an inflaton field. In that setting the matter/radiation density produced by Hawking radiation obeys
\[
\dot{\rho}=g_*\frac{3\sigma H^5}{32\pi^4}-4H\rho,
\]
with
\[
T=\left(\frac{\rho}{2\sigma g_*}\right)^{1/4},\qquad
s=\frac{8\sigma}{3}g_*T^3.
\]
The curvature is
\[
\mathcal R=-6pH^2,\qquad \dot{\mathcal R}=12(p+2)pH^3,
\]
and the observationally favored value \(p+2=0.013\) implies \(\dot{\mathcal R}\neq 0\) during inflation. Because entropy continues to increase until the end of inflation, the generated asymmetry is later diluted; the paper finds a numerical dilution factor of approximately \(0.77\), leading to
\[
\frac{n_B}{s}\approx 0.77 \frac{15g_b}{4\pi^2g_*}\frac{\dot{\mathcal R}}{M_*^2 T}\bigg|_{T_D},
\]
and, for the reduced Planck scale choice of \(M_*\), obtains
\[
\frac{n_B}{s}\approx 1.88\times 10^{-10},
\]
close to the quoted observed value \(0.864_{-0.015}^{+0.016}\times 10^{-10}\) [1708.09137].

A distinct FRW realization uses an exact solution with cosmological constant \(\Lambda\),
\[
a(t)=m\,\sinh^{\frac{2}{3(1+\omega)}}\!\left[\Lambda_0(1+\omega)t+t_0\right],
\qquad
\Lambda_0=\sqrt{\frac{3\Lambda}{4}},
\]
which interpolates from decelerated to accelerated expansion. Its Ricci derivative is
\[
\dot R =-\frac{8\Lambda_0^3}{3}(1+\omega)(3\omega-1)\,
\coth\!\left[\Lambda_0(1+\omega)t+t_0\right]\,
\mathrm{csch}^2\!\left[\Lambda_0(1+\omega)t+t_0\right].
\]
This immediately shows that \(\dot R=0\) for \(\omega=1/3\) and \(\omega=-1\), so neither pure radiation nor pure de Sitter generates asymmetry in that model. For the benchmark choice
\[
g_{\star s}=106,\qquad g_B\sim 1,\qquad T_D=2\times10^{12}\,\mathrm{GeV},\qquad M_\star=2\times10^{16}\,\mathrm{GeV},
\]
the derived \(n_B/s\) is reported to agree with the quoted observational level \(\eta_B\sim 9.0\times10^{-11}\) over an appropriate time range [2006.06534].

A third inflationary construction attributes the departure from \(R=0\) to gravitationally induced particle production. In that scenario
\[
\frac{\Gamma}{3H}=\left(\frac{H}{H_I}\right)^p,
\qquad
H=\frac{H_I}{(1+Da^{2p})^{1/p}},
\]
and the Ricci scalar becomes
\[
R=-4\frac{H^2}{H_I^p}\left(1-\frac{H^p}{H_I^p}\right)M_{\rm Pl}^2,
\]
with
\[
\dot R = 24(2+p)\frac{H^3}{H_I^p}\left(1-\frac{H^p}{H_I^p}\right).
\]
The source vanishes both at the exact de Sitter start and in the late radiation limit, so baryogenesis is localized in the transition epoch. For \(p=1\), \(M_*=M_{\rm Pl}\), and \(H_I\sim10^{-4}M_{\rm Pl}\), the paper finds \(\eta\sim 4.5\times10^{-10}\) and infers an inflationary scale \(H_I\sim10^{15}\,\mathrm{GeV}\) from the observed asymmetry [1610.01591].

## 3. The radiation-era obstruction and mechanisms that evade it

The central phenomenological question is how to maintain a nonzero source when the background is close to conformal radiation. A large fraction of the literature can be read as a sequence of attempts to evade the standard GR result \(R=\dot R=0\) for \(w=1/3\).

Several distinct strategies recur. Some change the geometric invariant in the baryogenesis operator, replacing \(R\) by \(\mathcal G\), torsional scalars, or boundary terms. Others keep the Ricci-scalar operator but modify the cosmological dynamics so that \(R\) or \(\dot R\) remains nonzero even for radiation-like matter. A further strategy uses anisotropy or particle-creation effects to introduce additional sources of curvature evolution [1607.00545, 1610.01453, 2012.12965, 1010.4966].

| Framework | Effective source | Radiation/conformal behavior |
|---|---|---|
| Standard GR GBG | \(\dot R\) | vanishes for \(w=1/3\) |
| Gauss-Bonnet GBG | \(\dot{\mathcal G}\) | generically nonzero in radiation |
| Running vacuum models | \(\dot R=-4\dot\Lambda(H)\) for \(w=1/3\) | nonzero if \(\Lambda(H)\) runs |
| DGP brane cosmology | brane-modified \(R,\dot R\) | nonzero for \(w=1/3\) |
| Anisotropic Bianchi I | shear-dependent \(R,\dot R\) | anisotropy can sustain a source |

In Bianchi I cosmology the Ricci scalar can be written as
\[
R=3H^2(1-3w)+\sigma^2(3w-1),
\]
so shear enters explicitly. The corresponding \(\dot R\) then depends on \(H\), \(\rho\), \(w\), \(\dot w\), and \(\sigma^2\). This means that anisotropy can enhance the baryon asymmetry and, in the simplest Bianchi-I analysis, can keep the source nonzero in situations where the isotropic contribution would be suppressed [1010.4966]. However, when anisotropy is generated by anisotropic inflation, the enhancement is constrained by the CMB bound
\[
1-\frac1c \lesssim 10^{-7}\times \left(\frac{|h_\ast|}{10^{-2}}\right)\left(\frac{N_k}{60}\right)^{-2},
\]
and by the gravitino bound \(T_{RD}\lesssim10^9\,\mathrm{GeV}\). The resulting conclusion is that anisotropy does enhance \(\dot R\), but generally not enough to explain the observed asymmetry without fine tuning [1603.02403].

Running vacuum models evade the conformal obstruction by promoting the vacuum energy to \(\Lambda(H)\). For \(w=\tfrac13\),
\[
\dot R=-32\pi G\,\dot\rho_\Lambda=-4\dot\Lambda(H),
\]
so gravitational baryogenesis becomes directly sensitive to the running vacuum sector. Two benchmark running-vacuum models yield
\[
\frac{n_B}{s}\sim 1.79316\times10^{-10}
\]
and
\[
\frac{n_B}{s}\sim 7.29947\times10^{-11},
\]
respectively, for representative choices of \(M_*\), \(T_D\), and the running parameters \(\nu\) or \(v\) [1610.01453]. DGP brane cosmology achieves a similar effect through braneworld corrections to the Friedmann equation. In the \(\epsilon=-1\) branch, the paper reports
\[
\frac{n_B}{s}\approx 8.38\times10^{-11}
\]
or
\[
\frac{n_B}{s}\approx 9.53\times10^{-11}
\]
for \(w=1/3\) on two decoupling branches, while standard GR would give zero [2012.12965].

Replacing \(R\) by the Gauss-Bonnet invariant \(\mathcal G\) is the most direct geometric workaround. Since
\[
\mathcal G=24H^2(\dot H+H^2),
\]
one generally has \(\dot{\mathcal G}\neq0\) even during radiation domination. In \(R+F(\mathcal G)\) gravity with
\[
F(\mathcal G)=f_0\mathcal G^\gamma,
\]
the baryon asymmetry scales as \(n_B/s\sim T_D^{5/\gamma-1}\), and for
\[
M_*=10^{12}\,\mathrm{GeV},\qquad T_D=2\times10^{16}\,\mathrm{GeV},\qquad \gamma=0.49,
\]
the bound \(n_B/s<9\times10^{-11}\) implies
\[
f_0>7.44\times10^{14}.
\]
The same paper notes that derivative couplings to \(R^{\alpha\beta}R_{\alpha\beta}\) or \(R^{\alpha\beta\gamma\delta}R_{\alpha\beta\gamma\delta}\) can also remain active in the conformal limit [1607.00545].

## 4. Modified-gravity generalizations

Beyond these background modifications, the literature generalizes gravitational baryogenesis by replacing the Ricci scalar with alternative gravitational or matter-geometry quantities, or by coupling the baryon current to \(\partial_\mu f\) rather than directly to \(\partial_\mu R\). The result is a large model space in which the freeze-out formula is structurally similar, but the source term is computed from the modified dynamics [1907.13460, 1607.08561].

In curvature-based modified gravity, a representative example is nonminimal \(f(R,T)\) theory with
\[
f(R,T)=R+\zeta RT.
\]
In that framework the paper studies three derivative couplings: \(\partial_\mu R\), \(\partial_\mu T\), and \(\partial_\mu f(R,T)\). For the scale factor \(a(t)=e^{\alpha t}t^\beta\), it finds that the \(\partial_\mu R\) and \(\partial_\mu f(R,T)\) cases can remain nonzero even in a radiation-dominated universe and can yield
\[
\frac{n_B}{s}\sim 8.05\times10^{-11}
\quad\text{and}\quad
\frac{n_B}{s}\simeq 8.03\times10^{-11}
\]
for \(\beta=\tfrac12\), \(\alpha=2\times10^{-25}\), and \(\zeta=1\), whereas the \(\partial_\mu T\) coupling is either zero or phenomenologically unacceptable in the same model [1907.13460]. In anisotropic \(f(R)\) gravity on a Bianchi-I background, the asymmetry is enhanced by anisotropy; Model I rises roughly linearly with the parameter \(n\), whereas Model II begins at a higher positive value and then decreases or stabilizes with \(n\), both using the target level \(n_b/s\approx 9.2\times10^{-11}\) as the phenomenological benchmark [2103.13328].

Teleparallel generalizations replace curvature by torsion. In \(f(T)\) gravity,
\[
T=-6H^2,
\]
and one studies either \(\partial_\mu(-T)J_B^\mu\) or \(\partial_\mu f(-T)J_B^\mu\). The TEGR limit performs poorly, with
\[
\frac{n_B}{s}\simeq -1.02\times10^{-42}
\]
for a benchmark radiation-era choice, but power-law models such as
\[
f(T)=B(-T)^n
\]
can reach
\[
\frac{n_B}{s}\simeq 7.53\times10^{-11},
\]
while a generalized \(\partial_\mu f(T)\) coupling can give
\[
\frac{n_B}{s}\simeq 9.24\times10^{-11}
\]
for suitable parameter choices [1607.08561]. Further teleparallel extensions include \(f(T,\mathcal T)\) gravity, where parameter sets reproducing
\[
9.42878\times10^{-11},\qquad 9.42029\times10^{-11},\qquad 9.42825\times10^{-11}
\]
are reported, albeit sometimes with very large or very small couplings [2305.09707], and \(f(T,T_G)\) or \(f(T,B)\) models, where generalized couplings to \(\partial_\mu(T+T_G)\), \(\partial_\mu(T+B)\), or \(\partial_\mu f\) are compared directly with the benchmark \(\eta_B/S=9.42\times10^{-11}\) [2009.13293].

Matter-geometry extensions broaden the source space further. In energy-momentum squared gravity,
\[
f(R,\mathcal T^2)=R+\eta M_{Pl}^{2-8n}(\mathcal T^2)^n,
\qquad
\mathcal T^2=T_{\mu\nu}T^{\mu\nu},
\]
the operator can involve \(\partial_\mu(\mathcal T^2)\) or \(\partial_\mu f(R,\mathcal T^2)\). The \(n=\tfrac12\) case is identified as the most promising; the \(n=1\) case has \(R=\dot R=0\) in radiation domination, but \(\partial_\mu\mathcal T^2\) can still generate asymmetry. The same paper also notes that some parameter choices imply implausibly low baryon-violating scales, such as \(M_B\sim160.96\) GeV in one \(n=1\) example, which would already have been observed [2409.04623].

Quantum-gravity-inspired and vector-tensor realizations have also been explored. In Hořava–Lifshitz gravity, modified FRW evolution leads to viable baryogenesis for quintessence-like and matter-dominated fluids, while radiation still gives zero asymmetry in the flat case [1802.00313]. In Extended Proca-Nuevo gravity, where a massive vector field couples nonminimally to curvature, three cosmological histories—power-law, exponential, and modified exponential—are analyzed, and all are reported to yield asymmetry consistent with the benchmark \(\eta_B/s\sim9.42\times10^{-11}\) for suitable ranges of the interaction scale \(\Lambda\) [2504.06133].

## 5. Instability, stabilization, and alternative formulations

The most serious criticism of standard gravitational baryogenesis is that the same derivative-curvature coupling that generates the chemical potential also modifies the gravitational field equations. Because the interaction contains derivatives of \(R\), varying the action with respect to the metric produces higher-order terms, so the usual algebraic trace relation of GR is replaced by a differential equation for the curvature scalar itself [1712.04627, 2301.08322].

In the critical analysis of the standard \((\partial_\mu R)J^\mu\) operator, the trace equation becomes fourth order. In a bosonic model, after thermal averaging and neglecting subleading terms, one obtains
\[
\frac{d^4R}{dt^4}+\mu^4R=-\frac12 T^{(tot)},
\qquad
\mu^4=\frac{m_{Pl}^2m_0^4}{8\pi q^2T^2},
\]
or, in the \(R^2\)-stabilized formulation,
\[
\frac{d^4R}{dt^4}+K^2\frac{d^2R}{dt^2}+M_R^2K^2R=0.
\]
The homogeneous solutions have exponentially growing modes, \(R\sim e^{\lambda t}\), with positive real parts for some roots. The conclusion is that the instability can develop on timescales much shorter than the Hubble time and can destroy standard cosmology unless additional physics regulates it [1712.04627].

A proposed cure is to add a Starobinsky-type term,
\[
\frac{M_{\rm Pl}^2}{16\pi}\left(R-\frac{R^2}{6M_R^2}\right),
\]
which introduces a scalaron mass \(M_R\). In that case the stability condition becomes
\[
K_{\rm eff}^2<4M_R^2.
\]
For \(M_R=3\times10^{13}\,\mathrm{GeV}\), one estimate given is
\[
\frac{M^2}{q\,T}\gtrsim 3\times10^4\,{\rm GeV}^{1/2}.
\]
This stabilizes the runaway in both bosonic and fermionic realizations, but the paper emphasizes that the resulting cosmology is still noticeably modified and the stabilized curvature amplitude remains much larger than in ordinary Friedmann evolution [2301.08322].

Two important alternative directions attempt to sidestep these problems rather than repair them. One is “gravitational baryogenesis without CPT violation,” based on nonminimal curvature-matter couplings rather than a direct \(J^\mu\partial_\mu R\) term. There the baryon current obeys
\[
\nabla_\mu J_B^\mu \simeq -\partial_\mu R\,J_B^\mu,
\]
CP is violated by an \(imR^5\) term, and out-of-equilibrium evolution is supplied by gravitational particle creation. The mechanism can generate \(\Delta_B\sim10^{-10}\) for
\[
1\ \mathrm{GeV}^{-2}\lesssim \lambda \lesssim 10^2\ \mathrm{GeV}^{-2},
\]
while explicitly avoiding higher derivatives of \(R\) in the baryogenesis operator [1909.03034].

The second is a recent unified framework in which the chemical potential is tied to entropy production rather than to \(\dot R\),
\[
\mu_B(t)=\kappa\,T(t)\,\frac{d}{dt}\ln S(t),
\qquad S\equiv a^3 s,
\]
and CP violation is supplied by a gravitational theta term,
\[
S_\theta=\frac{\theta}{64\pi^2}\int d^4x\,\sqrt{-g}\,R\tilde R.
\]
The formalism highlights adiabatic cancellation of oscillatory sources through a universal low-pass factor
\[
F(x)=\frac{1}{\sqrt{1+x^2}},
\qquad x=\omega\tau_{\rm off},
\]
and derives
\[
Y_B \simeq c\,\kappa\,\Pi_{\rm eff}\,\kappa_{\rm inst}\,\theta.
\]
A dilaton UV completion with \(f_\sigma\sim10^{17}\!-\!10^{18}\,\mathrm{GeV}\) is proposed, and the same framework predicts stochastic-gravitational-wave circular polarization
\[
|\Pi|\simeq \pi|\theta|
\]
for small \(|\theta|\) in a Loop Quantum Cosmology bounce [2601.11690].

## 6. Phenomenology, parameter dependence, and present status

Across the literature, the target baryon asymmetry is consistently of order \(10^{-10}\), but specific papers quote somewhat different observational reference values:
\[
0.864_{-0.015}^{+0.016}\times10^{-10},\qquad
9.0\times10^{-11},\qquad
9.2\times10^{-11},\qquad
9.42\times10^{-11}.
\]
This spread reflects differing conventions and observational inputs rather than a substantive disagreement about the required scale [1708.09137, 2006.06534, 2103.13328, 2305.09707].

The parametric structure is simple in most models. Larger \(\dot R\) or its generalized analog enhances the asymmetry, larger \(M_*^2\) suppresses it, and the decoupling temperature enters through the freeze-out relation. Some backgrounds add further dilution or enhancement factors. In vacuum inflation, entropy continues to grow and the frozen asymmetry is diluted by a factor \(\sim0.77\) before the end of inflation [1708.09137]. In anisotropic models, shear enhances \(\dot R\), but CMB and reheating constraints severely limit how much that enhancement can matter in practice [1603.02403]. In running-vacuum, DGP, and Gauss-Bonnet settings, the key advantage is precisely that the source does not vanish for radiation-like matter [1610.01453, 2012.12965, 1607.00545].

The phenomenological status is correspondingly mixed. Many constructions can be tuned to produce the correct order of magnitude: vacuum inflation, running vacuum models, DGP braneworld cosmology, several \(f(T)\), \(f(R,T)\), \(f(T,B)\), \(f(T,T_G)\), and EPN scenarios all report acceptable benchmark values [1708.09137, 1907.13460, 1607.08561, 2009.13293, 2504.06133]. Other cases are disfavored because the sign is wrong, the asymmetry is far too small, or the required parameters are implausible. Examples include the TEGR benchmark with \(n_B/s\simeq-1.02\times10^{-42}\), the \(\partial_\mu T\) coupling in the specific \(f(R,T)=R+\zeta RT\) model, the \(n=1\) energy-momentum-squared example requiring \(M_B\sim160.96\) GeV, and several \(f(T,\mathcal T)\) fits that use extreme parameter magnitudes [1607.08561, 1907.13460, 2409.04623, 2305.09707].

The current research landscape therefore supports two simultaneous conclusions. First, gravitational baryogenesis remains a productive effective framework for linking baryon asymmetry to background geometry, and modified gravity can provide many mechanisms for keeping the source active when standard GR would suppress it. Second, the canonical \((\partial_\mu R)J^\mu\) formulation is not automatically consistent as a dynamical theory of gravity, and a credible model must address higher-derivative backreaction, cosmological stability, and the plausibility of the required parameter region. Recent work has accordingly shifted part of the emphasis from merely obtaining \(n_B/s\sim10^{-10}\) to formulating versions that evade adiabatic cancellation, survive stability tests, and offer external observables such as stochastic-gravitational-wave polarization [1712.04627, 2301.08322, 2601.11690].

Source: https://www.emergentmind.com/topics/gravitational-baryogenesis