---
title: 'Plebanski''s Formulation: Constrained BF Gravity'
url: https://www.emergentmind.com/topics/plebanski-s-formulation
type: topic
---

# Plebanski's Formulation: Constrained BF Gravity

Plebanski’s formulation is a reformulation of four-dimensional general relativity as a constrained \(BF\) theory. In its non-chiral form the basic variables are a bivector-valued 2-form \(B^{IJ}\) and an \(SO(1,3)\) or \(Spin(4)\) connection; in its chiral form they are a triple of self-dual 2-forms \(\Sigma^i\) and an \(SO(3,\mathbb{C})\) or \(SU(2)\) connection. The defining feature is the simplicity constraint, which forces the 2-form data to arise from a tetrad, after which the metric can be reconstructed by the Urbantke map and the field equations become Einstein or Einstein–Cartan equations in the gravitational sector. The same framework underlies spin-foam models, pure-connection actions, heavenly equations for self-dual metrics, and several matter-coupled and unimodular variants [1009.5595][1204.4225][2405.15408].

## 1. Constrained \(BF\) structure and simplicity constraints

In a standard non-chiral Lorentzian normalization, the Plebanski variables are an antisymmetric 2-form \(B^{IJ}\), an \(SO(1,3)\) connection one-form \(A^{IJ}\), and a Lagrange multiplier \(\phi_{IJKL}\). The action with cosmological constant is
\[
S[B,A,\phi] \;=\; \int \left( B^{IJ} \wedge F_{IJ}(A) \;-\; \frac{1}{2}\,\phi_{IJKL}\,B^{IJ}\wedge B^{KL} \;+\; \frac{\Lambda}{12}\,\epsilon_{IJKL}\,B^{IJ}\wedge B^{KL} \right),
\]
with \(F^{IJ}(A)=dA^{IJ}+A^{I}{}_{K}\wedge A^{KJ}\). In the chiral \(SU(2)\) version one writes instead
\[
S[B^i,A^i,\Psi_{ij}] \;=\; \int \left( B^i \wedge F^i(A) \;-\; \frac{1}{2}\,(\Psi_{ij}-\tfrac{\Lambda}{3}\delta_{ij})\,B^i\wedge B^j \right),
\]
with symmetric traceless \(\Psi_{ij}\) [1009.5595].

The simplicity constraints are the algebraic conditions that distinguish gravity from unconstrained \(BF\) theory. In non-chiral form they may be written as
\[
B^{IJ}\wedge B^{KL} \;=\; \frac{1}{6}\,\epsilon^{IJKL}\, B^{MN}\wedge B_{MN},
\]
while in spinor language they appear as
\[
\Sigma^{(AB}\wedge \Sigma^{CD)} \;=\; 0,
\qquad
\bar\Sigma^{(A'B'}\wedge \bar\Sigma^{C'D')} \;=\; 0.
\]
These equations force the 2-forms to be simple, meaning built from tetrads rather than arbitrary bivector data [1009.5595].

Classically, the solutions separate into distinct Plebanski sectors. In the quadratic classification one has
\[
B^{IJ}=\pm e^I\wedge e^J,
\qquad
B^{IJ}=\pm \star(e^I\wedge e^J),
\]
together with a degenerate sector characterized by vanishing 4-volume density. The sectors \( (I\pm)\) and \((II\pm)\) differ by the internal Hodge star, whereas the degenerate sector does not admit a nondegenerate metric interpretation [1107.0709].

A closely related Euclidean \(Spin(4)\) form used in spin-foam work is the Holst-\(BF\) action
\[
S = -\frac{1}{2\kappa} \int \mathrm{tr}\big[(B + \frac{1}{\gamma} *B) \wedge F\big],
\]
with momentum
\[
J := \frac{1}{\kappa} \left(B + \frac{1}{\gamma} *B\right).
\]
This makes the Immirzi parameter \(\gamma\) explicit at the \(BF\)-theory level and is the starting point for the discrete \(Spin(4)\) analysis of the EPRL/FK vertex [1107.0709].

## 2. Metric reconstruction and equivalence to Einstein gravity

Once the simplicity constraints hold, the 2-forms determine a tetrad and hence a metric. In chiral notation, given a tetrad \(E^I\) with \(\eta_{IJ}=\mathrm{diag}(-1,+1,+1,+1)\), the self-dual 2-forms are
\[
\Sigma^i = i\,E^0 \wedge E^i - \frac{1}{2}\,\epsilon^{i}{}_{jk}\,E^j \wedge E^k.
\]
The metric is recovered from the 2-forms by the Urbantke formula; one representative form is
\[
g(\xi,\eta)\,\mathrm{vol}_g
=
-\frac{1}{12}\,\epsilon_{ijk}\,(\iota_\xi B^i)\wedge(\iota_\eta B^j)\wedge B^k.
\]
Thus the metric is not fundamental but reconstructed from a triple of 2-forms [2309.06333][2412.15779].

