---
title: Conformal Killing Gravity Cosmology
url: https://www.emergentmind.com/topics/conformal-killing-gravity-cosmology
type: topic
---

# Conformal Killing Gravity Cosmology

Conformal Killing Gravity cosmology denotes the cosmological sector of Harada’s third-order gravitational theory, later reformulated as Einstein gravity supplemented by a divergence-free conformal Killing tensor. In this framework, the cosmological constant can appear as an integration constant, the field equations can be reduced from third to second order on Robertson–Walker and generalized Robertson–Walker backgrounds, and the additional conformal-Killing contribution behaves as a geometric dark sector rather than as an independently postulated dark-energy fluid. The resulting cosmologies admit modified Friedmann equations with an extra term proportional to \(a^{2}\) or, equivalently in the usual redshift parametrization used in several works, \((1+z)^{-2}\); they can interpolate between decelerating matter eras and late-time acceleration, but they also raise selection and consistency problems that remain active subjects of discussion [2308.06803], [2404.11468], [2409.14353], [2502.06262].

## 1. Foundational equations and Einstein-type reformulation

Harada’s starting point is a third-order tensor equation
\[
H_{jkl}=8\pi G\,T_{jkl},
\]
with
\[
H_{jkl} = \nabla_j R_{kl} +\nabla_k R_{lj} +\nabla_l R_{jk} -\frac13\left(g_{kl}\nabla_jR+g_{lj}\nabla_kR+g_{jk}\nabla_lR\right).
\]
Because \(H_{jkl}\) contains derivatives of the Ricci tensor, the theory is third order in the metric. The central result of Mantica and Molinari is that, in the cosmological setting considered there, this equation can be integrated once and rewritten as
\[
R_{jk}-\frac12 g_{jk}R = T_{jk}+K_{jk},
\]
where \(K_{jk}\) is an arbitrary symmetric tensor required to be both divergence-free and conformal Killing. In that notation,
\[
K_{jk}=R_{jk}-\frac12 g_{jk}R-T_{jk}.
\]
The Bianchi identity then gives
\[
\nabla^j K_{jk}=0,
\]
and therefore
\[
\nabla^j T_{jk}=0,
\]
so energy-momentum conservation remains a consequence of the field equations in the integrated formulation [2308.06803].

This reduction is the decisive conceptual simplification in conformal Killing gravity cosmology. It turns Harada’s equations from third order to second order and places the theory within a standard Einstein-type cosmological workflow, with \(K_{jk}\) acting as an additional geometric source. In parallel presentations of the theory, Harada’s field equations are also written in trace-modified form as
\[
H_{abc}=G_{(ab;c)}=\tilde G_{(ab;c)}, \qquad T_{abc}=T_{(ab;c)}=\tilde T_{(ab;c)},
\]
with \(\tilde G_{ab}=R_{ab}-Rg_{ab}\) and \(\tilde T_{ab}=T_{ab}-\tfrac12 Tg_{ab}\). In this language, any GR solution of \(G_{ab}=T_{ab}\) and any Einstein-\(\Lambda\) solution of \(G_{ab}+\Lambda g_{ab}=T_{ab}\) also solves conformal Killing gravity, which explains why the theory inherits the standard homogeneous and isotropic sector while extending it by additional integration constants and source ambiguities [2404.09310].

## 2. Conformal Killing tensors, perfect-fluid form, and Robertson–Walker geometry

A symmetric tensor \(K_{ij}\) is conformal Killing when it satisfies
\[
\nabla_i K_{jk}+\nabla_j K_{ki}+\nabla_k K_{ij} = \eta_i g_{jk}+\eta_j g_{ki}+\eta_k g_{ij},
\]
where \(\eta_i\) is the associated conformal vector. If \(\eta_i=0\), one recovers an ordinary Killing tensor; if \(\eta_i=\nabla_i\phi\), one has a gradient conformal Killing tensor. In the divergence-free case, the contraction formula used in the literature gives
\[
\nabla^j K_{jk}=0 \quad \Longleftrightarrow \quad \eta_k=\nabla_k K,
\]
with \(K=g^{ij}K_{ij}\) [2308.06803].

For cosmology, the relevant tensors have perfect-fluid form,
\[
K_{ij}=A\,g_{ij}+B\,u_i u_j,\qquad u^i u_i=-1,
\]
or, in the explicit cosmological notation used in the integration theorem,
\[
K_{ij}=\lambda g_{ij}+\frac{4}{3}\lambda\,u_i u_j.
\]
In generalized Robertson–Walker spacetime, such a tensor exists precisely when the spacetime has the GRW structure and the velocity field is torse-forming,
\[
\nabla_i u_j = H\,(u_i u_j + g_{ij}),
\]
with \(H=\dot a/a\) in comoving coordinates. This geometrizes the Hubble function directly at the level of the conformal Killing tensor [2308.06803].

