Hostname: page-component-7bb8b95d7b-5mhkq Total loading time: 0 Render date: 2024-10-04T22:34:40.606Z Has data issue: false hasContentIssue false

KP governs random growth off a 1-dimensional substrate

Published online by Cambridge University Press:  21 April 2022

Jeremy Quastel*
Affiliation:
Department of Mathematics, University of Toronto, 40 St. George Street, Toronto, Ontario, Canada M5S 2E4
Daniel Remenik
Affiliation:
Departamento de Ingeniería Matemática and Centro de Modelamiento Matemático (UMI-CNRS 2807), Universidad de Chile, Av. Beauchef 851, Torre Norte, Piso 5, Santiago, Chile; E-mail: dremenik@dim.uchile.cl

Abstract

The logarithmic derivative of the marginal distributions of randomly fluctuating interfaces in one dimension on a large scale evolve according to the Kadomtsev–Petviashvili (KP) equation. This is derived algebraically from a Fredholm determinant obtained in [MQR17] for the Kardar–Parisi–Zhang (KPZ) fixed point as the limit of the transition probabilities of TASEP, a special solvable model in the KPZ universality class. The Tracy–Widom distributions appear as special self-similar solutions of the KP and Korteweg–de Vries equations. In addition, it is noted that several known exact solutions of the KPZ equation also solve the KP equation.

Type
Probability
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press

1 Matrix Kadomtsev–Petviashvili equation for multidimensional distributions

The $1$ -dimensional Kardar–Parisi–Zhang (KPZ) universality class consists of random growth models, last-passage percolation and directed polymers and random stirred fluids. All models in the class have an analogue of the height function $h(t,x)$ (free energy, integrated velocity), whose long-time large-scale evolution is the principal object of study. The name of the class comes from the KPZ equation,

(1.1) $$ \begin{align} \partial_t h = \lambda(\partial_xh)^2 + \nu \partial_x^2h + \sigma \xi \end{align} $$

(with $\xi $ a space-time white noise), a canonical continuum equation for random growth introduced in [Reference Kardar, Parisi and ZhangKPZ86]. However, the main interest is in the universal features which are found only at large space-time scales, under the 1:2:3 scaling corresponding to $\varepsilon \to 0$ in

(1.2) $$ \begin{align} \varepsilon^{1/2} h\left(\varepsilon^{-3/2} t, \varepsilon^{-1} x\right) - C_\varepsilon t. \end{align} $$

The KPZ equation is not invariant under this scaling, which sends $(\lambda ,\nu ,\sigma )$ to $\left (\lambda ,\varepsilon ^{1/2}\nu ,\varepsilon ^{1/4}\sigma \right )$ . A key problem is to find the true, scaling-invariant equation for random interface growth.

Since the early 2000s [Reference JohanssonJoh03; Reference SasamotoSas05; Reference Borodin, Ferrari, Prähofer and SasamotoBor+07] it has been known, for a number of models in the class and special scaling-invariant initial data narrow wedge and flat, that the distributional limits of formula (1.2) are the Tracy–Widom distributions of random matrix theory. In an earlier article [Reference Matetski, Quastel and RemenikMQR17] it was shown that, at least for one model in the class, TASEP, (1.2) converges to a Markov process $\mathfrak {h}(t,x)$ which is invariant under 1:2:3 scaling:

(1.3) $$ \begin{align} \alpha\mathfrak{h}\left(\alpha^{-3}t,\alpha^{-2}x;\alpha^{-1}\mathfrak{h}_0\left(\alpha^2x\right)\right)\stackrel{\text{dist}}{=}\mathfrak{h}(t,x;\mathfrak{h}_0), \end{align} $$

where $\mathfrak {h}_0$ after the semicolon denotes the initial data. It is widely believed that this KPZ fixed point governs the limiting fluctuation for all models in the class.

The KPZ fixed point does not satisfy a stochastic differential equation. In place of that, it inherits a variational formulation from TASEP, a Hopf–Lax-type formula involving a nontrivial input noise called the Airy sheet $\mathcal {A}(x,y)$ : for the KPZ fixed point starting from $\mathfrak {h}(0,x)=\mathfrak {h}_0(x)$ ,

(1.4) $$ \begin{align} \mathfrak{h}(t,x) \stackrel{\text{dist}}{=} \mathop{\mathrm{sup}}\limits_{y\in\mathbb{R}}\left\{ t^{1/3}\mathcal{A}\left(t^{-2/3} x,t^{-2/3} y\right)- \tfrac1{t}(x-y)^2 + \mathfrak{h}_0(y)\right\}. \end{align} $$

The Airy sheet $\mathcal {A}(x,y)$ can be thought of as the height function at x at time $1$ , starting from a narrow wedge at y at time $0$ , and therefore involves coupling different initial conditions. As far as we know at the present time, the problem of the coupled initial conditions is not integrable, and therefore the distribution of the Airy sheet is unknown. This led to a problem, that it was unclear that equation (1.4) even involved a unique object on the right-hand side. An important advance is in [Reference Dauvergne, Ortmann and VirágDOV19], where it is shown that the Airy sheet is a functional of the Airy line ensemble. This puts the variational formula (1.4) on solid footing, as it obviates the need for uniqueness of the Airy sheet. However, the functional is nonexplicit. In this sense, equation (1.4) is not satisfying as a universal scaling-invariant equation.

Instead of a universal stochastic equation, one can study the n-space point distribution functions

(1.5) $$ \begin{align} F(t,x_1,\dotsc,x_n,r_1,\dotsc, r_n) = \mathbb{P}_{\mathfrak{h}_0}\left(\mathfrak{h}(t,x_1)\leq r_1,\dotsc,\mathfrak{h}(t,x_n)\leq r_n\right), \end{align} $$

where the subscript in $\mathbb {P}_{\mathfrak {h}_0}$ denotes the initial data of the KPZ fixed point. In the cases of narrow wedge and flat initial data, it is known [Reference JohanssonJoh00; Reference SasamotoSas05; Reference Borodin, Ferrari, Prähofer and SasamotoBor+07] that the $1$ -dimensional distributions $F(1, x,r)$ are, respectively, the Tracy–Widom Gaussian unitary ensemble (GUE) and Gaussian orthogonal ensemble (GOE) random matrix distributions (but except in the particular case of narrow wedge initial data, the connection between random growth and random matrices has remained tangential and murky). The multidimensional distributions in these cases are given by Fredholm determinants and define the Airy $_2$ and Airy $_1$ processes. The $1$ -dimensional distributions can be written in terms of the Hastings–McLeod solution of Painlevé II; a longstanding open question is whether the distributions satisfy an equation in the more general setting.

1.1 Main results

In [Reference Matetski, Quastel and RemenikMQR17], it is shown that for initial data $\mathfrak {h}_0$ in $\operatorname {\mathrm {\text {UC}}}$ , a space of upper semicontinuous functions with a linear growth condition (see equation (1.24)), the n-space point distribution functions from equation (1.5) are given by Fredholm determinants

(1.6) $$ \begin{align} F(t,x_1,\dotsc,x_n,r_1,\dotsc, r_n) = \det(\mathbf{I} - \mathbf{K}), \end{align} $$

where $\mathbf {K} = \mathbf {K}( t, x_1,\dotsc ,x_n,r_1,\dotsc , r_n,\mathfrak {h}_0)$ is an operator on the n-fold direct sum of $L^2(\mathbb {R}^+)$ , which is given by an explicit $n\times n$ matrix kernel $\mathbf {K}_{ij}(u_1,u_2)=\mathbf {K}_{ij}(t, x_1,\dotsc ,x_n,r_1,\dotsc , r_n,\mathfrak {h}_0,u_1,u_2)$ , $i,j=1,\dotsc ,n$ . The determinant is nonzero, so we can define $\mathbf {R}= (\mathbf {I}-\mathbf {K})^{-1}$ . Furthermore, $\mathbf {K}\mathbf {R}= \mathbf {R}\mathbf {K}=\mathbf {I}-\mathbf {R}$ , and we can let

$$ \begin{align*} Q=\mathbf{R}\mathbf{K}(0,0), \end{align*} $$

which is an $n\times n$ matrix-valued function of $t,x_1,\dotsc ,x_n,r_1,\dotsc , r_n$ and the initial height profile $\mathfrak {h}_0$ .

Let

$$ \begin{align*} \mathcal{D}_r= \partial_{r_1}+\dotsb + \partial_{r_n}, \qquad \mathcal{D}_x= \partial_{x_1}+\dotsb + \partial_{x_n}. \end{align*} $$

In Section 4.1 we will show that

(1.7) $$ \begin{align} \mathcal{D}_r\log F = \operatorname{{\mathop{\mathrm{{tr}}}}} Q. \end{align} $$

Our main result is the following:

Theorem 1.1. For any $\mathfrak {h}_0\in \operatorname {\mathrm {\text {UC}}}$ , Q and its derivative $q= \mathcal {D}_rQ$ solve the matrix Kadomtsev–Petviashvili (KP) equation

(1.8) $$ \begin{align} \partial_tq+\tfrac12\hspace{0.05em}\mathcal{D}_rq^2+\tfrac1{12}\mathcal{D}_r^3q+\tfrac14\mathcal{D}_x^2Q+\tfrac12[q,\hspace{0.05em}\mathcal{D}_xQ]=0, \end{align} $$

where $[A,B] = AB-BA$ . In particular, for the one-point marginals of the KPZ fixed point (i.e., equation (1.5) in the case $n=1$ ), $\phi =\partial _r^2 \log F$ satisfies the scalar KP-II equation

(1.9) $$ \begin{align} \partial_t \phi + \tfrac12\partial_r \phi^2 + \tfrac1{12}\partial_r^3 \phi + \tfrac14\partial_r^{-1} \partial_x^2 \phi = 0. \end{align} $$

The KP equation (1.9) was originally derived from studies of long waves in shallow water [Reference Ablowitz and SegurAS79]. It has come to be accepted as the natural $2$ -dimensional extension of the Korteweg–de Vries (KdV) equation; when $\phi $ is independent of x, corresponding in our case to flat initial data, it reduces to the KdV equation,

(1.10) $$ \begin{align} \partial_t \phi + \tfrac12\partial_r \phi^2+ \tfrac1{12}\partial_r^3 \phi=0. \end{align} $$

The KP equation is completely integrable and plays an important role in the Sato theory as the first equation in the KP hierarchy [Reference Miwa, Jimbo and DateMJD00]. The matrix KP equation (1.8) exists in the literature (see, e.g., [Reference KonopelchenkoKon82; Reference SakhnovichSak03]). None of the previous physical derivations of the KP equation seem to be related to the problem at hand, and it could well be that our evolution is through a class of functions where the equation is formally the same, because of similarities in the weakly nonlinear asymptotics, but the physics is completely different. Note that the equations are usually written with with coefficients $3$ , $1$ and $3$ replacing our $\tfrac 12$ , $\tfrac 1{12}$ and $\tfrac 14$ , which is achieved by $t\mapsto 12t$ and $\phi \mapsto \tfrac 12 \phi $ .

Remark 1.2.

  1. 1. The fact that random interface growth is governed by the KP equation was not anticipated.Footnote 1 We do not have a physical intuition why it is true; it follows by, essentially, algebra from the form of the kernel in the Fredholm determinant for equation (1.5), and we believe it is the first example of a physical law having been obtained in such a fashion. In retrospect, in the scalar case it was known that the evolution equations for the kernel (Theorem 1.3(1)–(3)) lead to equation (1.9), a fact that seems to have been rediscovered many times, particularly by Pöppe [Reference PöppePo89] (see also [Reference Zaharov and ŠabatZS74; Reference PöppePo84; Reference Pöppe and SattingerPS88; Reference McKeanMcK11]), before we rediscovered it. The fact that the kernel evolves by Theorem 1.3(1)–(3), and the importance of this fact, was not recognised earlier and is one of the key contributions of this article.

  2. 2. There are not so many natural partial differential equations with the necessary (in view of equation (1.3)) invariance under

    $$ \begin{align*} \phi(t,x,r) \mapsto \alpha^{-2} \phi\left(\alpha^{-3} t,\alpha^{-2} x,\alpha^{-1} r\right),\qquad \mathfrak{h}_0(x)\mapsto \alpha^{-1} \mathfrak{h}_0\left(\alpha^2 x\right). \end{align*} $$
    So in retrospect, once one knows that the finite-dimensional marginals come from a closed equation and that special examples are connected to Painlevé transcendents, one might expect to look for integrable equations such as the KP equation. The question is why the finite-dimensional marginals should come from a closed equation at all.
  3. 3. The Lax pair formulation of equation (1.9) (and also equation (1.8), with the matrix q replacing $\phi $ ) is

    (1.11) $$ \begin{align} \partial_t L = [L,A], \qquad L= \partial_x + \partial^2_{r} + 2\phi. \end{align} $$
    Here $A= \tfrac 13\partial _r^3 + \tfrac 23 \phi \partial _r + \tfrac 12 \partial _r\phi - \tfrac 12 \partial _r^{-1} \partial _x \phi $ , which tells us that the ‘spectrum’ of L is conserved. A very interesting question (a version of which A. Borodin asked us) is to make this precise and understand its physical meaning.
  4. 4. The $1$ -dimensional distribution functions themselves therefore satisfy the equivalent Hirota bilinear equation,

    $$ \begin{align*} F\partial^2_{tr} F -\partial_t F\partial_r F+ \tfrac1{12}F\partial^4_r F - \tfrac13\partial_r F\partial^3_r F+ \tfrac14\left(\partial_r^2 F\right)^2+\tfrac14F\partial^2_x F-\tfrac14(\partial_x F)^2=0, \end{align*} $$
    which again has the necessary 1:2:3 invariance, now under
    $$ \begin{align*} F(t,x,r) \mapsto F\left(\alpha^{-3} t,\alpha^{-2} x,\alpha^{-1} r\right),\qquad \mathfrak{h}_0(x)\mapsto \alpha^{-1} \mathfrak{h}_0\left(\alpha^2 x\right). \end{align*} $$
  5. 5. Unlike other limit points for fluctuation universality classes in probability, the Tracy–Widom distributions themselves lack any invariance. Theorem 1.1 recovers the invariance of the scaling limit under the 1:2:3 scaling, and the Tracy–Widom distributions then appear in the context of the KPZ universality class as special self-similar solutions of the KP equation (see Section 2.1).

Theorem 1.1 follows from the form of the kernel in the determinantal formula (1.6) and the following result. In the scalar case this appeared earlier in [Reference PöppePo89].

Theorem 1.3. Let $\mathbf {K}=\mathbf {K}(t,x,r)$ be an operator on the n-fold direct sum of $L^2\left (\mathbb {R}_{\ge 0}\right )$ which is trace class uniformly in compact sets of $t>0$ , $x,r\in \mathbb {R}$ , with matrix kernel $\mathbf {K}(u,v)=\mathbf {K}(t,x,r,u,v))$ which satisfies the following:

  1. (1) $\partial _r\mathbf {K}(u,v)=(\partial _{u}+\partial _{v})\mathbf {K}(u,v)$ ,

  2. (2) $\partial _t\mathbf {K}(u,v)=-\tfrac 13\left (\partial _{u}^3+\partial _{v}^3\right )\mathbf {K}(u,v)$ ,

  3. (3) $\partial _x\mathbf {K}(u,v)=\left (\partial _{v}^2-\partial _{u}^2\right )\mathbf {K}(u,v)$ .

Suppose in addition that $\det (\mathbf {I}-\mathbf {K})>0$ for all finite $t,x,r$ ; that $\mathbf {K}$ is real analytic in t and in each $x_i$ and $r_i$ ; and that the trace norm $\lVert \mathbf {K}\rVert _1<1$ for r in some open real interval. Then the $n\times n$ matrix $Q=(\mathbf {I}-\mathbf {K})^{-1}\mathbf {K}(0,0)$ and its derivative $q= \partial _rQ$ satisfy the matrix KP equation

(1.12) $$ \begin{align} \partial_tq+\tfrac12\hspace{0.05em}\partial_rq^2+\tfrac1{12}\partial_r^3q+\tfrac14\partial_x^2Q+\tfrac12[q,\hspace{0.05em}\partial_xQ]=0.\end{align} $$

The kernel in equation (1.6) is given explicitly (see [Reference Matetski, Quastel and RemenikMQR17, Eqn. (4.2)]) as

(1.13) $$ \begin{align} &\mathbf{K}_{ij}(t, x_1,\dotsc,x_n,r_1,\dotsc, r_n,\mathfrak{h}_0, u,v)= -e^{\left(x_j-x_i\right)\partial^2}\left(u+r_i,v+r_j\right)\mathbf{1}_{x_i<x_j}\nonumber\\ &~+ \int_{s\in\mathbb{R}^+\atop z,b\in\mathbb{R} } p_-(z,ds,db) \mathbf{S}_{t,-x_i-s} (b,u+r_i)\mathbf{S}_{t,x_j}\left(z,v+r_j\right) \nonumber\\ &~ + \int_{s\in\mathbb{R}^+\atop z,b\in\mathbb{R} } p_+(z,ds,db) \mathbf{S}_{t,-x_i} (z,u+r_i) \mathbf{S}_{t,x_j-s} \left(b, v+r_j\right)\nonumber\\ &~ - \int_{s_-,s_+\in\mathbb{R}^+\atop z,b_-,b_+\in\mathbb{R} } p_-(z,ds_-,db_-)p_+(z,ds_+,db_+) \mathbf{S}_{t,-x_i-s_-} (b_-,u+r_i) \mathbf{S}_{t,x_j-s_+} \left(b_+, v+r_j\right), \end{align} $$