Varying the action with respect to the connection gives the compatibility equation
\[
D_A B^i = 0
\qquad\text{or}\qquad
D_\omega \Sigma^{AB}=0.
\]
In the gravitational sector this identifies the gauge connection with the self-dual part of the spin connection compatible with the tetrad. In the complex \(SO(3,\mathbb{C})\) treatment this relation is
\[
A^i = i\,\omega^{0i}(e) + \frac{1}{2}\,\epsilon^{i}{}_{jk}\,\omega^{jk}(e)
      \equiv \frac{1}{2}\,\omega^{IJ}(e)\,\Sigma^i{}_{IJ}.
\]
A central point of the gauge-versus-spacetime connection analysis is that this identification is not automatic from the Plebanski equations alone; it follows once the Urbantke metric and Cartan structure equations are used to relate the internal bundle to spacetime geometry. Under those hypotheses, the internal connection depends only on the Levi-Civita part, independently of torsion and non-metricity that may be present in an arbitrary spacetime connection used in the intermediate derivation [1204.4225].

Variation with respect to the 2-forms gives the curvature equation
\[
F^i = M^{ij}\,\Sigma^j
\]
or, in Einstein form,
\[
F^i_{\mu\nu}
=
\Psi_{ij}\,\Sigma^j_{\mu\nu}
+
\frac{\Lambda}{3}\,\delta_{ij}\,\Sigma^j_{\mu\nu}.
\]
Here the symmetric tracefree field \(\Psi_{ij}\) becomes the self-dual Weyl curvature on shell. In the \(SU(2)\)-structure formulation one may write
\[
R_{\mu\nu} = -\,\Sigma^i{}_\mu{}^\alpha\,F^i_{\alpha\nu},
\qquad
R = \Sigma^{i\,\mu\nu}\,F^i_{\mu\nu},
\]
so the Einstein condition is equivalent to the statement that the anti-self-dual part of \(F^i\) vanishes and the scalar part is fixed [2405.15408].

The canonical \(3+1\) decomposition recovers the complex Ashtekar phase space. With
\[
\tilde E^a{}_i = \frac{1}{2}\,\epsilon^{abc}\,\Sigma^i_{bc},
\]
the fundamental Poisson bracket is
\[
\{A^i_a(x),\tilde E^b{}_j(y)\} = 8\pi G\,\delta^b_a\,\delta^i{}_j\,\delta^{(3)}(x-y),
\]
and the first-class constraints are the Gauss, diffeomorphism, and Hamiltonian constraints. This is the standard canonical bridge from the covariant Plebanski action to Ashtekar’s variables [2411.09748].

## 3. Plebanski sectors, discretization, and spin-foam quantization

In spin-foam applications one discretizes the continuum theory on a triangulation, often focusing on a single oriented 4-simplex. The continuum 2-form is integrated over triangles to give bivectors \(B_{ab}^{IJ}\), or \(J_{ab}^{IJ}\) in Holst-\(BF\) variables, and two discrete conditions are fundamental: orientation,
\[
B_{ab}(S) = -B_{ba}(S),
\]
and closure,
\[
\sum_{b\neq a} B_{ab} = 0.
\]
For such data there exists a unique constant continuum 2-form \(B_{\mu\nu}^{IJ}(\sigma)\) on an embedded simplex \(\sigma\subset \mathbb{R}^4\) reproducing the discrete bivectors, which makes it meaningful to speak of the Plebanski sector of discrete data [1107.0709].

A crucial result for the EPRL/FK models is that the linear simplicity constraints do not isolate a single gravitational sector. In a time gauge they imply
\[
B^+_{ab}=B^-_{ab},
\qquad
u^+_{ab}=-u^-_{ab},
\]
so the discrete bivectors can be parameterized by areas and 3D normals. However, the stationary-point analysis shows that linear simplicity admits precisely three sectors: \((II+)\), \((II-)\), and the degenerate sector. For nondegenerate critical data one has
\[
B^{\mathrm{phys}}_{ab}(A,n,X^\pm)=\mu\,B^{\mathrm{geom}}_{ab}(\sigma),
\qquad
\mu=\pm 1,
\]
so \(\mu=+1\) gives sector \((II+)\) and \(\mu=-1\) gives sector \((II-)\); if the two chiral transport sets are equivalent, the data are degenerate. The resulting asymptotics attribute the two Regge terms \(e^{\pm i\sum A\Theta}\) to \((II\pm)\), while the extra \(\gamma\)-dependent terms \(e^{\pm (i/\gamma)\sum A\Theta}\) and vector-geometry contributions come from the degenerate sector. The paper explicitly argues that these extra terms arise from sector mixing rather than a sum over manifold orientations [1107.0709].

