---
title: Parabolic Airy Process
url: https://www.emergentmind.com/topics/parabolic-airy-process
type: topic
---

# Parabolic Airy Process

The parabolic Airy process is the process
\[
A(\alpha)=A_2(\alpha)-\alpha^2,\qquad \alpha\in\mathbb R,
\]
where \(A_2\) is the Airy\(_2\) process. It is a canonical KPZ scaling object whose finite-dimensional distributions are given by a Fredholm determinant with the extended Airy kernel, and it is the one-time spatial marginal of the KPZ fixed point with narrow wedge initial condition. In the multiline formulation, it is the top curve of the parabolic Airy line ensemble, obtained by subtracting the deterministic parabola from the stationary Airy line ensemble [2507.20487] [2302.00097].

## 1. Definition and exact finite-dimensional law

A standard exact definition takes the parabolic Airy process to be \(A(\alpha)=A_2(\alpha)-\alpha^2\). For ordered points \(\alpha_1<\cdots<\alpha_m\) and thresholds \(\beta_1,\dots,\beta_m\),
\[
P\left(\bigcap_{i=1}^m \{ A (\alpha_i) \le \beta_i \}\right)
=
P\left(\bigcap_{i=1}^m \{ A_2 (\alpha_i) \le \beta_i+\alpha_i^2 \}\right)
=
\det\!\left(I-\chi K_{\mathrm{ext}}\chi\right),
\]
with cutoff \(\chi(\alpha_i,x)=\mathbf 1_{(\beta_i+\alpha_i^2,\infty)}(x)\) and extended Airy kernel
\[
K_{\mathrm{ext}}(\alpha_i,x;\alpha_j,y)
=
\begin{cases}
\displaystyle \int_0^\infty e^{-z(\alpha_i-\alpha_j)}\,\Ai(x+z)\Ai(y+z)\,dz, & \alpha_i\ge \alpha_j,\\[1.2em]
\displaystyle -\int_{-\infty}^0 e^{-z(\alpha_i-\alpha_j)}\,\Ai(x+z)\Ai(y+z)\,dz, & \alpha_i< \alpha_j.
\end{cases}
\]
Its one-point law is therefore
\[
P(A(\alpha)\le \beta)=F_{\mathrm{GUE}}(\beta+\alpha^2),
\]
so the parabola appears as a deterministic shift of the GUE Tracy–Widom variable [2507.20487].

This Fredholm-determinantal description is equivalent to the usual extended-kernel definition of the Airy\(_2\) process. A direct equivalence proof between the classical extended-Airy formula and an alternative equal-time formula inherited from KPZ fixed-point formulas was given in 2025. In that formulation, the same multipoint distribution is expressed as a contour Fredholm determinant \(\det(I+K)\), and the proof identifies \(\det(I+K)\) with \(\det(I-\chi K_{\mathrm{ext}}\chi)\) by explicit kernel manipulations, contour deformations, and a generalized Andreief identity [2507.20487].

A complementary boundary-value formulation replaces the extended kernel by the Airy Hamiltonian \(H=-\Delta+x\). For \(u_1<\dots<u_n\),
\[
\mathbb{P}\!\left(\mathcal A_2(u_1)\le x_1,\dots,\mathcal A_2(u_n)\le x_n\right)
=
\det\!\left(I-K_{\Ai}+\bar P_{x_1}e^{(t_1-t_2)H}\bar P_{x_2}\cdots \bar P_{x_n}e^{(t_n-t_1)H}K_{\Ai}\right),
\]
and this formulation is the basis for continuum barrier problems and parabolic-variational identities [1301.0750].

## 2. KPZ fixed point, narrow wedge scaling, and last-passage origins

