---
title: Weylian Boundary Terms in Gravity and Field Theory
url: https://www.emergentmind.com/topics/weylian-boundary-terms
type: topic
---

# Weylian Boundary Terms in Gravity and Field Theory

Searching arXiv for recent and foundational papers on Weylian boundary terms, boundary Weyl anomaly, and Weyl-transverse gravity.
Weylian boundary terms are boundary contributions in gravitational or conformal field-theoretic constructions in which Weyl geometry, Weyl-covariant structures, or Weyl-anomalous boundary invariants are promoted from auxiliary boundary data to dynamical or variationally essential ingredients. In the recent cosmological framework of Harko, Lobo, and Shahidi, the Einstein–Hilbert bulk action is kept Riemannian while the boundary is treated as Weylian, so that the boundary variation is evaluated with a Weyl connection and Weyl covariant derivative rather than being canceled by a Gibbons–Hawking–York term [2603.21766]. In quantum field theory and holography, the phrase also refers to Gibbons–Hawking-type boundary companions of bulk Weyl invariants in conformal anomalies, as well as purely boundary conformal invariants built from trace-free extrinsic curvature [1510.04566], [1702.00566], [2111.14713]. A distinct but related usage appears in Weyl-transverse gravity, where Weyl-invariant boundary structures organize the variational principle, covariant phase space, and conserved charges in the presence of a fixed background volume form [2601.15976]. Across these settings, the common theme is that boundary geometry is no longer a passive bookkeeping device: it contributes independent Weyl-covariant data, modifies the variational problem, and may generate effective dynamics or anomaly coefficients.

## 1. Weylian boundary geometry and variational re-interpretation

In the cosmological construction of a Universe with a Weylian boundary, the bulk geometry remains Riemannian and metric-compatible, with Levi–Civita connection \(\Gamma^\lambda_{\mu\nu}\) and \(\nabla_\mu g_{\alpha\beta}=0\), while the boundary and its neighborhood are assumed to have Weyl-type non-metricity [2603.21766]. The defining condition is
\[
\tilde{\nabla}_\mu g_{\alpha\beta} = -\alpha\,\omega_\mu\, g_{\alpha\beta},
\]
where \(\omega_\mu\) is the Weyl gauge vector and \(\alpha\) is a dimensionless Weyl gauge coupling. The corresponding Weyl connection is
\[
\tilde{\Gamma}^\lambda_{\mu\nu} = \Gamma^\lambda_{\mu\nu} + \Psi^\lambda_{\mu\nu},
\]
with
\[
\Psi^\lambda_{\mu\nu} = \frac{\alpha}{2}\left(\delta^\lambda_\mu\,\omega_\nu + \delta^\lambda_\nu\,\omega_\mu - g_{\mu\nu}\,\omega^\lambda\right).
\]

This construction retains the Einstein–Hilbert bulk action
\[
S_g = -\frac{1}{2\kappa^2}\int_\Omega R(g)\sqrt{-g}\,d^4x + S_m,
\]
but departs from standard GR at the level of the boundary variation. Instead of canceling the variation of the Ricci tensor by adding a Gibbons–Hawking–York term or by imposing \(\delta g_{\mu\nu}|_{\partial\Omega}=0\), the boundary part is reinterpreted in Weyl geometry through the replacements
\[
\Gamma \to \tilde{\Gamma},\qquad \nabla \to \tilde{\nabla},\qquad \delta\Gamma \to \delta\tilde{\Gamma}.
\]
This yields a modified variation in which the boundary term becomes a term proportional to \(\delta g^{\mu\nu}\) with coefficients built from \(\omega_\mu\) and its Weyl-covariant derivatives [2603.21766]. The resulting expression,
\[
g^{\mu\nu}\tilde{\nabla}_\lambda\delta\tilde{\Gamma}^\lambda_{\mu\nu} - g^{\mu\lambda}\tilde{\nabla}_\lambda\delta\tilde{\Gamma}^\sigma_{\mu\sigma},
\]
is treated as a genuine contribution to the field equations. In this formulation, a “Weylian boundary” therefore means that the metric \(g_{\mu\nu}\) is the same as in the bulk, the non-metricity is confined to boundary-related structures, and the extra field \(\omega_\mu\) enters only through boundary-related terms [2603.21766].