The degenerate sector is independently important. In the \(Spin(4)\) theory with the usual trace condition on the Lagrange multiplier removed, the degenerate simplicity constraints are exactly solvable and split into “deg-gravitational” and “deg-topological” subsectors,
\[
B_{\mu\nu}^{IJ} = \frac{1}{2}\,\epsilon^{IJ}{}_{KL} x^K b^L_{\mu\nu},
\qquad
B_{\mu\nu}^{IJ} = x^{[I} b^{J]}_{\mu\nu},
\]
with \(x_I b^I_{\mu\nu}=0\). The deg-gravitational subsector reduces to covariantly embedded \(SU(2)\) \(BF\) theory, and its spin-foam quantization is the \(SU(2)\) Crane–Yetter state sum [1202.5039].

This solvable sector also clarifies a quantization issue. Restricting \(Spin(4)\) representations and intertwiners, the usual strategy in EPRL/FK, is not sufficient to recover the correct vertex amplitude. The proposed remedy is to impose the secondary second-class constraints in the vertex measure on an equal footing with the primary simplicity constraints. In the degenerate model this produces the \(SU(2)\) \(\{15j\}\) vertex and the full \(SU(2)\) Crane–Yetter/Ooguri state sum, in agreement with canonical path-integral quantization [1202.5039].

## 4. Pure-connection, \(BF+\)potential, instanton, and unimodular variants

A major development of the Plebanski framework is the elimination of auxiliary fields to obtain pure-connection or \(BF+\)potential descriptions. For \(\Lambda\neq 0\), one pure-connection action for self-dual gravity is
\[
S[A,\Psi]
=
\frac{1}{16\pi G}
\int
\left(\Psi_{ij}+\frac{\Lambda}{3}\delta_{ij}\right) F^i[A]\wedge F^j[A],
\]
and the self-dual gravity truncation takes the form
\[
S_{\mathrm{SDGR}}[A,\Psi] = \int \Psi_{ij} F^i\wedge F^j.
\]
Variation with respect to \(\Psi_{ij}\) imposes \(F^i\wedge F^j \propto \delta^{ij}\), so the metric can again be reconstructed by Urbantke, now directly from curvature rather than from an independent \(B\)-field [2412.15779].

A related line of work shows that there is a family of field redefinitions acting on \(BF\)-type Lagrangians of gravity. In the chiral case this is a one-parameter family; in the non-chiral case it is a two-parameter family. Lagrangians related by these transformations are classically equivalent, differing by topological terms and reparameterizations of auxiliary matrices. In particular, the chiral transformation gives an alternative derivation of the \(BF+\)potential formulation of general relativity, while the non-chiral analysis yields a new \(BF+\)potential form of GR and clarifies the status of non-chiral pure-connection actions [1708.07694].

Another Plebanski-like action modifies the way the symmetric auxiliary matrix enters:
\[
S[A,B,\rho,\Psi]
=
\int
\Bigl[
B_i\wedge F^i
+
\frac{1}{2}(\Psi_{ij}-\lambda\delta_{ij})\,B^i\wedge B^j
+
(\beta\,\mathrm{Tr}(\Psi^{-1})-\gamma)\rho
\Bigr].
\]
Its advantages are that the symmetric matrix can be integrated out, leading to a pure \(BF\)-type action for GR, and that the special choice \(\lambda=0\), \(\beta\neq 0\) yields conformally anti-self-dual gravity. The canonical analysis shows that the resulting phase space matches the Ashtekar formalism up to a canonical transformation induced by a topological term [1605.06212].

The instanton representation of Plebanski gravity is another reformulation adapted to the self-dual sector. Its basic fields are a gauge connection and a nondegenerate \(3\times 3\) matrix \(\Psi_{ae}\), with \(\Sigma_{\mu\nu}^a=\Psi_{ae}F_{\mu\nu}^e\). In the nondegenerate sector, together with the symmetry condition on \(\Psi\), the algebraic constraint \(\Lambda+\mathrm{tr}\,\Psi^{-1}=0\), and the Gauss law, this is equivalent to Plebanski’s equations. The formulation makes Hodge self-duality of \(F\) explicit and is tailored to gravitational instantons [1207.4534].