where $p_\pm (z,ds,db) = p_z\left (\tau _{\mathfrak {h}_0^\pm }\in ds, B\left (\tau _{\mathfrak {h}_0^\pm }\right )\in db\right )$ are hitting measures of the hypographs of $\mathfrak {h}_0^-(x):=\mathfrak {h}_0(-x)$ and $\mathfrak {h}_0^+(x):=\mathfrak {h}_0(x)$ , $x\ge 0$ , by a Brownian motion B with diffusion coefficient $2$ , and where $\mathbf {S}_{t,x}(z_1,z_2)=\mathbf {S}_{t,x}(z_1-z_2)$ with

(1.14) $$ \begin{align} \mathbf{S}_{t,x}(z)= t^{-1/3} e^{\frac{2x^3}{3t^2}-\frac{zx}{t}}\hspace{-0.1em}\operatorname{\mathrm{\text{Ai}}}\left(-t^{-1/3}z+t^{-4/3}x^2\right). \end{align} $$

This is real analytic, and one checks directly (using $\operatorname {\mathrm {\text {Ai}}}''(z)=z\hspace {-0.1em}\operatorname {\mathrm {\text {Ai}}}(z)$ ) that

(1.15) $$ \begin{align} \partial_t\mathbf{S}_{t,x}(z) = \tfrac13 \partial_z^3\mathbf{S}_{t,x}(z),\qquad \partial_x\mathbf{S}_{t,x}(z) = \partial_z^2\mathbf{S}_{t,x}(z). \end{align} $$

Now fix $(x_1,\dots c,x_n),(r_1,\dots c,r_n)\in \mathbb {R}^n$ , introduce auxiliary variables $x,r\in \mathbb {R}$ and let

(1.16) $$ \begin{align} \widehat{\mathbf{K}}(t,x,r)=\mathbf{K}(t,x_1+x,\dots c,x_n+x,r_1+r,\dots c, r_n+r,\mathfrak{h}_0,\cdot,\cdot). \end{align} $$

The differential relations in equation (1.15) together with the fact that r enters the kernel just as a shift give Theorem 1.3(1)–(3) for $\widehat {\mathbf {K}}$ . On the other hand, it is shown in [Reference Matetski, Quastel and RemenikMQR17, Appx. A.1] that $ \Vert \widehat {\mathbf {K}}\Vert _1\longrightarrow 0$ as $r\to \infty $ and that the integrals converge absolutely and uniformly in compact balls around $t>0$ , $x\in \mathbb {R}^n$ and $r\in \mathbb {R}^n$ , and thus $\widehat {\mathbf {K}}$ is real analytic in those variables as well. Thus we have the following:

Theorem 1.4. The kernel $\widehat {\mathbf {K}}$ in equation (1.16) satisfies the conditions of Theorem 1.3.

Translating the partial derivatives in r and x in equation (1.12) back to the $r_i$ s and $x_i$ s gives equation (1.8) for the kernel $\mathbf {K}$ .

There is an informal way [Reference Matetski, Quastel and RemenikMQR17] to write the kernel of the KPZ fixed point which makes Theorem 1.3(2) and (3) completely apparent (equation (1) comes simply again from the way the $r_i$ s come in the formula). For $\mathfrak {h}\in \operatorname {\mathrm {\text {UC}}}$ and $\ell _1<\ell _2$ , let

$$ \begin{align*} {\textbf{P}}_{\ell_1,\ell_2}^{\text{No hit}\,\mathfrak{h}}(u_1,u_2)\hspace{0.05em} \mathrm{d} u_2&=\mathbb{P}_{\mathbf{B}\left(\ell_1\right)=u_1}\!\left(\mathbf{B}(y)>\mathfrak{h}(y)\text{ on }[\ell_1,\ell_2],\,\mathbf{B}(\ell_2)\in\mathrm{d} u_2\right),\\ {\textbf{P}}_{\ell_1,\ell_2}^{\text{Hit}\,\mathfrak{h}}&=e^{\left(\ell_2-\ell_1\right)\partial^2}- {\textbf{P}}_{\ell_1,\ell_2}^{\text{No hit}\,\mathfrak{h}}, \end{align*} $$

where $\mathbf {B}$ is again a Brownian motion with diffusion coefficient $2$ and $\partial $ denotes the derivative operator. The Brownian scattering transform of $\mathfrak {h}$ is the formal object

(1.17) $$ \begin{align} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})} =\lim_{\ell_1\to -\infty\atop \ell_2\to \infty}e^{\ell_1\partial^2}{\textbf{P}}_{\ell_1,\ell_2}^{\text{Hit}\,\mathfrak{h}}e^{-\ell_2\partial^2} =\mathbf{I}-\lim_{\ell_1\to -\infty\atop \ell_2\to \infty} e^{\ell_1\partial^2}{\textbf{P}}_{\ell_1,\ell_2}^{\text{No hit}\,\mathfrak{h}}e^{-\ell_2\partial^2}, \end{align} $$

where ${\textbf {P}}_{\ell _1,\ell _2}^{\text {Hit/No hit}\,\mathfrak {h}}$ are thought of as operators with the given integral kernels. This doesn’t make sense, since the backward heat operator is asked to act on nonanalytic functions. However, in the formula for the KPZ fixed point, $\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})}$ is never actually used by itself, but only after conjugation by the Airy unitary group

(1.18) $$ \begin{align} \mathbf{U}_t=e^{ -\frac13t \partial^3}, \end{align} $$

with $t\neq 0$ ; and as we explain next, in our setting the conjugated kernel is in fact well defined.

For $a\in \mathbb {R}$ , define

(1.19) $$ \begin{align} \chi_a(x)=\mathbf{1}_{x>a}\qquad\text{and}\qquad\bar\chi_a(x)=\mathbf{1}_{x\leq a}, \end{align} $$

which we also regard as multiplication operators acting on $L^2(\mathbb {R})$ . For $t>0$ the Airy semigroup acts by convolution with Airy functions – that is, $\mathbf {U}_t$ has integral kernel $\mathbf {U}_t(x,y)=t^{-1/3}\operatorname {\mathrm {\text {Ai}}}\left (t^{-1/3}(x-y)\right )$ ; its inverse $\mathbf {U}_t^{-1}=\mathbf {U}_{-t}$ equals $\mathbf {U}_t^*$ . The Airy functions are not themselves in $L^2(\mathbb {R})$ ; however, for $t>0$ and $r>-\infty $ , $ \mathbf {U}_t^{-1}\chi _{r}$ maps $L^2(\mathbb {R})$ into the domain of $e^{x \partial ^2}$ for any $x\in \mathbb {R}$ . So for $t>0$ and $r>-\infty $ , we define on $L^2([r,\infty ))$

(1.20) $$ \begin{align} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_t= \lim_{\ell_1\to -\infty\atop \ell_2\to \infty} \mathbf{U}_t e^{\ell_1\partial^2}{\textbf{P}}_{\ell_1,\ell_2}^{\text{Hit}\,\mathfrak{h}}e^{-\ell_2\partial^2}\mathbf{U}_t^{-1}. \end{align} $$

For any $t>0$ and $r>-\infty $ , the limit on the right-hand side of equation (1.20) exists in trace class on $L^2([r,\infty ))$ [Reference Quastel and RemenikQR19; Reference Matetski, Quastel and RemenikMQR17] and defines the left-hand side as a trace class operator in this space. The limit is given by the right-hand side of equation (1.13) with $x_i=x_j=x$ and $r_i=r_j=0$ ; the convergence was proved first in [Reference Quastel and RemenikQR19]. It satisfies the semigroup property

$$ \begin{align*} \mathbf{U}_s\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_t\mathbf{U}_s^{-1}=\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_{\mathbf{t}+\mathbf{s}}, \end{align*} $$

so we can write (at least informally)

(1.21) $$ \begin{align} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_t=\mathbf{U}_t\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}\mathbf{U}_t^{-1}. \end{align} $$

Note that we avoid the problem of domains by defining the left-hand side of equation (1.20) not as a product of three operators but just as one operator with the semigroup property. In this sense the Brownian scattering operator is the germ of the semigroup. Alternatively, one can think of the Brownian scattering operator as the entire semigroup $\left (\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})}_t\right )_{t>0}$ from equation (1.21). The fact that equation (1.17) is formal is important, though. We will see in formula (5) that the limit of equation (1.21) as $t\searrow 0$ is not $\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})}$ .

From $\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})}$ we build an extended Brownian scattering operator acting on $L^2(\{x_1,\dotsc ,x_m\}\times \mathbb {R})$ ,

(1.22) $$ \begin{align} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_{\text{ext}}\left(x_i,\cdot;x_j,\cdot\right)=-e^{\left(x_j-x_i\right)\partial^2}\mathbf{1}_{x_i<x_j}+e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})} e^{x_j\partial^2}, \end{align} $$

with the analogous caveat that in order to make sense of it, each of the $\left (x_i,x_j\right )$ entries should be conjugated by $\mathbf {U}_t$ and surrounded by $\chi _{r_i}$ on the left and $\chi _{r_j}$ on the right. Then the kernel of the KPZ fixed point appearing on the right-hand side of equation (1.6) and given in equation (1.13) can be written as

(1.23) $$ \begin{align} \mathbf{K}_{ij}(u_1,u_2)=\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_{t,\text{ext}}\left(x_i,u_1+r_i;x_j,u_2+r_j\right)\quad\text{with}\quad \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{h}_0\right)}_{t,\text{ext}}=\mathbf{U}_t\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{h}_0\right)}_{\text{ext}}\mathbf{U}_t^{-1}, \end{align} $$

where the Airy operators act on the left and right of each entry of $\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {h}_0\right )}_{\text {ext}}$ . From this expression one sees that $\partial _t\mathbf {K}_{ij}=-\frac 13\left (\partial ^3\mathbf {K}_{ij}-\mathbf {K}_{ij}\partial ^3\right )$ and $\partial _{x_\ell }\mathbf {K}_{ij}=-\frac 12\left (\partial ^2\mathbf {K}_{ij}\mathbf {1}_{i=\ell }-\mathbf {K}^2_{ij}\partial ^2\mathbf {1}_{j=\ell }\right )$ , which, after integration by parts and summing over the $x_\ell $ s as before, become Theorem 1.3(2) and (3).

Remark 1.5. For compactly supported initial data, which in our context means that $\mathfrak {h}(y)=-\infty $ for y outside some compact interval, the definition of the Brownian scattering operator in equation (1.20) becomes explicit, without any limits. In fact, consider $\mathfrak {h}$ with support on some interval $[-a,a]$ and let $\ell _1<-a$ and $\ell _2>a$ . Clearly ${\textbf {P}}_{\ell _1,\ell _2}^{\text {Hit}\,\mathfrak {h}}=e^{-\left (a+\ell _1\right )\partial ^2}{\textbf {P}}_{-a,a}^{\text {Hit}\,\mathfrak {h}}e^{\left (\ell _2-a\right )\partial ^2}$ , so $\mathbf {U}_te^{\ell _1\partial ^2}{\textbf {P}}_{\ell _1,\ell _2}^{\text {Hit}\,\mathfrak {h}}e^{-\ell _2\partial ^2}\mathbf {U}_t^{-1}=\mathbf {U}_te^{-a\partial ^2}{\textbf {P}}_{-a,a}^{\text {Hit}\,\mathfrak {h}}e^{-a\partial ^2}\mathbf {U}_t^{-1}$ . Since $\mathbf {U}_te^{x\partial ^2}=\mathbf {S}_{t,x}$ and $\mathbf {U}_t^{-1}e^{x\partial ^2}=\mathbf {S}_{-t,x}$ (see equation (1.15)), this shows that

$$ \begin{align*} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_t=\mathbf{S}_{t,-a}{\textbf{P}}_{-a,a}^{\text{Hit}\,\mathfrak{h}}\mathbf{S}_{-t,-a}. \end{align*} $$

The right-hand side coincides with the right-hand side of equation (1.13) with $i=j$ and $x_i=r_i=0$ (see [Reference Matetski, Quastel and RemenikMQR17, Sec. 4.1]).

1.2 Initial data

The natural class of initial data for our problem (the ‘ $1$ -dimensional substrate’) corresponds to functions in the spaceFootnote 2

(1.24) $$ \begin{align} &\operatorname{\mathrm{\text{UC}}}=\big\{\mathfrak{h}\!:\mathbb{R}\longrightarrow[-\infty,\infty)\!:\mathfrak{h}\text{ is upper semicontinuous,}\nonumber\\ &\qquad\qquad\qquad\quad\mathfrak{h}(x)\leq A+B\lvert x\rvert\text{ for some}\ A,B>0\ \text{and }\mathfrak{h}\not\equiv-\infty\big\}. \end{align} $$

A function is upper semicontinuous if and only if its hypograph $\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h}) = \{(\mathbf {x},y): y\le \mathfrak {h}(\mathbf {x})\}$ is closed in $[-\infty ,\infty )\times \mathbb {R}$ . We endow $[-\infty ,\infty )$ with the distance $d_{\left [-\infty ,\infty \right )}(y_1,y_2) = \lvert e^{y_1} - e^{y_2}\rvert $ , and use the topology of local Hausdorff convergence, which means Hausdorff convergence of the restrictions to $-L\le x\le L$ of $\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h}_n)$ to $\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})$ for each $L>0$ .

Example 1.6 Finite collection of narrow wedges

Set $a_1<a_2<\dotsb <a_k$ , $b_1,b_2,\dotsc ,b_k\in \mathbb {R}$ . Then $\mathfrak {h}=\mathfrak {d}_{\vec a}^{\vec b}$ is in $\operatorname {\mathrm {\text {UC}}}$ , with $\mathfrak {d}_{\vec a}^{\vec b}$ defined by

$$ \begin{align*}\mathfrak{d}_{\vec a}^{\vec b}(x)=b_i\quad\text{if}\ x=a_i\ \text{for some}\ i,\qquad\mathfrak{d}_{\vec a}^{\vec b}(x)=-\infty\quad\text{otherwise}.\end{align*} $$

The initial data for the one-point case (equation (1.9)) is the escarpment

(1.25) $$ \begin{align} \phi(0,x,r) = 0\quad\text{for }r\ge \mathfrak{h}_0(x),\qquad\phi(0,x,r)=-\infty\quad\text{for }r<\mathfrak{h}_0(x). \end{align} $$

These are unusual and do not fit into any well-posedness schemes known for the KP equation.Footnote 3 Although the infinity looks formal, we believe the solutions to the equations with such initial data are well posed, but we leave the proofs for future work. They also appear not to develop solitons. Since F is given by a Fredholm determinant, these initial conditions represent an entirely new class of integrable initial data for the KP equation.

The initial data for the matrix KP equation (1.8) are formally

(1.26)

where

(1.27) $$ \begin{align} {\textbf{P}}^{{\geq\mathfrak{h}_0,\leq -\mathfrak{d}_{\vec x}^{-\vec r}}}_{x_i,x_j}\!\left(r_i,r_j\right)dr_j&=\mathbb{P}_{\mathbf{B}\left(x_i\right)=r_i}\Big(\mathbf{B}(y)\geq \mathfrak{h}_0(y)~\forall\,y\in\left[x_i,x_j\right],\nonumber\\ &\mathbf{B}(x_n)\leq r_n\text{ for each }x_n\in\left(x_i,x_j\right), \mathbf{B}\left(x_j\right)\in dr_j\Big). \end{align} $$

The probability is with respect to a Brownian motion $\mathbf {B}$ with diffusivity $2$ starting at $r_i$ at time $x_i$ . This is derived in Appendix A for finite collections of narrow wedges. Unlike the scalar case, one can see immediately that the initial data as written are insufficient, because in the matrix product the $0$ and $\infty $ interact. Wrongly interpreted, it appears to produce anomalous solutions, so the initial data written would have to be augmented by at least some description of the rate of convergence to $0$ and $\infty $ in the $t\downarrow 0$ limit. We leave this also for future work.

Remark 1.7. When $t=0$ , the escarpment initial data equation (1.25) just mean (sending x to $-x$ ) that L from equation (1.11) is the heat operator with Dirichlet boundary data on the hypograph of $\mathfrak {h}_0$ . This corresponds to the Brownian scattering transform which computes transition probabilities of Brownian motions killed when passing through that hypograph.

Outline

In Section 2 we show how the famous Painlevé expressions for the Tracy–Widom distributions just arise as special examples of the KP equation, and write equations for a few other special initial data, in particular the Airy processes. In Section 3 we show that some of the explicit formulas for the $1$ -dimensional distributions of the KPZ equation (as opposed to the KPZ fixed point) also satisfy the KP equation. Section 4 contains the proof of Theorem 1.3.

2 Examples

2.1 Tracy–Widom distributions

A key observation is that the GUE and GOE Tracy–Widom distributions are now seen to arise simply as special similarity solutions of the KP equation (1.9):

Example 2.1 Tracy–Widom GUE distribution

Consider $\mathfrak {h}_0=\mathfrak {d}_0$ , the narrow wedge initial condition defined as $\mathfrak {d}_0(0)=0$ and $\mathfrak {d}_0(x)=-\infty $ , $x\neq 0$ . With this choice of initial data one has $\mathfrak {h}(t,x)+x^2/t\stackrel {\text {dist}}{=}t^{1/3}\mathcal {A}\left (t^{-2/3}x\right )$ , where $\mathcal {A}$ is the Airy $_2$ process (see Section 2.3), which is stationary (in x). From this and the 1:2:3 scaling invariance of equation (1.9), it is natural to look for a self-similar solution of the form

