---
title: Toeplitz Determinants with FH Singularities
url: https://www.emergentmind.com/topics/toeplitz-determinants-with-fisher-hartwig-singularities
type: topic
---

# Toeplitz Determinants with FH Singularities

A Toeplitz determinant with Fisher-Hartwig (FH) singularities is a determinant of a Toeplitz matrix whose underlying symbol exhibits root-type and/or jump-type singularities at a finite set of points on the unit circle. These structures arise naturally in statistical mechanics, random matrix theory, fermionic physics, and integrable probability, where the determinants encode correlation functions, gap probabilities, or multiplicative statistics. The Fisher-Hartwig conjecture (and subsequent theorem) provides their precise asymptotic expansion as the matrix size grows, capturing leading powers, amplitudes, oscillatory terms, and the subtle effects of merging or interacting singularities.

## 1. Fisher-Hartwig Symbols and Toeplitz Determinants

A Toeplitz matrix of dimension $n$ is specified by a function (“symbol”) $f(z)$ defined on the unit circle $|z|=1$, with entries $[f_{j-k}]_{0\leq j,k\leq n-1}$, where $f_{m}$ is the $m$-th Fourier coefficient of $f$. The Fisher-Hartwig symbol is a product of a smooth (nonvanishing, analytic) factor and finitely many singular terms located at points $z_j=e^{i\varphi_j}$:
\[
f(e^{i\theta})=e^{V(e^{i\theta})} \prod_{j=1}^m |e^{i\theta}-z_j|^{2\alpha_j} e^{i\beta_j \,\mathrm{sgn}(\theta-\varphi_j)}
\]
where $V(z)$ is real-analytic near $|z|=1$, and $(\alpha_j, \beta_j)$ are the root and jump exponents at $z_j$. The Toeplitz determinant is defined as $D_n[f] = \det(f_{j-k})_{0\leq j,k\leq n-1}$ or equivalently via the Heine integral as a partition function with log-gas interactions on the circle [2512.09123][1108.1355].

## 2. Asymptotic Formulas: The Fisher-Hartwig Expansion

For large $n$, the Fisher-Hartwig formula provides the leading behavior of $D_n[f]$:
\[
D_n[f] = E[f]^n\, n^{\sum_{j=1}^m (\alpha_j^2 - \beta_j^2)}\, \prod_{j=1}^m \frac{G(1+\alpha_j+\beta_j)G(1+\alpha_j-\beta_j)}{G(1+2\alpha_j)} \prod_{1\leq j<k\leq m} |z_j-z_k|^{2(\beta_j\beta_k - \alpha_j\alpha_k)} (1+o(1))
\]
where $E[f]=\exp(V_{0})$ with $V_0$ the zero Fourier mode of $V$ and $G$ is the Barnes G-function [2512.09123][1108.1355][1909.07362]. The powers of $n$ encode local singular behavior at each $z_j$. Further corrections involve oscillatory and polynomial-in-$1/n$ terms arising from “branches” in the FH representations, corresponding to shifts in the integer part of the jumps:
\[
D_n[f] \sim \sum_{\{k_j\}} \exp[n V_0 + i n\sum_j k_j \varphi_j - \sum_j (\beta_j - k_j)^2 \ln n]\, \mathcal{A}(\{\alpha_j, \beta_j; k_j\}) \exp\left(\sum_{\ell=1}^{\infty} \frac{c_\ell(\{\alpha_j,\beta_j; k_j\})}{n^{\ell}}\right)
\]
The sum is over all relevant “branches” (integer shifts $k_j$ such that $\sum_j k_j = 0$) [1112.2530][1309.4867]. At special values of the exponents, competing branches contribute equally, producing so-called “switching” and crossover behavior [1112.2530]. When singularities merge, the asymptotics are given by a uniform expression involving Painlevé V transcendent [1004.3696][1403.3639][2311.06616].

## 3. Merging Singularities and Painlevé Transitions

As two Fisher-Hartwig singularities at $z_1, z_2$ coalesce ($z_1\to z_2$), the simple product formula is no longer sufficient. Instead, the determinant exhibits a crossover, governed by a universal Painlevé V equation (“double-scaling” coalescence regime):
\[
\log D_n(f_t) = \log D_n(f_0) + i n t (\beta_2 - \beta_1) + i\int_{0}^{-2i n t} [\sigma(s) - c ] \frac{ds}{s} + O(1)
\]
Here $f_t$ is the symbol with two singularities at separation $2t$, $f_0$ the merged (single singularity) symbol, and $\sigma(s)$ is a solution to the Jimbo–Miwa–Okamoto $\sigma$-form of Painlevé V with explicit exponents determined by $\alpha_j$, $\beta_j$ [1403.3639][1004.3696][2311.06616]. This formula interpolates between regimes of two isolated singularities ($t>0$ fixed, $n\to\infty$) and the merged-case ($t\to0$), capturing the relevant logarithmic corrections and oscillatory behavior.