A more specialized construction introduces Sinyukov-like tensors. In the perfect-fluid ansatz
\[
K_{ij}=A g_{ij}+B u_i u_j,
\]
\(K_{ij}\) is conformal Killing iff the velocity field is shear-free. If it is also divergence-free, then
\[
H=\frac{\dot B}{2B},\qquad \nabla_j B=-\dot B\,u_j+2B\,\dot u_j,
\]
and
\[
\eta_j=\nabla_j A+\dot B\,u_j.
\]
If, in addition, the flow is acceleration-free, the tensor becomes a special divergence-free CKT of Sinyukov-like type, and the spacetime is generalized RW; with conformal flatness it is RW. In RW spacetime the tensor is fixed by the scale factor \(a(t)\) and two constants \(C\) and \(\Lambda\):
\[
B(t)=\frac{C}{n-1}a^2(t),\qquad A(t)=\frac{n+1}{2}B(t)-\Lambda.
\]
For \(n=4\),
\[
K_{kl} = g_{kl}\left[\frac{5}{6}Ca^2-\Lambda\right] + u_k u_l \frac{C a^2}{3}.
\]
This is the tensorial origin of the cosmological dark term in the Robertson–Walker sector [2404.11468].

## 3. Modified Friedmann equations and the geometric dark sector

With RW geometry and ordinary matter
\[
T_{jk}=p\,g_{jk}+(p+\rho)u_j u_k,
\]
the conformal-Killing contribution converts the Einstein equation into modified Friedmann relations. In the notation with \(8\pi G\) restored, Mantica and Molinari derive
\[
8\pi G\,p = -\frac13 R-\xi-\lambda-\frac{2}{3}\Lambda = -\frac{R^\star}{6a^2}-\frac{\ddot a}{a}-\frac{\Lambda}{3}+\frac{Ca^2}{6},
\]
and
\[
8\pi G\,\rho = -\frac13 R-\xi+\lambda = \frac{R^\star}{2a^2}+3\frac{\dot a^2}{a^2}-\Lambda+\frac{Ca^2}{2}.
\]
Here \(R^\star\) is the scalar curvature of the spatial slices, \(k=R^\star/6\), \(\xi=3(\dot H+H^2)=3\ddot a/a\), and \(C\) is the integration constant generated by the conformal-Killing integration. Eliminating the conformal-tensor parameter gives
\[
\frac{5\rho+3p}{4\pi G} = \frac{R^\star}{3a^2} +2\frac{\ddot a}{a} -\frac{4\Lambda}{3},
\]
while the continuity equation retains the standard form
\[
\rho(a)=\rho_0\left(\frac{a_0}{a}\right)^{3(w+1)}
\]
for \(p=w\rho\) [2308.06803].

In the equivalent Sinyukov-like description, the same cosmology is written as
\[
\mu=\frac{R^\star}{2a^2}+3H^2+\frac12 Ca^2-\Lambda,
\]
\[
p=-\frac{R^\star}{6a^2}-3H^2-2\dot H-\frac{5}{6}Ca^2+\Lambda.
\]
For dust plus radiation,
\[
\mu=\mu_m+\mu_r,\qquad p=p_r=\frac13\mu_r,
\]
with
\[
\mu_m=\mu_{m0}\left(\frac{a}{a_0}\right)^{-3},\qquad \mu_r=\mu_{r0}\left(\frac{a}{a_0}\right)^{-4},
\]
the generalized Friedmann equation becomes
\[
\left(\frac{H}{H_0}\right)^2 = \Omega_r\left(\frac{a}{a_0}\right)^{-4} + \Omega_m\left(\frac{a}{a_0}\right)^{-3} + \Omega_k\left(\frac{a}{a_0}\right)^{-2} + \Omega_\Lambda + \Omega_D\left(\frac{a}{a_0}\right)^2,
\]
with
\[
\Omega_m+\Omega_r+\Omega_k+\Omega_\Lambda+\Omega_D=1.
\]
In redshift space,
\[
\left(\frac{H(z)}{H_0}\right)^2 = \Omega_r(1+z)^4+\Omega_m(1+z)^3+\Omega_k(1+z)^2 +\Omega_\Lambda+\frac{\Omega_D}{(1+z)^2}.
\]
The associated effective dark density and pressure are
\[
\mu_D=-\frac12 Ca^2+\Lambda,\qquad p_D=\frac56 Ca^2-\Lambda.
\]
Thus the conformal Killing sector is represented as a geometric dark fluid whose density grows as \(a^2\) [2404.11468].