$$ \begin{align*} \phi^{\text{nw}} (t,x,r)=t^{-2/3} \psi^{\text{nw}} \left(t^{-1/3} r + t^{-4/3}x^2\right). \end{align*} $$

This turns equation (1.9) into

(2.1) $$ \begin{align} (\psi^{\text{nw}})'''+12\psi^{\text{nw}}(\psi^{\text{nw}})'-4r(\psi^{\text{nw}})'-2\psi^{\text{nw}}=0. \end{align} $$

The transformation $\psi ^{\text {nw}}=-\mathfrak {q}^2$ takes equation (2.1) into Painlevé II:

(2.2) $$ \begin{align} \mathfrak{q}''=r\mathfrak{q}+2\mathfrak{q}^3. \end{align} $$

As $r\to -\infty $ , the solution is approximately $\phi ^{\text {nw}}(t,x,r)\sim -\left (\tfrac {r}{2t} +\tfrac {x^2}{2t^2}\right )$ , picking out the Hastings–McLeod solution $\mathfrak {q}(r)\sim -\operatorname {\mathrm {\text {Ai}}}(r)$ as $r\to \infty $ . Thus we recover

$$ \begin{align*} F(t,x,r)= F_{\text{GUE}}\left(t^{-1/3} r+t^{-4/3}x^2\right), \end{align*} $$

where $F_{\text {GUE}}$ is the GUE Tracy–Widom distribution [Reference Tracy and WidomTW94], usually written in the equivalent form

$$ \begin{align*} F_{\text{GUE}}(s) = \exp\left\{ -\int_s^\infty\!\mathrm{d} u\, (u-s)\mathfrak{q}^2(u)\right\}. \end{align*} $$

Example 2.2 Tracy–Widom GOE distribution

If $\mathfrak {h}_0(x)\equiv 0$ , the flat initial condition, there is no x dependence and equation (1.9) reduces to the KdV equation. Now we look for a self-similar solution of equation (1.10) in the form

$$ \begin{align*} \phi^{\text{fl}}(t,r) = (t/4)^{-2/3} \psi^{\text{fl}}\left((t/4)^{-1/3} r\right) \end{align*} $$

(the extra factor of $1/4$ is just to coordinate conventions with random matrix theory), obtaining the ordinary differential equation

$$ \begin{align*} \left(\psi^{\text{fl}}\right)''' +12\left(\psi^{\text{fl}}\right)'\psi^{\text{fl}}-r\left(\psi^{\text{fl}}\right)' -2\psi^{\text{fl}} =0. \end{align*} $$

Miura’s transform

$$ \begin{align*} \psi^{\text{fl}} = \tfrac12 \left( \mathfrak{q}' -\mathfrak{q}^2\right) \end{align*} $$

brings this to Painlevé II (equation (2.2)), with the same behaviour as $r\to -\infty $ . So we recover

$$ \begin{align*} F(t,x,r)= F_{\text{GOE}}\left(4^{1/3}t^{-1/3} r\right), \end{align*} $$

where $F_{\text {GOE}}$ is the GOE Tracy–Widom distribution [Reference Tracy and WidomTW96], usually written in the equivalent form

$$ \begin{align*} F_{\text{GOE}}(r) = \exp\left\{ -\frac12\int_r^\infty\!\mathrm{d} u\,\mathfrak{q}(u)\right\}F_{\text{GUE}}(r)^{1/2}. \end{align*} $$

These two examples also have the following interpretation: Let $\lambda ^{\text {max,GUE}}_N$ and $\lambda ^{\text {max,GOE}}_N$ be the largest eigenvalues of $N\times N$ matrices chosen from the Gaussian unitary and Gaussian orthogonal ensembles multiplied by $\sqrt {N}$ , so that $\lambda ^{\text {max,GUE}}_N\sim 2N + N^{1/3}\zeta _{\text {GUE}}$ and $\lambda ^{\text {max,GOE}}_N\sim 2N + N^{1/3}\zeta _{\text {GOE}}$ , with $\zeta _{\text {GUE}}$ and $\zeta _{\text {GOE}}$ the standard Tracy–Widom GUE and GOE random variables. Let

$$ \begin{align*} F_1(t,r)& = {\lim_{N\to \infty}} \mathbb{P}\left(N^{-1/3}\left( \lambda^{\text{max,GOE}}_{Nt}-2Nt\right) \le 4^{1/3}r\right)=F_{\text{GOE}}\left(4^{1/3}t^{-1/3}r\right),\\ F_2(t,x,r) &= {\lim_{N\to \infty}} \mathbb{P}\left(N^{-1/3}\left( \lambda^{\text{max,GUE}}_{Nt}-2 Nt\right) \le r+x^2/t\right)=F_{\text{GUE}}\left(t^{-1/3} r+t^{-4/3}x^2\right). \end{align*} $$

As we have seen, $\partial _r^2\log F_1$ and $\partial _r^2\log F_2$ satisfy the KP equation (1.9). In the former case, there is no dependence on x, and the KP equation reduces to the KdV equation (1.10).

2.2 Partial differential equations for other initial data

Another question is whether there are analogues of Painlevé II for other self-similar solutions. It is natural to observe $\phi $ in the frame of the inviscid solution $\tfrac 14\left (\partial _x\bar {\mathfrak {h}}\right )^2-\partial _t\bar {\mathfrak {h}}=0$ of Burgers’ equation,

$$ \begin{align*} \phi(t,x,r) := \bar\phi\left(t,x,r-\bar{\mathfrak{h}}(t,x)\right); \end{align*} $$

one obtains

(2.3) $$ \begin{align}\partial_t \bar\phi +\bar\phi\hspace{0.1em}\partial_r \bar\phi + \tfrac1{12}\partial_r^3 \bar\phi + \tfrac14\partial_r^{-1} \partial_x^2 \bar\phi -\tfrac12\partial_x\bar{\mathfrak{h}}\hspace{0.1em}\partial_x\bar\phi+ V \bar\phi= 0,\end{align} $$

with $V=-\tfrac {1}4\partial ^2_x\bar {\mathfrak {h}}$ and with initial data $\bar \phi (0,x,r) = 0$ for $r\ge 0$ and $-\infty $ for $r<0$ .

In order to get a solution for the rescaled spatial process, write

$$ \begin{align*} \phi(t,x,r) = t^{-2/3}\psi\left(t,t^{-2/3}x,t^{-1/3}\left(r- \bar{\mathfrak{h}}(t,x)\right)\right). \end{align*} $$

Then

$$ \begin{align*}t\hspace{0.05em}\partial_t\psi-\tfrac23\psi-\tfrac13r\partial_r\psi+\tfrac1{12}\partial_r^3\psi+\psi\partial_r\psi+\tfrac14\partial_r^{-1}\partial_x^2\psi -\tfrac{1}4\partial_x^2\tilde{\mathfrak{h}}\hspace{0.05em}\psi -\left(\tfrac23x +\tfrac12\partial_x\tilde{\mathfrak{h}}\right)\partial_x\psi =0,\end{align*} $$

with $\tilde {\mathfrak {h}}(t,x)=t^{-1/3}\bar {\mathfrak {h}}\left (t,t^{2/3}x\right )$ .

Example 2.3 Half-flat initial data

Consider $\mathfrak {h}_0(x)=0$ , $x\le 0$ and $\mathfrak {h}_0(x)=-\infty $ , $x>0$ . Now $\bar {\mathfrak {h}}(t,x) = -x^2/t\hspace {0.05em}\mathbf {1}_{x\geq 0}$ . There is dependence on x, though not on t. This gives rise to a partial differential equation for $\psi (x,r)$ ,

$$ \begin{align*}-\tfrac13r\partial_r\psi+\tfrac1{12}\partial_r^3\psi-\left(\tfrac16\mathbf{1}_{x\ge0}+\tfrac23\mathbf{1}_{x<0}\right)\psi+\psi\partial_r\psi+\tfrac14\partial_r^{-1}\partial_x^2\psi +\left(\tfrac13\mathbf{1}_{x\ge 0}-\tfrac23\mathbf{1}_{x<0}\right)x\partial_x\psi=0.\end{align*} $$

Remark 2.4 Lower-tail heuristics

Typically the equation is controlled on large scales by the equation with the third derivative dropped, and Burgers’ equation makes sense for such wedge-type initial data. Let $\bar \phi $ be as in equation (2.3) and

$$ \begin{align*} \check\phi = \bar\phi -\eta, \end{align*} $$

where $\eta =\tfrac {c }{t} r\mathbf {1}_{r<0}$ . Then $\bar \phi (0,x,r)=0$ and

$$ \begin{align*} \partial_t \check\phi + \tfrac12\partial_r \check\phi^2+ \tfrac1{12}\partial_r^3\check\phi+ \tfrac14\partial_r^{-1} \partial_x^2\check\phi + \partial_r\left(\eta\check\phi\right) &-\tfrac12 \partial_x \bar{\mathfrak{h}}\partial_x\check\phi-\tfrac14 \partial^2_x \bar{\mathfrak{h}}~ \check\phi \\ &+\left(\tfrac{c-1}{t} -\tfrac14 \partial^2_x \bar{\mathfrak{h}}\right) \eta - \tfrac{c }{12t} \delta'_0(r) = 0. \end{align*} $$

One hopes to set things up so that $\check \phi $ has good decay at $\pm \infty $ . Consider our two basic examples. In the flat case, $\bar {\mathfrak {h}}\equiv 0$ , and we want to take $c=1$ to make the second-to-last term drop out, which will lead to the conclusion that $\phi (t,r)\sim r/t$ , or $\log F(t,r)\sim - \tfrac 1{6t} r^3$ as $r\to -\infty $ , recovering the Tracy–Widom GOE lower tail. In the narrow wedge, case $\bar {\mathfrak {h}} = -\tfrac {x^2}{t}$ , so $\tfrac 14 \partial ^2_x \bar {\mathfrak {h}}=-\tfrac {1}{2t}$ , and we want to take $c=1/2$ leading to $\log F(t,x,r)\sim -\tfrac 1{12 t} \left (r+ \tfrac {x^2}{t}\right )^3 $ , recovering the Tracy–Widom GUE lower tail.

The conclusion is that the lower tail of the distributions can be seen directly from the ‘Burgers’ part of the KP equation, which dominates in that region.

Note that this is only intended to be a quick heuristic. If one were to try to make it rigorous, say, for the flat and narrow wedge case, one would be led through self-similar solutions back to the Painlevé equation and the standard tools there to get the left tails (see, e.g., [Reference Baik, Buckingham and DiFrancoBBD08]). Such Riemann–Hilbert methods give much finer information, for example, subleading terms. It is an interesting question whether such results can be obtained from the partial differential equation (2.3) so as to get finer information on the left tail for more general initial data.

2.3 Airy process

The Airy process $\mathcal {A}(x)$ is defined as

$$ \begin{align*} \mathcal{A}(x) := \mathfrak{h}(1,x;\mathfrak{d}_0)+x^2, \end{align*} $$

where $\mathfrak {h}(t,x;\mathfrak {d}_0)$ is the KPZ fixed point starting from a narrow wedge $\mathfrak {d}_0$ at the origin. The Airy process is stationary and the one-point distribution is the Tracy–Widom GUE distribution.

K. Johansson famously asked whether there is an equation for the multipoint distribution. Partial differential equations (PDEs) were given by [Reference Adler and van MoerbekeAM05; Reference Tracy and WidomTW03]. Starting with narrow wedge initial data, the matrix KP equation (1.8) gives us another PDE. This equation lives at the same level as the system of PDEs in [Reference Tracy and WidomTW03]; it could well be that they are equivalent, but we have not succeeded in confirming that yet, and leave it for future work.

Note that the KP-II equation is written in variables t, $r=r_1+\dotsb +r_n$ and $x=x_1+\dotsb +x_n$ , the rest of the variables entering only from the boundary condition. Exploiting skew time-reversal invariance, we find an extra symmetry in the narrow wedge case. Let $\mathfrak {d}^{\vec r}_{\vec x}$ denote a multiple narrow wedge (see Example 1.6). For $\vec x\in \mathbb {R}^m$ and $z\in \mathbb {R}$ , write $\vec x+z=(x_1+z,\dotsc ,x_m+z)$ . By skew time reversibility and translation and affine invariance of the KPZ fixed point [Reference Matetski, Quastel and RemenikMQR17, Thm. 4.5],

(2.4) $$ \begin{align} F\left(t,\vec x+z,\vec r+a\right)&=\mathbb{P}\big(\mathfrak{h}(t,x_i+z;\mathfrak{d}_0)\leq r_i+a\text{ for all }i\big) =\mathbb{P}\left(\mathfrak{h}(t,\cdot;\mathfrak{d}_0)\leq-\mathfrak{d}^{-\vec r-a}_{\vec x+z}\right)\nonumber\\ &=\mathbb{P}\left(\mathfrak{h}\left(t,\cdot;\mathfrak{d}^{-\vec r-a}_{\vec x+z}\right)\leq-\mathfrak{d}_0\right) =\mathbb{P}\left(\mathfrak{h}\left(t,0;\mathfrak{d}^{-\vec r}_{\vec x+z}\right)\leq a\right)\nonumber\\ &=\mathbb{P}\left(\mathfrak{h}\left(t,-z;\mathfrak{d}^{-\vec r}_{\vec x}\right)\leq a\right). \end{align} $$

Now the right-hand side is just the one-point distribution at $-z$ with a given, fixed initial condition. So if we let $G(t,z,a)=F\left (t,\vec x+z,\vec r+a\right )$ , we see that $\partial _a^2\log G$ satisfies equation (1.9) in $(t,z,a)$ (in terms of equation (1.25), the initial data are $G(0,z,a)=-\infty $ if $z=-x_i$ for some i and $a<-r_i$ and $G(0,z,a)=0$ otherwise). But $\partial _zG(t,z,a)=\mathcal {D}_xF\left (t,\vec x+z,\vec r+a\right )$ , and similarly, $\partial _aG(t,z,a)=\mathcal {D}_rF\left (t,\vec x+z,\vec r+a\right )$ . So now setting $a=z=0$ , we deduce that $\phi =\mathcal {D}_r^2\log F$ satisfies

(2.5) $$ \begin{align} \partial_t \phi + \tfrac12 \mathcal{D}_r \phi^2 + \tfrac1{12}\mathcal{D}_r^3 \phi + \tfrac14\mathcal{D}_r^{-1} \mathcal{D}_x^2 \phi = 0. \end{align} $$

The initial data are now similarly $\phi \left (0,\vec x,\vec r\right )=-\infty $ if $x_i=0$ for some i and $r_i<0$ and $\phi \left (0,\vec x,\vec r\right )=0$ otherwise.

Note that from Theorem 1.1, we know that $\mathcal {D}_r^2\log F=\operatorname {{\mathop {\mathrm {{tr}}}}} q$ and $\partial _t\operatorname {{\mathop {\mathrm {{tr}}}}} q+\operatorname {{\mathop {\mathrm {{tr}}}}}(q\mathcal {D}_rq)+\tfrac 1{12}\mathcal {D}_r^3\operatorname {{\mathop {\mathrm {{tr}}}}} q+\tfrac 14\mathcal {D}_x^2\operatorname {{\mathop {\mathrm {{tr}}}}} Q=0$ (using $\operatorname {{\mathop {\mathrm {{tr}}}}}(AB)=\operatorname {{\mathop {\mathrm {{tr}}}}}(BA)$ ). But the foregoing argument implies then that $\operatorname {{\mathop {\mathrm {{tr}}}}} q$ itself solves the KP equation (2.5), and as a consequence, we deduce in the narrow wedge case that $\operatorname {{\mathop {\mathrm {{tr}}}}}(q\mathcal {D}_rq)=\operatorname {{\mathop {\mathrm {{tr}}}}}(q)\operatorname {{\mathop {\mathrm {{tr}}}}}(\mathcal {D}_rq)$ . This can also be proved directly using the fact that, in this case, $(\partial _u+\partial _v)\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {h}_0\right )}_t(u,v)$ is a rank $1$ kernel, which implies that q is a rank $1$ matrix.

An alternative derivation of equation (2.5) using the path integral formula for the KPZ fixed point can be found in Appendix B.

Next we use the 1:2:3 scaling invariance of the KPZ fixed point and stationarity of the Airy process, as in Example 2.1. Let H, $\Psi $ , $\psi $ denote the $(m+1)$ -point distribution function of the Airy process and its logarithmic derivatives,

(2.6) $$ \begin{align}H(y_1,\dotsc,y_{m},r_0,\dotsc,r_m) &= \mathbb{P}( \mathcal{A}(0)\le r_0,\mathcal{A}(y_1)\le r_1,\dotsc,\mathcal{A}(y_m)\leq r_m), \nonumber\\ \Psi&= \mathcal{D}_r \log H,\qquad \psi=\mathcal{D}_r \Psi. \end{align} $$

Since $F(t,x_0,\dotsc ,x_m,r_0,\dotsc ,r_m)=H\left (t^{-2/3}\bar x_1,\dotsc ,t^{-2/3}\bar x_m, t^{-1/3} \hat r_0, \dotsc ,t^{-1/3}\hat r_m\right )$ , with $\bar x_i=x_i-x_0$ and $\hat r_i=r_i+\tfrac {x^2_i}{t}$ , equation (2.5) leads to the following:

Theorem 2.5 Airy $_2$ -process multipoint function