A closely related warm-inflation model adopts the same conceptual move, but specializes to an integrable Weyl geometry \(\omega_\mu=\nabla_\mu\psi\), thereby trading the Weyl vector for a scalar \(\psi\) that encodes the boundary effect [2407.10262]. In that setting, the boundary term in the Einstein–Hilbert variation is re-evaluated in Weyl geometry and produces effective scalar contributions in the field equations. This suggests a broader usage of “Weylian boundary terms” for geometric boundary corrections induced by Weyl non-metricity rather than by explicit bulk Weyl actions.

## 2. Field equations generated by Weylian boundary terms

When the Weylian boundary variation is inserted into the modified action principle and rewritten in Riemannian covariant form, the generalized Einstein equations become [2603.21766]
\[
G_{\mu\nu} + \frac{1}{2}\alpha\left( \nabla_\mu\omega_\nu + \nabla_\nu\omega_\mu - 2g_{\mu\nu}\nabla_\lambda\omega^\lambda \right) - \alpha^2\left( \omega_\mu\omega_\nu + \frac{1}{2}\omega^2 g_{\mu\nu} \right) = \kappa^2 T_{\mu\nu},
\]
with \(\omega^2=\omega_\lambda\omega^\lambda\). The additional terms are purely geometric and involve only the Weyl vector and its first derivatives.

These equations may be read as Einstein equations sourced by an effective geometric sector. The formulation
\[
G_{\mu\nu} = \kappa^2 T_{\mu\nu} + T^{(\omega)}_{\mu\nu}
\]
isolates the Weyl contribution as an effective stress-energy built from derivative terms and quadratic condensate-like terms in \(\omega_\mu\) [2603.21766]. The paper does not introduce \(T^{(\omega)}_{\mu\nu}\) as a primary tensorial object, but the later cosmological expressions for effective density and pressure are direct projections of these boundary-induced terms.

A significant feature of the construction is that ordinary matter is not conserved. Taking the divergence of the generalized field equations and using the Bianchi identity yields [2603.21766]
\[
\begin{aligned}
\kappa^2\nabla^\nu T_{\mu\nu} &= \frac{1}{2}\alpha\left( \omega^\alpha R_{\mu\alpha} + \Box\omega_\mu - \nabla_\mu\nabla_\alpha\omega^\alpha \right) \\
&\quad - \alpha^2\left( \omega_\mu\nabla_\alpha\omega^\alpha + \omega^\alpha\nabla_\alpha\omega_\mu + \omega^\alpha\nabla_\mu\omega_\alpha \right).
\end{aligned}
\]
The matter and Weyl sectors therefore exchange energy-momentum. The warm-inflation realization reaches a parallel conclusion: the total effective fluid satisfies a global conservation law, but the radiation, inflaton, and Weyl scalar sectors exchange energy according to how the total conservation equation is split phenomenologically [2407.10262]. In both constructions, non-conservation is not an incidental artifact but a direct consequence of promoting boundary Weyl contributions to dynamical relevance without supplementing them by an independent bulk action for \(\omega_\mu\) or \(\psi\).

## 3. Cosmological realizations: dark energy, inflation, and BBN