A closely related exact solution writes the cosmological evolution in Friedmann form as
\[
\left(\frac{H}{H_0}\right)^2 = \Omega_{\rm m}\left(\frac{a}{a_0}\right)^{-3} +\Omega_{\rm r}\left(\frac{a}{a_0}\right)^{-4} +\Omega_k\left(\frac{a}{a_0}\right)^{-2} +\Omega_\Lambda +\Omega_{\rm eff}\left(\frac{a}{a_0}\right)^2,
\]
with
\[
\Omega_{\rm eff}\equiv 1-\Omega_{\rm m}-\Omega_{\rm r}-\Omega_k-\Omega_\Lambda.
\]
Since \(\rho_{\rm eff}\propto a^2\), the continuity equation fixes
\[
\omega=-\frac53.
\]
This is the origin of the phantom-like effective equation of state often associated with the conformal Killing dark term [2308.07634].

## 4. Cosmological solution space and late-time evolution

The flat, dust-dominated FLRW sector is the canonical example. Setting
\[
R^\star=0,\qquad \Lambda=0,\qquad p=0,
\]
one obtains
\[
H^2(a) = -\frac{C}{6}a^2+\frac{D}{3a^3},
\]
equivalently
\[
\frac{\dot a^2}{a^2} = -\frac{C}{6}a^2+\frac{D}{3a^3},
\]
with \(D=8\pi G\rho_0 a_0^3\), and
\[
\frac{\ddot a}{a} = -\frac{C}{3}a^2-\frac{D}{6a^3}.
\]
This is precisely the form that yields decelerating early evolution and late-time acceleration when the \(Ca^2\) term becomes dominant; in the notation used there, \(C<0\) corresponds to accelerated expansion at late times. This is the standard conformal Killing gravity realization of “accelerating expansion without dark energy” in a matter-dominated universe [2308.06803].

A second exact sector is the constant-curvature case. The scale factor is written as
\[
a(t)=A\sinh\theta,
\]
with
\[
\theta=\frac{t}{A}\sqrt{\frac{R^\star}{3}},
\]
and the effective equation of state tends to
\[
w\to \frac13 \quad (t\sqrt{R^\star}\ll 1),
\]
and
\[
w\to -1 \quad (t\sqrt{R^\star}\gg 1),
\]
or, in the Harada extension discussed there, to a phantom-like regime with \(w=-5/3\) in the far future when radiation is included. This provides a radiation-to-dark-energy interpolation generated by geometry [2308.06803].

Vacuum or sourceless FLRW solutions further enlarge the cosmological phase space. With
\[
\frac{\dot a^2}{a^2}=\frac{8\pi G\rho}{3}+\frac{\Lambda}{3}-\frac{k}{a^2}+\alpha a^2,
\]
the new integration constant \(\alpha\) measures the deviation from GR; \(\alpha=0\) reproduces the standard Friedmann equation. In vacuum,
\[
\dot a^2+V(a)=0,\qquad V(a)=-\alpha a^4-\frac{\Lambda}{3}a^2+k.
\]
The classification by \(\alpha\), \(k\), and \(\Delta=\Lambda+4k\alpha\) yields singularity-free eternal cosmologies, closed oscillating universes, and universes evolving symmetrically from a big bang to a big crunch within a finite lapse of time. For example, with \(k>0\), \(\alpha<0\), and \(\Lambda>6\sqrt{-k\alpha}\), the scale factor oscillates periodically between two finite radii; with \(k<0\) and \(\alpha<0\), the universe evolves symmetrically from a big bang to a big crunch in finite time [2404.00328].

Late-time fate depends sharply on the sign of the dark parameter in data-driven RW analyses. If \(\Omega_D>0\), then as \(z\to -1\),
\[
H^2(z)\to +\infty,
\]
and the model predicts a big-rip singularity; one estimate quoted for the future time interval is
\[
\tau \approx 21.0~\text{Gy}.
\]
If \(\Omega_D<0\), then \(H^2\) reaches zero at finite scale factor, and in the CC+BAO fit this occurs around
\[
z_c\simeq -0.407.
\]
This suggests that conformal Killing cosmology naturally supports both phantom-ending and finite-turnaround futures, depending on the sign of the geometric dark term [2404.11468].

## 5. Observational fits, growth of structure, and early-universe consistency