Unimodular Plebański gravity replaces the fixed cosmological constant parameter by an integration constant. In the fixed-volume version the action is
\[
S[A,\Sigma,M;\omega_0]
=
\frac{-i}{8\pi G}
\int
\left(
\Sigma^i\wedge F^i
-
\frac{1}{2}M^{ij}\Sigma^i\wedge\Sigma^j
+
\frac{1}{2}\omega_0\,\mathrm{tr}M
\right),
\]
while the parametrized Henneaux–Teitelboim version uses an exact 4-form \(dT\) instead of \(\omega_0\). In both cases \(\mathrm{tr}M\) is not fixed a priori: in the fixed-volume theory it becomes constant on shell by the Bianchi identity, and in the parametrized theory \(d\,\mathrm{tr}M=0\) is a field equation. These theories also admit pure-connection reductions, such as
\[
S_{\mathrm{UPC}}[A;\omega_0]
=
\frac{-i}{8\pi G}\int \omega_0\,\mathrm{tr}\sqrt{X},
\]
with \(F^i\wedge F^j = X^{ij}\omega_0\) [2309.06333].

The canonical analysis of the unimodular variants retains the Ashtekar pair \((A^i_a,\tilde E^a{}_i)\) but changes the status of the Hamiltonian. In the preferred-volume theory the Hamiltonian density is constrained to be spatially constant rather than to vanish, while in the parametrized theory a new canonical pair \((\tilde\rho,\lambda)\) appears, with \(\lambda\) constrained to be spatially and temporally constant. This yields one additional global degree of freedom interpreted as the cosmological constant, together with a natural “volume time” built from the exact 4-form \(dT\) [2411.09748].

## 5. Matter couplings, source terms, and torsion

Matter coupling in the Plebanski formalism is less direct than in metric GR because the basic gravitational variables are 2-forms and connections rather than a metric and its Levi-Civita connection. For tensor matter, a general prescription is to take an ordinary matter action \(S_m[g,Q]\) and replace the metric by the Urbantke metrics built from the self-dual and anti-self-dual 2-forms:
\[
S_m[\Sigma,\bar\Sigma,Q]
=
S_m(g(\Sigma),Q)+S_m(g(\bar\Sigma),Q).
\]
For Dirac spinors one reconstructs a tetrad from the Urbantke metric and inserts it into the Einstein–Cartan Dirac action. In the gravitational sector these prescriptions reproduce exactly the Einstein–Cartan equations with matter; for bosonic matter they give the usual torsion-free coupling, while for fermions the torsion equation is algebraic and sourced by the spin current [1009.5595].

The same paper verifies concrete scalar, Yang–Mills, and Dirac couplings. For Yang–Mills, an auxiliary-field action based on \(\Sigma^{AB}\), \(\bar\Sigma^{A'B'}\), and symmetric spinors \(\phi_{AB}\), \(\bar\phi_{A'B'}\) reproduces the standard Yang–Mills stress tensor after eliminating auxiliaries. For spinors, the \(\omega\)-equation becomes the Einstein–Cartan torsion equation and the curvature equations reduce to the Einstein–Cartan relations once the translation formulas between \(\Sigma\)-based and tetrad-based sources are used [1009.5595].

A more direct translation of metric matter sources into chiral Plebanski variables starts from the decomposition
\[
T_{\mu\nu}=\hat T_{\mu\nu}+\frac{T}{4}g_{\mu\nu},
\qquad
\hat T_{\mu\nu}=T_{\mu\nu}-\frac{T}{4}g_{\mu\nu},
\]
and lifts \(\hat T_{\mu\nu}\) into the algebraic curvature space by the Kulkarni–Nomizu product
\[
(\hat T \odot g)_{\mu\nu\rho\sigma}
=
\frac{1}{2}\left(g_{\mu[\rho}\hat T_{\sigma]\nu}-g_{\nu[\rho}\hat T_{\sigma]\mu}\right).
\]
The chiral matter source is then
\[
T^i
=
4i\,\big[P^+(\hat T\odot g)\big]^{0i}{}_{\mu\nu}\,dx^\mu\wedge dx^\nu
=
\Sigma^{i\,\rho}{}_{[\mu}\hat T_{\nu]\rho}\,dx^\mu\wedge dx^\nu,
\]
and the matter-coupled field equation becomes
\[
F_i
=
\left(
\Psi_{ij}
+
\frac{1}{3}(\Lambda-2\pi G T)\delta_{ij}
\right)\Sigma^j
+
8\pi G\,T_i.
\]
Applying the chiral Bianchi identity \(d^A F_i=0\) then reproduces the standard conservation law \(\nabla_\mu T^{\mu\nu}=0\). For a spherically symmetric Maxwell source, the anti-self-dual part of the equations yields the Reissner–Nordström–de Sitter metric [2606.20042].