The logarithmic derivatives (2.6) of the multipoint function of the Airy process satisfy the PDE

$$ \begin{align*} &\textstyle-\tfrac1{6}\psi + \tfrac12\mathcal{D}_r\psi^2 + \tfrac1{12} \mathcal{D}_r^3 \psi -\tfrac13\sum_{a=0}^mr_a\partial_{r_a}\psi-\tfrac23\sum_{a=1}^my_a\partial_{y_a}\psi\\ &\qquad\qquad\qquad\qquad\textstyle-\sum_{a=1}^m\sum_{b=0,\,b\neq a}^my_a^2\partial_{r_a}\hspace{-0.1em}\partial_{r_b}\Psi+\sum_{a,b=1,\,a\neq b}^my_ay_b\hspace{0.05em}\partial_{r_a}\hspace{-0.1em}\partial_{r_b}\Psi =0, \end{align*} $$

with boundary condition $\psi (0,\dotsc ,0,r_1,\dotsc ,r_m) =\psi ^{\text {nw}}( \min (r_1,\dotsc ,r_m))$ from equation (2.1).

In particular, for the logarithmic derivative of the Airy-process two-point function $\Psi (y,r_1,r_2)=\mathcal {D}_r\log \mathbb {P}(\mathcal {A}(0)\leq r_1,\,\mathcal {A}(y)\leq r_2)$ , we get

(2.7) $$ \begin{align} -\tfrac1{6}\psi+\tfrac12\mathcal{D}_r\psi^2+\tfrac1{12}\mathcal{D}_r^3\psi-\tfrac13(r_1\partial_{r_1}\psi+r_2\partial_{r_2}\psi)-\tfrac23y\hspace{0.05em}\partial_y\psi-y^2\partial_{r_1}\partial_{r_2}\Psi =0. \end{align} $$

In our notation, the equation derived in [Reference Adler and van MoerbekeAM05] reads

(2.8) $$ \begin{align} (r_2-r_1)\partial_{r_1}\partial_{r_2}\psi+y\hspace{0.05em}\partial_y\left(\partial_{r_1}-\partial_{r_2}\right)\psi+y^2(\partial_{r_1}-\partial_{r_2})\partial_{r_1}\partial_{r_2}\Psi +\partial_{r_1}\!\Psi\hspace{0.1em}\partial_{r_2}\hspace{-0.1em}\mathcal{D}_r\psi-\partial_{r_2}\!\Psi\hspace{0.1em}\partial_{r_1}\hspace{-0.1em}\mathcal{D}_r\psi=0. \end{align} $$

Although equations (2.7) and (2.8) have similarities, they do not appear to be equivalent. We have not succeeded in reconciling them, and we leave this for future work. Of course the Airy process has many symmetries, and it is plausible that the equations are not equivalent yet both hold.

Example 2.6 Airy $_1$ -process multipoint function

Consider now the KPZ fixed point with flat initial data, $\mathfrak {h}_0\equiv 0$ . The multipoint function satisfies

$$ \begin{align*}F\left(t,\vec x,\vec r\right)=\mathbb{P}\left(\mathcal{A}_1(0)\leq t^{-1/3}r_1,\,\mathcal{A}_1\left(t^{-2/3}\bar x_2\right)\leq t^{-1/3}r_2,\dotsc,\mathcal{A}_1\left(t^{-2/3}\bar x_m\right)\leq t^{-1/3}r_m\right) \end{align*} $$

with $\bar x_i=x_i-x_1$ , where $\mathcal {A}_1(x):=\mathfrak {h}(1,x;0)$ is the Airy $_1$ process and we have used the fact that it is stationary. We know that the left-hand side equals $\operatorname {{\mathop {\mathrm {{tr}}}}} Q$ where the $m\times m$ matrix Q and its derivative q solve equation (1.8). But the right-hand side depends only on the differences $x_i-x_1$ , so Q actually solves the (integrated) matrix KdV equation

$$ \begin{align*}\partial_tQ+\tfrac12(\mathcal{D}_rQ)^2+\tfrac1{12}\mathcal{D}_r^3Q=0.\end{align*} $$

Take now $m=2$ for simplicity. The (formal) initial data $Q(0,x_1,x_2,r_1,r_2)$ (equation (1.26)) are also invariant under $r_i,x_i\mapsto t^{-1/3}r_i,t^{-2/3}x_i$ followed by $Q\mapsto t^{-1/3}Q$ , and depend only on $\lvert x_2-x_1\rvert \,(\text {in fact, }Q_{12}(0,r_1,r_2,x_1,x_2) = \frac {\mathbf {1}_{r_1\ge 0,r_2\ge 0}}{\sqrt {4\pi \left \lvert x_2-x_1\right \rvert }}(e^{- \left (r_1+r_2\right )^2/4\left \lvert x_2-x_1\right \rvert } -e^{-\left (r_1-r_2\right )^2/4\left \lvert x_2-x_1\right \rvert } )\text { by the reflection principle} )$ . In this case we look for a matrix solution of equation (1.8) with the scaling

$$ \begin{align*}Q(t,r_1,r_2,x_1,x_2)= t^{-1/3}\bar{Q}\left(t^{-2/3}(x_2-x_1),t^{-1/3}r_1,t^{-1/3}r_2\right).\end{align*} $$

The conclusion is that

(2.9) $$ \begin{align} \mathcal{D}_r\log\mathbb{P}(\mathcal{A}_1(0)\leq r_1,\,\mathcal{A}_1(y)\leq r_2)=\operatorname{{\mathop{\mathrm{{tr}}}}}\!\left(\bar Q(y,r_1,r_2)\right), \end{align} $$

with $\bar Q(y,r_1,r_2)$ solving the matrix PDE

$$ \begin{align*}-\tfrac13\bar Q+\tfrac12\left(\mathcal{D}_r\bar Q\right)^2+\tfrac1{12}\mathcal{D}_r^3\bar Q-\tfrac13\left(r_1\partial_{r_1}\bar Q+r_2\partial_{r_2}\bar Q\right)-\tfrac23y\partial_y\bar Q=0.\end{align*} $$

It does not seem to be possible to turn this $4\times 4$ system into a closed equation for the left-hand side of equation (2.9); the fact that this works for the Airy $_2$ process is very particular to narrow wedge initial data (and follows from skew time reversibility as used in equation (2.4)).

3 The KP-II equation in special solutions of the KPZ equation

Theorem 1.3 also implies that some of the special explicit solutions for $1$ -dimensional distributions of the KPZ equation (1.1) satisfy the (scalar) KP-II equation (1.9) as well. At this point we do not know whether this is part of a more general fact or whether the KP-II equation only holds in these special cases because of some symmetry. All we have is examples.

3.1 Narrow wedge solution of the KPZ equation

Let ${\textrm {h}}_{\text {nw}}$ be the narrow wedge solution of equation (1.1) with $\lambda =\nu =\tfrac 14$ and $\sigma =1$ . In other words, ${\textrm {h}}_{\text {nw}}=\log Z$ , where Z is the fundamental solution of the stochastic heat equation with multiplicative noise:

(3.1) $$ \begin{align} \partial_t Z = \tfrac14 \partial^2_xZ + \xi Z,\qquad Z(0,x)=\delta_0(x). \end{align} $$

The KPZ generating function is

(3.2) $$ \begin{align} G_{\text{nw}}(t,x,r) = \mathbb{E}\left[\exp\!\left\{-e^{{\textrm{h}}_{\text{nw}}\left(t,x\right)+\frac{t}{12}-r}\right\}\right]. \end{align} $$

The distribution of ${\textrm {h}}_{\text {nw}}(t,x)$ was computed in 2010 in [Reference Amir, Corwin and QuastelACQ11; Reference Sasamoto and SpohnSS10; Reference DotsenkoDot10; Reference Calabrese, Doussal and RossoCDR10], with the resultFootnote 4 that $G_{\text {nw}}(t,x,r)=\det (\mathbf {I} - \mathbf {K})_{L^2(\mathbb {R}^+)}$ with

$$ \begin{align*}\mathbf{K}(u,v)=\int_{-\infty}^\infty dy\,t^{-2/3}\frac{1}{1+e^{y}}\operatorname{\mathrm{\text{Ai}}}\left(t^{-1/3}(u+r-y)+t^{-4/3}x^2\right)\operatorname{\mathrm{\text{Ai}}}\left(t^{-1/3}(v+r-y)+t^{-4/3}x^2\right).\end{align*} $$

If we conjugate the operator by multiplying the kernel by $e^{\left (v-u\right )x/t}$ , the value of the Fredholm determinant doesn’t change – that is, we also have

$$ \begin{align*}G_{\text{nw}}(t,x,r)=\det\left(\mathbf{I} - \widetilde{\mathbf{K}}\right)_{L^2(\mathbb{R}^+)}\quad\text{with}\quad\widetilde{\mathbf{K}}\left(u,v\right)=e^{(v-u)x/t}\mathbf{K}(u,v).\end{align*} $$

One checks directly that $\widetilde {\mathbf {K}}$ satisfies the differential relations Theorem 1.3(1)–(3). This can also be seen by writing $\widetilde {\mathbf {K}}(u,v) = \mathbf {U}_te^{-x\partial ^2}\hspace {0.05em}\mathbf {M}\hspace {0.05em} e^{x\partial ^2}\mathbf {U}_t^{-1}(u+r,v+r)$ , with $\mathbf {M}$ the multiplication operator $\mathbf {M} f(u) = (1+e^{u})^{-1}f(u)$ and $\mathbf {U}_t$ the Airy unitary operator defined in equation (1.18). We then have the following:

Theorem 3.1. $\phi _{\text {nw}}:=\partial _r^2\log G_{\text {nw}}$ satisfies the KP-II equation (1.9).

The initial condition is $\lim _{t\searrow 0}\phi _{\text {nw}}\left (t,x,r-\tfrac {x^2}{t} -\log \sqrt {\pi t} \right ) = -e^{-r}$ . This suggests defining the x-independent shifted variable $\hat \phi _{\text {nw}}(t,r )= \phi _{\text {nw}}\left (t,x,r-\tfrac {x^2}{t} -\log \sqrt {\pi t} \right )$ , which now satisfies the cylindrical KdV equation

$$ \begin{align*} \partial_t \hat \phi_{\text{nw}} +\tfrac1{2t}\partial_r \hat \phi_{\text{nw}}+ \hat \phi_{\text{nw}}\partial_r \hat \phi_{\text{nw}}+ \tfrac1{12}\partial_r^3 \hat \phi_{\text{nw}} +\tfrac1{2t} \hat \phi_{\text{nw}}= 0,\qquad \hat \phi_{\text{nw}}(0,r) =-e^{-r}. \end{align*} $$

3.2 KPZ equation with spiked/half-Brownian initial data

Consider now the solution ${\textrm {h}}_{b}$ of equation (1.1) with $\lambda =\nu = \tfrac 14$ , $\sigma =1$ and m-spiked initial data, where $b=(b_1,\dotsc ,b_m)\in \mathbb {R}^m$ are the spike parameters. When $m=1$ , this corresponds to half-Brownian initial data (more precisely, at the level of the stochastic heat equation one sets $Z(0,x)=e^{B(x)+b_1x}\mathbf {1}_{x\geq 0}$ , where $B(x)$ is a Brownian motion with diffusivity 2); for the definition in the general case $m\geq 1$ , we refer to [Reference Borodin, Corwin and FerrariBCF14, Defn. 1.9]. Define $G_b$ as in equation (3.2) (with ${\textrm {h}}_b$ in place of ${\textrm {h}}_{\text {nw}}$ ). Then from [Reference Borodin, Corwin and FerrariBCF14, Thm. 1.10] we get that $G_b(t,0,r)=\det (\mathbf {I} - \mathbf {K}_0)_{L^2\left [0,\infty \right )}$ (see footnote 4 again), with

