Exact calculations beyond charge neutrality in timelike Liouville field theory

Timelike Liouville field theory (also known as imaginary Liouville theory or imaginary Gaussian multiplicative chaos) is expected to describe two-dimensional quantum gravity in a positive-curvature regime, but its path integral is not a probability measure and rigorous exact computations are currently available only in the charge-neutral (integer screening) case. In this paper we show that at the special coupling $b=1/\sqrt{2}$, the Coulomb-gas expansion of the timelike path integral becomes explicitly computable beyond charge neutrality. The reason is that the $n$-fold integrals generated by the interaction acquire a Vandermonde/determinantal structure at $b=1/\sqrt{2}$, which allows exact evaluation in terms of classical special functions. We derive Mellin-Barnes type representations (involving the Barnes $G$-function and, in a three-point case, Gauss hypergeometric functions) for the zero- and one-point functions, for an antipodal two-point function, and for a three-point function with a resonant insertion $α_2=b$. We then address the subtle zero-mode integration: after a Gaussian regularization we obtain an explicit renormalized partition function $C(1/\sqrt{2},μ)=e(4π\sqrt2 μ)^{-1}$, identify distributional limits in the physically relevant regime $α_j=\frac{1}{2}Q+\mathrm{i} P_j$, and compare with the Hankel-contour prescription recently proposed in the physics literature. These results provide the first rigorously controlled family of exact calculations in timelike Liouville theory outside charge neutrality.

💡 Research Summary

The paper addresses a long‑standing gap in the rigorous understanding of timelike (or imaginary) Liouville field theory, which is expected to describe two‑dimensional quantum gravity in a positive‑curvature regime. While the spacelike version admits a probabilistic construction via Gaussian multiplicative chaos, the timelike theory suffers from a “wrong‑sign” kinetic term, making its path integral non‑probabilistic and limiting exact results to the charge‑neutral (integer screening) sector. The authors overcome this limitation by focusing on the special coupling constant (b=1/\sqrt{2}). At this value the Coulomb‑gas expansion of the timelike path integral acquires a Vandermonde (determinantal) structure after stereographic projection, allowing the multi‑fold integrals to be evaluated exactly using Selberg‑type techniques and orthogonal polynomial methods.

The paper first decomposes the Liouville field (\phi=c+X) into a zero mode (c) and a mean‑zero Gaussian free field (X). The fixed‑zero‑mode functional (C(\alpha,x,b,\mu,c)) is defined via an imaginary Gaussian free field and expanded into a Coulomb‑gas series. For (b=1/\sqrt{2}) the coefficients (a_n) become determinants, leading to closed‑form expressions in terms of the Barnes (G) function and Gamma functions. The authors then resum the series by a Mellin–Barnes transform, obtaining explicit integral representations for several correlation functions:

• Zero‑point function (Theorem 2.1) – expressed as an integral over (y) involving (\Gamma(1-iy)) and a function (f) built from Barnes (G) and Gamma factors.
• Renormalized partition function (Theorem 2.2) – after regularizing the divergent zero‑mode integral with a Gaussian cutoff (e^{-\epsilon c^{2}}) and taking (\epsilon\to0), the result is (C(1/\sqrt{2},\mu)=e,(4\pi\sqrt{2},\mu)^{-1}).
• One‑point function – similarly obtained, featuring ratios of Barnes (G) functions.
• Antipodal two‑point function – for insertions (\alpha_j=Q/2+iP_j) (the physically relevant regime), the regularized correlator converges to a distribution implementing momentum conservation, i.e. a Dirac delta in the continuous momentum variables.
• Three‑point function with a resonant insertion (\alpha_2=b) – exhibits a genuine pole matching the singularity structure predicted by the timelike DOZZ formula.

A major focus is the treatment of the zero‑mode integration. Besides the real‑line Gaussian regularization, the authors implement the Hankel‑type contour proposed by Usciati et al. (2025). They compute the corresponding results for the partition function and the two‑point function, showing that the two prescriptions lead to different numerical prefactors and phases, thereby highlighting the ambiguity in defining timelike Liouville correlators beyond charge neutrality.

The paper situates its findings within the broader literature: it reproduces the scaling (\mu^{\frac{1}{b^2}}) expected from Liouville theory, contrasts its non‑zero result with the vanishing outcome of analytic continuation of the spacelike DOZZ formula, and notes the connection of the special coupling (b=1/\sqrt{2}) to the free‑fermion point in the massless sine‑Gordon model. While the determinantal simplification appears unique to (b^2=1/2), the authors discuss the possibility of similar integrable structures at other rational values of (b^2), as suggested by recent bootstrap results, but leave this as an open problem.

In summary, the work provides the first rigorously controlled family of exact calculations in timelike Liouville theory outside the charge‑neutral sector. By exploiting a special coupling that renders the Coulomb‑gas integrals determinantal, the authors obtain closed‑form expressions for several low‑point functions, resolve the zero‑mode ambiguity via explicit regularizations, and compare competing contour prescriptions. These results open new avenues for studying non‑neutral timelike Liouville correlators, modular and fusion kernels at rational couplings, and potential supersymmetric extensions.