In a spatially flat FLRW geometry,
\[
ds^2 = dt^2 - a^2(t)\big(dx^2+dy^2+dz^2\big),
\]
homogeneity and isotropy force the Weyl vector in the dark-energy model to take the purely temporal form \(\omega_\mu=(\omega_0(t),0,0,0)\) [2603.21766]. The generalized Friedmann equations then become
\[
3H^2 = \kappa^2\rho + \frac{3}{2}\alpha\omega_0(2H + \alpha\omega_0),
\]
\[
2\dot{H} + 3H^2 = -\kappa^2 p + \frac{1}{2}\alpha\omega_0(4H + \alpha\omega_0) + \alpha\dot{\omega}_0.
\]
The effective geometric dark-energy density and pressure are
\[
\rho_{\text{eff}} = \frac{3}{2}\alpha\omega_0(2H + \alpha\omega_0),
\qquad
p_{\text{eff}} = -\frac{1}{2}\alpha\omega_0(4H + \alpha\omega_0) - \alpha\dot{\omega}_0.
\]
Late-time acceleration is thus attributed to Weylian boundary terms of purely geometric origin, with no scalar potential, no explicit cosmological constant, and no exotic matter in the action [2603.21766]. To close the background system, the effective dark-energy equation of state is imposed phenomenologically through the Barboza–Alcaniz parametrization,
\[
\gamma(z) = \gamma_0 + \gamma_a \frac{z(1+z)}{1+z^2},
\]
so that \(w_{\text{eff}}(z)=\gamma(z)\) [2603.21766].

The warm-inflation model uses an integrable Weyl boundary, \(\omega_\mu=\partial_\mu\psi\), and introduces a dissipative scalar inflaton \(\phi\) with Lagrangian
\[
\mathcal L_\phi = e^{\beta\phi}\left[\frac12 g^{\mu\nu}\partial_\mu\phi\,\partial_\nu\phi - V(\phi)\right]\sqrt{-g}
\]
[2407.10262]. The FLRW equations acquire Weylian boundary contributions through \(\psi\),
\[
3H^2 = \frac{8\pi G}{c^2}\rho_m + e^{\beta\phi}\left(\frac12\dot\phi^2 + V(\phi)\right) + \frac{\alpha^2}{2}\left(\ddot\psi + 3H\dot\psi + \dot\psi^2\right),
\]
\[
2\dot H + 3H^2 = -\frac{8\pi G}{c^2}p_m - e^{\beta\phi}\left(\frac12\dot\phi^2 - V(\phi)\right) + \frac{\alpha^2}{2}\left(\ddot\psi + 3H\dot\psi\right).
\]
Three distinct models are then defined by different splittings of the total conservation equation, corresponding respectively to mixed inflaton–Weyl radiation production, radiation driven only by the inflaton, and radiation driven only by the Weyl boundary field [2407.10262]. Numerical solutions show that radiation is created during the very early phases of expansion, reaches a maximum, and is followed by a transition from an accelerating inflationary phase to a decelerating one, with the Weyl vector playing a significant role during radiation creation [2407.10262].

The post–warm-inflationary consequences of this framework have been constrained at the BBN epoch. In the BBN analysis, the same generalized Friedmann system is specialized to the radiation era and the extra energy density \(\rho_w\) coming from the dissipative scalar and the Weyl boundary is treated as a small correction to standard radiation [2509.01162]. The analytic Helium-4 constraint yields
\[
\left|\frac{\delta T_f}{T_f}\right| < 4.7\times 10^{-4},
\]
which implies
\[
\rho_w < 1.35\times 10^{-15}\ \text{GeV}^4 = 1.35\times 10^{-3}\ \text{MeV}^4
\]
at freeze-out [2509.01162]. A numerical implementation based on \texttt{PRyMordial}, together with a Genetic Algorithm and MCMC pipeline implemented in the python code \texttt{genesys}, finds that the Weylian boundary contribution at BBN corresponds to \(\Delta N_{\rm eff}\sim \mathcal O(10^{-3})\), so it remains effectively indistinguishable from standard BBN at current observational precision [2509.01162]. This constrains any strong cosmological role of Weylian boundary terms to epochs prior to the MeV era.

## 4. Observational status and comparison with \(\Lambda\)CDM

The late-time FLRW realization of Weylian boundary terms has been confronted with cosmic chronometers, Pantheon\(^+\), and DESI DR2 BAO data [2603.21766]. In the dimensionless redshift formulation, the Hubble rate is written as
\[
h(z) = \omega_1 + \omega_2,
\qquad
\omega_2 = \sqrt{3\omega_1^2 + \bar{\rho}_m},
\]
with evolution equations for \(\bar{\rho}_m(z)\) and \(\omega_1(z)\) closed by the Barboza–Alcaniz equation of state [2603.21766]. The model is then fit jointly in the parameter space \(H_0,\Omega_{m0},\gamma_0,\gamma_a\).

