Operator Product Expansion (OPE)

• Operator Product Expansion (OPE) is a formalism in quantum field theory that rewrites products of nearby fields as an asymptotic sum of local operators with singular coefficient functions.
• It clearly separates leading short-distance singularities from sub-leading terms using scaling analysis, and applies in both flat and curved spacetimes.
• OPE underpins renormalization group flow and conformal bootstrap methods, providing a unifying structure for symmetry constraints and non-perturbative quantum field analyses.

The operator product expansion (OPE) is a foundational structure in quantum field theory (QFT) that encodes the short-distance singular behavior of products of local quantum fields in a universal, model-independent expansion. Originally introduced by Wilson and Zimmermann, the OPE expresses the product of fields at nearby spacetime points as an asymptotic sum of local fields at a single point, each accompanied by a singular coefficient function (“OPE coefficient”). The OPE formalism provides a unifying language for renormalization, symmetry constraints, and non-perturbative analysis in quantum and conformal field theory, with deep connections to RG flows, algebraic QFT in curved backgrounds, integrability, the conformal bootstrap, and modern approaches to quantum gravity.

1. Formal Definition and Mathematical Structure

Let $\{O_i(x)\}$ be a basis of composite local fields on a (possibly curved) spacetime $M$. The OPE posits that for any finite set of points $x_1, \ldots, x_n$ and a reference point $y$, the product

$$O_1(x_1) \cdots O_n(x_n) \sim \sum_{k} C_k(x_1,\ldots,x_n; y) O_k(y)$$

as $x_1,\ldots,x_n \to y$, where $\sim$ denotes equality in the sense of scaling degree (asymptotic expansion), and the $C_k$ are c-number distributions (OPE coefficients) capturing all singular structure. The basis $\{O_k(y)\}$ comprises local fields (composite operators), and the expansion is organized by engineering or canonical dimension $d_k$, defined via the scaling degree of the leading OPE coefficient with the identity operator (Hollands et al., 2023).

The precise nature of the asymptotic relation uses the scaling degree $\mathrm{sd}[\cdot]$ of distributions at the total diagonal $x_1 = \dots = x_n = y$. The singularity structure is thus governed by the lowest-dimension fields appearing on the right-hand side.

2. Singular Structure and Scaling Analysis

The OPE separates leading short-distance singularities from finite or sub-leading trends. In the canonical case of the free massless Klein–Gordon field in four-dimensional Minkowski space,

$$\phi(x_1)\,\phi(x_2) \sim \frac{1}{2\pi^2 (x_1 - x_2)^2} + \phi^2(y) + \ldots$$

as $x_1, x_2 \to y$, with the first term proportional to the identity operator and encoding universal, state-independent divergences. The next term involves a new local field (e.g., the renormalized composite $\phi^2(y)$), which is finite in the coincidence limit. Sub-leading terms correspond to operators of increasing dimension, whose coefficients vanish more rapidly as the points coalesce.

In curved spacetime, the Minkowski interval is replaced by the squared geodesic distance and a Van Vleck determinant factor, with subleading terms involving local curvature invariants such as the Ricci scalar $R$ (Hollands et al., 2023). The classification of singularities is via the scaling degree: the lowest-dimensional operators contribute the most divergent terms as $x_i \to y$.

3. Renormalization Group Flow and Flow Relations

The OPE coefficients satisfy exact differential “flow equations” under deformations by coupling parameters, which generalize renormalization group concepts. For a theory depending on a parameter $\lambda$ (e.g., a mass squared or a $\phi^4$ coupling), one can derive differential relations:

$$\partial_\lambda C^{(n)}_k = \sum_p \int d^dy\, X(y;z)\, C^{(n+1)}_p (y,x_1,\ldots,x_n; z) - \text{subtractions}$$

where $X(y;z)$ is an IR cutoff, and the subtraction terms ensure covariance and remove UV divergences (Hollands et al., 2023). These relations, together with initial data from the free theory, permit the perturbative construction of interacting OPE coefficients order by order in the coupling.

This structure underlies the renormalization group interpretation of the OPE, making the entire collection of coefficients a dynamical object that flows with the theory’s parameters. In CFT, this viewpoint leads directly to RG equations for the scaling dimensions and structure constants (Hollands, 2017).

4. OPEs in Curved Spacetime and Axiomatic Constraints