## 4. Riemann–Hilbert Approach and Subleading Expansions

The technical machinery underpinning the FH asymptotics exploits the Riemann–Hilbert (RH) problem for orthogonal polynomials on the unit circle (OPUC) with respect to the symbol $f$. The Deift–Zhou nonlinear steepest-descent method is applied to construct global (“outer”) and local (“inner,” singular) parametrices that asymptotically solve the RH problem [1206.1292][1909.07362][2512.09123]. Local parametrices near singularities are built from Bessel, confluent hypergeometric, or Painlevé functions, depending on the nature of the singularity and the scaling regime. Matching and recursions yield not only the leading order but also subleading corrections (in inverse powers of $n$), with polynomial or oscillatory structure; periodicity of the amplitudes under shifts of the jump exponents has been established [1112.2530][1309.4867].

## 5. Physical and Probabilistic Applications

Toeplitz determinants with FH singularities are central in quantum many-body physics and random matrix theory. In non-equilibrium Fermi-edge physics and Luttinger liquid models, multiple-step distribution functions yield Toeplitz symbols with multiple jump-type singularities; the resulting determinants encode the asymptotics of Green functions, tunneling currents, and full counting statistics, with exponents and oscillations dictated by the FH parameters and branch structure [1203.6418][1010.5645][1212.0708]. In random matrix theory, the moments of characteristic polynomials (and their distributions) are written as Toeplitz determinants with FH singularities; in the scaling limits they interpolate between regimes of Gaussian multiplicative chaos and classical log-correlated fields [2311.06616][1509.03446].

In gap probabilities for CUE, symbols with jump/root singularities describe hard-edge or soft-edge universality, with the exact gap probabilities expressed via Fredholm determinants of hypergeometric or sine kernels, and their asymptotics reduced to a Painlevé or RHP analysis [1907.11509][1004.3696]. In the Ising model, the diagonal susceptibility and spin-spin correlation functions are governed by sums over Toeplitz determinants with FH symbols deformed by external parameters, exhibiting phase transitions and singularities tracked by the FH structure [1502.04922][1004.3696].

## 6. Discrete and Non-Hermitian Generalizations

In discrete settings (e.g., log-gas on the roots of unity, characteristic polynomials of random permutation matrices), discrete Toeplitz determinants with FH singularities exhibit the same asymptotic formulas as in the continuum as long as the system is sufficiently dilute; deviations and additional correction factors arise for dense occupation [1509.03446]. The Fisher-Hartwig asymptotics have also been generalized to two-dimensional settings, notably for non-Hermitian random matrix ensembles, where the partition function with boundary singularities is governed by analogous, but spatially extended, expressions and universality classes [2512.09123].

## 7. Summary Table: Fisher-Hartwig Asymptotics

| Context                      | Symbol Structure                     | Leading Asymptotic Behavior                                          |
|------------------------------|--------------------------------------|----------------------------------------------------------------------|
| Isolated singularities       | $f(z) = e^{V(z)}\prod |z-z_j|^{2\alpha_j} e^{i\beta_j\arg(z-z_j)}$     | $E[f]^n\, n^{\Sigma_j(\alpha_j^2-\beta_j^2)} \times$ $G$-factors [2512.09123]             |
| Multiple-step/jump scenario  | Discontinuities in $n(\epsilon)$     | Sum over FH branches; exponents determined by $\beta_j$ [1203.6418]  |
| Merging singularities        | $f_t$ with $z_1\!\to\! z_2$          | Uniform expansion involving Painlevé V $\sigma$-transcendent [1403.3639]|
| Discrete log-gas             | $f$ sampled at $M$-th roots of unity | FH formula holds for $N\ll M$, corrections otherwise [1509.03446]     |
| Non-Hermitian ensembles      | Edge singularities on $|z|=1$        | 2d analog of FH; partition function asymptotics [2512.09123]            |

This synthesis of the Fisher-Hartwig theory and its extensions reveals the central role of Toeplitz determinants with singular symbols in capturing universal scaling, transition phenomena, and fine structure across mathematical physics, with methods and results grounded in Riemann–Hilbert analysis, asymptotic expansions including both universal constants and highly nontrivial oscillatory and merging corrections, and rich connections to integrable systems and random matrices.

Source: https://www.emergentmind.com/topics/toeplitz-determinants-with-fisher-hartwig-singularities