For the combined CC + Pantheon\(^+\) + BAO dataset, the best-fit parameters are reported as
\[
H_0 = 67.31^{+1.61}_{-1.67}\,\mathrm{km\,s^{-1}\,Mpc^{-1}},
\quad
\Omega_{m0} = 0.111^{+0.017}_{-0.017},
\]
\[
\gamma_0 \simeq -0.609^{+0.031}_{-0.030},
\quad
\gamma_a \simeq -0.597^{+0.073}_{-0.082},
\]
with \(\omega_1(0)\approx 0.333^{+0.0051}_{-0.0053}\), \(\mathcal M\approx -19.416\pm0.051\), and \(r_d\approx147.3\pm3.4\) Mpc [2603.21766]. For \(\Lambda\)CDM on the same dataset, the fit is
\[
H_0 = 68.65^{+1.63}_{-1.63}\,\mathrm{km\,s^{-1}\,Mpc^{-1}},
\qquad
\Omega_{m0} = 0.307^{+0.008}_{-0.008}.
\]
The reduced chi-squared values are essentially identical,
\[
\chi_\text{red}^2 \approx 1.048\quad\text{for }\Lambda\text{CDM},
\qquad
\chi_\text{red}^2 \approx 1.045\quad\text{for Weyl boundary},
\]
while the Bayesian evidence gives
\[
\ln\mathcal Z_{\Lambda\text{CDM}} = -929.38\pm0.28,
\qquad
\ln\mathcal Z_\text{Weyl} = -924.49\pm0.28,
\]
hence
\[
\ln B_{W\Lambda} = 4.887\pm0.335,
\]
interpreted on Jeffreys’ scale as strong (near “moderate–strong”) evidence in favor of the Weyl boundary model over \(\Lambda\)CDM [2603.21766]. The best-fit \(H_0\) is also consistent with Planck 2018 at the level of \(\sim 0.02\sigma\), although no SH0ES calibration is included [2603.21766].

A recurrent misconception is that these models simply repackage a cosmological constant. The data analysis does show that the background expansion can reproduce almost exactly the predictions of the \(\Lambda\)CDM paradigm, but the mechanism is structurally different: the effective dark energy arises from Weylian boundary terms rather than from an explicit \(\Lambda\) term in the action [2603.21766]. A second misconception is that the framework is a fully specified vector-tensor theory. It is not: there is no dynamical equation for \(\omega_\mu\) derived from a bulk action, and the closure of the cosmological system relies on phenomenological assumptions such as the effective equation of state [2603.21766].

## 5. Weylian boundary terms in conformal anomaly and boundary CFT

A separate and older line of research uses “Weylian boundary terms” to denote boundary contributions associated with Weyl anomalies on manifolds with boundary. In this setting, the key statement is that each type B bulk invariant \(I_k(W)\), polynomial in the Weyl tensor, must be supplemented by a Gibbons–Hawking-type boundary companion \(J_k(W,K)\) to ensure a well-posed variational principle [1510.04566]. Solodukhin’s construction introduces a Weylian GH tensor \(P_{\alpha\beta\mu\nu}\), built from extrinsic geometry and the metric, with the same algebraic symmetries as the bulk Weyl tensor. This tensor is trace-free and depends only on the trace-free extrinsic curvature \(\hat K_{ij}\), making it a conformal tensor suitable for boundary anomaly functionals [1510.04566].