In globally hyperbolic curved spacetimes, the OPE remains well-defined when the coefficients $C_k[M](x_1,\ldots,x_n; y)$ are required to depend locally and covariantly on the metric, orientation, spin structure, and other background data (Axiom C1). Additional axioms include:

• Curvature ambiguities: Allowed in the definition of composite fields, e.g., $\phi^2 \to \phi^2 + c R$ for a scalar field.
• Microlocal spectrum condition: A wave-front set bound reminiscent of the Hadamard condition ensures correct singularity structure (Axiom C7).
• Associativity: Multiple ways of fusing subsets of points and then approaching the diagonal must yield identical total OPE coefficients (Axiom C8).
• Hermiticity, symmetry, and smooth background dependence, among others (Hollands et al., 2023).

Under these axioms, the OPE delivers the full local content of an interacting quantum field theory: specifying all $C_k$ consistent with (C1)–(C9) is equivalent to fully defining the QFT on $M$.

5. OPE in Conformal Field Theory: Structure and Bootstrap

In conformal field theories (CFTs), the OPE is highly constrained:

• Local fields decompose into conformal multiplets labeled by primary operators ($\mathcal{O}$) of scaling dimension $\Delta$ and spin $\ell$, with descendants generated by derivatives.
• The OPE for primaries is

$$\mathcal{O}_i(x_1)\mathcal{O}_j(x_2) = \sum_k f_{ijk}\, |x_{12}|^{\Delta_k - \Delta_i - \Delta_j}\, P_{ij}^k(x_{12},\partial) \mathcal{O}_k(x_2)$$

where $f_{ijk}$ are constant structure constants, and $P_{ij}^k$ are fixed differential operators determined by the conformal group (Hollands et al., 2023).

• The “conformal data” $\{\Delta_i, \ell_i, f_{ijk}\}$ fully encode the CFT. The OPE, when substituted into four-point functions, leads to crossing relations

$$\sum_{\mathcal{O}} f_{12\mathcal{O}} f_{34\mathcal{O}} G^s_{\Delta,\ell}(u,v) = \sum_{\mathcal{O}} f_{14\mathcal{O}} f_{23\mathcal{O}} G^t_{\Delta,\ell}(u,v)$$

with $G^{s,t}_{\Delta,\ell}(u,v)$ the conformal blocks and $(u,v)$ the standard cross ratios. This constraint underlies the conformal bootstrap.

• In unitary CFTs on $\mathbb{R}^d$, the OPE converges (not just asymptotes) for sufficiently small separations (Hollands et al., 2023).

6. Structural Theorems: PCT via OPE and Fundamental Examples

The OPE axioms permit the formulation and proof of structural field-theoretic properties, including the PCT (parity, charge conjugation, time reversal) theorem, even in curved spacetime without global discrete symmetries. Specifically, for $M$ and its time-reversed partner $\bar M$,

$$C_{j_1 \ldots j_n}^i[M](x_1,\ldots,x_n;y) = (-1)^{\sum U_j} (-1)^{\sum F_j} \overline{C_{j_1 \ldots j_n}^{\,i}[\bar{M}](x_1,\ldots,x_n;y)}$$

where $U_j$ counts spinor indices and $F_j$ fermionic parity. In flat Minkowski space, additional symmetry identifies this formula with the standard PCT theorem.

Illustrative examples include:

• Free fields: The OPE for the free massless Klein–Gordon field and its curved spacetime generalization.
• CFT two- and three-point functions: Standard power-law dependence dictated by conformal invariance, with explicit OPE structure constants.
• Flow equations: Explicit $\phi^4$-theory OPE coefficient relations and conformal bootstrap equations (Hollands et al., 2023).

7. Significance and Universality

The OPE has become a universal organizing principle in quantum and conformal field theory due to its ability to:

• Encode singularities, operator content, and short-distance expansions into finitely many (but generally infinite) sets of coefficient functions.
• Systematize renormalization and RG flow, both perturbative and non-perturbative.
• Provide a foundation for axiomatic and algebraic formulations of QFT, including in curved backgrounds and for theories lacking standard Lagrangian descriptions.
• Directly constrain physical observables through bootstrap and associativity equations, giving rise to rigorous non-perturbative results in CFT (Hollands et al., 2023).

The OPE thus serves not only as an expansion or computational tool but as a candidate for the “defining structure” of full quantum field theories, bridging perturbative, geometric, and algebraic formulations.