$$ \begin{align*} \mathbf{K}_0(u,v) =\frac{1}{(2\pi{\textrm{i}})^{2}}\int_{\mathcal{C}_t} d\eta\int_{\mathcal{C}_t'} d\xi\,\frac{t^{-1/3}\pi}{\sin\left(t^{-1/3}\pi(\xi-\eta)\right)}\frac{e^{\xi^3/3-\left(u+t^{-1/3}r\right)\xi}}{e^{\eta^3/3-\left(v+t^{-1/3}r\right)\eta}}\prod_{k=1}^m\frac{\Gamma\left(t^{-1/3}\eta-b_k\right)}{\Gamma\left(t^{-1/3}\xi-b_k\right)}, \end{align*} $$

where $\mathcal {C}_t$ goes from $-\frac 14t^{1/3}-{\textrm {i}}\infty $ to $-\tfrac 14t^{1/3}+{\textrm {i}}\infty $ crossing the real axis to the right of $t^{1/3}b_1,\dotsc ,t^{1/3}b_m$ and $\mathcal {C}_t'=\mathcal {C}_t+\frac 12t^{1/3}$ . We scale $(\eta ,\xi )\mapsto \left (t^{1/3}\eta ,t^{1/3}\xi \right )$ and $(u,v)\mapsto \left (t^{-1/3}u,t^{-1/3}v\right )$ (in the Fredholm determinant) so that $\mathbf {K}_0(u,v)$ is now given as

$$ \begin{align*}\int_{\mathcal{C}_1} d\eta\int_{\mathcal{C}_1'} d\xi\,\frac{\pi}{\sin(\pi(\xi-\eta))}\frac{e^{t\xi^3/3-(u+r)xi}}{e^{t\eta^3/3-(v+r)\eta}}\prod_{k=1}^m\!\frac{\Gamma(\eta-b_k)}{\Gamma(\xi-b_k)}.\end{align*} $$

Since ${\textrm {h}}_{b+x/t}(t,x)+x^2/t$ is stationary in x (see [Reference Borodin, Corwin and FerrariBCF14, Rem. 1.14]), we may write $G_{b}(t,x,r)=\det (\mathbf {I}-\mathbf {K})_{L^2\left [0,\infty \right )}$ , with

$$ \begin{align*} \mathbf{K}(u,v)=\frac{1}{(2\pi{\textrm{i}})^{2}}\int_{\widetilde{\mathcal{C}}_1} d\xi\int_{\widetilde{\mathcal{C}}_1'} d\eta\,\frac{\pi}{\sin(\pi(\xi-\eta))}\frac{e^{t\xi^3/3-\left(u+r+x^2/t\right)\xi}}{e^{t\eta^3/3-\left(v+r+x^2/t\right)\eta}}\prod_{k=1}^m\frac{\Gamma(\eta-b_k-x/t)}{\Gamma(\xi-b_k-x/t)}. \end{align*} $$

Here the contour $\widetilde {\mathcal {C}}_1$ has to cross the real axis to the right of $b_i+x/t$ for all i. Conjugating the kernel by $e^{ux/t}$ as in the previous case (that is, multiplying the kernel by $e^{(u-v)x/t}$ ) – and changing variables $\eta \longmapsto \eta +x/t$ , $\xi \longmapsto \xi +x/t$ , we get

$$ \begin{align*} G_{b}(t,x,r)=\det\left(\mathbf{I}-\widetilde{\mathbf{K}}\right)_{L^2\left[0,\infty\right)}, \end{align*} $$

with

$$ \begin{align*}\widetilde{\mathbf{K}}(u,v)=\frac{1}{(2\pi{\textrm{i}})^{2}}\int_{\mathcal{C}_1} d\xi\int_{\mathcal{C}_1'} d\eta\,\frac{\pi}{\sin(\pi(\xi-\eta))}\frac{e^{t\xi^3/3+x\xi^2-(u+r)\xi}}{e^{t\eta^3/3+x\eta^2-(v+r)\eta}}\prod_{k=1}^m\frac{\Gamma(\eta-b_k)}{\Gamma(\xi-b_k)}.\end{align*} $$

As in the previous case, $\tilde {\mathbf {K}}$ satisfies the necessary differential relations, so we have the following:

Theorem 3.2. $\phi _b:= \partial _r^2\log G_b$ satisfies the KP-II equation (1.9).

3.3 KPZ equation with two-sided Brownian initial data

In [Reference Le DoussalLD20] (which appeared about a month after the first version of the present article), Le Doussal suggests that a certain modified generating function for the solution ${\textrm {h}}_{w_\pm }$ of equation (1.1) with the same scaling as in the of regoing and initial data of the form ${\textrm {h}}_{w_\pm }(0,x)=\mathbf {B}(x)+w_-x\mathbf {1}_{x<0}+w_+x\mathbf {1}_{x\geq 0}$ with $\mathbf {B}$ a double-sided Brownian motion with $\mathbf {B}(0)=0$ and $w_->w_+$ should also satisfy the KP equation. This is actually true, as we explain next. Define

$$ \begin{align*}\widetilde G_{w_\pm}(t,x,r)=\Gamma(w_--w_+)^{-1}\mathbb{E}\left[2\hspace{0.05em} e^{\frac12\left(w_--w_+\right)\left({\textrm{h}}_{w_\pm}\left(t,x\right)+\frac{t}{12}-r\right)}K_{w_+-w_-}\left(2e^{\frac12\left({\textrm{h}}_{w_\pm}\left(t,x\right)+\frac{t}{12}-r\right)}\right)\right],\end{align*} $$

where $K_\nu $ is the modified Bessel function of order $\nu $ . This modified generating function $\widetilde G_{w_\pm }$ can alternatively be expressed as the analog of equation (3.2) where the KPZ height function ${\textrm {h}}_{w_\pm }$ is replaced by a randomly shifted height function ${\textrm {h}}_{w_\pm }(t,x)+\Upsilon $ :

$$ \begin{align*}\widetilde G_{w_\pm}(t,x,r)=\mathbb{E}\left[\exp\!\left\{-e^{{\textrm{h}}_{w_\pm}\left(t,x\right)+\Upsilon+\frac{t}{12}-r}\right\}\right],\end{align*} $$

with $\Upsilon $ an independent log-gamma random variable with parameter $w_--w_+$ ; that is, $e^{-\Upsilon }$ has density $\Gamma (w_--w_+)^{-1}x^{w_--w_+-1}e^{-x}$ . Explicit formulas for the distribution of this shifted height function were obtained in [Reference Imamura and SasamotoIS12; Reference Imamura and SasamotoIS13] using the nonrigorous replica method. A similar formula, which is the one we will use, was obtained rigorously later on [Reference Borodin, Corwin, Ferrari and VetőBor+15]; the equality between the two expressions for $\widetilde G_{w_\pm }$ is essentially [Reference Borodin, Corwin, Ferrari and VetőBor+15, Corollary 2.6] (see also [Reference Borodin, Corwin, Ferrari and VetőBor+15, Remark 2.10]).

From [Reference Borodin, Corwin, Ferrari and VetőBor+15, Thm. 2.9] (and footnote 4 again), we have $\widetilde G_{w_\pm }(t,x,r)=\det (\mathbf {I}-\mathbf {K})_{L^2\left [0,\infty \right )}$ , with

$$ \begin{align*} &\mathbf{K}(u,v)=\frac{1}{(2\pi{\textrm{i}})^{2}}\int_{\mathcal{C}_t} d\eta\int_{\tilde{\mathcal{C}}_t} d\xi\,\frac{t^{-1/3}\pi}{\sin\left(t^{-1/3}\pi(\xi-\eta)\right)}\frac{e^{\xi^3/3-\left(u+t^{-1/3}\left(r+\frac{x^2}{t}\right)\right)\xi}}{e^{\eta^3/3-\left(v+t^{-1/3}\left(r+\frac{x^2}{t}\right)\right)\eta}}\\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times\frac{\Gamma\left(w_-+\frac{x}{t}-t^{-1/3}\xi\right)}{\Gamma\left(t^{-1/3}\xi-w_+-\frac{x}{t}\right)}\frac{\Gamma\left(t^{-1/3}\eta-w_+-\frac{x}{t}\right)}{\Gamma\left(w_-+\frac{x}{t}-t^{-1/3}\eta\right)}, \end{align*} $$

where $\mathcal {C}_t$ goes from $-\frac 14t^{1/3}-{\textrm {i}}\infty $ to $-\frac 14t^{1/3}+{\textrm {i}}\infty $ crossing the real axis between $t^{1/3}w_+$ and $t^{1/3}w_-$ , and $\tilde {\mathcal {C}}_t$ goes from $\frac 14t^{1/3}-{\textrm {i}}\infty $ to $\frac 14t^{1/3}+{\textrm {i}}\infty $ staying to the right of $\mathcal {C}_t$ and also crossing the real axis between $t^{1/3}w_+$ and $t^{1/3}w_-$ . Changing variables as in the last example leads to

$$ \begin{align*} \widetilde G_{w_\pm}(t,x,r)=\det\left(\mathbf{I}-\widetilde{\mathbf{K}}\right)_{L^2\left[0,\infty\right)}, \end{align*} $$

with

$$ \begin{align*}\widetilde{\mathbf{K}}(u,v)=\frac{1}{(2\pi{\textrm{i}})^{2}}\int_{\mathcal{C}_1} d\eta\int_{\tilde{\mathcal{C}}_1} d\xi\,\frac{\pi}{\sin(\pi(\xi-\eta))}\frac{e^{t\xi^3/3+x\xi^2-(u+r)\xi}}{e^{t\eta^3/3+x\eta^2-(v+r)\eta}}\frac{\Gamma(w_--\xi)}{\Gamma(\xi-w_+)}\frac{\Gamma(\eta-w_+)}{\Gamma(w_--\eta)}.\end{align*} $$

As in the previous examples, $\widetilde {\mathbf {K}}$ satisfies the hypotheses of Theorem 1.3, so we have the following:

Theorem 3.3. $\widetilde \phi _{w_\pm }(t,x,r):=\partial _r^2\log \widetilde G_{w_\pm }(t,x,r)$ satisfies the KP-II equation (1.9).

Our next two examples are related to the setting of Theorem 3.2. Here we abandon the setting of the KPZ equation and go back to random matrix distributions and the KPZ fixed point.

3.4 The Baik–Ben Arous–Péché distribution for spiked random matrices

Consider $G_b$ as in the case of spiked initial data. It is known (see [Reference Borodin, Corwin and FerrariBCF14, Cor. 1.15] for the case $t=1,x=0$ ; the general case follows in the same way or by scaling and shift invariance) that

$$ \begin{align*}\lim_{\varepsilon\to0}G_{\varepsilon^{1/2}b}\left(\varepsilon^{-3/2}t,\varepsilon^{-1}x,\varepsilon^{-1/2}r\right)=F_{\text{BBP},t^{1/3}b-t^{-2/3}x}\left(t^{-1/3}r+t^{-4/3}x^2\right),\end{align*} $$

where $F_{\text {BBP},b}$ is the Baik–Ben Arous–Péché (BBP) distribution arising from spiked (unitarily invariant) random matrices [Reference Baik, Ben Arous and PéchéBBP05]. On the other hand, by the scaling invariance of the KP-II equation (see Remark 1.2), for each fixed $\varepsilon>0$ , $\phi ^\varepsilon _b(t,x,r):=\partial _r^2\log G_{\varepsilon ^{1/2}b}\left (\varepsilon ^{-3/2}t,\varepsilon ^{-1}x,\varepsilon ^{-1/2}r\right )$ satisfies the KP-II equation as well. As a consequence, one expects the following to hold:

Theorem 3.4. Define

$$ \begin{align*}\tilde F_{\text{BBP},b}(t,x,r)=F_{\text{BBP},t^{1/3}b-t^{-2/3}x}\left(t^{-1/3}r+t^{-4/3}x^2\right).\end{align*} $$

Then $\phi _{\text {BBP},b}:=\partial _r^2\log \tilde F_{\text {BBP},b}$ satisfies the KP-II equation (1.9).

This is indeed the case, which can be checked in a similar way as before using Theorem 1.3 and the explicit Fredholm determinant formula for $F_{\text {BBP},b}$ (see, e.g., [Reference Borodin, Corwin and FerrariBCF14, Eqn. (1.4)]).

3.5 KPZ fixed point with half-Brownian initial data

Consider the KPZ fixed point $\mathfrak {h}$ with half-Brownian initial data $\mathfrak {h}_0(x)=\mathbf {B}(-x)$ for $x\leq 0$ and $\mathfrak {h}_0(x)=-\infty $ for $x>0$ , where $\mathbf {B}$ is a Brownian motion started at $0$ with diffusivity $2$ . In the setting of [Reference Matetski, Quastel and RemenikMQR17], this corresponds to starting TASEP with a product measure with density $\frac 12$ on the negative integers and no particles on the positive integers, and known results [Reference Corwin, Ferrari and PéchéCFP10] in this case give $\mathbb {P}(\mathfrak {h}(1,x)\leq r)=F_{\text {BBP},x}\left (r+x^2\right )$ . By the scaling invariance of the KPZ fixed point, in this case for $F_{\text {half-BM}}(t,x,r)=\mathbb {P}_{\mathfrak {h}_0}(\mathfrak {h}(t,x)\leq r)$ this gives

$$ \begin{align*}F_{\text{half-BM}}(t,r,x)=F_{\text{BBP},t^{-2/3}x}\left(t^{-1/3}r+t^{-4/3}x^2\right).\end{align*} $$

Comparing with Theorem 3.4 we deduce the following:

Theorem 3.5. $\phi _{\text {half-BM}}=\partial _r^2\log F_{\text {half-BM}}$ satisfies the KP-II equation (1.9).

Similar statements can be written for more general spiked initial conditions in this setting as well as for multipoint distributions (now in terms of the matrix KP equation as in Theorem 1.1).

4 Derivation of the KP equation

This section contains the proof of Theorem 1.3. After we performed the complicated computation, we discovered that a very similar argument was actually known in the $1$ -dimensional case [Reference PöppePo89]. It is shown there that the Fredholm determinant of a kernel satisfying suitable differential relations solves the Hirota equations. The differential relations turn out to be equivalent to the way the kernel depends on t, x and r in Theorem 1.3. It seems to actually go back to [Reference Zaharov and ŠabatZS74; Reference Zaharov and ŠabatZS79], though it is not explicit there, and has been rediscovered in the literature multiple times.

4.1 The logarithmic derivative

Let

$$ \begin{align*}\Phi(t,r,x)=\partial_r\log\det(\mathbf{I}-\mathbf{K})_{L^2\left(\mathbb{R}^+\right)\oplus\dotsm\oplus L^2\left(\mathbb{R}^+\right)},\end{align*} $$

where $\mathbf {K}$ is a matrix kernel satisfying the assumptions of Theorem 1.3. From now on we omit the subscript on the Fredholm determinant and traces.

Given an operator $\mathbf {A}$ acting on the n-fold direct sum of $L^2(\mathbb {R}^+)$ with matrix kernel $(\mathbf {A}_{ab}(u,v))_{a,b=1}^n$ , we will write $\mathrm {d}_1\mathbf {A}$ and $\mathrm {d}_2\mathbf {A}$ for the operators with kernels, given by

$$ \begin{align*}\mathrm{d}_1\mathbf{A}(u,v)=\big(\partial_u\mathbf{K}_{ab}(u,v)\big)_{a,b=1,\dotsc,n}\quad\text{and}\quad\mathrm{d}_2\mathbf{A}(u,v)=\big(\partial_v\mathbf{K}_{ab}(u,v)\big)_{a,b=1,\dotsc,n}.\end{align*} $$

Note that this is just a notational device; $\mathrm {d}_i$ is not meant to denote an operator. By Theorem 1.3(1),

$$ \begin{align*} \Phi(t,r,x)&=\partial_{r}\hspace{-0.1em}\log(\det(\mathbf{I}-\mathbf{K}))=-\operatorname{{\mathop{\mathrm{{tr}}}}}\left((\mathbf{I}-\mathbf{K})^{-1}\partial_{r}\mathbf{K}\right)\\ &=-\operatorname{{\mathop{\mathrm{{tr}}}}}\left((\mathbf{I}-\mathbf{K})^{-1}(\mathrm{d}_1+\mathrm{d}_2)\mathbf{K}\right)=-\operatorname{{\mathop{\mathrm{{tr}}}}}\left((\mathbf{I}-\mathbf{K})^{-1}\mathrm{d}_1\mathbf{K}+\mathrm{d}_2\left((\mathbf{I}-\mathbf{K})^{-1}\mathbf{K}\right)\right)\\ &=-\operatorname{{\mathop{\mathrm{{tr}}}}}\left(\mathrm{d}_1\mathbf{K}(\mathbf{I}-\mathbf{K})^{-1}+\mathrm{d}_2\left((\mathbf{I}-\mathbf{K})^{-1}\mathbf{K}\right)\right)=-\operatorname{{\mathop{\mathrm{{tr}}}}}\left((\mathrm{d}_1+\mathrm{d}_2)\left((\mathbf{I}-\mathbf{K})^{-1}\mathbf{K}\right)\right)\\ &=-\sum_a\int_0^\infty d\xi\,\partial_\xi\left((\mathbf{I}-\mathbf{K})^{-1}\mathbf{K}\right)_{aa}(\xi,\xi)=\sum_a\left((\mathbf{I}-\mathbf{K})^{-1}\mathbf{K}\right)_{aa}(0,0), \end{align*} $$

where we used the cyclicity of the trace. Introducing the notation

$$ \begin{align*}[\mathbf{A}]=\big(\mathbf{A}_{a,b}(0,0)\big)_{a,b=1,\dotsc,n},\end{align*} $$

this tells us that $\Phi $ can be expressed as an n-dimensional trace,

$$ \begin{align*}\Phi(t,r,x)=\operatorname{{\mathop{\mathrm{{tr}}}}}[\mathbf{R}\mathbf{K}]\qquad\text{with}\quad\mathbf{R}=(\mathbf{I}-\mathbf{K})^{-1}.\end{align*} $$

Note here that $\mathbf {I}-\mathbf {K}$ is invertible because the determinant is nonzero. This proves equation (1.7) (after identifying $\partial _r$ with $\mathcal {D}_r$ as in Section 1.1) with, in the present notation,

$$ \begin{align*} Q=[\mathbf{R}\mathbf{K}]. \end{align*} $$

4.2 Formulas for the partial derivatives

The goal is now is to use Theorem 1.3(1)–(3) to show that algebraically, Q satisfies the matrix KP equation (1.8).

Write $\mathbf {K}'=\partial _r\mathbf {K}$ . Using the general formula $\partial _a(\mathbf {I}-\mathbf {A}(a))^{-1}=(\mathbf {I}-\mathbf {A}(a))^{-1}\partial _a\mathbf {A}(a)(\mathbf {I}-\mathbf {A}(a))^{-1}$ for an operator $\mathbf {A}(a)$ depending smoothly on a parameter a, together with the identity

$$ \begin{align*}\mathbf{K}\mathbf{R}=\mathbf{R}\mathbf{K}=\mathbf{R}-\mathbf{I},\end{align*} $$

which we will use repeatedly, we have

$$ \begin{align*} \partial_rQ&=[\partial_r(\mathbf{R}\mathbf{K})]=[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}+\mathbf{R}\mathbf{K}']=[\mathbf{R}\mathbf{K}'\mathbf{R}]\\ \text{and similarly}\\ \partial_r^2Q&=2[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]+[\mathbf{R}\mathbf{K}''\mathbf{R}],\\ \partial_r^3Q&=6[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]+3[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}'\mathbf{R}]+3[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}''\mathbf{R}]+[\mathbf{R}\mathbf{K}'''\mathbf{R}],\\ \partial_tQ&=[\mathbf{R}\partial_t\mathbf{K}\mathbf{R}]. \end{align*} $$

Next we want to compute $(\partial _rQ)^2$ . Note that, in general,

$$ \begin{align*} ([\mathbf{A}][\mathbf{B}])_{a,b}&=-\sum_{c}\int_0^{\infty}d\eta\,\partial_\eta(\mathbf{A}_{ac}(0,\eta)\mathbf{B}_{cb}(\eta,0))\\ &=-\sum_{c}\int_0^{\infty}d\eta\,\big(\partial_\eta\mathbf{A}_{ac}(0,\eta)\mathbf{B}_{cb}(\eta,0)+\mathbf{A}_{ac}(0,\eta)\partial_\eta\mathbf{B}_{cb}(\eta,0)\big)\\ &=-(\mathrm{d}_2\mathbf{A}\mathbf{B})_{ab}(0,0)-(\mathbf{A}\mathrm{d}_1\mathbf{B})_{ab}(0,0), \end{align*} $$

so that the following integration-by-parts formula holds:

(4.1) $$ \begin{align} [\mathbf{A}][\mathbf{B}]=-[\mathbf{A}\mathrm{d}_1\mathbf{B}+\mathrm{d}_2\mathbf{A}\mathbf{B}]. \end{align} $$

We will use this in the formula

$$ \begin{align*} (\partial_rQ)^2&=([\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]+[\mathbf{R}\mathbf{K}'])^2 \\ &=[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]^2+[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}][\mathbf{R}\mathbf{K}']+[\mathbf{R}\mathbf{K}'][\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]+[\mathbf{R}\mathbf{K}']^2. \end{align*} $$

Using $\mathrm {d}_2(\mathbf {K}_a\mathbf {K}_b)=\mathbf {K}_a\mathrm {d}_2\mathbf {K}_b$ , the first term equals

$$ \begin{align*} -[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_2\mathbf{K}\mathbf{R}&+\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}\mathbf{K}]\\ &=-[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]+[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}\mathbf{K}]. \end{align*} $$

Similarly, the fourth term equals

$$ \begin{align*}-[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}']+[\mathbf{R}(\mathrm{d}_1\mathbf{K}'\mathbf{R}-\mathbf{K}'\mathrm{d}_1\mathbf{R})\mathbf{K}']\end{align*} $$

and the two middle ones equal

$$ \begin{align*} &-[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}']-[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]\\ &\qquad\qquad+[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}']+[\mathbf{R}(\mathrm{d}_1\mathbf{K}'\mathbf{R}-\mathbf{K}'\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}\mathbf{K}]. \end{align*} $$

Using this together with Theorem 1.3(1) and (2), which in the current notation read $\partial _r\mathbf {K}=(\mathrm {d}_1+\mathrm {d}_2)\mathbf {K}$ and $\partial _t\mathbf {K}=-\frac 13\left (\mathrm {d}_1^3+\mathrm {d}_2^3\right )\mathbf {K}$ , we get