In four dimensions, the integrated anomaly contains the bulk invariant \(\int_{M_4}\sqrt{g}\,\mathrm{Tr}W^2\) together with the boundary term
\[
-8\int_{\partial M_4}\sqrt{\gamma}\,W_{i n j n}\,\hat K^{ij},
\]
and a purely extrinsic boundary invariant \(\int_{\partial M_4}\sqrt{\gamma}\,\mathrm{Tr}\hat K^3\) [1510.04566]. In six dimensions, the bulk type B invariants \(\mathrm{Tr}_1(W^3)\), \(\mathrm{Tr}_2(W^3)\), and \(\mathrm{Tr}(W\Delta W)+\mathrm{Tr}_2(W X W)\) acquire analogous boundary companions built by replacing one occurrence of \(W\) with \(P\), together with purely extrinsic terms such as \(\mathrm{Tr}\hat K^4\) and \((\mathrm{Tr}\hat K^2)^2\) [1510.04566]. In odd dimensions, there is no local bulk Weyl anomaly density, but the integrated anomaly can still be non-vanishing because of boundary conformal invariants; in \(d=3\), for example, it is entirely a boundary effect [1510.04566].

The holographic analysis of the integrated conformal anomaly in \(d=3\) and \(d=4\) confirms the same general structure [1702.00566]. In \(d=4\), the integrated anomaly contains boundary terms involving \(W_{ninj}\hat k^{ij}\) and \(\mathrm{Tr}\,\hat k^3\), with coefficients interpreted as boundary central charges. For \(\mathcal N=4\) SYM, the minimal-surface holographic prescription reproduces the free-field boundary anomaly for supersymmetric boundary conditions, while the Takayanagi prescription introduces a parameter \(m\) that changes the boundary charges \(b_1\) and \(c\) [1702.00566].

The classification problem becomes considerably richer for four-dimensional conformal boundaries and defects in ambient dimension \(d\ge 5\). A systematic Wess–Zumino-consistent construction finds, for a 4D defect of codimension \(q\ge2\), one A-type invariant \(\overline E_4\), twenty-two parity-even B-type invariants, and six universal parity-odd B-type invariants, with additional parity-odd structures for special codimensions \(q=2,4\) [2111.14713]. For codimension one, corresponding to a 4D boundary or defect in \(d=5\), the anomaly reduces to the structure previously identified by Astaneh and Solodukhin, together with three parity-odd terms [2111.14713]. Several of these central charges enter the displacement operator two-point function, the stress-tensor one-point function, and the universal part of entanglement entropy [2111.14713]. This establishes that Weylian boundary terms in the anomaly context are not merely variational counterterms; they define intrinsic boundary or defect data of the quantum field theory.

## 6. Boundary conditions, holographic Weyl structures, and Weyl-transverse gravity

The technical treatment of Weylian boundary terms is closely tied to the broader problem of variational principles with boundary. In higher-order theories, adding total divergences is harmless only when the correct higher-order boundary data are fixed. If only the fields are fixed and derivatives are left unconstrained, boundary terms can change the solution space and may even render it empty [1106.3738]. This warning is especially relevant for Weyl or metric-affine gravity, because any attempt to treat boundary terms as dynamical must specify boundary conditions consistent with the formal order of the action [1106.3738]. A plausible implication is that the cosmological Weylian-boundary construction, being defined through a modified variation rather than a complete covariant bulk action for \(\omega_\mu\), should be interpreted as an effective field-equation-level model rather than a standard higher-order variational theory.

Weyl-transverse gravity provides a more systematic framework in which boundary terms are constructed directly from Weyl-invariant geometric data [2601.15976]. WTG uses a dynamical metric \(g_{\mu\nu}\), a fixed non-dynamical volume form \(\boldsymbol{\omega}\), and an auxiliary Weyl-invariant metric \(\tilde g_{\mu\nu}\) satisfying \(\sqrt{-\tilde g}=\omega\) [2601.15976]. The natural boundary Lagrangians are built from the extrinsic curvature \(\tilde K\) of the auxiliary geometry, for example
\[
\ell_b = \frac{b}{16\pi G}\,\tilde K\,\bm\omega_\Gamma,
\]
and admit auxiliary Dirichlet, auxiliary Neumann, dynamical Dirichlet, and York boundary conditions [2601.15976]. In particular, York boundary conditions take the form
\[
\delta q_{\langle\mu\nu\rangle}|_\Gamma = 0,
\qquad
\delta\tilde K|_\Gamma = 0,
\]
which fix the conformal class of the induced dynamical metric and the mean extrinsic curvature defined with \(\tilde g\) [2601.15976]. The resulting covariant phase-space formalism shows that Weyl transformations lie in the kernel of the presymplectic form, while conserved charges receive a distinctive contribution proportional to \((\Lambda-\lambda)\xi\cdot\bm\omega\) due to the fixed background volume form [2601.15976].