In KPZ theory the parabolic Airy process arises from curved, or droplet, geometry. For the one-dimensional KPZ equation with sharp wedge initial data,
\[
\partial_t h(x,t)=\frac{\lambda}{2}\big(\partial_x h(x,t)\big)^2+\nu\,\partial_x^2 h(x,t)+\sqrt{D}\,\eta(x,t),
\]
the exact narrow-wedge analysis isolates a deterministic growth term, a deterministic spatial parabola, and a fluctuation field \(\xi(x,t)\). Under the scaling \(\alpha=(t/2)^{1/3}\), \(x_j=2\alpha^2 y_j\), \(s_j=\alpha a_j\), the centered field satisfies
\[
\lim_{t\to\infty}\Pr\big(\xi(2\alpha^2 y_1,t)<\alpha a_1,\ldots,\xi(2\alpha^2 y_n,t)<\alpha a_n\big)=F_n(a,y),
\]
where \(F_n\) is the finite-dimensional distribution of the Airy process. Because the deterministic term \(-x^2/(2t)\) has already been subtracted, the limit is \(\mathcal A_2(y)\); without that subtraction, the same statement is equivalently written as convergence to \(\mathcal A_2(y)-y^2\) [1101.4622].

At the level of the KPZ fixed point, the relation becomes exact and scale-covariant. If \(H(\alpha,\tau)\) denotes the KPZ fixed point with narrow wedge initial condition, then
\[
H(\alpha,\tau)\stackrel{d}{=}\tau^{1/3}A(\tau^{-2/3}\alpha),
\qquad
A(\alpha)\stackrel{d}{=}H(\alpha,1).
\]
Thus the parabolic Airy process is precisely the equal-time spatial field of the narrow-wedge KPZ fixed point [2507.20487].

In directed polymer and last-passage formulations, the same object is the curved-profile scaling limit. A canonical statement is
\[
H^{\rm point}_N(u)\xrightarrow[N\to\infty]{}\mathcal A_2(u)-u^2
\]
in the topology of uniform convergence on compact sets. This identifies \(\mathcal A_2-u^2\) as the genuine spatial fluctuation field in point-to-point geometry, rather than as an auxiliary recentering [1301.0750]. In exponential LPP coupling arguments, the narrow-wedge profile likewise converges to \(\mathcal A_2(u)-u^2\), and local information about the parabolically shifted process transfers to \(\mathcal A_2\) because the quadratic term is negligible under \(\varepsilon^{-1/2}\) zooming [1704.01903].

The replica-Bethe-ansatz derivation of KPZ multipoint statistics for sharp wedge also enters the history of the subject. For \(n\ge2\), the finite-time \(n\)-point generating function in the 2011 KPZ paper is obtained under a natural factorization assumption that is not exact at finite time, but its long-time limit reproduces the standard Airy-process Fredholm determinant formula [1101.4622]. This methodological caveat is important in tracing how exact Airy-process formulas emerged from KPZ calculations.

## 3. Parabolic Airy line ensembles, Brownian Gibbs structure, and Wiener densities

The parabolic Airy process is only the top curve of a richer multiline object. If \(\widetilde A\) denotes the stationary Airy line ensemble, the parabolic Airy line ensemble is
\[
A_i(t)=\widetilde A_i(t)-t^2,
\]
a random element of \(C^{\mathbb N}(\mathbb R)\) with almost sure ordering
\[
A_1>A_2>\cdots.
\]
The stationary ensemble is recovered by adding back \(t^2\), and \(A(\cdot)+(\cdot)^2\) is stationary under spatial shifts and time reversal. In this notation the parabolic Airy process is simply the top line \(A_1\) [2302.00097].

This parabolic shift is the normalization in which the Brownian Gibbs property takes its cleanest form. In Hammond’s diffusion-\(1\) convention,
\[
\mathcal L(i,x)=2^{-1/2}\big(\mathcal A(i,x)-x^2\big),
\]
and the factor \(2^{-1/2}\) is chosen so that the curves are locally Brownian with diffusion constant \(1\). The Brownian Gibbs property says that after conditioning on the ensemble outside a strip, the conditional law inside is that of independent Brownian bridges conditioned on mutual avoidance and on staying above the next lower curve. This is the structural reason the parabolically shifted ensemble, rather than the unshifted Airy line ensemble, is central in Gibbsian and polymeric formulations [1609.02971].