$$ \begin{align*} \begin{array}{llr} &12\partial_tQ(t,r,x)+\partial_r^3Q(t,r,x)+6(\partial_rQ(t,r,x))^2&\\ &=-\,4\left[\mathbf{R}\left(\mathrm{d}_1^3+\mathrm{d}_2^3\right)\mathbf{K}\mathbf{R}\right]+6[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]+3[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}'\mathbf{R}] &\quad{-(1)+(2)+(3)}\\ &\qquad+\,3[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}''\mathbf{R}]+\left[\mathbf{R}\left(\mathrm{d}_1^3+\mathrm{d}_2^3+3\mathrm{d}_1^2\mathrm{d}_2+3\mathrm{d}_1\mathrm{d}_2^2\right)\mathbf{K}\mathbf{R}\right]\quad &{+(4)+(5)}\\ &\quad-\,6[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}']+6[\mathbf{R}(\mathrm{d}_1\mathbf{K}'\mathbf{R}-\mathbf{K}'\mathrm{d}_1\mathbf{R})\mathbf{K}']&\quad{-(6)+(7)}\\ &\quad-\,6[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}']-6[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]&\quad{-(8)-(9)}\\ &\quad+\,6[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}']&\quad{+(10)}\\ &\quad+\,6[\mathbf{R}(\mathrm{d}_1\mathbf{K}'\mathbf{R}-\mathbf{K}'\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}\mathbf{K}]&\quad{+(11)}\\ &\quad-\,6[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}]+6[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}\mathbf{K}{}].&\quad{-(12)+(13)} \end{array} \end{align*} $$

The element we have marked (12) equals $6[\mathbf {R}\mathbf {K}'\mathbf {R}\mathbf {K}'\mathbf {R}\mathbf {K}'\mathbf {R}]-6[\mathbf {R}\mathbf {K}'\mathbf {R}\mathbf {K}'\mathbf {R}\mathbf {K}']$ , so

$$ \begin{align*} S_{\text{I}}&:=-\text{(1)}+\text{(2)}+\text{(5)}-\text{(8)}-\text{(12)}=-3\left[\mathbf{R}\left(\mathrm{d}_1^3+\mathrm{d}_2^3-\mathrm{d}_1^2\mathrm{d}_2-\mathrm{d}_1\mathrm{d}_2^2\right)\mathbf{K}\mathbf{R}\right]\\ &=-3\left[\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)^2(\mathrm{d}_1+\mathrm{d}_2)\mathbf{K}\mathbf{R}\right]=-3\left[\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)^2\mathbf{K}'\mathbf{R}\right]. \end{align*} $$

Similarly,

$$ \begin{align*} S_{\text{II}}:=+\text{(3)}+\text{(4)}-\text{(6)}-\text{(9)}=-3[\mathbf{R}\mathbf{K}''\mathbf{R}\mathbf{K}'\mathbf{R}]+3[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}''\mathbf{R}] \end{align*} $$

and

$$ \begin{align*} S_{\text{III}}&:=+\text{(7)}+\text{(10)}+\text{(11)}+\text{(13)}\\ &=+6[\mathbf{R}(\mathrm{d}_1\mathbf{K}'\mathbf{R}-\mathbf{K}'\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}]+6[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}'\mathbf{R}]. \end{align*} $$

So

$$ \begin{align*} \begin{aligned} &4\partial_tQ(t,r,x)+\tfrac13\partial_r^3Q(t,r,x)+2(\partial_rQ(t,r,x))^2=\tfrac13(S_{\text{I}}+S_{\text{II}}+S_{\text{III}})\\ &\quad=-[\mathbf{R}\mathbf{K}''{\mathbf{R}}\mathbf{K}'{\mathbf{R}}]+[\mathbf{R}\mathbf{K}'{\mathbf{R}}\mathbf{K}''{\mathbf{R}}]-\left[\mathbf{R}\left(\mathrm{d}_1^3+\mathrm{d}_2^3-\mathrm{d}_1^2\mathrm{d}_2-\mathrm{d}_1\mathrm{d}_2^2\right)\mathbf{K}\mathbf{R}\right]\\ &\qquad+2[\mathbf{R}(\mathrm{d}_1\mathbf{K}'{\mathbf{R}}-\mathbf{K}'{\mathrm{d}}_1\mathbf{R})\mathbf{K}'{\mathbf{R}}]+2[\mathbf{R}\mathbf{K}'{\mathbf{R}}(\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathbf{K}\mathrm{d}_1\mathbf{R})\mathbf{K}'{\mathbf{R}}]\\ &\quad=-[\mathbf{R}(\mathrm{d}_2-\mathrm{d}_1)\mathbf{K}'{\mathbf{R}}\mathbf{K}'{\mathbf{R}}]+[\mathbf{R}\mathbf{K}'{\mathbf{R}}\mathbf{K}''{\mathbf{R}}]-\left[\mathbf{R}\left(\mathrm{d}_1^3+\mathrm{d}_2^3-\mathrm{d}_1^2\mathrm{d}_2-\mathrm{d}_1\mathrm{d}_2^2\right)\mathbf{K}\mathbf{R}\right]\\ &\qquad-2[\mathbf{R}\mathbf{K}'{\mathbf{R}}\mathrm{d}_1\mathbf{R}\mathbf{K}'{\mathbf{R}}]+2[\mathbf{R}\mathbf{K}'{\mathbf{R}}\mathrm{d}_1\mathbf{K}\mathbf{R}\mathbf{K}'{\mathbf{R}}]. \end{aligned} \end{align*} $$

Now using $-[\mathbf {R}\mathbf {K}'\mathbf {R}\mathrm {d}_1\mathbf {R}\mathbf {K}'\mathbf {R}]+[\mathbf {R}\mathbf {K}'\mathbf {R}\mathrm {d}_1\mathbf {K}\mathbf {R}\mathbf {K}'\mathbf {R}]=-[\mathbf {R}\mathbf {K}'\mathbf {R}\mathrm {d}_1(\mathbf {I}-\mathbf {K})\mathbf {R}\mathbf {K}'\mathbf {R}]$ , which equals $-[\mathbf {R}\mathbf {K}'\mathbf {R}\mathrm {d}_1\mathbf {K}'\mathbf {R}]$ , yields

(4.2) $$ \begin{align} &4\partial_tQ(t,r,x)+\tfrac13\partial_r^3Q(t,r,x)+2(\partial_rQ(t,r,x))^2\nonumber\\ &=-[\mathbf{R}(\mathrm{d}_2-\mathrm{d}_1)\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)\mathbf{K}'\mathbf{R}]-\left[\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)^2\mathbf{K}'\mathbf{R}\right]. \end{align} $$

Remark 4.1. At this stage we can already see that the one-point distribution in the flat case $\mathfrak {h}_0\equiv 0$ satisfies the (integrated) KdV equation. In fact, the arguments in [Reference Matetski, Quastel and RemenikMQR17, Sec. 4.4] lead in this case (in the language of the Brownian scattering operator (1.21)) to $\mathbf {K}^{\operatorname {\mathrm {\text {{hypo}}}}(\mathfrak {h})}_t=\mathbf {I}-\mathbf {U}_t(\mathbf {I}-\varrho )\bar \chi _0(\mathbf {I}-\varrho )\mathbf {U}_{t}^{-1}=\mathbf {U}_t\varrho \hspace {0.05em}\mathbf {U}_{t}^{-1}$ , with $\varrho f(x)=f(-x)$ , which yields $\mathbf {K}(u,v)=(2t)^{-1/3}\operatorname {\mathrm {\text {{Ai}}}}\left ((2t)^{-1/3}(u+v+2r)\right )$ . Hence $\mathbf {K}$ is a Hankel kernel, and thus the right-hand side in equation (4.2) vanishes.

Next we add the derivatives in the $x_i$ variables. As for $\partial _r^2Q$ , we have

(4.3) $$ \begin{align} \partial_x^2Q=2[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]+\left[\mathbf{R}\partial_x^2\mathbf{K}\mathbf{R}\right]. \end{align} $$

On the other hand, if we apply $\partial _r$ to equation (4.2) and use Theorem 1.3(3) to write $(\mathrm {d}_2-\mathrm {d}_1)\mathbf {K}'=\partial _x\mathbf {K}$ , we get

$$ \begin{align*} \begin{split} &-[\mathbf{R}\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\partial_x\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]-2[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}''\mathbf{R}]\hskip-9pt\\ &\quad+2[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]+[\mathbf{R}\mathbf{K}''\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]+[\mathbf{R}\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}'\mathbf{R}]+[\mathbf{R}\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]\hskip-9pt\\ &\quad+[\mathbf{R}\mathbf{K}'\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)\partial_x\mathbf{K}\mathbf{R}]+[\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)\partial_x\mathbf{K}'\mathbf{R}]+[\mathbf{R}(\mathrm{d}_1-\mathrm{d}_2)\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]. \end{split} \end{align*} $$

Note that the first and eighth terms cancel. We want to add $\partial _x^2Q(t,r,x)$ . Since $(\mathrm {d}_1-\mathrm {d}_2)\partial _x\mathbf {K}'=\partial _x\left (\mathrm {d}_1^2-\mathrm {d}_2^2\right )\mathbf {K}=-\partial _x^2\mathbf {K}$ , the next-to-last term in the last expression cancels the second bracket on the right-hand side of equation (4.3). Using additionally $\partial _x\mathbf {K}'+(\mathrm {d}_1-\mathrm {d}_2)\partial _x\mathbf {K}=2\mathrm {d}_1\partial _x\mathbf {K}$ and $-\partial _x\mathbf {K}'+(\mathrm {d}_1-\mathrm {d}_2)\partial _x\mathbf {K}=-2\mathrm {d}_2\partial _x\mathbf {K}$ and writing

$$ \begin{align*} q=\partial_rQ, \end{align*} $$

we deduce that

(4.4) $$ \begin{align} &4\partial_tq+\tfrac13\partial_r^3q+2(q\hspace{0.05em}\partial_rq+\partial_rq\hspace{0.1em} q)+\partial_x^2Q(t,r,x)\nonumber\\ &\quad=2\hspace{0.1em}\Big(\hspace{-0.1em}\!-[\mathbf{R}\mathrm{d}_2\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}]-\tfrac12[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}''\mathbf{R}]\nonumber\\ &\qquad+[\mathbf{R}\mathbf{K}'\mathbf{R}\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]+\tfrac12[\mathbf{R}\mathbf{K}''\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]+[\mathbf{R}\mathbf{K}'\mathbf{R}\mathrm{d}_1\partial_x\mathbf{K}\mathbf{R}]\nonumber\\ &\qquad+[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]\Big). \end{align} $$

We claim that the right-hand side equals two times

(4.5) $$ \begin{align} &-\big( [\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathrm{d}_1\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}] + [\mathbf{R}\partial_x\mathrm{d}_2\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]\big)\nonumber\\ &+ \big( [\mathbf{R}\mathbf{K}'\mathbf{R}\mathrm{d}_2\mathbf{K}\mathbf{R}\partial_x\mathbf{K}\mathbf{R}] + [\mathbf{R}\mathbf{K}'\mathbf{R}\mathrm{d}_1\partial_x\mathbf{K}\mathbf{R}]\big)\nonumber\\ &-\big( [\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathrm{d}_2\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}] + [\mathbf{R}\partial_x\mathbf{K}\mathbf{R} \mathrm{d}_1\mathbf{K}'\mathbf{R}]\big)\nonumber\\ &+ \big( [\mathbf{R}\mathbf{K}'\mathbf{R}\mathrm{d}_1\mathbf{K}\mathbf{R}\partial_x\mathbf{K}\mathbf{R}] + [\mathbf{R}\mathrm{d}_2\mathbf{K}'\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]\big). \end{align} $$

To see this, express the right-hand side of equation (4.4) as $2(r_1+r_2+\dotsb +r_7)$ , express formula (4.5) as $q_1+q_2+\dotsb +q_8$ and note first that $r_1=q_2$ , $r_6=q_4$ , $r_2=q_1+q_5$ and $r_4=q_3+q_7$ . On the other hand, we have $r_3=-\frac 12[\mathbf {R}\partial _x\mathbf {K}\mathbf {R}(\mathrm {d}_1+\mathrm {d}_2)\mathbf {K}'\mathbf {R}]=q_6+\frac 12[\mathbf {R}\partial _x\mathbf {K}\mathbf {R}(\mathrm {d}_1-\mathrm {d}_2)\mathbf {K}'\mathbf {R}]=q_6-\frac 12r_7$ and similarly $r_5=\frac 12[\mathbf {R}(\mathrm {d}_1+\mathrm {d}_2)\mathbf {K}'\mathbf {R}\partial _x\mathbf {K}\mathbf {R}]=q_8+\frac 12[\mathbf {R}(\mathrm {d}_1-\mathrm {d}_2)\mathbf {K}'\mathbf {R}\partial _x\mathbf {K}\mathbf {R}]=q_8-\frac 12r_7$ . This gives $r_3+r_5+r_7=q_6+q_8$ , and finishes proving the claim.

Integrating by parts (i.e., using equation (4.1)) within each parenthesis in formula (4.5), we get

$$ \begin{align*} 2\partial_tq+\tfrac16\partial_r^3q&+(q\hspace{0.05em}\partial_rq+\partial_rq\hspace{0.1em} q)+\tfrac12\partial_x^2Q(t,r,x)\nonumber\\ =&-[\mathbf{R}\partial_x\mathbf{K}(\mathbf{R}\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathrm{d}_1\mathbf{R}) \mathbf{K}'\mathbf{R}] + [\mathbf{R}\partial_x\mathbf{K}][\mathbf{R}\mathbf{K}'\mathbf{R}]\nonumber\\ &\quad + [\mathbf{R}\mathbf{K}'(\mathbf{R}\mathrm{d}_2\mathbf{K}\mathbf{R}-\mathrm{d}_2\mathbf{R})\partial_x\mathbf{K}\mathbf{R}] - [\mathbf{R}\mathbf{K}'\mathbf{R}][\partial_x\mathbf{K}\mathbf{R}]\nonumber\\ &\quad - [\mathbf{R}\partial_x\mathbf{K}(\mathbf{R}\mathrm{d}_2\mathbf{K}\mathbf{R}-\mathrm{d}_2\mathbf{R})\mathbf{K}'\mathbf{R}] + [\mathbf{R}\partial_x\mathbf{K}\mathbf{R}][\mathbf{K}'\mathbf{R}] \nonumber\\ &\quad + [\mathbf{R}\mathbf{K}'(\mathbf{R}\mathrm{d}_1\mathbf{K}\mathbf{R}-\mathrm{d}_1\mathbf{R})\partial_x\mathbf{K}\mathbf{R}] - [\mathbf{R}\mathbf{K}'][\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]. \end{align*} $$

Write this as $s_1+\dotsb +s_8$ . Notice that $s_1+s_5$ yields a term involving

$$ \begin{align*}( \mathbf{R}\mathrm{d}_1\mathbf{K} \mathbf{R}-\mathrm{d}_1\mathbf{R})+ (\mathbf{R}\mathrm{d}_2\mathbf{K}\mathbf{R}-\mathrm{d}_2\mathbf{R}) = \mathbf{R}(\mathbf{K}\mathrm{d}_1\mathbf{K}+\mathrm{d}_2\mathbf{K}\mathbf{K})\mathbf{R}-(\mathrm{d}_1+\mathrm{d}_2)\mathbf{I}, \end{align*} $$

where we have used $\mathbf {R}\mathbf {K}=\mathbf {K}\mathbf {R}=\mathbf {R}-\mathbf {I}$ again, and thus integrating by parts one more time we get

$$ \begin{align*} s_1+ &s_5=-[\mathbf{R}\partial_x\mathbf{K}(\mathbf{R}(\mathbf{K}\mathrm{d}_1\mathbf{K}+\mathrm{d}_2\mathbf{K}\mathbf{K})\mathbf{R}-(\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{K}'\mathbf{R}]\\ &=[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}\mathbf{K}][\mathbf{K}\mathbf{R}\mathbf{K}'\mathbf{R}]+[\mathbf{R}\partial_x\mathbf{K}((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{K}'\mathbf{R}]\\ &=([\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]-[\mathbf{R}\partial_x\mathbf{K}])([\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{K}'\mathbf{R}])+[\mathbf{R}\partial_x\mathbf{K}((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{K}'\mathbf{R}]\\ &=[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}][\mathbf{R}\mathbf{K}'\mathbf{R}]-s_2-s_6+[\mathbf{R}\partial_x\mathbf{K}][\mathbf{K}'\mathbf{R}]+[\mathbf{R}\partial_x\mathbf{K}((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{K}'\mathbf{R}]. \end{align*} $$

In a similar fashion, we get

$$ \begin{align*}s_3+s_7=-[\mathbf{R}\mathbf{K}'\mathbf{R}][\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]-s_4-s_8-[\mathbf{R}\mathbf{K}'][\partial_x\mathbf{K}\mathbf{R}]-[\mathbf{R}\mathbf{K}'((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\partial_x\mathbf{K}\mathbf{R}].\end{align*} $$

Therefore,

$$ \begin{align*} 2\partial_tq&+\tfrac16\partial_r^3q+(q\hspace{0.05em}\partial_rq+\partial_rq\hspace{0.1em} q)+\tfrac12\partial_x^2Q(t,r,x)\\ &\quad=[\mathbf{R}\partial_x\mathbf{K}][\mathbf{K}'\mathbf{R}]+[\mathbf{R}\partial_x\mathbf{K}((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{K}'\mathbf{R}]-[\mathbf{R}\mathbf{K}'][\partial_x\mathbf{K}\mathbf{R}]\\ &\qquad-[\mathbf{R}\mathbf{K}'((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\partial_x\mathbf{K}\mathbf{R}]+[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}][\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\mathbf{K}'\mathbf{R}][\mathbf{R}\partial_x\mathbf{K}\mathbf{R}]. \end{align*} $$

In order to complete the proof we note that if $\mathbf {A}$ and $\mathbf {B}$ are nice kernels, then integrating by parts we get

$$ \begin{align*} [\mathbf{A}((\mathrm{d}_1+\mathrm{d}_2)\mathbf{I})\mathbf{B}]=[\mathbf{A}\mathrm{d}_1\mathbf{B}]+[\mathbf{A}\mathrm{d}_2\mathbf{I}\hspace{0.05em}\mathbf{B}]=-[\mathbf{A}][\mathbf{B}], \end{align*} $$

which immediately yields

$$ \begin{align*} 2\partial_tq+\tfrac16\partial_r^3q+(q\hspace{0.05em}\partial_rq+\partial_rq\hspace{0.1em} q)+\tfrac12\partial_x^2Q(t,r,x)=[\mathbf{R}\partial_x\mathbf{K}\mathbf{R}][\mathbf{R}\mathbf{K}'\mathbf{R}]-[\mathbf{R}\mathbf{K}'\mathbf{R}][\mathbf{R}\partial_x\mathbf{K}\mathbf{R}].\hskip-8pt \end{align*} $$

The right-hand side equals $\partial _xQ\hspace {0.1em} q-q\hspace {0.05em}\partial _xQ$ , so the equation becomes equation (1.8), as needed.

4.3 Proof of Theorem 1.3

For r inside the interval where the trace norm of $\mathbf {K}$ is strictly less than $1$ , everything in the previous section is well defined, because that condition ensures that $(\mathbf {I}-\mathbf {K})^{-1}-\mathbf {I}$ is trace class, and therefore the series for the brackets are convergent. So the algebraic computations hold pointwise for such r.

Since the kernel $\mathbf {K}(u,v)$ is real analytic in $t>0$ , x and r, the Fredholm determinant $\det (\mathbf {I}-\mathbf {K}(t,x,r))$ is as well, since it is given by its Fredholm series, each of whose terms is real analytic and which converges uniformly because $\mathbf {K}$ is trace class uniformly in compact sets of $t,x,r$ . Since the determinant never vanishes, $\partial ^2_r \log \det (\mathbf {I}-\mathbf {K})$ is also real analytic, as are all the terms in the KP equation (either the scalar or matrix version). Therefore the left-hand side of the KP equation (1.8) (or equation (1.9)) is a real analytic function. We have proved that this function vanishes for r in an open interval, and therefore it vanishes everywhere.

Appendix A Multipoint initial data

A.1 $t\to 0$ limit of the Brownian scattering operator

Let the initial data for the KPZ fixed point be a finite collection of narrow wedges $\mathfrak {d}^{\vec b}_{\vec a}$ as in Example 1.6. Fix $x_1<\dotsb <x_m$ . Our goal is to compute the limit

$$ \begin{align*} \lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{d}_{\vec a}^{\vec b}\right)}_te^{x_j\partial^2} \end{align*} $$

in operator norm in $L^2([r,\infty ))$ for any fixed r, with $\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}(\mathfrak {h})}_t$ defined in equation (1.21). The convergence could be upgraded to trace norm in most cases, but it does not hold, for example, in the important case $m=1$ , $x_1=a_1$ , where it is easy to see that the limit itself is not trace class (see equation (A.2)). Nor should it be; one should not expect convergence to the initial data in such a strong sense. What we are aiming for here is just to understand the $t\to 0$ behaviour of the matrix kernel.

Recall the operators $\mathbf {S}_{t,x}$ introduced in equation (1.14), which satisfy the differential relations (1.15). They can be expressed as $\mathbf {S}_{t,x}=e^{-\frac 13t\partial ^3+x\partial ^2}=e^{x\partial ^2}\mathbf {U}_t$ , which means that

$$ \begin{align*}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}(\mathfrak{h})}_t=\lim_{\ell\to\infty}\mathbf{S}_{-t,-\ell}{\textbf{P}}_{-\ell,\ell}^{\text{Hit}\,\mathfrak{h}}\mathbf{S}_{t,-\ell}.\end{align*} $$

We begin by studying the single narrow wedge case, $k=1$ , writing $a=a_1$ , $b=b_1$ . In this case we have, for $\ell>a$ , ${\textbf {P}}_{-\ell ,\ell }^{\text {Hit}\,\mathfrak {d}_a^b}=e^{(a+\ell )\partial ^2}\bar \chi _be^{(\ell -a)\partial ^2}$ by definition of the left-hand side and where $\bar \chi _b$ was defined in equation (1.19). Then

(A.1) $$ \begin{align} \mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{d}_{a}^{b}\right)}_t =\lim_{\ell\to\infty}\mathbf{S}_{-t,-\ell}e^{(a+\ell)\partial^2}\bar\chi_be^{(\ell-a)\partial^2}\mathbf{S}_{t,-\ell}=\mathbf{S}_{-t,a}\bar\chi_b\mathbf{S}_{t,-a}. \end{align} $$

Consider first the case $x_i\leq a\leq x_j$ . Since $\mathbf {S}_{\pm t,x}\longrightarrow e^{x\partial ^2}$ as $t\to 0$ for all $x\geq 0$ in operator norm, by equation (A.1) we get

(A.2) $$ \begin{align} \lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}\left(\mathfrak{d}_{a}^{b}\right)}_te^{x_j\partial^2}=\lim_{t\to0}\mathbf{S}_{-t,a-x_i}\bar\chi_b\mathbf{S}_{t,x_j-a} =e^{\left(a-x_i\right)\partial^2}\bar\chi_be^{\left(x_j-a\right)\partial^2}={\textbf{P}}^{\text{Hit }}\mathfrak{d}_{a}^{b}}_{x_i,x_j} \end{align} $$

