⚛️ high-energy theory

Conformal blocks and braiding matrices for RCFTs I: Virasoro minimal models

This paper presents a systematic method for computing braiding matrices in Virasoro minimal models by directly calculating conformal block series via the Shapovalov form to fix accessory parameters in Fuchsian ODEs, thereby deriving unitary F-matrices that determine structure constants and conjecturing their entries lie in cyclotomic extensions to obtain exact results.

Original authors: Suresh Govindarajan, Aditya Jain, Akhila Sadanandan, Prasanta K. Tripathy

Published 2026-09-18

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Suresh Govindarajan, Aditya Jain, Akhila Sadanandan, Prasanta K. Tripathy

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast landscape of theoretical physics, there exists a special class of theories known as rational conformal field theories. These are mathematical frameworks used to describe systems that look the same at every scale, a property called scale invariance. They are particularly useful for understanding the behavior of materials at critical points, such as when a magnet loses its magnetism or when a liquid turns into a gas. Within these theories, the fundamental building blocks are not particles in the traditional sense, but rather fields that carry specific weights, determining how they interact and transform. The rules governing how these fields can combine are called fusion rules, and the way they exchange positions in space is described by braiding. Together, these rules form a complex algebraic structure known as a modular tensor category, which acts like a blueprint for the entire theory. While physicists have long been able to map out the basic fusion rules and the global symmetries of these systems, a crucial piece of the puzzle has remained elusive: the precise numerical values that dictate how these fields braid around one another. Without these specific numbers, the blueprint is incomplete, and it is impossible to fully distinguish between different theories that might otherwise look identical from a distance.

A team of researchers has now developed a systematic method to calculate these missing braiding numbers, specifically for a family of models known as Virasoro minimal models. These models are among the simplest and most well-understood examples of rational conformal field theories, yet even here, determining the exact braiding values has been a tedious and difficult task, often requiring complex algebraic machinery that becomes unwieldy as the systems grow more complicated. The researchers approached this problem by treating the mathematical equations that describe these fields not as abstract symbols, but as solutions to a specific type of differential equation. They realized that if they could compute just the first few terms of the interaction between four fields, they could reverse-engineer the entire equation that governs their behavior. This equation, which has three special points where the behavior changes, acts like a map. By solving this map numerically, the team could trace how the solutions at one point connect to the solutions at another, effectively revealing the braiding matrices that were previously hidden.

The process began with a direct calculation of the interaction between four fields, focusing on the very first few layers of complexity. Using a method that involves summing up the contributions of all possible intermediate states, the team generated a short series of numbers that describe how the fields interact at close range. They then used these numbers to fix the parameters of the differential equation, a step that usually requires knowing the answer in advance. Once the equation was fully determined, they solved it to find the connection between the starting point and the endpoint of the interaction. This connection is the braiding matrix, a table of numbers that tells us exactly how the system changes when two fields swap places. The team tested their method on numerous examples, ranging from simple cases where all four fields were identical to more complex scenarios where every field was different. In every instance, their numerical results matched the known theoretical values with high precision, confirming that their approach was sound.

A significant part of their work involved refining these numerical results into exact mathematical forms. The researchers noticed a pattern: when they adjusted their calculations to ensure the system behaved in a physically consistent way, the resulting numbers were not random decimals but belonged to a very specific family of numbers related to roots of unity. This observation allowed them to convert their high-precision numerical estimates into exact formulas. For example, in one of their test cases, they were able to derive the exact values for the three-point interaction strengths, which are fundamental constants of the theory. They also found that their method could be extended to theories with more complex symmetries, such as the three-state Potts model, which describes a system with three possible states for each site. By applying their technique to this model, they successfully calculated the braiding matrices for a system that possesses a symmetry beyond the standard rules, demonstrating the versatility of their approach.

The researchers emphasize that their method bypasses the need for the traditional, laborious derivation of these equations from first principles, which often involves checking for "null states" that should not exist in a consistent theory. Instead, by computing the initial terms of the interaction directly, they let the data dictate the form of the governing equation. This bypasses the algebraic tedium that has historically slowed down progress in this area. The team also verified that their calculated matrices satisfy a complex consistency condition known as the hexagon identity, a rule that ensures the braiding operations are logically coherent no matter how the fields are rearranged. This verification serves as a strong check on the validity of their results. While their current work focuses on the Virasoro minimal models, the authors suggest that this technique could be applied to a much wider class of theories, including those based on more exotic symmetries and even to hypothetical theories that have not yet been realized in nature. By providing a reliable way to compute these braiding matrices, the researchers have offered a new tool for classifying and understanding the deep algebraic structures that underpin the quantum world.

Technical Summary: Conformal Blocks and Braiding Matrices for RCFTs I: Virasoro Minimal Models