The most refined local comparison presently available is the Wiener-density analysis on compact spacetime patches \(K=\{1,\dots,k\}\times[a,a+t]\). Writing \(\nu_{k,t}\) for the law of the recentered top \(k\) curves and \(\mu_{k,t}\) for \(k\) independent Brownian motions of variance \(2\), one has a Radon–Nikodym derivative \(X_{k,t}=d\nu_{k,t}/d\mu_{k,t}\) satisfying
\[
X_{k,t}(f)=\exp\Big(-S(f)+O_k\big(t^3+\sqrt t\,S(f)^{5/6}\big)\Big),
\]
and, in the abstract form,
\[
X_K(f)=\exp\big(-\textsf{S}(f)+o(\textsf{S}(f))\big).
\]
The action is explicit: if \(Tf\) is the “Tetris” stacking map that minimally shifts a \(k\)-tuple into ordered position above \(0\) while preserving increments, then
\[
S(f)=\frac{2}{3}\Big(\sum_{i=1}^k [Tf_i(0)]^{3/2}+[Tf_i(t)]^{3/2}\Big).
\]
For the top curve \(k=1\), this reduces to the cost of the minimal upward shift needed to keep the path above the parabolic background. The same analysis yields a uniform bound
\[
X_{k,t}(f)\le e^{O_k(t^3)}
\]
for \(\mu_{k,t}\)-a.e. \(f\), sharp lower-tail control for the density, and a small-amplitude large deviation principle with Schilder rate function
\[
I(f)=\frac14\int_0^t \sum_{i=1}^k |f_i'(s)|^2\,ds
\]
on absolutely continuous paths [2302.00097].

Brownian regularity results complement these density formulas. After affine recentering on a fixed interval \([K,K+d]\), the \(k\)-th parabolically shifted Airy curve has Radon–Nikodym derivative \(f_k\) with respect to Brownian bridge such that
\[
j_k:=\int \exp\big\{\alpha_k(\log f_k)^{6/5}\big\}\,dB<\infty,
\]
and hence \(f_k\in L^p\) for every \(p<\infty\). This yields Brownian-bridge comparison estimates, Gaussian-type suprema bounds up to lower-order corrections, and Brownian \(1/2\)-Hölder modulus behavior for the top curve [1609.02971]. A separate LPP-coupling approach shows the local scaling limit
\[
\lim_{\varepsilon\downarrow 0}\varepsilon^{-1/2}\bigl(\mathcal A_2(\varepsilon x)-\mathcal A_2(0)\bigr)\stackrel{dist.}=2B(x),
\]
and the same local limit applies to \(\mathcal A_2(x)-x^2\), since the parabola contributes only \(O(\varepsilon^{3/2})\) under this scaling [1704.01903].

## 4. Variational identities, extremal laws, and deformations

A major part of the subject concerns functionals of the parabolic Airy process rather than its finite-dimensional distributions. The most classical identity is the full-line supremum law
\[
\mathbb P\!\left(\sup_{x\in\mathbb R}\{\mathcal A_2(x)-x^2\}\le m\right)=F_{\rm GOE}(4^{1/3}m),
\]
which identifies point-to-line curved KPZ geometry with GOE Tracy–Widom fluctuations. The same framework yields a continuum barrier formula
\[
\mathbb P\!\left(\mathcal A_2(t)\le g(t)\text{ for }t\in[\ell,r]\right)
=
\det\!\left(I-K_{\Ai}+\Theta^g_{[\ell,r]}e^{(r-\ell)H}K_{\Ai}\right),
\]
and for the parabolic barrier \(g(t)=t^2+m\) the Brownian-bridge representation of \(\Theta^g\) becomes explicitly tractable. The joint law of the maximizer
\[
\mathcal T=\argmax_{t\in\mathbb R}\{\mathcal A_2(t)-t^2\}
\]
and the maximum
\[
M=\max_{t\in\mathbb R}\{\mathcal A_2(t)-t^2\}
\]
also admits an explicit determinant formula [1301.0750].