in operator norm and in all of $L^2(\mathbb {R})$ .

Next consider the case $x_i\leq x_j$ and $a\notin \left [x_i,x_j\right ]$ . We will show that in this case our operator goes to $0$ as $t\to 0$ . We will assume for simplicity that $a>x_j$ ; the case $a<x_i$ works in the same way. We have $\chi _re^{-x_i\partial ^2}\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {d}_{a}^{b}\right )}_te^{x_j\partial ^2}\chi _r=\left (\chi _r\mathbf {S}_{-t,a-x_i}\bar \chi _b\right )\left (\bar \chi _b\mathbf {S}_{t,x_j-a}\chi _r\right )$ , and the first factor goes to $\chi _r e^{\left (a-x_i\right )\partial ^2}\bar \chi _b$ as $t\to 0$ in operator norm, so it is enough to show that $\bar \chi _b\mathbf {S}_{t,x_j-a}\chi _r$ goes to $0$ . We estimate its Hilbert–Schmidt norm,

(A.3) $$ \begin{align} \left\lVert\bar\chi_b\mathbf{S}_{t,x_j-a}\chi_r\right\rVert_2^2=\int_r^\infty dv\int_{-\infty}^b d\eta\,t^{-2/3}e^{4\bar x_j^3/3t^2-2(\eta-v)\bar x_j/t}\hspace{-0.1em}\hspace{-0.05em}\operatorname{\mathrm{\text{{Ai}}}}\left(t^{-1/3}(v-\eta)+t^{-4/3}\bar x_j^2\right)^2, \end{align} $$

where $\bar x_j=x_j-a$ . Split the $\eta $ integral according to whether $\eta \leq v\wedge b$ or $v\wedge b<\eta \leq b$ . On the first piece we may use the classical bound on the Airy function, $\lvert Ai(s)\rvert \leq C\hspace {0.05em} e^{-\frac 23(s\vee 0)^{3/2}}$ , to see that the integral is bounded by $C\hspace {0.05em} t^{-2/3}e^{\frac 43\bar x_j^3/t^2}\int _r^\infty dv\int _{-\infty }^{v\wedge b}d\eta \,e^{-\frac 43(v-\eta )^{3/2}/t^{1/2}-\frac 43\left \lvert \bar x_j\right \rvert ^3/t^2-2(\eta -v)\bar x_j/t}$ . The exponent in the $\eta $ integral is maximised at the edge of the integration, $\eta =v\wedge b$ , so applying Laplace’s method we deduce that the same integral is bounded by $C\hspace {0.05em} t^{-c}\hspace {0.05em} e^{-\frac 83\left \lvert \bar x_j\right \rvert ^3/t^2}\hspace {-0.1em}\left [\lvert r-b\rvert +\int _b^\infty dv\,e^{-\frac 43(v-b)^{3/2}/t^{1/2}-2(b-v)\bar x_j/t}\right ]$ , for some $c,C>0$ , which clearly goes to $0$ as $t\to 0$ . On the second piece the integration region is bounded, so we get directly (since $\bar x_j<0$ ) that the integral goes to $0$ . This shows that the left-hand side of equation (A.3) goes to $0$ as $t\to 0$ , as desired.

The last possibility is that $x_i>x_j$ . In this case the operator goes to $0$ again as $t\to 0$ . Now one has to estimate the Hilbert–Schmidt norm of the whole operator $e^{-x_i\partial ^2}\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {d}_a^b\right )}_t e^{x_j\partial ^2}=\mathbf {S}_{-t,a-x_i}\bar \chi _b\mathbf {S}_{t,x_j-a}$ on $L^2([r,\infty ))$ ; the estimates are a bit more tedious but very similar to the ones we used in the previous case, so we skip the details.

The conclusion of all this is that, in the case of narrow wedge initial data $\mathfrak {d}_{a}^{b}$ ,

(A.4) $$ \begin{align} \lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{d}_a^b\right)}_t e^{x_j\partial^2}&=\lim_{t\to0}\mathbf{S}_{-t,a-x_i}\bar\chi_b\mathbf{S}_{t,x_j-a}\nonumber\\ &=e^{\left(a-x_i\right)\partial^2}\bar\chi_b e^{\left(x_j-a\right)\partial^2}\mathbf{1}_{x_i\leq a\leq x_j}={\textbf{P}}^{\text{Hit }\mathfrak{d}_a^b}_{x_i,x_j}\mathbf{1}_{x_i\leq x_j} \end{align} $$

in operator norm in $L^2([r,\infty ))$ .

Now we turn to the general case $\mathfrak {h}=\mathfrak {d}_{\vec a}^{\vec b}$ . For $\ell>\lvert a_1\rvert \vee \lvert a_k\rvert $ we have, by inclusion-exclusion, that ${\textbf {P}}_{-\ell ,\ell }^{\text {Hit}\,\mathfrak {d}_{\vec a}^{\vec b}}$ equals

$$ \begin{align*} \sum_{n=1}^k(-1)^{n+1}\!\!\sum_{1\leq p_1<\dotsb<p_n\leq k}\!\!e^{\left(a_{p_1}-\ell\right)\partial^2}\bar\chi_{b_{p_1}}e^{\left(a_{p_2}-a_{p_1}\right)\partial^2}\bar\chi_{b_{p_2}}\dotsm e^{\left(a_{p_n}-a_{p_{n-1}}\right)\partial^2}\bar\chi_{b_{p_n}}e^{\left(\ell-a_{p_n}\right)\partial^2}, \end{align*} $$

so proceeding as in equation (A.1) we get that $e^{-x_i\partial ^2}\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {d}_{\vec a}^{\vec b}\right )}_te^{x_j\partial ^2}$ equals

$$ \begin{align*} \sum_{n=1}^k(-1)^{n+1}\sum_{1\leq p_1<\dotsb<p_n\leq k}\mathbf{S}_{-t,a_{p_1}-x_i}\bar\chi_{b_{p_1}}e^{\left(a_{p_2}-a_{p_1}\right)\partial^2}\bar\chi_{b_{p_2}}\dotsm e^{\left(a_{p_n}-a_{p_{n-1}}\right)\partial^2}\bar\chi_{b_{p_n}}\mathbf{S}_{t,x_j-a_{p_n}}. \end{align*} $$

Each summand can be factored as

$$ \begin{align*} \left(\mathbf{S}_{-t,a_{p_1}-x_i}\bar\chi_{b_{p_1}}\mathbf{S}_{t,0}\right)\left(\mathbf{S}_{-t,a_{p_2}-a_{p_1}}\bar\chi_{b_{p_2}}\mathbf{S}_{t,0}\right)\dotsm\left(\mathbf{S}_{-t,a_{p_n}-a_{p_{n-1}}}\bar\chi_{b_{p_n}}\mathbf{S}_{t,x_j-a_{p_n}}\right). \end{align*} $$

By equation (A.4), as $t\to 0$ the first factor goes to ${\textbf {P}}^{\text {Hit }\mathfrak {d}_{a_{p_1}}^{b_{p_1}}}_{x_i,a_{p_1}}\mathbf {1}_{x_i\leq a_{p_1}}$ , the last factor goes to ${\textbf {P}}^{\text {Hit }\mathfrak {d}_{a_{p_n}}^{b_{p_n}}}_{a_{p_{n-1}},x_j}\mathbf {1}_{a_{p_n}\leq x_j}$ and each of the inner factors goes to ${\textbf {P}}^{\text {Hit }\mathfrak {d}_{a_{p_s}}^{b_{p_s}}}_{a_{p_{s-1}},a_{p_s}}$ , $2\leq s\leq n-1$ . Therefore

$$ \begin{align*} &\lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}}\left(\mathfrak{d}_{\vec a}^{\vec b}\right)}_te^{x_j\partial^2}\\ &\qquad\qquad=\sum_{n=1}^k(-1)^{n+1}\sum_{1\leq p_1<\dotsb<p_n\leq k}{\textbf{P}}^{\text{Hit }\mathfrak{d}_{a_{p_1}}^{b_{p_1}}}_{x_i,a_{p_1}}{\textbf{P}}^{\text{Hit }\mathfrak{d}_{a_{p_2}}^{b_{p_2}}}_{a_{p_1},a_{p_2}}\dotsm{\textbf{P}}^{\text{Hit }\mathfrak{d}_{a_{p_n}}^{b_{p_n}}}_{a_{p_{n-1}},x_j}\mathbf{1}_{x_i\leq a_{p_1},\,x_j\geq a_{p_n}}, \end{align*} $$

and then, using inclusion-exclusion again, we deduce finally that

$$ \begin{align*} \lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}\left(\mathfrak{d}_{\vec a}^{\vec b}\right)}_te^{x_j\partial^2}={\textbf{P}}^{\text{Hit }}\mathfrak{d}_{\vec a}^{\vec b}}_{x_i,x_j}\mathbf{1}_{x_i\leq x_j} \end{align*} $$

in operator norm in $L^2([r,\infty ))$ .

Appendix A.2 Matrix KP initial data

Now we proceed formally. Consider compactly supported initial data $\mathfrak {h}\in \operatorname {\mathrm {\text {UC}}}$ , meaning that $\mathfrak {h}(y)=-\infty $ for y outside some compact interval. Approximating $\mathfrak {h}$ by initial data of the form $\mathfrak {d}_{\vec a}^{\vec b}$ , we obtain

$$ \begin{align*} \lim_{t\to0}e^{-x_i\partial^2}\mathbf{K}^{\operatorname{\mathrm{\text{hypo}}(\mathfrak{h})}_te^{x_j\partial^2}={\textbf{P}}^{\text{Hit }}\mathfrak{h}}_{x_i,x_j}\mathbf{1}_{x_i\leq x_j}. \end{align*} $$

In terms of the extended Brownian scattering operator (1.22), this gives

(A.5)

Remark A.1. This formula formally recovers the correct initial data for the KPZ fixed point: the kernel $\mathbf {K}_{ij}(u,v)$ from equation (1.23) as $t\to 0$ becomes $(\mathbf {K}_0)_{ij}\left (u+r_i,v+r_j\right )$ and, in particular, it is upper triangular, so $\det (\mathbf {I}-\mathbf {K})_{L^2\left (\mathbb {R}^+\right )\oplus \dotsm \oplus L^2\left (\mathbb {R}^+\right )}$ as $t\to 0$ becomes

$$ \begin{align*}\det\!\big(\mathbf{I}-\mathbf{K}_0\big)_{L^2\left(\left[r_1,\infty\right)\right)\oplus\dotsm\oplus L^2\left(\left[r_n,\infty\right)\right)}=\prod_{i=1}^n\det\left(\mathbf{I}-\bar\chi_{\mathfrak{h}\left(x_i\right)}\right)_{L^2\left(\left[r_i,\infty\right)\right)}=\prod_{i=1}^n\mathbf{1}_{r_i\geq\mathfrak{h}\left(x_i\right)},\end{align*} $$

as desired.

Now we compute $[\mathbf {R}_0\mathbf {K}_0]$ with $\mathbf {K}_0$ as in formula (A.5). In Section 4 our kernels acted on the n-fold direct sum of $L^2(\mathbb {R}^+)$ , but here $\mathbf {K}_0$ acts on $L^2([r_1,\infty ))\oplus \dotsm \oplus L^2([r_n,\infty ))$ ; this corresponds in the current setting to evaluating the $(i,j)$ entry of $(\mathbf {I}-\mathbf {K}_0)^{-1}\mathbf {K}_0$ at $\left (r_i,r_j\right )$ . Additionally, it is on this last space that we need to compute compositions of operators; to make this explicit it is convenient to replace the kernel entries $(\mathbf {K}_0)_{ij}$ by $\chi _{r_i}(\mathbf {K}_0)_{ij}\chi _{r_j}$ and simply compute on the n-fold direct sum of $L^2(\mathbb {R})$ . Doing this, and since $\mathbf {K}_0$ is upper triangular, we can expand (formally) the entries of $(\mathbf {I}-\mathbf {K}_0)^{-1}\mathbf {K}_0$ as