Problem Statement
The classification of Rational Conformal Field Theories (RCFTs) via the holomorphic modular bootstrap relies on identifying "tenable" solutions to Modular Linear Differential Equations (MLDEs). While the modular data (S- and T-matrices) of these solutions can be classified, this data alone is insufficient to uniquely determine the underlying Modular Tensor Category (MTC). Specifically, the braiding matrices (F-matrices and R-matrices), which encode the associativity and exchange properties of the theory, are not fixed by modular data alone. Existing methods to compute these matrices—such as the Coulomb-gas formalism, boundary CFT, or quantum group approaches—often involve complex algebraic derivations or rely on specific constructions that become tedious for higher-level null states. Furthermore, deriving the Fuchsian Ordinary Differential Equations (ODEs) (BPZ equations) directly from null-state conditions is computationally intensive for orders greater than two. This paper addresses the need for a systematic, direct method to compute braiding matrices for Virasoro minimal models, serving as a proof of concept for extending these methods to tenable solutions with unknown RCFT realizations.

Methodology
The authors develop a hybrid formalism combining direct power-series computation with the theory of Fuchsian ODEs. The methodology proceeds in five distinct steps:

Decomposition and Direct Computation: Using the characters of a tenable solution, the associated modules are decomposed into irreducible Virasoro Verma modules. For a chosen four-point function, the authors compute the first few terms of the "core" conformal block series directly using the Shapovalov form. This involves calculating the Shapovalov matrix (Gram matrix of descendants) and the coefficients $\alpha$ and $\beta$ derived from the Operator Product Expansion (OPE) up to a specific level (typically level 10, sufficient for the cases studied).
ODE Reconstruction: Instead of deriving the BPZ equation from the null-state condition (which is tedious for high levels), the authors use the computed series terms and the known local exponents (determined by the conformal weights of the intermediate primaries) to fix the Fuchsian ODE. The ODE is of order $N$ , determined by the lowest null level of the external operators. The "accessory parameters" of the ODE, which are not fixed by local exponents, are determined by matching the series expansion of the ODE solution to the directly computed conformal block series.
Numerical F-Matrix Extraction: Once the ODE is fully determined, the authors compute the Frobenius solutions at the regular singular points $z=0$ and $z=1$ . The F-matrices (braiding matrices) are then obtained numerically as the connection matrices relating the bases of solutions at these two points.
Unitarization and Exact Reconstruction: The numerically computed F-matrices are transformed into unitary (or orthogonal) matrices by a change of basis involving the square roots of the products of three-point structure constants. The authors propose Conjecture 3.1, stating that the squares of the entries of the unitary F-matrix lie in a cyclotomic extension of the rational numbers, $\mathbb{Q}(\zeta_{2\tilde{N}})$ . This conjecture allows the conversion of high-precision numerical estimates into exact algebraic expressions.
Handling Exceptional Cases: The method is adapted for cases where the ODE order $N$ exceeds the number of physical intermediate primaries ( $N > n$ ) or where exponents differ by integers (resonance). In these cases, "spurious" solutions or logarithmic ambiguities arise, which are resolved by the direct knowledge of the physical conformal block series.

Key Contributions and Results

Systematic ODE Determination: The paper provides a robust algorithm to determine Fuchsian ODEs for conformal blocks without explicitly solving the null-state decoupling equations for high levels. This bypasses the algebraic complexity of deriving BPZ equations for order-3 and order-4 cases.
Numerical-to-Exact Conversion: By leveraging the conjecture on the cyclotomic nature of F-matrix entries, the authors successfully convert numerical connection matrices into exact analytic forms.
Verification: The method is applied to numerous examples in Virasoro minimal models $M(p, q)$ , covering cases with identical, pairwise identical, and distinct operators, as well as orders 2, 3, and 4. The results are cross-verified against the Dotsenko–Fateev Coulomb-gas formalism, showing agreement to arbitrary precision.
Extended Symmetries: The formalism is extended to RCFTs with symmetries extending the Virasoro algebra (e.g., the three-state Potts model with $W_3$ symmetry). The authors demonstrate how to construct blocks for these theories by decomposing extended modules into Virasoro modules and applying gluing coefficients, successfully computing braiding matrices for the Potts model.
Structure Constants: As a by-product of the unitarization process, the method determines products of three-point structure constants (OPE coefficients).

Significance and Claims
The paper claims that this formalism offers a practical and general route to computing braiding matrices, which are essential for fully characterizing the MTC of a theory. The authors emphasize that while the R-matrices are trivial (determined solely by conformal weights), the F-matrices carry the non-trivial calculational complexity.

The significance of this work lies in its potential application to "tenable" solutions of MLDEs that do not yet have a known RCFT realization. By computing the F-matrices for these solutions, one can determine if they correspond to a valid MTC, thereby aiding in the classification of RCFTs. The authors note that modular data alone cannot distinguish between different MTCs with the same fusion rules; the computation of braiding matrices is the necessary step to break this circularity.

The paper concludes by stating that the ultimate goal is to apply this method to tenable solutions with no known realization, potentially identifying their MTC class and verifying their existence as RCFTs. The authors also plan to extend this formalism to RCFTs based on affine Kac-Moody algebras, where the F-matrices relate to quantum group 6j-symbols.

More like this