Torsion occupies a special place in the non-chiral formulation. In the ordinary gravitational sector, if \(B^{IJ}=e^I\wedge e^J\) or its dual, then the connection equation
\[
D_A B^{IJ}=0
\]
implies
\[
D_A(e^I\wedge e^J)=T^I\wedge e^J-e^I\wedge T^J=0,
\]
and for nondegenerate tetrads this gives \(T^I=0\). However, through the Nieh–Yan identity
\[
d(e^I\wedge T_I)=T^I\wedge T_I - e^I\wedge e^J\wedge F_{IJ}(A),
\]
a torsion-squared action can be related to the \(BF\) term up to a boundary contribution. In the self-dual setting this leads to the statement that a torsionful phase is equivalent to ordinary gravity or to a topological \(BF\) phase, so torsion need not be treated as an independent propagating variable in these formulations [1209.2155].

## 6. Self-dual sector, heavenly equations, and later developments

The self-dual sector of Plebanski theory admits a remarkable scalar reduction. In the second heavenly formulation, the problem of finding self-dual Einstein metrics is reduced to a nonlinear second-order PDE for a single potential. In one standard coordinate form the equation is
\[
\Omega_{xz} + \Omega_{yw} + \Omega_{xx}\Omega_{yy} - (\Omega_{xy})^2 = 0.
\]
A recent pure-connection derivation obtains both the flat-space version and the constant-curvature analogue by using a complex basis of self-dual 2-forms and a covariant light-cone ansatz for the chiral connection, showing that the heavenly equation emerges directly from the self-dual gravity equations written in pure-connection form [2412.15779].

The same pure-connection analysis emphasizes structural parallels between self-dual gravity and self-dual Yang–Mills. In particular, the vector fields
\[
X_f^\mu := \Omega^{\mu\nu}\partial_\nu f
\]
satisfy
\[
[X_{f_1},X_{f_2}] = X_{[f_1,f_2]},
\qquad
[f_1,f_2]:=\Omega^{\mu\nu}(\partial_\mu f_1)(\partial_\nu f_2),
\]
which gives a concrete realization of the self-dual Yang–Mills kinematic algebra as the Lie algebra of \((0,1)\) vector fields on \(\mathbb{R}^4\) endowed with a complex structure. This provides a direct bridge between Plebanski’s self-dual gravity, pure-connection methods, and color–kinematics structures [2412.15779].

The second heavenly equation also supports a perturbiner expansion for self-dual spacetimes. In the marked-tree formulation, a self-dual background with arbitrarily many positive-helicity gravitons is constructed as a sum over tree graphs, while negative-helicity gravitons are added as linear anti-self-dual perturbations. Evaluating the self-dual Plebanski action supplemented by a boundary term on this background produces the graviton MHV tree amplitude. A distinctive feature of that derivation is that the amplitude comes entirely from the boundary term and reproduces the NSVW tree formula without using BCFW recursion or twistor theory [2408.11139].

More recent work extends the chiral differential-form viewpoint beyond formal equivalence statements. The \(SO(3,\mathbb{C})\)-bundle interpretation treats the basic 2-forms as soldering forms on an associated bundle, develops linearized and nonlinear gauge fixings for the Einstein equations in chiral variables, and uses the self-dual/anti-self-dual decomposition to analyze type \(D\) spacetimes and evolution schemes for numerical relativity. A plausible implication is that Plebanski’s formulation is not only a reformulation of classical GR but also a flexible organizational framework for canonical, covariant, integrable, and numerical approaches [2604.20772].

Taken together, these developments show that Plebanski’s formulation is simultaneously a constrained \(BF\) theory, a route to Ashtekar variables, a foundation for spin-foam amplitudes, a source of pure-connection and \(BF+\)potential actions, and a natural language for self-dual geometry. Its central operation remains unchanged across these contexts: replacing the metric by 2-form data, and recovering Einstein geometry through simplicity, compatibility, and curvature equations.

Source: https://www.emergentmind.com/topics/plebanski-s-formulation