(A.6) $$ \begin{align} \left((\mathbf{I}-\mathbf{K}_0)^{-1}\mathbf{K}_0\right)_{ij}=\mathbf{1}_{i\leq j}\sum_{\substack{\pi:i\to j\\\pi\text{ incr.}}}\prod_{\ell=1}^{\lvert\pi\rvert-1}\chi_{r_{\pi(\ell)}}(\mathbf{K}_0)_{\pi(\ell),\pi(\ell+1)}\chi_{r_{\pi(\ell+1)}}, \end{align} $$

where the sum is over nondecreasing paths $\pi $ going from i to j along integers and $\lvert \pi \rvert $ denotes the length of the path. Fix $i\leq j$ and assume first that $r_\ell \geq \mathfrak {h}(x_\ell )$ for each $i\leq \ell \leq j$ . Consider a fixed path $\pi $ from i to j. If $\pi (\ell )=\pi (\ell +1)$ for some $\ell $ , then the corresponding factor in the product inside the sum will be $\chi _{r_{\pi (\ell )}}\bar \chi _{\mathfrak {h}\left (x_{\pi (\ell )}\right )}\chi _{r_{\pi (\ell )}}=0$ , so only strictly increasing paths contribute to the sum and we get

$$ \begin{align*}\left((\mathbf{I}-\mathbf{K}_0)^{-1}\mathbf{K}_0\right)_{ij}=\mathbf{1}_{i<j}\sum_{\substack{\pi:i\to j\\\pi\text{ str. incr.}}}(-1)^{\lvert\pi\rvert-1}\prod_{\ell=1}^{\lvert\pi\rvert-1}\chi_{r_{\pi(\ell)}}{\textbf{P}}^{\text{No hit } \mathfrak{h}}_{x_{\pi(\ell)},x_{\pi(\ell+1)}}\chi_{r_{\pi(\ell+1)}}\end{align*} $$

(note that this sum is now finite). Evaluating at $\left (r_i,r_j\right )$ and applying inclusion-exclusion again, we deduce that, as desired (compare with equation (1.27)),

(A.7) $$ \begin{align} [\mathbf{R}_0\mathbf{K}_0]_{ij} &=-\mathbf{1}_{i<j}\mathbb{P}_{\mathbf{B}\left(x_i\right)=r_i}\Big(\mathbf{B}(y)\geq \mathfrak{h}(y)~\forall \ y\in\left[x_i,x_j\right],\nonumber\\ &\qquad\mathbf{B}(x_\ell)\leq r_\ell\text{ for each }x_\ell\in\left(x_i,x_j\right),\,\mathbf{B}\left(x_j\right)\in dr_j\Big)/dr_j\nonumber\\ &=-\mathbf{1}_{i<j}{\textbf{P}}^{\geq\mathfrak{h},\leq-\mathfrak{d}^{-\vec r}_{\vec x}}_{x_i,x_j}\left(r_i,r_j\right). \end{align} $$

Suppose next that $r_m<\mathfrak {h}(x_m)$ for some $i\leq m\leq j$ , and for simplicity assume that this is the only such index satisfying the condition (the argument can be generalised easily). Assume also that $i<j$ . From the argument in the previous case, we know that if $\pi \!:i\to j$ has a constant piece which stays at any index other than m, then $\pi $ does not contribute to the sum in equation (A.6). Hence any path $\pi $ from i to j which does contribute to the sum can be decomposed as $\pi _1\circ \upsilon \circ \pi _2$ , with $\pi _1\!:i\to m$ and $\pi _2\!:m\to j$ strictly increasing (we allow for $\pi _1$ or $\pi _2$ to be empty if $m=i$ or $m=j$ ), and $\upsilon $ staying at m for a given number of steps (which could be $0$ ). The product inside the sum in equation (A.6) splits between factors coming from the three pieces of the path, and from the middle part we get a factor $\left (\chi _{r_\ell }\bar \chi _{\mathfrak {h}\left (x_m\right )}\chi _{r_m}\right )^{\lvert \upsilon \rvert }=\mathbf {I}\cdot \mathbf {1}_{\lvert \upsilon \rvert =0}+\chi _{r_m}\bar \chi _{\mathfrak {h}\left (x_m\right )}\mathbf {1}_{\lvert \upsilon \rvert>0}$ . In other words, and repeating the previous argument,

$$ \begin{align*} &\left((\mathbf{I}-\mathbf{K}_0)^{-1}\mathbf{K}_0\right)_{ij}=-\!\!\!\!\sum_{\substack{\pi_1:i\to m\\\pi_1\text{ str. incr.}}}\sum_{\substack{\pi_2:m\to j\\\pi_2\text{ str. incr.}}}\!\!\!\sum_{\nu\geq0}\left((-1)^{\lvert\pi_1\rvert-1}\prod_{\ell=1}^{\lvert\pi_1\rvert-1}\chi_{r_{\pi_1(\ell)}}{\textbf{P}}^{\text{No hit } \mathfrak{h}}_{x_{\pi_1(\ell)},x_{\pi_1(\ell+1)}}\chi_{r_{\pi_1(\ell+1)}}\right)\\ &\qquad\times\left(\mathbf{I}\cdot\mathbf{1}_{\nu=0}+\chi_{r_m}\bar\chi_{\mathfrak{h}\left(x_m\right)}\mathbf{1}_{\nu>0}\right)\left((-1)^{\lvert\pi_2\rvert-1}\prod_{\ell=1}^{\lvert\pi_2\rvert-1}\chi_{r_{\pi_2(\ell)}}{\textbf{P}}^{\text{No hit } \mathfrak{h}}_{x_{\pi_2(\ell)},x_{\pi_2(\ell+1)}}\chi_{r_{\pi_2(\ell+1)}}\right)\\ &\quad=-\sum_{\nu\geq0}\chi_{r_i}{\textbf{P}}^{{\geq\mathfrak{h},\leq\mathfrak{d}_{\vec x}^{\vec r}}}_{x_i,x_m}\left(\mathbf{I}\cdot\mathbf{1}_{\nu=0}+\chi_{r_m}\bar\chi_{\mathfrak{h}\left(x_m\right)}\mathbf{1}_{\nu>0}\right){\textbf{P}}^{{\geq\mathfrak{h},\leq\mathfrak{d}_{\vec x}^{\vec r}}}_{x_m,x_j}\chi_{r_j}. \end{align*} $$

But ${\textbf {P}}^{\geq \mathfrak {h},\leq \mathfrak {d}_{\vec x}^{\vec r}}_{x_i,x_m}(u,v)=0$ , because at the endpoint v it requires $\mathfrak {h}(x_m)\leq v\leq r_m$ (the analogous statement holds for the other factor). Hence we conclude that, in this case, $[\mathbf {R}_0\mathbf {K}_0]_{i,j}=0$ , which for the same reason means that equation (A.7) still holds.

Suppose finally that $i=j$ and $r_i<\mathfrak {h}(x_i)$ . Now the only possible paths in equation (A.6) are constant paths $\pi $ of arbitrary length $\lvert \pi \rvert \geq 1$ . Each such $\lvert \pi \rvert $ contributes a term of the form $\chi _{r_i}\bar \chi _{\mathfrak {h}\left (x_i\right )}$ , which evaluated at $(r_i,r_i)$ is taken to be $1$ , and hence $[\mathbf {R}_0\mathbf {K}_0]_{i,i}$ diverges to $\infty $ in this case (which coincides with the physical meaning of this quantity, namely $\partial _{r_i}\hspace {-0.1em}\log F(t,x_1,\dotsc ,x_n,r_1,\dotsc ,r_n)$ ).

The conclusion is then that $Q(0,x_1,\dotsc ,x_n,r_1,\dotsc ,r_n)=[\mathbf {R}_0\mathbf {K}_0]$ satisfies equation (1.26).

Appendix B Alternative derivation of the KP-II equation for narrow wedge multipoint distributions

In this section we will derive the KP-II equation (2.5) for the Airy $_2$ process directly using the path-integral formula for the KPZ fixed point [Reference Matetski, Quastel and RemenikMQR17, Prop. 4.3]. Define $F\left (t,\vec x+y,\vec r+a\right )$ as in the first equality of equation (2.4). Then letting $\mathbf {K}_{t,x}=\mathbf {K}^{\operatorname {\mathrm {\text {hypo}}}\left (\mathfrak {h}_0\right )}_{t}(x,\cdot ;x,\cdot )$ , we have

$$ \begin{align*} F\left(t,\vec x+y,\vec r+a\right)&=\det\left(\mathbf{I}-\mathbf{K}_{t,x_1+y}+\bar\chi_{r_1+a}e^{\left(x_2-x_1\right)\partial^2}\bar\chi_{r_2+a}\dotsm\bar\chi_{r_m+a}e^{\left(x_1-x_m\right)\partial^2}\mathbf{K}_{t,x_1+y}\right)\\ &=\det\left(\mathbf{I}-\mathbf{K}+\bar\chi_{r_1}e^{\left(x_2-x_1\right)\partial^2}\bar\chi_{r_2}\dotsm\bar\chi_{r_m}e^{\left(x_1-x_m\right)\partial^2}\mathbf{K}\right), \end{align*} $$

with $\mathbf {K}=\mathbf {K}_{t,x_1+y}(a+\cdot ,a+\cdot )=e^{a\partial }\mathbf {K}_{t,x_1+y}e^{-a\partial }$ . Note that the product of operators preceding $\mathbf {K}$ in the last term does not depend on t, y or a; call it $\mathbf {I}-\mathbf {P}$ , so that $F=\det (\mathbf {I}-\mathbf {P}\mathbf {K})$ . Up to here this is general, but now we specialise to the narrow wedge case, for which $\mathbf {K}=e^{a\partial }\left (\mathbf {S}_{t,-x_1-y}\right )^*\bar \chi _0\mathbf {S}_{t,x_1+y}e^{-a\partial }$ (see equation (A.1)). Using the cyclic property of the determinant, we get

$$ \begin{align*} F&=\det\left(\mathbf{I}-\bar\chi_0\mathbf{S}_{t,x_1+y}e^{-a\partial}\mathbf{P} e^{a\partial}\left(\mathbf{S}_{t,-x_1-y}\right)^*\bar\chi_0\right)\\ &=\det\left(\mathbf{I}-\bar\chi_0e^{-a\partial}\mathbf{S}_{t,x_1+y}\mathbf{P} \left(\mathbf{S}_{t,-x_1-y}\right)^*e^{a\partial}\bar\chi_0\right)\\ &=\det\left(\mathbf{I}-\chi_0\varrho e^{-a\partial}\mathbf{S}_{t,x_1+y}\mathbf{P} \left(\mathbf{S}_{t,-x_1-y}\right)^*e^{a\partial}\varrho\chi_0\right), \end{align*} $$

with $\varrho $ the reflection operator $\varrho f(x)=f(-x)$ . So letting

$$ \begin{align*}\mathbf{L}=\varrho e^{-a\partial}\mathbf{S}_{t,x_1+y}\mathbf{P} \left(\mathbf{S}_{t,-x_1-y}\right)^*e^{a\partial}\varrho=e^{a\partial}\left(\mathbf{S}_{t,x_1+y}\right)^*(\varrho\mathbf{P}\varrho)\mathbf{S}_{t,-x_1-y}e^{-a\partial}\end{align*} $$

(the second equality is a simple computation), we get that $F=\det (\mathbf {I}-\mathbf {L})$ . Now $\partial _a\mathbf {L}(u,v)=(\partial _u+\partial _v)\mathbf {L}(u,v)$ , $\partial _t\mathbf {L}(u,v)=-\frac 13\left (\partial _u^3+\partial _v^3\right )\mathbf {L}(u,v)$ and $\partial _y\mathbf {L}(u,v)=\left (\partial _u^2-\partial _v^2\right )\mathbf {L}(u,v)$ , which correspond to Theorem 1.3(1), (2) and (3) (except for the change $y\mapsto -y$ , which as in Section 2.3 is irrelevant), so the theorem implies that $\phi =\partial _a^2\log (F)$ solves the KP-II equation in $(t,y,a)$ , and translating back to the $\mathcal {D}_r$ , $\mathcal {D}_x$ derivatives yields equation (2.5).

Acknowledgements

The authors thank P. Le Doussal for pointing out that a claim made in the first version of this paper about the generating function of the KPZ equation with flat initial condition solving the KdV equation was inconsistent with known moment formulas [Reference Le Doussal and CalabreseLDC12]. The first author was supported by the Natural Sciences and Engineering Research Council of Canada and by a Killam research fellowship. The second author was supported by CMM ANID PIA AFB170001, by Programa Iniciativa Científica Milenio grant NC120062 through Nucleus Millenium Stochastic Models of Complex and Disordered Systems, and by Fondecyt Grant 1201914.

Conflict of Interest

None.

Footnotes

Dedicated to the memory of Harold Widom

1 However, see [Reference ProlhacPro20], which appeared on arXiv two days before this article was first posted, where it is shown that particular finite-volume solutions [Reference Baik and LiuBL19] can be written as superpositions of solitons. [Reference Baik, Liu and SilvaBLS20] treats other finite-volume initial conditions.

2 With some work, the growth condition on the initial data can be relaxed to $\mathfrak {h}(x) \le A x^2 + B$ for the problem up to a finite time $t=t(A)$ .

3 In [Reference Kenig, Ponce and VegaKPV97], an odd polynomial with positive leading coefficient is considered as initial data for the KdV equation, which is somewhat in the same spirit.

4 See [Reference Borodin, Corwin and FerrariBCF14, Thm. 1.10], which computes the generating function directly. Compared to that formula, we are changing variables $(t,x)\longmapsto (2t,2x)$ to match the two different scaling conventions for equation (3.1), and using the fact that ${\textrm {h}}_{\text {nw}}(t,x)+x^2/t$ is stationary.

References

Ablowitz, M. J. and Segur, H., ‘On the evolution of packets of water waves’, J. Fluid Mech. 92(4) (1979), 691715.CrossRefGoogle Scholar
Adler, M. and van Moerbeke, P., ‘PDEs for the joint distributions of the Dyson, Airy and sine processes’, Ann. Probab. 33(4) (2005), 13261361.CrossRefGoogle Scholar
Amir, G., Corwin, I. and Quastel, J., ‘Probability distribution of the free energy of the continuum directed random polymer in 1 + 1 dimensions’, Comm. Pure Appl. Math. 64(4) (2011), 466537.CrossRefGoogle Scholar
Baik, J., Ben Arous, G. and Péché, S., ‘Phase transition of the largest eigenvalue for nonnull complex sample covariance matrices’, Ann. Probab. 33(5) (2005), 16431697.CrossRefGoogle Scholar
Baik, J., Buckingham, R. and DiFranco, J., ‘Asymptotics of Tracy-Widom distributions and the total integral of a Painlevé II function’, Comm. Math. Phys. 280(2) (2008), 463497.CrossRefGoogle Scholar
Baik, J. and Liu, Z., ‘Multipoint distribution of periodic TASEP’, J. Amer. Math. Soc. 32(3) (2019), 609674.CrossRefGoogle Scholar
Baik, J., Liu, Z. and Silva, G., ‘Limiting one-point distribution of periodic TASEP’, Ann. Inst. Henri Poincaré Probab. Stat., forthcoming; Preprint, 2020, arXiv:2008.07024.Google Scholar
Borodin, A., Corwin, I. and Ferrari, P., ‘Free energy fluctuations for directed polymers in random media in $1+1$ dimension’, Comm. Pure Appl. Math. 67(7) (2014), 11291214.CrossRefGoogle Scholar
Borodin, A., Corwin, I., Ferrari, P. and Vető, B., ‘Height fluctuations for the stationary KPZ equation’, Math. Phys. Anal. Geom. 18(1) (2015), 20.CrossRefGoogle Scholar
Borodin, A., Ferrari, P. L., Prähofer, M. and Sasamoto, T., ‘Fluctuation properties of the TASEP with periodic initial configuration’, J. Stat. Phys. 129(5-6) (2007), 10551080.CrossRefGoogle Scholar
Calabrese, P., Doussal, P. L. and Rosso, A., ‘Free-energy distribution of the directed polymer at high temperature’, EPL (Europhys. Lett.) 90(2) (2010), 20002.CrossRefGoogle Scholar
Corwin, I., Ferrari, P. L. and Péché, S., ‘Limit processes for TASEP with shocks and rarefaction fans’, J. Stat. Phys. 140(2) (2010), 232267.CrossRefGoogle Scholar
Dauvergne, D., Ortmann, J. and Virág, B., ‘The directed landscape’, Acta Math. forthcoming; Preprint, 2019, arXiv:1812.00309.Google Scholar
Dotsenko, V., ‘Bethe ansatz derivation of the Tracy-Widom distribution for one-dimensional directed polymers’, EPL (Europhys. Lett.) 90(2) (2010), 20003.CrossRefGoogle Scholar
Imamura, T. and Sasamoto, T., ‘Exact solution for the stationary Kardar-Parisi-Zhang equation’, Phys. Rev. Lett. 108(19) (2012), 190603.CrossRefGoogle ScholarPubMed
Imamura, T. and Sasamoto, T., ‘Stationary correlations for the 1D KPZ equation’, J. Stat. Phys. 150(5) (2013), 908939.CrossRefGoogle Scholar
Johansson, K., ‘Shape fluctuations and random matrices’, Comm. Math. Phys. 209(2) (2000), 437476.CrossRefGoogle Scholar
Johansson, K., ‘Discrete polynuclear growth and determinantal processes’, Comm. Math. Phys. 242(1-2) (2003), 277329.CrossRefGoogle Scholar
Kardar, M., Parisi