The Robertson–Walker version with
\[
\left(\frac{H(z)}{H_0}\right)^2 = \Omega_m(1+z)^3+\Omega_\Lambda+\frac{\Omega_D}{(1+z)^2}
\]
was tested on cosmic chronometers and BAO. In the flat, late-time approximation \((\Omega_k=\Omega_r=0)\), the CC-only fit gave approximately
\[
\Omega_m \approx 0.311,\qquad \Omega_\Lambda \approx 0.368,\qquad \Omega_D \approx +0.321,
\]
while the CC+BAO fit gave approximately
\[
\Omega_m \approx 0.306,\qquad \Omega_\Lambda \approx 1.103,\qquad \Omega_D \approx -0.410.
\]
The matter density remains close to the \(\Lambda\)CDM value, whereas \(\Omega_\Lambda\) and \(\Omega_D\) have large uncertainties and can change sign; their sum remains broadly close to the \(\Lambda\)CDM dark-energy fraction. In the same framework, the linear matter contrast in the matter-dominated era obeys
\[
\ddot\delta_m+2H\dot\delta_m - \left[ \frac32H^2-\frac{5C}{12}a^2-\frac{\Lambda}{2} \right]\delta_m=0,
\]
and the numerical conclusion is that the conformal Killing dark sector gives no significant deviation from \(\Lambda\)CDM and GR results during matter domination [2404.11468].

A later Bayesian analysis used COBAYA with DESI DR2 or SDSS DR16 BAO data and Pantheon+ or Union3 supernova compilations in a spatially flat FRW background. In that presentation, the late-time Hubble law was written as
\[
\left(\frac{H(z)}{H_{0}}\right)^{2} =\Omega_M(1+z)^{3}+\Omega_{\Lambda}+\Omega_{D}(1+z)^{2},
\]
with \(\Omega_\Lambda = 1-\Omega_M-\Omega_D\). The reported constraints were
\[
H_0 = 66.44 \pm 0.93,\quad \Omega_M = 0.310 \pm 0.009,\quad \Omega_D = -0.235 \pm 0.095
\]
for DESI+Union3,
\[
H_0 = 66.61 \pm 0.58,\quad \Omega_M = 0.309 \pm 0.008,\quad \Omega_D = -0.221 \pm 0.057
\]
for DESI+Pantheon+,
\[
H_0 = 65.77 \pm 1.04,\quad \Omega_M = 0.311 \pm 0.037,\quad \Omega_D = -0.25 \pm 0.19
\]
for DR16+Union3, and
\[
H_0 = 66.11 \pm 0.74,\quad \Omega_M = 0.301 \pm 0.036,\quad \Omega_D = -0.243 \pm 0.14
\]
for DR16+Pantheon+. The same study reported \(q_0\) values near \(-0.3\), present dark-energy parameters \(w_{D0}>-1\) in the quintessence regime, and first-acoustic-peak estimates \(100\theta_s\) extraordinarily near the Planck best value \(1.04097\pm0.00046\) [2508.02603].

Taken together, these studies indicate that conformal Killing gravity can reproduce a background expansion history close to \(\Lambda\)CDM, with the main empirical novelty concentrated in the split between \(\Omega_\Lambda\) and the geometric dark parameter \(\Omega_D\), rather than in large deviations of \(\Omega_m\) or in large departures of linear growth. A plausible implication is that the model’s observational distinctiveness lies primarily in late-time background evolution and future asymptotics rather than in matter-era structure formation [2404.11468], [2508.02603].

## 6. Symmetry-based extensions and non-conservative generalizations