Restricting the supremum to a half-line produces the Airy\(_{2\to1}\) crossover. For every \(\alpha,m\in\mathbb R\),
\[
\mathbb P\!\left(\sup_{t\le \alpha}\big(\mathcal A_2(t)-t^2\big)\le m-\min\{0,\alpha\}^2\right)=G^{2\to1}_\alpha(m),
\]
where \(G^{2\to1}_\alpha(m)=\mathbb P(\mathcal A_{2\to1}(\alpha)\le m)\). As \(\alpha\to-\infty\) this family tends to \(F_{\rm GUE}\), and as \(\alpha\to+\infty\) it tends to \(F_{\rm GOE}(4^{1/3}m)\); for each fixed \(\alpha\), its right tail has the same \(\exp(-(4/3)m^{3/2})\) decay rate at the level of the paper’s upper bound [1111.2565].

Finite-rank deformations of the parabolic Airy functional arise from the Airy process with wanderers. The central law is
\[
F^{(\vec a,\vec b)}_m(r)
=
\mathbb P\!\left(\sup_{t\in\mathbb R}\big(A^{(\vec a,\vec b)}_m(t)-t^2\big)\le r\right),
\]
which reduces to \(F_1(2^{2/3}r)\) when \(m=0\). This law admits equivalent descriptions by a Fredholm determinant, a Painlevé II/Lax-pair formula, a KPZ fixed-point marginal, a Bloemendal–Virág-type PDE, and a KdV characterization. It also appears as the scaling limit of the maximal height of non-intersecting Brownian bridges with finitely many outlier endpoints [2009.07781].

Stationary KPZ introduces a Brownian deformation of the parabolic Airy variational problem. The paper on two-time KPZ memory studies
\[
\max_{\hat y}\left({\cal A}_2(\hat y)-\hat y^2+\sqrt{2}\,B(\hat y)\right),
\]
identifying it with \({\cal A}_{\rm stat}(0)\), and more generally
\[
{\cal A}_{\rm stat}(\hat x)
=
\max_{\hat y}\left({\cal A}_2(\hat y-\hat x)-(\hat y-\hat x)^2+\sqrt{2}B(\hat y)\right)
\]
at the one-point level. The corresponding maximum has the Baik–Rains or extended Baik–Rains distribution, and the argmax distribution is exactly the KPZ scaling function \(f_{\rm KPZ}\). The same paper derives a joint law with \({\cal A}_2(0)\) and, from it, the large-time-separation KPZ correlation constant
\[
C_\infty \approx 0.623.
\]
Its derivation uses replica Bethe ansatz together with a decoupling assumption in the large-time limit, a qualification explicitly stated in the paper [1709.06264].

## 5. Coupled parabolic Airy processes and the Airy sheet

The parabolic Airy process acquires a genuinely multipoint geometry when coupled through the Airy sheet. In Brownian last-passage scaling, the time-one Airy sheet \(W\) is the limit of the centered weights \(W_n(x,y)\), and for fixed \(x\) the parabolically adjusted profile
\[
z\mapsto W(x,z)+Q(z-x),
\qquad
Q(z)=2^{-1/2}z^2,
\]
is an Airy\(_2\)-type process in the paper’s normalization. Equivalently, \(W(x,\cdot)\) itself is an Airy\(_2\) profile with the parabola already subtracted [1904.01717].

The simplest nontrivial coupling is the difference
\[
Z(z)=W(1,z)-W(-1,z).
\]
This process is almost surely continuous and non-decreasing, and it is constant in a random neighbourhood of almost every \(z\in\mathbb R\). The exceptional set where it fails to be locally constant has Hausdorff dimension \(1/2\). Geometrically, this set corresponds to endpoints for which the maximizing paths from \((-1,0)\) and \((1,0)\) to \((z,1)\) fail to have already coalesced before reaching the endpoint [1904.01717].

A sharper description uses the parabolic Airy sheet \(\mathcal S\) and the difference profile
\[
\mathcal D(x)=\mathcal S(1,x)-\mathcal S(-1,x).
\]
On every compact interval \([c,d]\), the increment process
\[
\mathcal D(\cdot)-\mathcal D(c)\big|_{[c,d]}
\]
is absolutely continuous, in law, with respect to Brownian local time at zero for Brownian motion of rate four, again in the sense of increments. At random points of increase—such as the first point of increase after a deterministic location, or a point sampled from the Stieltjes measure induced by \(\mathcal D\)—one has the local limit
\[
\varepsilon^{-1/2}\bigl(\mathcal D(\tau+\varepsilon t)-\mathcal D(\tau)\bigr)\Longrightarrow L(t),
\]
where \(L\) is Brownian local time of rate four. This identifies the singular monotone coupling between two parabolic Airy\(_2\) profiles with Brownian local time, not merely at the level of Hausdorff dimension but at the level of tangent-process structure [2103.12029].