A holographic realization of boundary Weyl structure appears in the generalized Fefferman–Graham gauge for asymptotically AdS\(_3\) spacetimes [2411.12513]. In this gauge, the bulk metric is written as
\[
ds^2 = \ell^2\left(\frac{dz}{z} - a_i(z,x)\,dx^i\right)^2 + h_{ij}(z,x)\,dx^i dx^j,
\]
and the leading coefficient \(a_i^{(0)}\) transforms inhomogeneously under radial Weyl diffeomorphisms,
\[
a^{(0)}_i \to a^{(0)}_i - \partial_i\log \mathcal B,
\]
exactly like a boundary Weyl connection [2411.12513]. The pair \((\gamma^{(0)}_{ij},a^{(0)}_i)\) therefore defines a boundary Weyl structure. Holographic renormalization in this framework uncovers new subleading divergences generated by the cross term \(g_{zi}\), which are canceled by a codimension-two counterterm,
\[
\Upsilon = \gamma^{ij}\,\bar\nabla_i n_j,
\]
so that the renormalized action remains finite [2411.12513]. The Weyl anomaly then takes the Weyl-covariant form
\[
\langle T^i{}_i + \nabla_i J^i\rangle = -\frac{\ell}{16\pi G}\,\hat{\mathcal R}^{(0)},
\]
with \(J^i\) the current conjugate to the Weyl connection [2411.12513]. This explicitly realizes Weylian boundary terms as boundary and corner counterterms required by a Weyl-covariant holographic variational principle.

In Palatini gravity and its dual formulations, Weyl invariance also constrains admissible boundary terms. For the unique Palatini \(f(R)\) theory admitting an Eddington reduction in \(D\neq2\), namely \(f(R)=cR^{D/2}\), the bulk action is Weyl invariant and the connection assumes the form of a Weyl connection relative to \(g_{\mu\nu}\) [2210.10383]. The appropriate generalized Gibbons–Hawking–York term is
\[
S_{\partial}^{f(R)} = \int_{\partial\mathcal{M}} d^{D-1}y\, \sqrt{|\gamma|}\, f(R)\, K,
\]
with \(\delta R|_{\partial\mathcal{M}}=0\) [2210.10383]. The paper shows that this boundary term maps consistently under the duality between Weyl-invariant Palatini \(f(R)\), Born–Infeld–Einstein gravity, and the Eddington affine action [2210.10383]. In two dimensions, the corresponding first-order bosonic string formulation requires the boundary term
\[
S_{\partial}^{\text{string}} = \int_{\partial\mathcal{M}} dy\, \big(K + 2 U^\mu n_\mu\big),
\]
which is invariant under the conformal symmetry of the 2D theory [2210.10383].

Taken together, these developments show that “Weylian boundary terms” is not a single concept but a family of related constructions. In cosmology, they are dynamical geometric sources induced by a Weylian reinterpretation of the boundary variation [2603.21766], [2407.10262], [2509.01162]. In conformal anomaly theory, they are boundary companions of bulk Weyl invariants and independent conformal invariants of extrinsic geometry [1510.04566], [1702.00566], [2111.14713]. In Weyl-transverse and holographic frameworks, they are Weyl-invariant boundary and corner terms that make the action differentiable and render the charge algebra and anomaly structure manifestly Weyl covariant [2601.15976], [2411.12513]. The unifying idea is that the boundary carries genuine Weyl-covariant information, and once that information is retained rather than discarded, it can alter field equations, anomaly functionals, conserved charges, and cosmological evolution.

Source: https://www.emergentmind.com/topics/weylian-boundary-terms