A symmetry-oriented extension studies scalar-field flat FLRW cosmology through the Eisenhart lift. The lifted minisuperspace metric is
\[
ds^2 = -6a\,da^2 + a^3\,d\phi^2 + \frac{1}{2a^3V(\phi)}\,d\chi^2,
\]
with null Hamiltonian constraint
\[
G_{AB}\dot\varphi^A\dot\varphi^B=0.
\]
A non-trivial conformal Killing vector of this lifted geometry supplies an extra conserved quantity and complete integrability. The determinant condition for the prolonged conformal Killing equations reduces to a nonlinear second-order ODE for
\[
h(\phi)=\frac{V'(\phi)}{V(\phi)},
\]
and its regular branch yields the most general local potential in the \(\chi\)-independent sector,
\[
V(\phi)=V_0 \left(\cos\beta\,e^{\alpha\phi}+\sin\beta\,e^{-\alpha\phi}\right)^{-2+\frac{\sqrt6}{\alpha}}.
\]
The singular branch solves the determinant equation but is incompatible with the full conformal Killing equations. This establishes a precise integrable class of conformal-Killing-invariant flat FLRW scalar cosmologies [2604.22247].

A different extension abandons separate matter conservation. In non-conservative conformal Killing gravity,
\[
\nabla^{p}T_{pl}=\nabla_{l}(\sigma R-\tau T),
\]
and the Einstein-like form becomes
\[
R_{kl}-\frac{R}{2}g_{kl}=T_{kl}-\sigma Rg_{kl}+\tau Tg_{kl}+K_{kl}.
\]
Here
\[
\Theta_{kl}=K_{kl}-\sigma Rg_{kl}+\tau Tg_{kl}
\]
is the effective dark-sector tensor. In FRW, CMB considerations are used to impose \(\sigma=0\), so the dark sector couples only with the trace of the stress-energy tensor. Dust then satisfies
\[
(1+\tau)\dot{\mu}_{M}+3H\mu_{M}=0,
\]
with solution
\[
\mu_{M}=\mu_{M0}\left(\frac{a}{a_{0}}\right)^{-\frac{3}{1+\tau}},
\]
and the Hubble function becomes
\[
\left(\frac{H(z)}{H_{0}}\right)^{2}=\Omega_{R}(1+z)^{4}+\Omega_{M}(1+z)^{\frac{3}{1+\tau}}+\Omega_{k}(1+z)^{2}+\Omega_{\Lambda}+\frac{\Omega_{D}}{(1+z)^{2}}.
\]
The corresponding effective dark fluid has \(w_D(a)\approx -1\) at early times and \(w_D(a)\approx -5/3\) at late times, so the interacting variant preserves the phantom-like asymptotics of the conservative model while changing matter dilution [2606.08213].

## 7. Open problems, selection issues, and formal criticisms

The cosmological promise of conformal Killing gravity is counterbalanced by substantial open problems. In vacuum, assuming constant scalar curvature reduces the field equations to
\[
\nabla_{(\alpha}R_{\beta\gamma)}=0,
\]
so the Ricci tensor becomes a Killing tensor. This admits Einstein metrics
\[
R_{\mu\nu}=\Lambda g_{\mu\nu},
\]
symmetric spaces, Ricci-symmetric spaces, and many non-Einstein Kundt solutions. In the Kundt sector, the metric function satisfies a linear third-order PDE, and examples include Siklos and Defrise geometries, \(\mathrm{AdS}_3\times\mathbb{R}\), \(\mathbb{M}_2\times\mathbb{H}_2\), \(\mathbb{M}_2\times S^2\), and higher-dimensional \(\mathrm{AdS}_m\times\mathbb{R}^n\). The resulting “Pandora’s box” problem is that the theory contains a myriad of vacuum solutions, but does not by itself provide a criterion selecting the observed isotropic and homogeneous FLRW branch with positive effective cosmological term [2409.14353].

A second issue is consistency. One detailed criticism is that a gravitational theory based on a rank-3 field equation does not arise from a diffeomorphism-invariant action with the metric as dynamical field, so the usual route from Bianchi identities and Noether charges to conserved quantities is unavailable. In this analysis, Schwarzschild remains a solution, but the usual mass parameter \(m\) has no clear interpretation as a conserved charge of the theory. Attempts to build conserved currents yield quantities that vanish identically on Einstein manifolds, including Schwarzschild and Kerr-type solutions. This undermines standard notions of energy, angular momentum, and gravitational-wave flux, and therefore raises doubts about the physical completeness of the cosmological framework [2502.06262].

A related structural caveat is matter-source non-uniqueness. In the pp-wave sector, the general vacuum solution is
\[
H=f(z,u)+\bar f(\bar z,u)+c\,z\bar z,
\]
so the GR wave profile is supplemented only by a non-propagating quadratic term. Yet the same metric can correspond to different matter sources, because if two trace-modified tensors differ by a Killing tensor, they generate the same \(H_{abc}\). Even without additional symmetries, one may add
\[
T_{ab}=\lambda g_{ab}
\]
without changing the three-tensor source, and with Killing vectors the ambiguity increases further. This suggests that the “dark sector” of conformal Killing cosmology is geometrically natural but not uniquely tied to a single matter interpretation [2404.09310].

These objections do not eliminate the cosmological constructions summarized above, but they fix the present status of the subject. Conformal Killing gravity cosmology is a mathematically explicit program in which accelerated expansion can arise from a divergence-free conformal Killing tensor, \(\Lambda\) can emerge as an integration constant, and RW dynamics can be treated within a second-order Einstein-like system. At the same time, vacuum multiplicity, source ambiguity, and the absence of a fully satisfactory conserved-charge structure remain central obstacles to a definitive physical interpretation [2308.06803], [2409.14353], [2502.06262].

Source: https://www.emergentmind.com/topics/conformal-killing-gravity-cosmology