These results suggest a precise geometric principle: the coupling of parabolic Airy processes through the Airy sheet is governed by geodesic coalescence, and the residual increments of the difference profile arise from local competition between the top two Airy lines. In the local theorem this is made explicit through a running-maximum representation involving \(P_2-P_1\), where \(P\) is the parabolic Airy line ensemble [2103.12029].

## 6. Terminology, model variants, and boundary analogues

The phrase “parabolic Airy process” is not completely uniform across the literature. Some papers define the process directly as \(A_2(\alpha)-\alpha^2\), while others say “Airy process” for \(A_2\) after the finite-time deterministic parabola has already been subtracted from the KPZ height field. In particular, the 2011 KPZ paper studies the already centered field \(\xi(x,t)\), so its scaling limit is written as \(\mathcal A_2\); without that centering, the same statement becomes convergence to \(\mathcal A_2(y)-y^2\) [1101.4622]. Likewise, Hammond’s terminology is “Airy line ensemble after the subtraction of a parabola” or “parabolically shifted ensemble,” with the preferred Brownian normalization
\[
2^{-1/2}\big(\mathcal A_2(x)-x^2\big)
\]
for the top curve [1609.02971].

The same parabolic shift appears outside KPZ growth and LPP in other edge-scaling problems. In the two-periodic Aztec diamond with a thin mesoscopic rough region, the last path at the rough–smooth boundary satisfies
\[
\mathbb P\big(\Gamma(t_1)\le \xi_1,\dots,\Gamma(t_j)\le \xi_j\big)
\longrightarrow
\mathbb P\big(A(t_1)\le \xi_1+t_1^2,\dots,A(t_j)\le \xi_j+t_j^2\big),
\]
equivalently \(\Gamma(t)\Rightarrow A(t)-t^2\) in finite-dimensional distributions, under the scaling \(a\sim n^{-1+\gamma}\) with \(\gamma\in(0,1/3)\) [2302.04663].

At the level of random-matrix cusp-to-edge crossover, the Pearcey process converges to the unshifted Airy process after centering at the deterministic edge location
\[
\frac{2}{27}(3\tau)^{3/2}
\]
and scaling by \((3\tau)^{1/6}\). That centering is distinct from subtracting \(t^2\), but after matching conventions it leads back to the same Airy\(_2\)/parabolic-Airy dichotomy familiar in KPZ theory [1009.0683].

A boundary analogue now exists in half-space KPZ. The pinned half-space Airy line ensemble \(\mathcal A^{\mathrm{hs};\infty}\) on \([0,\infty)\) is the \(\varpi\to\infty\) limit of critical half-space Airy line ensembles, and after parabolic shift
\[
\sqrt{2}\cdot L_i(t)+t^2=A_i(t)
\]
it satisfies a one-sided Brownian Gibbs property described by pairwise pinned Brownian motions. Far from the origin it converges to the standard Airy line ensemble, while at the origin its law is that of the ordered eigenvalues of the stochastic Airy operator at \(\beta=4\), each with doubled multiplicity. This places the parabolic Airy formalism within a broader family of boundary-sensitive KPZ scaling objects [2601.04546].

The modern picture is therefore two-layered. As a single process, the parabolic Airy process is \(A_2-\alpha^2\), with exact multipoint laws, variational identities, and KPZ fixed-point interpretations. As part of a line or sheet ensemble, it is the top curve of a more rigid Brownian-Gibbs object whose multiline interactions govern endpoint repulsion, geodesic coalescence, and boundary effects [2507.20487] [2302.00097].

Source: https://www.emergentmind.com/topics/parabolic-airy-process