# Half-Space Stationary Kardar-Parisi-Zhang Equation

Guillaume Barraquand, Alexandre Krajenbrink, Pierre Le Doussal

► To cite this version:

Guillaume Barraquand, Alexandre Krajenbrink, Pierre Le Doussal. Half-Space Stationary Kardar-Parisi-Zhang Equation. Journal of Statistical Physics, 2020, 181, pp.1149 - 1203. 10.1007/s10955-020-02622-z. hal-04974636

HAL Id: hal-04974636

https://hal.science/hal-04974636v1

Submitted on 3 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Half-Space Stationary Kardar–Paris–Zhang Equation

Guillaume Barraquand<fnref n="1" /> · Alexandre Krajenbrink<fnref n="2" /> ID · Pierre Le Doussal<fnref n="1" />

Received: 24 March 2020 / Accepted: 31 July 2020 / Published online: 7 August 2020
© The Author(s) 2020

## Abstract

We study the solution of the Kardar–Paris–Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x, t)$ on the positive half line, equivalently the free energy of the continuum directed polymer in a half space with a wall at $x = 0$. The boundary condition $\partial_x h(x, t)|_{x=0} = A$ corresponds to an attractive wall for $A < 0$, and leads to the binding of the polymer to the wall below the critical value $A = -1/2$. Here we choose the initial condition $h(x, 0)$ to be a Brownian motion in $x > 0$ with drift $-(B + 1/2)$. When $A + B \to -1$, the solution is stationary, i.e. $h(\cdot, t)$ remains at all times a Brownian motion with the same drift, up to a global height shift $h(0, t)$. We show that the distribution of this height shift is invariant under the exchange of parameters $A$ and $B$. For any $A, B > -1/2$, we provide an exact formula characterizing the distribution of $h(0, t)$ at any time $t$, using two methods: the replica Bethe ansatz and a discretization called the log-gamma polymer, for which moment formulae were obtained. We analyze its large time asymptotics for various ranges of parameters $A, B$. In particular, when $(A, B) \to (-1/2, -1/2)$, the critical stationary case, the fluctuations of the interface are governed by a universal distribution akin to the Baik–Rains distribution arising in stationary growth on the full-line. It can be expressed in terms of a simple Fredholm determinant, or equivalently in terms of the Painlevé II transcendent. This provides an analog for the KPZ equation, of some of the results recently obtained by Betea–Ferrari–Occelli in the context of stationary half-space last-passage-percolation. From universality, we expect that limiting distributions found in both models can be shown to coincide.

**Keywords** Kardar–Paris–Zhang · Random matrix theory · Growth process · Stationary measure · Halfspace · Lieb–Liniger

Communicated by Ivan Corwin.

<figure><img src="image_2.png" /></figure>

Alexandre Krajenbrink
akrajenb@sissa.it

<footnote n="1">1 Laboratoire de Physique de l'École Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, 75005 Paris, France</footnote>

<footnote n="2">2 SISSA and INFN, via Bonomea 265, 34136 Trieste, Italy</footnote>
# 1 Introduction

The Kardar–Parisier–Zhang equation [1] describes the stochastic dynamics of the height field, $h(x, t)$, of a growing interface in the continuum, as a function of time, or equivalently, of the free energy of a continuum directed polymer in a random potential as a function of its length. In one space dimension, $x \in \mathbb{R}$, it is at the center of a vast universality class, the KPZ class, which contains numerous solvable discrete models with the same large scale behavior. Examples of such solvable models include random permutations and the associated PNG growth model, interacting particle systems (TASEP, ASEP, and variants), the stochastic six-vertex model and other stochastic vertex models, and several discrete models of directed polymers (DP). There has been many papers on the subject and we refer the reader to the reviews and lecture notes [2–14]. Exact results have also been obtained for the KPZ equation itself, and its equivalent system, the continuous directed polymer [15–32]. These results have been obtained on the full line, $x \in \mathbb{R}$, and have allowed proving the existing conjectures for the scaling exponents of the height fluctuations, $\delta h \sim t^{1/3} \sim x^{1/2}$, and to predict and classify the universal probability distributions which arise for various initial conditions. Most notably, the droplet and flat initial conditions were shown to lead, at large time, to the Tracy–Widom distributions [33,34] for the height at one point (centered and scaled by $t^{1/3}$) associated respectively to the Gaussian unitary and orthogonal ensembles, GUE and GOE, of random matrix theory. Some of these predictions have been successfully tested in experiments [35–40].

It is also interesting for applications [41] to study models in the KPZ class restricted to the half line $x \in \mathbb{R}_+$, for which fewer results are available at present. For some specific coupling to the boundary (at $x = 0$) the solvability properties can sometimes be preserved. One can usually define a parameter, that we will call $A$, and which will be defined below in Eq. (2.3) for the KPZ equation, characterizing this coupling. For the KPZ equation, $A = 0$ and $A = +\infty$ correspond respectively to Neumann and to Dirichlet boundary conditions. The half-line KPZ equation is equivalent to a continuum directed polymer on a half-space with a wall at $x = 0$, which is repulsive for $A > 0$, and attractive for $A < 0$. A remarkable feature of the half-space problem is the existence of a critical value $A_c$ of the parameter $A$ at which a phase transition occurs. In the polymer language it corresponds to a binding of the polymer to the wall if the attraction is strong enough $A < A_c$. The existence of this transition was predicted by Kardar in [42] for the continuum directed polymer, using the replica Bethe ansatz. No prediction for the height distribution was obtained however.

In mathematics, in a pioneering paper in 1999, Baik and Rains [2] proved the existence of a similar transition in the context of the longest increasing sub-sequences (LIS) of symmetrized random permutations. There it was shown that, in the unbound phase, for the droplet initial condition, and near $x = 0$, the (scaled) height fluctuations obey the Tracy–Widom distribution associated to the Gaussian symplectic ensemble (GSE). In the bound phase $A < A_c$, the fluctuations are simply Gaussian. Exactly at the transition, $A = A_c$, fluctuations also obey the Tracy–Widom distribution but the one associated to the Gaussian orthogonal ensemble (GOE). These results were extended in other models, e.g. a GSE-GUE crossover was shown as the endpoint position is varied towards the bulk in the PNG model [43]. For the TASEP in a half-space, equivalent to last passage percolation in a half-quadrant [3], similar results were obtained in [44,45] using Pfaffian-Schur processes [46], in particular concerning the crossover as the parameter controlling the boundary is varied simultaneously with the distance to the boundary. All these half-space models discussed above can be studied via the framework of Pfaffian point processes and random matrix theoretic techniques (these mod-
els are free-fermionic). However these models do not converge to the KPZ equation. The models converging to the KPZ equation such as the asymmetric simple exclusion process (ASEP) or directed polymers are not directly related to Pfaffian point processes (they are non-free-fermionic). Among those, ASEP was studied in [47–49] and the half-space log-gamma polymer was studied in [50–52].

For the KPZ equation in the half-space, a solution for the one-point height distribution near the wall valid at any time, was obtained for $A = +\infty$, for the droplet initial condition using the replica Bethe ansatz in [53] (see also [54]). Another solution (also non-rigorous) was obtained for $A = 0$ in [55], with related but different methods using nested contour integral representations. In both cases the large time limit is found to be GSE Tracy–Widom, consistent with the general picture obtained by Baik and Rains as discussed above. Next, taking the limit from a half-space ASEP, a rigorous solution for the critical case $A = -\frac{1}{2}$ was obtained [49], leading to GOE Tracy–Widom fluctuations. More recently, in [56], a solution valid for any time was found using the replica Bethe ansatz for any $A > -1/2$, and which leads to the GOE Tracy–Widom distribution at the critical point $A = -1/2$. In [57] a solution from the RBA, taking into account bound states, was obtained, in agreement with the results of [56].

A remarkable property of the KPZ equation on the full line is that the stationary measure is the Brownian motion in the sense that if the initial condition $h(x, 0)$ is a two sided Brownian motion (with the appropriate amplitude) the PDF of the height difference between two space points is time independent. This leaves a uniform shift $h(0, t)$ whose fluctuations, scaled by $t^{1/3}$, was shown to follow the so-called Baik–Rains distribution, which is universal over the KPZ class. For the KPZ equation, the solution for all time with Brownian initial condition was found using the RBA in [30] and proved rigorously in [31]. Early investigations on stationary models in the KPZ universality class started with [58] in the context of the polynuclear growth model, which introduced the Baik–Rains distribution as the limiting distribution of height fluctuations. For a very similar model (TASEP), the spatial correlations were investigated in [59]. Outside the class of free fermionic models, besides [30,31] that we have already discussed, let us also mention [60] which proved the one point convergence of ASEP height function towards the Baik–Rains distribution.

The aim of the present paper is to address the same problem, but for the half-line. We study the KPZ equation on the half-line with a boundary parameter $A$ and an initial condition chosen as a unit one-sided Brownian motion with a drift, which we denote for later convenience as $-(B + \frac{1}{2})$. The problem is thus determined by two parameters, $A$ and $B$ and we will study the phase diagram in the $(A, B)$ plane. As we show, on the line $A + B + 1 = 0$ the initial condition is stationary, in the same sense as above, i.e. the PDF of the height difference between two space points remains at all time the one of the same Brownian motion. We show furthermore that the distribution of the height shift, $h(0, t)$, is invariant under the exchange of parameters $A$ and $B$. We will use two methods to obtain the exact generating function which characterizes the distribution of $h(0, t)$, at any time $t$. The first one is the replica Bethe ansatz and is a generalization of the calculation presented in [56]. The second one starts from a known moment formula for the so-called log-gamma polymer [61], and takes the continuum limit to the KPZ equation. The obtained formula is valid for any point in the quadrant $A, B > -1/2$. We then study its large time limit, leading to the phase diagram of Fig. 1. In the $(A, B)$ plane, the point $(A, B) = (-1/2, -1/2)$ plays a special role, as the system is critical with respect to the boundary, and at the same time stationary, and will be called critical stationary. At this point, we show that in the large time limit the fluctuations of the interface are governed by a universal distribution, that we express in terms of a simple Fredholm determinant, equivalently in terms of the Painlevé II transcendent. In a sense, it is
the analog for the half-line problem, of the Baik–Rains distribution for the full line. Inside the quadrant A, $B > -1/2$ the distribution is obtained as the GSE Tracy–Widom distribution and on scales $A + 1/2$, $B + 1/2 \sim t^{-1/3}$, there is a universal two-parameter crossover distribution that we obtain in the quadrant. In the so-called bound phase $A < -1/2$ or $B < -1/2$ away from the critical stationary point (including along the line $A + B + 1 = 0$), the fluctuations are expected to be Gaussian, except when $A = B < -1/2$ where we expect the fluctuations to be distributed as the maximum eigenvalue of a $2 \times 2$ GUE matrix. The crossover to this behavior is however beyond the scope of this paper. Finally, note that our results will be consistent in the limit $B \to +\infty$ with all previous results for the droplet initial condition, in particular with the ones in [56] for $A \ge -1/2$. The level of mathematical rigor of all these results is discussed in Sect. 2.4.

It is important to mention that very recently Betea et al. [62] studied stationary half-space KPZ growth for a discrete model, the last-passage-percolation with exponential weights (i.e. a zero-temperature polymer). They obtained a formula for the asymptotic height distribution, depending on several parameters controlling the distance to the boundary and the position on the line $A + B + 1$ near the critical stationary point. We expect, from the universality within the KPZ class, that our present result and theirs should match. The Pfaffian formula of [62, Theorem 2.7] and our formulae (2.25), (2.27) and (2.28) look different. Although we also have a Pfaffian representation, Eq. (2.20), the associated kernels are different and we do not have a general method to show the equivalence of the Fredholm Pfaffians. A similar issue is discussed in [56, Sect. 4.2, Eq. (74)]. Note that our formula allows for a very easy numerical evaluation of the CDF, given below in Fig. 3, and of the first moments, and allows us to determine tail estimates.

It would be interesting to study the distribution of $h(x, t)$ when $x$ is at a distance of order $t^{2/3}$ from 0, and $A, B$ are scaled close to $-1/2$. This would correspond to varying the parameter $\eta$ in [62]. However, while we can obtain some integral formula for the moments of $Z(x, t)$ in the case $x > 0$, see (4.19) below, we do not expect that it can be rewritten as a Pfaffian formula and the asymptotic analysis would require to develop other methods. We leave this for future consideration.

## Outline

First in Sect. 2 we define the models and make a summary of the main results and formulae obtained in this paper. In the two following sections we compute the moments of the polymer partition sum, i.e. the exponential moments of the KPZ field, by two methods. In Sect. 3 we present the derivation using the Bethe ansatz. In Sect. 4 we obtain the moments starting from the log-gamma polymer, and check that the two moment formulae coincide. In Sect. 5 we obtain, a Pfaffian formula for the Laplace transform generating function by summing up the moments. This leads to our first result, valid for all times and any $A, B > -1/2$, for the generating function as a Fredholm Pfaffian in terms of a matrix kernel. The large time limit of this formula, and of the matrix kernel, is studied in Sect. 6. In Sect. 7 we extend the method described in [54] to obtain a formula for the Laplace transform generating function in terms of a scalar kernel, valid for all times and $A, B > -1/2$. In Sect. 7.2 we perform the large time limit on this scalar kernel, which leads to a two parameter family of interpolating kernels near the point $(A, B) = (-1/2, -1/2)$. From it we obtain various limits, including our formula for the critical stationary distribution.
# 2 Model and Main Results

## 2.1 Model

In this paper we study the KPZ equation, which reads, in dimensionless units

$$ \partial_t h(x, t) = \partial_x^2 h(x, t) + (\partial_x h(x, t))^2 + \sqrt{2} \xi(x, t) \quad (2.1) $$

where $\xi(x, t)$ is the standard space-time white noise, with $E[\xi(x, t)\xi(x', t')] = \delta(x-x')\delta(t-t')$. One introduces, via the Cole-Hopf mapping, the directed polymer partition sum $Z(x, t) = e^{h(x,t)}$, where $h(x, t)$ is solution of the KPZ equation (2.1). It satisfies the multiplicative noise stochastic heat equation (SHE)

$$ \partial_t Z(x, t) = \partial_x^2 Z(x, t) + \sqrt{2} \xi(x, t) Z(x, t) \quad (2.2) $$

understood here with the Ito prescription. Equation (2.2) means that $Z(x, t)$ can be seen as a partition sum over continuum directed paths in the random potential $-\sqrt{2}\xi(x, t)$, with the endpoint at time $t$ fixed at position $x$.

**Definition 2.1** We consider the SHE on the half-line $x \ge 0$ with boundary parameter $A$ and $(x, t) \mapsto Z(x, t)$ to be the solution to (2.2) (it can be shown that the solution is unique, see [63,64] and references therein) with the boundary condition

$$ \partial_x Z(x, t)|_{x=0} = AZ(0, t). \quad (2.3) $$

and with the Brownian initial data, in presence of a drift $-1/2 - B$

$$ Z(x, 0) = e^{B(x)-(1/2+B)x} \quad (2.4) $$

We have shifted by $1/2$ the drift parameter to make more explicit a remarkable symmetry between parameters $A$ and $B$. Indeed, we show (see Claim 4.9) that for any $A, B \in \mathbb{R}$, we have the equality in distribution

$$ Z_A^B(x = 0, t) = Z_B^A(x = 0, t), \text{ for any } t > 0. $$

When $B$ goes to $+\infty$, we recover a result recently proved in [65, Theorem 1.1].

When $A + B + 1 = 0$, the model defined in Definition 2.1 is stationary in the sense that for any fixed time $t$, the spatial process $\{Z(x, t)/Z(0, t)\}_{x \ge 0}$ has the same distribution as $\{Z(x, 0)\}_{x \ge 0}$, that is the exponential of a standard Brownian motion with drift $-1/2 - B$. Equivalently, the distribution of the slope field $\partial_x h(x, t)$ is time-stationary. Let us explain where this condition $A + B + 1 = 0$ comes from. In Sect. 4, we consider a discretization of the KPZ equation, the log-gamma directed polymer. We identify in Sect. 4.2 initial and boundary conditions for the log-gamma polymer which make increments of the partition function stationary in time, using a result from [61] which deals with the full-space case. The log-gamma polymer partition function converges weakly to the stochastic heat equation from Definition 2.1 at high temperature, see details in Sect. 4.3 (note that we provide only a sketch of proof of this result based on the combination of results from [66] to [65]). Hence we may pass to the limit, and taking into account the precise scalings, we obtain that the spatial process $\{Z(x, t)/Z(0, t)\}_{x \ge 0}$ is stationary when $A + B + 1 = 0$.

There may exist other initial conditions for the half-space KPZ equation (or equivalently the stochastic heat equation) such that the slope field $\partial_x h(x, t)$ is time-stationary. Indeed, the KPZ equation arises as a scaling limit of the height function of particle systems such as ASEP for which other stationary distributions exist (see Refs. [67–71]).
Before presenting our main results, let us clarify the meaning of the boundary condition (2.3). As a process in $x$, $Z(x, t)$ has the same regularity as a Brownian motion, hence $\partial_x Z(x, t)$ cannot be associated to a real value. To make sense of (2.3), we say [63] that $Z(x, t)$ is a solution of (2.2) if it satisfies

$$Z(x, t) = \int_{0}^{\infty} p_{t}^{A}(x, y)Z(y, 0)dy + \int_{0}^{\infty} dy \int_{0}^{t} p_{t-s}^{A}(x, y)Z(y, s)\xi(y, s)ds, \quad (2.5)$$

where the last integral is an Ito integral and $p_t^A(x, y)$ is the heat kernel on the positive half-line (i.e. it solves the equation $\partial_t u = \Delta u$ with initial data $\delta_y$) that satisfies the boundary condition

$$\partial_x p_t^A(x, y)|_{x=0} = A p_t^A(0, y), \quad t > 0, y > 0. \quad (2.6)$$

The main consequence that we will use below is that

$$\partial_{x_i} \mathbb{E}[Z(x_1, t) \dots Z(x_n, t)]|_{x_i=0} = A \mathbb{E}[Z(x_1, t) \dots Z(x_n, t)]|_{x_i=0}, \quad 1 \le i \le n, \quad (2.7)$$

which can be obtained by replacing $Z(x_i, t)$ using (2.5) inside the expectation and differentiating with respect to $x_i$.

In terms of directed polymers, $Z(x, t)$ can be represented as a partition sum over directed paths

$$Z(x, t) = \mathbb{E}_{\mathcal{B}} \left[ \int_{0}^{+\infty} dy e^{\mathcal{B}(y)-(B+1/2)y} \int_{x(0)=y}^{x(t)=x} \mathcal{D}x(\tau) e^{-\int_{0}^{t} d\tau \left[ \frac{1}{4} \left(\frac{dx}{d\tau}\right)^2 - \sqrt{2\eta(x(\tau), \tau)} + 2A\delta(x(\tau)) \right]} \right], \quad (2.8)$$

where $\eta$ is a space-time white noise and $\mathcal{D}$ denotes the “measure on paths”; more precisely the path integral is defined as an expectation value over reflected Brownian bridges $x(\tau) \in \mathbb{R}_+$ (reflected at $x=0$) for a given realization of the Brownian initial condition $\mathcal{B}(y)$, followed by an expectation over the Brownian $\mathcal{B}$. The extra $\delta$ interaction ensures the proper boundary condition at $x=0$ for $Z(x, t)$, see [55, Sect. 3.2].

## 2.2 Presentation of the Main Results

Our main results concern the height at $x=0$, $h(0, t)$. In general its large time behavior is expected to be

$$h(0, t) \approx v_{\infty}^{A,B} t + t^{\beta} \chi \quad (2.9)$$

where $\chi$ is an $\mathcal{O}(1)$ random variable, and $\beta$ the growth fluctuation exponent. In the quadrant $A, B \ge -1/2$, to which our exact results are restricted, one has $\beta = 1/3$ and $v_{\infty}^{A,B} = -\frac{1}{12}$. Hence to present these results, we define everywhere the shifted variable

$$H(t) = h(0, t) + \frac{t}{12} \quad (2.10)$$

Note however that in the so-called bound phase, which will not be studied in great detail here, we expect a different value of $v_{\infty}^{A,B}$ with $\beta = 1/2$ and different distributions for $\chi$ (see Sect. 4.6). In the limit $B \to +\infty$, i.e. for the droplet initial condition, it was found [57] that $v_{\infty}^{A,+\infty} = -\frac{1}{12} + (A + \frac{1}{2})^2$ for $A \le -1/2$. In the general $A, B$ case, we expect that $v_{\infty}^{A,B} = -\frac{1}{12} + (\min\{A + \frac{1}{2}, B + \frac{1}{2}, 0\})^2$, based on a heuristic argument presented in Sect. 4.6.
### 2.2.1 Finite Time: Fredholm Pfaffian of Matrix Kernel

Our main result valid for all time $t \ge 0$ and all $A, B > -\frac{1}{2}$ is that the following generating function defined for $\zeta > 0$ can be written as a Fredholm Pfaffian

$$ \mathbb{E}[\exp(-\zeta We^{H(t)})] = 1 + \sum_{n_s=1}^{+\infty} \frac{(-1)^{n_s}}{n_s!} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-r_p}} \mathrm{Pf} [K(r_i, r_j)]_{n_s \times n_s}. \quad (2.11) $$

Here $W$ is a random variable, with an inverse gamma distribution of parameter $A+B+1$, see (3.22), independent from $H(t)$, which enters in the construction of the generating function. The kernel $K$ is matrix valued and represented by a $2 \times 2$ block matrix with elements

$$ \begin{aligned} K_{11}(r, r') &= \iint_{C_2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} G(w)G(z) \cos(\pi w) \cos(\pi z) e^{-rw-r'z+t \frac{w^3+z^3}{3}}, \\ K_{22}(r, r') &= \iint_{C_2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} G(w)G(z) \frac{\sin(\pi w)}{\pi} \frac{\sin(\pi z)}{\pi} e^{-rw-r'z+t \frac{w^3+z^3}{3}}, \quad (2.12) \\ K_{12}(r, r') &= \iint_{C_2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} G(w)G(z) \cos(\pi w) \frac{\sin(\pi z)}{\pi} e^{-rw-r'z+t \frac{w^3+z^3}{3}}, \\ K_{21}(r, r') &= -K_{12}(r', r). \end{aligned} $$

where the dependence in parameters $A, B$ only appears in the function

$$ G(z) = \frac{\Gamma(A + \frac{1}{2} - z) \Gamma(B + \frac{1}{2} - z)}{\Gamma(A + \frac{1}{2} + z) \Gamma(B + \frac{1}{2} + z)} \Gamma(2z) \quad (2.13) $$

and the contour $C$ is an upwardly oriented vertical line parallel to the imaginary axis with real part between 0 and $\min\{A + \frac{1}{2}, B + \frac{1}{2}, 1\}$. The series in (2.11) can also be interpreted as a Fredholm Pfaffian, see Eq. (5.17) and Appendix B. The kernel (2.12) has a similar structure as the kernel defining the GSE Tracy-Widom distribution [72]. It is not entirely obvious that the integrals over the $r_i$ in (2.11) are well-defined, but this is the case. Indeed, one can show that (i) all the entries of $K$ have exponential decay as $r, r'$ go too $+\infty$, using a standard contour shift argument, see e.g. [51, Lemma 6.4] (ii) all entries of $K$ grow at most polynomially with $|r|, |r'|$, which can be shown using a variant of [51, Lemma 7.11].

### 2.2.2 Finite Time: Result in Terms of a Scalar Kernel

The matrix kernel in (2.12) has the structure of a Schur Pfaffian. Following [54] and Appendix B, we are able to express the generating function in terms of the Fredholm determinant of a scalar kernel

$$ \mathbb{E}[\exp(-\zeta We^{H(t)})] = \sqrt{\mathrm{Det}(I - \bar{K}_{t,\zeta})_{\mathbb{R}^2(\mathbb{R}_+)}}. \quad (2.14) $$

were the kernel $\bar{K}_{t,\zeta}$ defined for all $(x, y) \in \mathbb{R}_+^2$ as

$$ \bar{K}_{t, \zeta}(x, y) = 2\partial_x \iint_{C_2} \frac{dw dz}{(2i\pi)^2} G(z)G(w) \frac{\sin(\pi(z-w))}{\sin(\pi(z+w))} e^{w+z} e^{-xz-yw+t \frac{w^3+z^3}{3}} \quad (2.15) $$

where the function $G(z)$ is defined in (2.13). Again this formula is valid for all time $t \ge 0$ and all $A, B > -\frac{1}{2}$. In principle from formula (2.11) or (2.14) the PDF of $H(t)$ for any time $t$ can be extracted, see e.g. [53] for the case $A, B = +\infty$. Here, we only extract the PDF's in the large time limit, as we now discuss.
**Remark 2.2** The kernel $\bar{K}_{t, \zeta}$ can be extended by adding a fictitious variable so that the r.h.s of (2.14) is a $\tau$-function of the Kadomtsev-Petviashvili (KP) equation, see Appendix E.

**Remark 2.3** Other cases where it is possible to transform a matrix-valued kernel into a scalar kernel have been considered in the random matrix literature, see Refs. [34, Sects. II–III] and [73, Sect. 3.1].

### 2.2.3 Large Time Limit and Phase Diagram

The phase diagram in the $(A, B)$ plane in the large time limit is shown in Fig. 1. Qualitatively there are three regions. In the region $A < -1/2$ with $A < B$, we expect, from the results of [57] for $B = +\infty$, that the polymer is “bound to the wall” and that the (scaled) height distribution at large time is Gaussian (see also analogous results for other models in [74], [44, Sect. 6], [51, Sect. 8.1]). By “bound to the wall”, we mean that the polymer path spends most of its time at the boundary and does not significantly venture into the bulk, this phenomenon was predicted by Kardar [42] who studied the depinning of the polymer by the random environment. By symmetry the same can be expected for the region $B < -1/2$ with $B < A$, which corresponds to a polymer “bound to the Brownian”. In the special case where $A = B < -1/2$, the nature of fluctuations is different, there is a competition between the boundary and the initial condition and we expect that the fluctuations have the same distribution as the largest eigenvalue of a $2 \times 2$ GUE matrix, based on heuristic arguments presented in Sect. 4.6. We have, however, no exact formula for the region $A < -1/2$ or $B < -1/2$ called the bound phase. Our exact results concern the third region, the quadrant $A, B \ge -1/2$.

The first result is that for any fixed $A, B > -1/2$, the distribution of the height $H(t)$ converges at large time to the GSE Tracy–Widom distribution $F_4$

$$ \lim_{t \to \infty} \mathbb{P} \left( \frac{H(t)}{t^{1/3}} \le s \right) = F_4(s), \quad A, B > -1/2. \qquad (2.16) $$

When $A = -1/2$, and for any $B > -1/2$ we find that the fluctuations are given the GOE Tracy–Widom distribution.

$$ \lim_{t \to \infty} \mathbb{P} \left( \frac{H(t)}{t^{1/3}} \le s \right) = F_1(s), \quad A = -1/2, B > -1/2. \qquad (2.17) $$

By symmetry the same holds for $B = -1/2$ and any $A > -1/2$. These results are natural to expect. Indeed, when $B > -1/2$, the initial condition has a drift so negative that the asymptotics of the height function should be the same as for the narrow wedge initial data. The limiting distribution then depends on the value of the boundary parameter $A$ according to the Baik–Rains transition discussed in the Introduction, hence the GSE and GOE Tracy–Widom asymptotics.

To study the “critical stationary” point $(A, B) = (-1/2, -1/2)$ we write<fnref n="1" />

$$ A + \frac{1}{2} = t^{-1/3}a, \quad B + \frac{1}{2} = t^{-1/3}b, \qquad (2.18) $$

and consider the large time limit at fixed values of $a, b$. This corresponds to zooming around the critical stationary point as represented in Fig. 2. It is natural to expect – and we indeed

<footnote n="1">¹ Note that the parameter $A$ plays the same role as in [49] and [56] but was denoted $a$ in [55] and $b$ in [57]. The parameter that we denote $a$ was denoted $\epsilon$ in [56].</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Phase diagram indicating the distribution of height fluctuations at large time, as a function of the parameters A, B. The nature of fluctuations in the dashed area around (A, B) = (-1/2, -1/2) is explained in Fig. 2</figcaption></figure>

show – that in the large time limit there is a two parameter family of CDFs $F^{(a,b)}(s)$, indexed by a, b such that

$$ \lim_{t \to \infty} \mathbb{P} \left( \frac{H(t)}{t^{1/3}} \le s \right) := F^{(a,b)}(s). \qquad (2.19) $$

Here, we first obtain a Fredholm Pfaffian formula for the CDF $F^{(a,b)}(s)$ with a, b > 0, by taking the large time limit of the matrix kernel formula (2.11). It reads

$$ F^{(a,b)}(s) = \left(1 + \frac{\partial_s}{a+b}\right) \mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)} \qquad (2.20) $$

where the large time matrix kernel $K^{(a,b)}$ is given in (6.5). Equivalently, a more convenient formula is obtained by taking the large time limit of the scalar kernel formula (2.14). We then obtain the CDF of the one-point KPZ height $H(t)$ in the critical region, in terms of a Fredholm determinant, for a, b > 0

$$ F^{(a,b)}(s) = \left(1 + \frac{\partial_s}{a+b}\right) \sqrt{\mathrm{Det}(I - \bar{K}^{(a,b)})_{\mathbb{L}^2(s, +\infty)}} \qquad (2.21) $$
<figure><img src="image_4.png" /><figcaption>Fig. 2 Zoom into the vicinity of (A, B) = (-1/2, -1/2). The distribution of height fluctuations at large time is indicated as a function of parameters a = t^{1/3}(A + 1/2), b = t^{1/3}(B + 1/2)</figcaption></figure>

where the scalar transition kernel $\bar{K}^{(a,b)}$ takes the form

$$ \bar{K}^{(a,b)}(x, y) = \int_{0}^{+\infty} d\lambda A^{(a,b)}(x+\lambda) A^{(a,b)}(y+\lambda) - \frac{1}{2} A^{(a,b)}(x) \int_{0}^{+\infty} d\lambda A^{(a,b)}(y+\lambda). \quad (2.22) $$

Here the function $A^{(a,b)}(x)$ is defined by the integral representation

$$ A^{(a,b)}(x) = \int \frac{dz}{2i\pi} \frac{a+z}{a-z} \frac{b+z}{b-z} e^{-xz+z_s^3}, \quad (2.23) $$

where the contour is a upwardly oriented vertical line with real part between 0 and $\min\{a, b\}$. Finally, introducing the operator $\hat{A}_s$ with kernel $\hat{A}_s(x, y) = A^{(a,b)}(x + y + s)$, the final and simplest expression for the cross-over CDF obtained from an algebraic manipulations of (2.21) is

$$ F^{(a,b)}(s) = \frac{1}{2} \left( 1 + \frac{\partial_s}{a+b} \right) \left( \mathrm{Det}(I - \hat{A}_s)_{\mathbb{L}^2(\mathbb{R}_+)} + \mathrm{Det}(I + \hat{A}_s)_{\mathbb{L}^2(\mathbb{R}_+)} \right) \quad (2.24) $$

It is clear on this formula that if $a, b \to +\infty$ simultaneously, then $A^{(a,b)}(x)$ converges to the standard Airy function, and $\bar{K}^{(a,b)}$ to the kernel associated to the GSE Tracy-Widom distribution (in the form found in [53]). This thus matches smoothly with the result (2.16) valid for any fixed A, B > -1/2. Another interesting limit, that we call $F^{(a)}(s) = \lim_{b \to +\infty} F^{(a,b)}(s)$, is the limit $b \to +\infty$ at fixed a, which corresponds to the droplet initial condition in the critical region for the wall parameter. A formula for that CDF was obtained, for $a \ge 0$, using the RBA in [56]. It was conjectured to coincide with the GSE-GOE-Gaussian crossover introduced by Baik and Rains [74], see also [44,45] in the context of last passage percolation. This crossover was also studied in the context of spiked models of random matrices from the GSE [75]. By the $A \leftrightarrow B$ symmetry, the case $A \to +\infty$ is similar, and we obtain the same distribution $F^{(b)}(s) = \lim_{a \to +\infty} F^{(a,b)}(s)$, which corresponds to
the model with Brownian initial data in presence of an infinite repulsive wall (see [64] for a more mathematical interpretation). This is consistent with Tracy–Widom GOE fluctuations for any fixed $A = -1/2$ and $B > -1/2$ ($a = 0, b = \infty$) or fixed $B = -1/2$ and $A > -1/2$ ($b = 0, a = \infty$).

A more difficult limit, which we discuss now, is the stationary critical point $(A, B) = (-1/2, -1/2)$ corresponding to both $a, b \to 0$.

## 2.3 Stationary Critical Distribution

At the point $(A, B) = (-1/2, -1/2)$ we have found a remarkable universal distribution, corresponding to the CDF $F(s) := F^{(0,0)}(s)$. Taking the limit $a, b \to 0$ is delicate (as it is also to obtain the Baik–Rains distribution) and we found that the representation of the kernel $\tilde{K}^{(a,b)}$ as in (2.22) was crucial. The limit is performed in Sect. 7.3. We have shown that the limit is well defined, i.e. independent of the ratio $r = b/a$. We have obtained the result in several equivalent forms. We recall that we are characterizing the CDF $F(s)$ such that

$$ \lim_{t \to \infty} \mathbb{P} \left( \frac{H(t)}{t^{1/3}} \le s \right) = F(s), \quad A = B = -1/2. \qquad (2.25) $$

The first form is in terms of the sum of two Fredholm determinants. Defining the following two kernels acting on functions in $L^2(\mathbb{R}_+)$,

$$ \mathrm{Ai}_s(x, y) = \mathrm{Ai}(s + x + y), \quad \tilde{\mathrm{Ai}}_s(x, y) = \mathrm{Ai}(s + x + y) + \int_0^{+\infty} d\lambda \, \mathrm{Ai}(s + x + \lambda), \qquad (2.26) $$

then

$$ F(s) = \partial_s \left[ 2\mathrm{Det}(I + \tilde{\mathrm{Ai}}_s) + (s-2)\mathrm{Det}(I + \mathrm{Ai}_s) \right]. \qquad (2.27) $$

The second form is expressed in terms of the CDF's of the GOE and GUE Tracy–Widom distributions $F_1$ and $F_2$ respectively, as

$$ F(s) = \partial_s \left[ \frac{F_2(s)}{F_1(s)} \int_{-\infty}^{s} dt \frac{F_1(t)^4}{F_2(t)^2} \right]. \qquad (2.28) $$

which is very reminiscent of the formula for the Baik–Rains distribution for the full space stationary problem [recalled in (C.12)]. The third form is expressed in terms of the Hastings–McLeod solution $q(s)$ to the Painlevé II equation as

$$ F(s) = \partial_s \left[ e^{-\frac{1}{2} \int_s^{+\infty} dr [(r-s)q^2(r) - q(r)]} \int_{-\infty}^{s} dr e^{-2 \int_r^{+\infty} dt q(t)} \right]. \qquad (2.29) $$

The first moments and cumulants are given in the Table 1 and we plot in Fig. 3, the CDF $F$ together with its derivative, the PDF. Finally, we computed and plotted in Appendix D the asymptotics of the CDF $F(s)$

* for large positive $s$

$$ \begin{aligned} 1 - F(s) = & \frac{s^{3/4} e^{-\frac{2s^{3/2}}{3}}}{4\sqrt{\pi}} \left[ 1 + \frac{139s^{-3/2}}{48} - \frac{11423s^{-3}}{4608} + \frac{3907027s^{-9/2}}{663552} \right. \\ & \left. - \frac{2886147455s^{-6}}{127401984} + o(s^{-6}) \right] \end{aligned} \qquad (2.30) $$
<table><caption>**Table 1** Mean, variance, skewness and excess kurtosis of the half-space critical stationary distribution and comparison with the Tracy–Widom and Baik–Rains distributions (see [76, Sect. 9.4.1] and [77])</caption><thead><tr><th>Distribution</th><th>Mean</th><th>Variance</th><th>Skewness</th><th>Excess kurtosis</th></tr></thead><tbody><tr><td>Half-space stationary</td><td>0</td><td>1.649</td><td>0.266</td><td>0.134</td></tr><tr><td>Tracy–Widom β = 1</td><td>-1.2065...</td><td>1.6078...</td><td>0.2935...</td><td>0.1652...</td></tr><tr><td>Tracy–Widom β = 2</td><td>-1.7711...</td><td>0.8132...</td><td>0.2241...</td><td>0.0934...</td></tr><tr><td>Tracy–Widom β = 4</td><td>-2.3069...</td><td>0.5177...</td><td>0.1655...</td><td>0.0492...</td></tr><tr><td>Baik–Rains</td><td>0</td><td>1.1504...</td><td>0.3594...</td><td>0.2892...</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Left: Critical stationary CDF *F*. Right: corresponding PDF. See Fig. 6 in Appendix. D for the comparison with the asymptotics ($s \to \pm\infty$) in true and logarithmic scales</figcaption></figure>

and for large negative *s*

$$F(s) = 2^{-203/48} e^{\xi'(-1)/2} \exp \left[ -\frac{|s|^3}{24} - \frac{|s|^{3/2}}{\sqrt{2}} + \frac{23}{16} \log |s| + \frac{91}{8\sqrt{2}|s|^{3/2}} \right. \\ \left. - \frac{3957}{128|s|^3} + \frac{28717}{128\sqrt{2}|s|^{9/2}} - \frac{469683}{512|s|^6} + o(s^{-6}) \right] \quad (2.31)$$

As we mentioned above, it remains to be shown that our formula is equivalent to the Fredholm Pfaffian formula obtained in [62, Theorem 2.7] (setting there $\delta = 0$ and $u = 0$) as expected from universality.

## 2.4 Mathematical Aspects

The results presented in this article rely on a combination of physics and mathematics methods, but we focus in this article on physics results and make clear here that most of our results are not proved according to the standards of rigor of the mathematics literature (in particular the results stated as Claims below). It remains a challenge to turn the arguments that we present here into mathematical theorems. Let us comment further on these aspects for the mathematically inclined reader.

The first difficulty from the mathematical point of view is that we cannot rigorously characterize the distribution of the KPZ equation through its moments, because they grow too fast to uniquely determine the distribution. This is why the moment generating series that we consider in Sect. 5 are actually divergent series, but the formal power series become convergent after certain manipulations and exchanges of series/integrals. For the full-space
KPZ equation, it has been proved (see e.g. [55]) that these manipulations lead to the correct answer. It could be possible to overcome this issue in our case by working on a model for which the moment problem is well-defined, and take a scaling limit to the KPZ equation. Such a strategy has been implemented for instance in [22,31,78,79] for the full-space KPZ equation and in [49] for the half-space KPZ equation with $A = -1/2$ and droplet initial data. Another possible approach is provided by the framework of half-space Macdonald processes [51] which allows to prove Laplace transform formulae despite the divergence of moments.

The second obstacle is that in order to prove the results of Sect. 3 below, one would need to prove the completeness of the Bethe ansatz eigenfunctions. Actually, we present in Sect. 4 another approach to obtain the same moment formulae. It relies on rigorous formulae for the log-gamma polymer from [51], and we take a scaling limit to the KPZ equation. We obtain a nested contour integral formulae for the moments of $Z_A^B(x, t)$ in Claim 4.7. Note that this formula allows to take $x > 0$. Then, for $x = 0$, we may move the contours together appealing to a combinatorial conjecture from Borodin, Bufetov and Corwin [55, Conjecture 5.2], and a Pfaffian structure appears. Hence we see that assuming completeness of Bethe ansatz eigenfunctions or using this conjecture leads to the same moment formula. The results from [80,81] suggest that the two problems are indeed related.

Finally, the asymptotic analysis of Fredholm Pfaffians such as (2.11) is delicate, especially in the critical stationary regime. In particular, we have first performed a large time limit for positive $a, b$, and then let the parameters $a, b$ go to 0. This allowed us to benefit from the structure of the kernel (2.22), which eventually lead to a very simple formula for the distribution $F(s)$. It would be interesting to find a generalization of the form (2.22) at finite time, and prove that the limits commute, i.e. one can take first $A, B \to -1/2$ and then study the large time limit.

# 3 Moments from the Replica Bethe Ansatz

## 3.1 Quantum Mechanics and Bethe Ansatz

In this Section we use the replica Bethe ansatz method to calculate the integer moments of the partition sum. The equal time multi-point moments of the solution of the SHE, $Z(x, t)$, over the KPZ noise can be expressed [82] as a matrix element of the quantum mechanical evolution operator in imaginary time of the Lieb–Liniger model [83]

$$ \mathbb{E}[Z(x_1, t) \dots Z(x_n, t)] = \langle x_1 \dots x_n | e^{-tH_n} | \Psi(t=0) \rangle \quad (3.1) $$

Here $H_n$ is the Hamiltonian of the Lieb Liniger model [83] for $n$ quantum particles with attractive delta function interactions of strength $c = -\bar{c} < 0$

$$ H_n = - \sum_{i=1}^{n} \partial_{x_i}^2 - 2\bar{c} \sum_{1 \le i < j \le n} \delta(x_i - x_j) \quad (3.2) $$

with here an below, in our units $\bar{c} = 1$. The initial state $|\Psi(t=0)\rangle$ is such that

$$ \mathbb{E}[Z(x_1, 0) \dots Z(x_n, 0)] = \langle x_1 \dots x_n | \Psi(t=0) \rangle \quad (3.3) $$
Since here we are considering the Brownian initial condition and interested in averages both over the Brownian and the KPZ noise we must take the initial state $|\Psi(t=0)\rangle$ as

$$ \langle x_1 \dots x_n | \Psi(t=0) \rangle = \Phi_0(x_1, \dots, x_n) := \mathbb{E}_B \left[ \exp \left( \sum_{j=1}^n \mathcal{B}(x_j) - (B+1/2)x_j \right) \right] \quad (3.4) $$

A simple calculation shows that $\Phi_0(x_1, \dots, x_n)$ is the fully symmetric function which in the sector $0 \le x_1 < \dots < x_n$ takes the form

$$ \Phi_0(x_1, \dots, x_n) = \exp \left( \sum_{j=1}^{n} \frac{1}{2} (2n - 2j + 1)x_j - (B + 1/2)x_j \right), \quad (3.5) $$

We can now rewrite (3.1) at coinciding points using the decomposition of the evolution operator $e^{-tH_n}$ in terms of the eigenstates of $H_n$ as

$$ \mathbb{E}[Z(x, t)^n] = \sum_{\mu} \Psi_{\mu}(x, \dots, x) \langle \Psi_{\mu} | \Phi_0 \rangle \frac{1}{||\mu||^2} e^{-tE_{\mu}} \quad (3.6) $$

Here the un-normalized eigenfunctions of $H_n$ are denoted $\Psi_\mu$ (of norm denoted $||\mu||$) with eigenenergies $E_\mu$. Here we used the fact that only symmetric (i.e. bosonic) eigenstates contribute since the initial and final states are fully symmetric in the $x_i$. Hence the $\sum_\mu$ denotes a sum over all bosonic eigenstates of the Lieb–Liniger model, also called delta Bose gas, and $\langle \Psi_\mu | \Phi_0 \rangle$ denotes the overlap, i.e. the Hermitian scalar product of the initial state (3.5) with the eigenstate $\Psi_\mu$.

We should remember now that $H_n$ is defined on the half-line $x \ge 0$. The boundary condition at the wall with parameter A translates into the same boundary condition for the wavefunctions (in each of their coordinate). This half-line quantum mechanical problem can be solved by the Bethe ansatz for $A = +\infty$, i.e. for Dirichlet boundary condition [47,84,85] (see also [86, Sect. 5.1]) and this fact was used in [87]. It can also be solved for arbitrary A, [55,86,88–93] which led to the moment formula in [94], [56] and [57].

From the Bethe ansatz the eigenstates $\Psi_\mu$ are thus Bethe states, i.e. superpositions of plane waves over all permutations $P$ of the $n$ rapidities $\lambda_j$ for $j \in [1, n]$ with an additional summation over opposite pairs $\pm \lambda_j$ due to the infinite hard wall. The bosonic (fully symmetric) eigenstates can be obtained everywhere from their expression in the sector $0 \le x_1 \le \dots \le x_n$, which reads

$$ \begin{aligned} \Psi_{\mu}(x_1, \dots, x_n) &= \frac{1}{(2\mathbf{i})^n} \sum_{P \in S_n} \prod_{p=1}^{n} \left( \sum_{\varepsilon_p = \pm 1} \varepsilon_p e^{i\varepsilon_p x_p \lambda_{P(p)}} A[\varepsilon_1 \lambda_{P(1)}, \varepsilon_2 \lambda_{P(2)}, \dots, \varepsilon_n \lambda_{P(n)}] \right) \\ A[\lambda_1, \dots, \lambda_n] &= \prod_{n \ge \ell > k \ge 1} \left( 1 + \frac{\mathbf{i}\bar{c}}{\lambda_\ell - \lambda_k} \right) \left( 1 + \frac{\mathbf{i}\bar{c}}{\lambda_\ell + \lambda_k} \right) \prod_{\ell=1}^{n} \left( 1 + \frac{\mathbf{i}\lambda_\ell}{A} \right) \end{aligned} \quad (3.7) $$

This wavefunction automatically satisfies both

1. The matching condition arising from the $\delta(x_i - x_j)$ interaction

$$ (\partial_{x_{i+1}} - \partial_{x_i} + \bar{c}) \Psi_{\mu}(x_1, \dots, x_n) \Big|_{x_{i+1}=x_i=0} = 0 \quad (3.8) $$

2. The boundary condition $\partial_{x_i} \Psi_\mu(x_1, \dots, x_n)|_{x_i=0} = A \Psi_\mu(x_1, \dots, x_n)|_{x_i=0}$ for all $i \in [0, n]$.
The allowed values for the rapidities $\lambda_i$, which parametrize the true physical eigenstates are determined by the Bethe equations arising from the boundary conditions at $x = L$ as discussed below. One will find that the normalized eigenstates $\psi_\mu = \Psi_\mu/||\mu||$ vanish as $(\lambda_i - \lambda_j)$ or $(\lambda_i + \lambda_j)$ when two rapidities become equal or opposite: hence the rapidities obey an exclusion principle.

The detailed Bethe equations, which determine the allowed values for the set of rapidities $\{\lambda_j\}$, depend on the choice of boundary condition at $x = L$. However, in the $L \to +\infty$ limit, these details do not matter. For simplicity we choose a hardwall at $x = L$. The Bethe equations then read

$$e^{2i\lambda_j L} = \frac{A - i\lambda_j}{A + i\lambda_j} \prod_{\ell \neq j} \frac{\lambda_j - \lambda_\ell - i\bar{c} \lambda_j + \lambda_\ell - i\bar{c}}{\lambda_j - \lambda_\ell + i\bar{c} \lambda_j + \lambda_\ell + i\bar{c}} \quad (3.9)$$

In the case of the infinite hardwall, these equations are also given in Ref. [85] and their solutions in the large $L$ limit were studied in Ref. [95]. The structure of the states for infinite $L$ is found similar to the standard case, i.e. the general eigenstates are built by partitioning the $n$ particles into a set of $n_s$ bound-states formed by $m_j \ge 1$ particles with $n = \sum_{j=1}^{n_s} m_j$. Each bound state $\mu$ is indexed by a set of $\{k_j, m_j\}_{j=1...n_s}$ where the $k_j$'s are real numbers. These states are perfect strings [96], i.e. a set of rapidities

$$\lambda^{j,a} = k_j + \frac{i\bar{c}}{2}(m_j + 1 - 2a) \quad (3.10)$$

where $a = 1, ..., m_j$ labels the rapidities within the string. Such eigenstates have momentum and energy

$$K_{\mu} = \sum_{j=1}^{n_s} m_j k_j, \quad E_{\mu} = \sum_{j=1}^{n_s} m_j k_j^2 - \frac{\bar{c}^2}{12} m_j (m_j^2 - 1). \quad (3.11)$$

The ground-state corresponds to a single $n$-string with $k_1 = 0$. The difference with the standard case is that the states are now invariant by a sign change of any of the momenta $\lambda_j \to -\lambda_j$, i.e. $k_j \to -k_j$. From now on, we will denote the wavefunctions of the string states as $\Psi_{\{k_\ell, m_\ell\}}$.

It is important to note that although for $A = +\infty$ the strings are the only solutions of the Bethe equations at large $L$, for finite $A$ there are other solutions which correspond to so-called boundary bound states. These solutions have been obtained and studied in details in [57]. As we will see below we will not need them in this work.

## 3.2 Moment Formula

To calculate the $n$-th moments of $Z(x, t)$ from formula (3.6), we need to perform a summation over the eigenstates. For $A < +\infty$ these eigenstates contain both the string states and the boundary bound states mentioned above. Our strategy here will be similar to the one in [56], i.e. we will calculate the moments for $n < 2A + 1$ which turn out to be sufficient to perform the analytic continuation in $n$ and obtain the generating function for any $A > -1/2$, using a method similar to the one in [30]. The nice feature is that when $n < 2A + 1$ there are no boundary bound states. To see that, consider the Table 1 in [57] which contains the classification of the boundary bound states for this problem. For $A > -1/2$, a $m$ particle boundary bound state must obey $m \ge [2A] + 2 \ge 2A + 1$. On the other hand from $n < 2A + 1$, one must have $m \le n < 2A + 1$, which excludes the bound state. Hence we need to consider only the string states.
A formula for the inverse of the squared norm of an arbitrary string state was obtained for $A < +\infty$ in [94] and [57], consistent with the results of [56], as

$$ 
\begin{align} 
\|\mu\|^2 &:= \int_0^L dx_1 \dots \int_0^L dx_n |\Psi_{\{k_\ell, m_\ell\}}(x_1, \dots, x_n)|^2 \nonumber \\
\frac{1}{||\mu||^2} &= \frac{1}{n!} \bar{c}^{n-n_s} 2^{n_s} \prod_{i=1}^{n_s} S_{k_i, m_i} H_{k_i, m_i} \prod_{1 \le i < j \le n_s} D_{k_i, m_i, k_j, m_j} L^{-n_s} \nonumber \\
D_{k_1, m_1, k_2, m_2} &= \left( \frac{4(k_1 - k_2)^2 + (m_1 - m_2)^2 c^2}{4(k_1 - k_2)^2 + (m_1 + m_2)^2 c^2} \right) \times \left( \frac{4(k_1 + k_2)^2 + (m_1 - m_2)^2 c^2}{4(k_1 + k_2)^2 + (m_1 + m_2)^2 c^2} \right) \nonumber \\
S_{k,m} &= \frac{2^{2m-2}}{m^2} \prod_{p=1}^{[m/2]} \frac{4k^2 + c^2(m-2p)^2}{4k^2 + c^2(m+1-2p)^2} \nonumber \\
H_{k,m} &= \prod_{a=1}^m \frac{A^2}{A^2 + (k + \frac{i\bar{c}}{2}(m+1-2a))^2} \tag{3.12}
\end{align} 
$$

with $S_{k,1} = 1$. Note that we have only kept the leading term in $L$ as $L \to +\infty$. Inserting the norm formula (3.12) into (3.6), we obtain the starting formula for the integer moments of the partition sum with Brownian weight on the endpoint in the limit $L \to +\infty$

$$ 
\begin{align} 
\mathbb{E}[Z(x, t)^n] &= \sum_{n_s=1}^n \frac{2^{n_s} \bar{c}^n}{n_s! \bar{c}^{n_s} n!} \prod_{p=1}^{n_s} \sum_{m_p \ge 1} \int_{\mathbb{R}} \frac{dk_p}{2\pi} m_p S_{k_p, m_p} H_{k_p, m_p} e^{(m_p^3 - m_p) \frac{\bar{c}^2 t}{12} - m_p k_p^2 t} \tag{3.13} \\
&\quad \times \delta_{n, \sum_{j=1}^{n_s} m_j} \prod_{i<j}^{n_s} D_{k_i, m_i, k_j, m_j} \Psi_{\{k_\ell, m_\ell\}}(x, \dots, x) \langle \Psi_{\{k_\ell, m_\ell\}} | \Phi_0 \rangle 
\end{align} 
$$

Here the Kronecker delta enforces the constraint $\sum_{j=1}^{n_s} m_j = n$ with $m_j \ge 1$ and in the summation over states we used $\sum_{k_j} \to m_j L \int_{\mathbb{R}} \frac{dk}{2\pi}$ which holds also here in the large $L$ limit: the momenta sums become continuous and one can use that the string momenta $m_j k_j$ correspond to free particles as in Refs. [15,24,25,57,87].

We can simplify the factor $\Psi_{\{k_\ell, m_\ell\}}(x, \dots, x)$ in (3.13). For the general Bethe state (3.7) (before insertion of the string solution), the $x=0$ limit then reads

$$ \Psi_{\mu}(0, \dots, 0) = \frac{n!}{A^n} \prod_{j=1}^{n} \lambda_j \qquad (3.14) $$

Inserting the string solution we see that we can replace in (3.13) at the wall $x=0$

$$ \Psi_{\{k_\ell, m_\ell\}}(0, \dots, 0) = \frac{n!}{A^n} \prod_{j=1}^{n_s} A_{k_j, m_j} \qquad (3.15) $$

$$ A_{k,m} = \prod_{a=1}^{m} \left( k + i \frac{\bar{c}}{2} (m + 1 - 2a) \right) = (-i \bar{c})^m \frac{\Gamma\left(\frac{1+m}{2} + \frac{ik}{c}\right)}{\Gamma\left(\frac{1-m}{2} + \frac{ik}{c}\right)} \qquad (3.16) $$

To obtain the $n$-th moment in (3.13) we still need to calculate the overlap $\langle \Psi_{\{k_\ell, m_\ell\}} | \Phi_0 \rangle$ where $\Phi_0$ is given in (3.5). In general it involves sums over permutations and leads to complicated expressions but in our case, a simple structure emerges akin to the one known
in full-space for a few initial conditions (droplet, half-flat, Brownian). Here, as we find in the Appendix A, the result in the half-space for Brownian initial conditions is quite simple

$$ \langle \Psi_{\mu} | \Phi_0 \rangle = \frac{n!}{A^n} \frac{\Gamma(A+B+1)}{\Gamma(A+B-n+1)} \prod_{j=1}^{n} \frac{\lambda_j}{B^2 + \lambda_j^2} \quad (3.17) $$

This holds under the condition that the integral converge, that is $\frac{n}{2} < B + \frac{1}{2}$, which we will also assume from now on. Inserting the rapidities $\lambda_j$ of the string state one see that the denominator in the product in the overlap (3.17) read

$$ \begin{aligned} E_{k,j} &= \prod_{a=1}^{m} \frac{1}{B^2 + (k + i \frac{\tilde{c}}{2}(m+1-2a))^2} \\ &= \frac{1}{\tilde{c}^{2m}} \frac{\Gamma(\frac{1-m}{2} + \frac{B+ik}{\tilde{c}}) \Gamma(\frac{1-m}{2} + \frac{B-ik}{\tilde{c}})}{\Gamma(\frac{1+m}{2} + \frac{B+ik}{\tilde{c}}) \Gamma(\frac{1+m}{2} + \frac{B-ik}{\tilde{c}})} \end{aligned} \quad (3.18) $$

while the numerator was already calculated in (3.15). We can thus define $C_{k,j} = A_{k,j}^2 E_{k,j}$ and putting all together we obtain the starting expression for the integer moments, denoting here and below $Z(0, t) = Z(t)$

$$ \begin{aligned} \mathbb{E}[Z(t)^n] &= \frac{\Gamma(A+B+1)}{\Gamma(A+B-n+1)} \sum_{n_s=1}^{n} \frac{2^{n_s} \tilde{c}^n n!}{n_s! \tilde{c}^{n_s}} \\ & \quad \times \prod_{p=1}^{n_s} \sum_{m_p \ge 1} \int_{\mathbb{R}} \frac{dk_p}{2\pi} m_p C_{k_p, m_p} S_{k_p, m_p} H_{k_p, m_p} e^{(m_p^3 - m_p) \frac{\tilde{c}^2 t}{12} - m_p k_p^2 t} \\ & \quad \delta_{n_s, \sum_{j=1}^{n_s} m_j} \prod_{i<j}^{n_s} D_{k_i, m_i, k_j, m_j} \end{aligned} \quad (3.19) $$

where we recall the constraint $\sum_{j=1}^{n_s} m_j = n$. Let us use $\tilde{c} = 1$ from now on. Denoting

$$ \begin{aligned} B_{k,m} &= 4m^2 C_{k,m} S_{k,m} H_{k,m} \\ &= \frac{2k}{\pi} \sinh(2\pi k) \Gamma(2ik + m) \Gamma(-2ik + m) \\ & \quad \times \frac{\Gamma(\frac{1-m}{2} + A + ik) \Gamma(\frac{1-m}{2} + A - ik) \Gamma(\frac{1-m}{2} + B + ik) \Gamma(\frac{1-m}{2} + B - ik)}{\Gamma(\frac{1+m}{2} + A + ik) \Gamma(\frac{1+m}{2} + A - ik) \Gamma(\frac{1+m}{2} + B + ik) \Gamma(\frac{1+m}{2} + B - ik)} \end{aligned} \quad (3.20) $$

The starting formula for the moments is then

$$ \begin{aligned} \mathbb{E}[Z(t)^n] &= \frac{\Gamma(A+B+1)}{\Gamma(A+B-n+1)} \\ & \quad \sum_{n_s=1}^{n} \frac{n!2^{n_s}}{n_s!} \prod_{p=1}^{n_s} \sum_{m_p \ge 1} \int_{\mathbb{R}} \frac{dk_p}{2\pi} \frac{B_{k_p, m_p}}{4m_p} e^{(m_p^3 - m_p) \frac{t}{12} - m_p k_p^2 t} \delta_{n_s, \sum_{j=1}^{n_s} m_j} \prod_{i<j}^{n_s} D_{k_i, m_i, k_j, m_j} \end{aligned} \quad (3.21) $$

where $B_{k,m}$ is given in (3.20) and $D_{k_i, m_i, k_j, m_j}$ is given in (3.12) and where we recall the constraint $\sum_{j=1}^{n_s} m_j = n$.
## 3.3 Decorated Moments

As in Refs. [30,97], it is useful to eliminate the Gamma factor $\frac{\Gamma(A+B+1)}{\Gamma(A+B-n+1)}$ in (3.21). To this aim we introduce a random variable $W \sim \text{Gamma}^{-1}(A + B + 1)$, independent of the KPZ height, which is inverse gamma distributed with parameter $A + B + 1$, in this case

$$p_W(x) = \frac{1}{\Gamma(A + B + 1)} x^{-A-B-2} e^{-\frac{1}{x}} \Theta(x) \quad (3.22)$$

The n-th moment of $W$ is given by

$$\mathbb{E}[W^n] = \frac{\Gamma(A + B - n + 1)}{\Gamma(A + B + 1)} \quad (3.23)$$

As we will see $\mathbb{E}[W^n Z(t)^n]$ will serve as the basis to form a Fredhom Pfaffian.

## 3.4 Moments in Terms of a Pfaffian

An important identity, which makes the problem solvable in the end, is that the inverse norms of the states can be expressed as a Schur Pfaffian. Introducing the reduced variables $X_{2p-1} = m_p + 2ik_p$ and $X_{2p} = m_p - 2ik_p$ for $p \in [1, n_s]$, the norm reads

$$\prod_{1 \le i < j \le n_s} D_{k_i, m_i, k_j, m_j} = \prod_{j=1}^{n_s} \frac{m_j}{2ik_j} \text{Pf} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \quad (3.24)$$

where we recall that the Pfaffian of an anti-symmetric matrix $A$ of size $N \times N$ is defined by

$$\text{Pf}(A) = \sqrt{\det(A)} = \sum_{\substack{\sigma \in S_N \\ \sigma(2p-1) < \sigma(2p)}} \text{sign}(\sigma) \prod_{p=1}^{N/2} A_{\sigma(2p-1), \sigma(2p)} \quad (3.25)$$

and that the Schur Pfaffian is given by (see Ref. [98])

$$\text{Pf} \left[ \frac{X_i - X_j}{X_i + X_j} \right] = \prod_{i < j} \frac{X_i - X_j}{X_i + X_j} \quad (3.26)$$

Hence the starting formula for the moments now becomes:

$$\mathbb{E}[W^n Z(t)^n] = \sum_{n_s=1}^{n} \frac{n!}{n_s!} \prod_{p=1}^{n_s} \sum_{m_p \ge 1} \int_{\mathbb{R}} \frac{dk_p}{2\pi} \frac{B_{k_p, m_p}}{4ik_p} e^{(m_p^3 - m_p) \frac{t}{12} - m_p k_p^2 t} \delta_{n, \sum_{j=1}^{n_s} m_j} \text{Pf} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \quad (3.27)$$

## 4 Moments from the Log-Gamma Polymer

In this section we compute again the moments of the solution of the SHE, $Z(x, t)$, taking a limit of a known formula for the moments of the partition function of the log-gamma polymer on the half-quadrant square lattice. This method uses the convergence of the log-gamma
<figure><img src="image_6.png" /><figcaption>Fig. 4 An admissible path in the half space log-gamma polymer model, that is a path proceeding by unit steps rightward and upward in the half quadrant $\{(i, j) \in \mathbb{Z}_{>0}^2 : i \ge j\}$</figcaption></figure>

polymer to the KPZ equation at high temperature and a combinatorial conjecture of Borodin-Bufetov-Corwin [55, Conjecture 5.2]. We also discuss in Sects. 4.5 and 4.6 useful identities in distribution coming from symmetries in so-called half-space Macdonald processes [51].

## 4.1 Moment Formula for the Log-Gamma Polymer

**Definition 4.1** (Half-space log-gamma polymer) Let $\alpha_0, \alpha_1, \alpha_2, \dots$ be real parameters such that $\alpha_i + \alpha_0 > 0$ for all $i \ge 1$ and $\alpha_i + \alpha_j > 0$ for all $i \ne j \ge 1$. The half-space log-gamma polymer is a probability measure on up-right paths confined in the half-quadrant $\{(i, j) \in \mathbb{Z}_{>0}^2 : i \ge j\}$ (see Fig. 4), where the probability of an admissible path $\pi$ between (1, 1) and $(n, m)$ is given by

$$ \frac{1}{\mathcal{Z}(n, m)} \prod_{(i,j) \in \pi} w_{i,j}, $$

and where $(w_{i,j})_{i \ge j}$ is a family of independent random variables such that for $i > j$, $w_{i,j} \sim \text{Gamma}^{-1}(\alpha_i + \alpha_j)$ and $w_{i,i} \sim \text{Gamma}^{-1}(\alpha_0 + \alpha_i)$. The notation $\text{Gamma}^{-1}(\theta)$ denotes the inverse of a Gamma distributed random variable with shape parameter $\theta$. The partition function $\mathcal{Z}(n, m)$ is given by

$$ \mathcal{Z}(n, m) = \sum_{\pi:(1,1) \to (n,m)} \prod_{(i,j) \in \pi} w_{i,j}. $$

The moments of the partition function were computed using half-space Macdonald processes in [51].

**Proposition 4.2** ([51, Corollary 6.36]) For $n \ge m$ and $k \in \mathbb{Z}_{>0}$ such that $k < \min\{2\alpha_i, \alpha_i + \alpha_0\}$,

$$ \mathbb{E}[\mathcal{Z}(n, m)^k] = \oint \frac{dz_1}{2i\pi} \cdots \oint \frac{dz_k}{2i\pi} \prod_{1 \le a < b \le k} \frac{z_a - z_b}{z_a - z_b - 1} \frac{z_a + z_b}{1 + z_a + z_b} \\ \times \prod_{i=1}^k \frac{2z_i}{z_i - \alpha_0 + 1/2} \prod_{j=1}^n \left( \frac{1}{\alpha_j - z_i - 1/2} \right) \prod_{j=1}^m \left( \frac{1}{z_i + \alpha_j - 1/2} \right) \quad (4.1) $$
where the contours are such that for all $1 \le c \le k$, the contour for $z_c$ encloses $\{-\alpha_j + 1/2\}_{1 \le j \le m}$ and $\{z_{c+1} + 1, \dots, z_k + 1\}$, and excludes the poles of the integrand at $\alpha_0 - 1/2$ and $\alpha_j - 1/2$ (for $1 \le j \le n$). Because the integrand decays at least quadratically at infinity, one may choose the contours to be all vertical lines such that the contour for the variable $z_i$ is $r_i + i\mathbb{R}$ where

$$ \max_j \{k - \alpha_j - 1/2\} < r_k + k - 1 < \dots < r_2 + 2 < r_1 < \min_j \{\alpha_0 - 1/2, \alpha_j - 1/2\}. $$

Note that if $k > \alpha_i + \alpha_j$ or $k > \alpha_i + \alpha_0$ for some $i < j$, the k-th moment of $Z(n, m)$ fails to exist.

**Remark 4.3** One may also compute mixed moments of the partition functions at several points along a down-right path.

## 4.2 Stationary Structure for the Log-Gamma Polymer

In this paragraph, we will need to assume $\alpha_0 + \alpha_1 = 0$. In order to do so, we need to consider a modified partition function where we have removed the weight $w_{1,1}$, i.e. we define

$$ Z^{\text{stat}}(n, m) = \frac{\mathcal{Z}(n, m)}{w_{1,1}}. \qquad (4.2) $$

Following [61], we define horizontal and vertical increments of the partition function as

$$ U_{n,m} = \frac{Z^{\text{stat}}(n, m)}{Z^{\text{stat}}(n-1, m)}, \quad V_{n,m} = \frac{Z^{\text{stat}}(n, m)}{Z^{\text{stat}}(n, m-1)}. \qquad (4.3) $$

The partition function satisfies the recurrence

$$ Z^{\text{stat}}(n, m) = w_{n,m}(Z^{\text{stat}}(n-1, m) + Z^{\text{stat}}(n, m-1)), \qquad (4.4) $$

where by convention we have assumed that $Z^{\text{stat}}(n, m) = 0$ if $(n, m)$ does not belong to the half quadrant $\{(i, j) \in \mathbb{Z}_{>0}^2 : i \ge j\}$. From there, one may deduce a recurrence for the increments

$$ U_{n,m} = w_{n,m} \left( 1 + \frac{U_{n,m-1}}{V_{n-1,m}} \right), \quad V_{n,m} = w_{n,m} \left( 1 + \frac{V_{n-1,m}}{U_{n,m-1}} \right). \qquad (4.5) $$

We will need the following lemma from [61] where the stationary structure for the full-space log-gamma polymer was introduced.

**Lemma 4.4** ([61, Lemma 3.2]) Let $U, V, w$ be independent random variables. Let

$$ U' = w \left( 1 + \frac{U}{V} \right), \quad V' = w \left( 1 + \frac{V}{U} \right), \quad w' = \left( \frac{1}{U} + \frac{1}{V} \right)^{-1}. \qquad (4.6) $$

If for some $\alpha > 0$ and $\theta \in (-\alpha, \alpha)$, $U \sim \text{Gamma}^{-1}(\alpha + \theta)$, $V \sim \text{Gamma}^{-1}(\alpha - \theta)$, $w \sim \text{Gamma}^{-1}(2\alpha)$, then the triples $(U, V, w)$ and $(U', V', w')$ have the same distribution.

Coming back to our model $Z^{\text{stat}}(n, m)$, when $\alpha_0 + \alpha_1 = 0$ the model is stationary in the following sense.

**Proposition 4.5** Let $k \in \mathbb{Z}_{\ge 1}$. Assume that $\alpha_2 = \alpha_3 = \dots = \alpha > 0$ and $\alpha_0 + \alpha_1 = 0$. Consider a down-right path in the lattice going through the points $\{(n_i, m_i)\}_{1 \le i \le k}$, such that
<figure><img src="image_7.png" /><figcaption>**Fig. 5** The two types of elementary local transformations of down-right paths considered in the proof of Proposition 4.5. The thick black path represents an arbitrary down-right path. The portions in red represent the local modifications of the path that we consider</figcaption></figure>

$(n_{i+1}, m_{i+1}) - (n_i, m_i)$ equals either $(0, -1)$ or $(1, 0)$ (see Fig. 5). We associate to this down-right path increments $\{I_j\}_{1\le j\le k-1}$ where

$$I_j = \begin{cases} U_{n_{j+1},m_{j+1}} & \text{when } n_{j+1} > n_j, \\ V_{n_j,m_j} & \text{when } m_j > m_{j+1}. \end{cases}$$

Then the increments $\{I_j\}_{1\le j\le k-1}$ are all independent and distributed as $I_j \sim \text{Gamma}^{-1}(\alpha_1 + \alpha)$ when $I_j$ is a horizontal U increment, and $I_j \sim \text{Gamma}^{-1}(\alpha_0 + \alpha)$ when $I_j$ is a vertical V increment. In particular, for any $m$, the increments $\{U_{n,m}\}_{n\ge m+1}$ are independent and distributed as $U_{n,m} \sim \text{Gamma}^{-1}(\alpha_1 + \alpha)$.

**Proof** The distribution of increments along the first row is completely constrained by the definition of the model. Indeed, we have that $\mathcal{Z}^{\text{stat}}(n, 1) = \prod_{i=2}^n w_{i,1}$, so that the increments along the first row are given by $U_{n,1} = w_{n,1}$ and the definition of the model implies that weights $w_{n,1} \sim \text{Gamma}^{-1}(\alpha_1 + \alpha)$ are independent. Hence, for $m = 1$, the increments $\{U_{n,m}\}_{n\ge 2}$ are independent and distributed as $\text{Gamma}^{-1}(\alpha_1 + \alpha)$ as claimed.

In other terms, we have seen that the statement of the Proposition is true for the infinite path going through the points $(n, 1)$ for all $n \ge 1$. We will show that the property is preserved under two types of local transformation of paths, depicted in Fig. 5, that consist in

1. (boundary update) Lifting one unit upwards the starting point of the path along the boundary;

2. (bulk update) Transforming a down-right step into a right-down step.

It is clear that any infinite down-right path can be obtained by iteration of these local transformations, starting from the path going through the points $(n, 1)$ for all $n \ge 1$. Furthermore, if the statement of the Proposition is true for any infinite down-right path, it is true as well for any subpath of the form $\{(n_i, m_i)\}_{1\le i\le k}$ as in the statement of the Proposition. Hence we only need to show that the distribution of increments is preserved under the two local moves.

The distribution of increments on the boundary is constrained by the definition of the model. We have that $\mathcal{Z}^{\text{stat}}(n, n) = w_{n,n}\mathcal{Z}^{\text{stat}}(n, n-1)$, so that $V_{n,n} = w_{n,n}$ and we recall that $w_{n,n} \sim \text{Gamma}^{-1}(\alpha_0 + \alpha_n)$ is independent from all other weights. Hence, for any $n = m$, $V_{n,m}$ is distributed as $\text{Gamma}^{-1}(\alpha_0 + \alpha_n)$ and is independent from the increments $\{U_{n',m'}, V_{n',m'}\}$ for $m' < m$ (since those increments are independent from $w_{n,n}$). Hence, after a boundary update, the distribution of increments is preserved.

In order to show that the property is preserved under bulk update, we use Lemma 4.4. After a bulk update, the increments are updated according to (4.5), where $w_{n,m}$ is independent from the increments on the earlier path and distributed as $\text{Gamma}^{-1}(2\alpha)$ (recall that we have
assumed that $\alpha_2 = \alpha_3 = \cdots = \alpha_0$). If $\alpha_0 + \alpha_1 = 0$, we may set $\theta = \alpha_1 = -\alpha_0$ and Lemma 4.4 implies that increments along the new path will be distributed as $U_{n,m} \sim \text{Gamma}^{-1}(\alpha + \alpha_1)$, $V_{n,m} \sim \text{Gamma}^{-1}(\alpha + \alpha_0)$. This shows that the distribution of increments is preserved under bulk update. Because before the bulk update, the variables $(U, V, w)$ are independent from the rest of the increments $I_j$ by induction, and the new random variables $(U', V')$ are just measurable functions of $(U, V, w)$, the new variables are also independent of the other increments $I_j$. This concludes the proof. $\square$

One consequence of the stationary structure is that we may compute the expectation of $\log \mathcal{Z}^{\text{stat}}(n, m)$. We assume that parameters $\alpha_i$ are chosen as in Proposition 4.5. Observe that $\log \mathcal{Z}^{\text{stat}}(n, m)$ is equal to the sum of the logarithms of increments of the partition function along any path from $(1, 1)$ to $(n, m)$. These increments are not independent, so that their sum is a highly non trivial random variable, but we know the expectation of each increment. Since the vertical increments are distributed as $\text{Gamma}^{-1}(\alpha_0 + \alpha)$ and the horizontal increments are distributed as $\text{Gamma}^{-1}(\alpha_1 + \alpha)$, we have that $(\text{for } \alpha_0 + \alpha_1 = 0)$

$$ \mathbb{E}[\log \mathcal{Z}^{\text{stat}}(n, m)] = -(n-1)\psi(\alpha_1 + \alpha) - (m-1)\psi(\alpha_0 + \alpha), \quad (4.7) $$

where we have used that $\mathbb{E}[\log(\text{Gamma}^{-1}(\theta))] = -\psi(\theta)$ and $\psi$ is the digamma function.

## 4.3 Convergence to the Half-Space KPZ Equation

At high temperature (when the parameters of inverse gamma random variables go to infinity and space-time coordinates are rescaled appropriately), the partition function $\mathcal{Z}(n, m)$ converges to the multiplicative noise stochastic heat equation on $\mathbb{R}_{\ge 0}$ with Robin type boundary condition [66].

Although the convergence of discrete directed polymers to half-space KPZ equation was proved rigorously in [66] (based on the full-space analogous result in [99]), we will rederive (heuristically) this convergence in order to adapt it to our units and initial condition (which is not covered in [66]). Let us change coordinates and use more natural time and space coordinates $\tau = n + m - 2$ and $\varkappa = n - m$. The partition function $Z_d(\varkappa, \tau) := \mathcal{Z}(n, m)$ satisfies the discrete version of the stochastic heat equation

$$ Z_d(\varkappa, \tau) = w_{\varkappa, \tau}(Z_d(\varkappa - 1, \tau - 1) + Z_d(\varkappa + 1, \tau - 1)), \quad \varkappa > 0 \quad (4.8) $$

where $w_{\varkappa, \tau} \sim \text{Gamma}^{-1}(\gamma_{\varkappa, \tau})$ with parameter $\gamma_{\varkappa, \tau} = \alpha_n + \beta_m$ (independent for each $\varkappa, \tau$). The boundary condition at $\varkappa = 0$ is given by

$$ Z_d(0, \tau) = w_{0, \tau} Z_d(1, \tau - 1). \quad (4.9) $$

Let us renormalize $Z_d$ and define $Z_r(\varkappa, \tau) = C^{-\tau} Z_d(\varkappa, \tau)$. The correct factor to use is such that $C^\tau$ behaves asymptotically as the point to line partition function where weights would be replaced by their average. Hence we set $C = 2\mathbb{E}[w_{\varkappa, \tau}]$. Note that in the following, we will choose parameters so that $\mathbb{E}[w_{\varkappa, \tau}]$ does not depend on $\tau, \varkappa$ (except along the lines $\tau = \varkappa$ or $\varkappa = 0$). We may rewrite (4.8) as

$$ \nabla_{\tau} Z_{r}(\varkappa, \tau) = \frac{1+b_{\varkappa, \tau}}{2} \Delta_{\varkappa} Z_{r}(\varkappa, \tau - 1) + b_{\varkappa, \tau} Z_{r}(\varkappa, \tau - 1), \quad (4.10) $$

where $b_{\varkappa, \tau} = \frac{2w_{\varkappa, \tau}}{C} - 1$, $\nabla_{\tau}$ is the discrete time derivative and $\Delta_{\varkappa}$ is the discrete Laplacian. Let us fix $\alpha_0 \in \mathbb{R}$, $\alpha_1 \in \mathbb{R}$ and set $\alpha_i = 1/2 + \sqrt{n}/2$ for all $i \ge 2$ and use the scalings

$$ \tau = nt/2, \quad \varkappa = \sqrt{n}x/2. \quad (4.11) $$
In this case, we may choose $C = 2/\sqrt{n}$ and the family of random variables $w_{\varkappa, \tau}$ rescales to a white noise in the sense that $n b_{\varkappa, \tau} \Rightarrow \sqrt{2\xi}(x, t)$.

At $\tau = \varkappa$, we have that for large $\varkappa$,

$$Z_r(\varkappa, \varkappa) = \text{Gamma}^{-1}(\alpha_0 + \alpha_1) \times e^{B(x)-\alpha_1 x} + o\left(\frac{1}{n}\right), \quad (4.12)$$

where the inverse Gamma random variable (coming from $w_{1,1}$) and the Brownian motion $B(x)$ are independent.

It is then natural to define the continuous limit

$$Z_{\infty}(x, t) = \lim_{n \to \infty} Z_{r}(\varkappa, \tau), \quad (4.13)$$

so that $Z_{\infty}$ has the initial data $Z_{\infty}(x, 0) = \text{Gamma}^{-1}(\alpha_0 + \alpha_1) \times e^{B(x)-\alpha_1 x}$. Under the scalings that we consider, the boundary condition (4.9) becomes

$$Z_{\infty}(0, t) \approx \frac{w}{C} Z_{\infty}\left(\frac{2}{\sqrt{n}}, t\right). \quad (4.14)$$

Let us take the average on both sides of (4.14). We use that $\mathbb{E}[\frac{w}{C}] = \frac{1}{1+\frac{2\alpha_0-1}{\sqrt{n}}}$, and consider that the weight $w$ is independent from $Z_{\infty}(\frac{2}{\sqrt{n}}, t)$, as this is true in (4.9). We obtain

$$\mathbb{E}[Z_{\infty}(0, t)] = \left(1 - \frac{2}{\sqrt{n}}(\alpha_0 - 1/2)\right) \mathbb{E}\left[Z_{\infty}\left(\frac{2}{\sqrt{n}}, t\right)\right] + o(1/\sqrt{n}), \quad (4.15)$$

which, by Taylor approximation, leads to

$$\partial_x \mathbb{E}[Z_{\infty}(x, t)] \Big|_{x=0} = (\alpha_0 - 1/2) \mathbb{E}[Z_{\infty}(0, t)]. \quad (4.16)$$

Note that one may also obtain the more general boundary condition (2.7) for mixed moments by multiplying both sides of (4.14) by $Z_{\infty}(x_2, t) \dots Z_{\infty}(x_n, t)$ before taking the average. Finally, multiplying (4.10) by $n$ we obtain, when taking formally the $n \to \infty$ limit, that $Z_{\infty}(x, t)$ should satisfy the SHE (2.2). Thus, we have arrived at the following.

**Claim 4.6** (Combining [66] and [65].) Let $\mathcal{Z}(n, m)$ be the partition function of the log-gamma polymer (see Definition 4.1) where $\alpha_2 = \alpha_3 = \dots = \frac{\sqrt{n}}{2} + \frac{1}{2}$, $\alpha_0 = A + \frac{1}{2}$ and $\alpha_1 = B + \frac{1}{2}$. Let $\mathcal{Z}(x, t)$ be the solution of the multiplicative noise stochastic heat equation from Definition 2.1 with boundary parameter $A$ and initial drift $-1/2 - B$. Fix $t > 0$, $x \ge 0$. Then the family of random variables

$$\left\{ \frac{\mathcal{Z}\left(\frac{nt+\sqrt{nx}}{4}, \frac{nt-\sqrt{nx}}{4}\right)}{\left(\frac{2}{\sqrt{n}}\right)^{nt/2-2}} \right\}_{t>0, x\ge0} \quad (4.17)$$

converges in distribution to $\text{Gamma}^{-1}(A + B + 1) \times \mathcal{Z}(x, t)$ (in the space of continuous space-time trajectories), where the inverse Gamma random variable is independent from the process $\mathcal{Z}(x, t)$. Moreover, the partition function $\mathcal{Z}^{\text{stat}}(n, m)$ from Sect. 4.2 converges, under the exact same scalings, to $\mathcal{Z}(x, t)$ (not multiplied by $\text{Gamma}^{-1}(A + B + 1)$).

Note that the derivation that we presented above is only heuristic. We will not provide a complete proof of this result, though we indicate where the needed arguments can be found. The convergence of the polymer partition function for the half-space log-gamma polymer to the multiplicative noise stochastic heat equation is proved in [66] using a chaos series
representation of the polymer partition function, see in particular Section 5 therein. However, the setting of [66] restricts to delta initial data $B = +\infty$ (and Robin type boundary with arbitrary parameter A). The convergence for Brownian initial data with arbitrary parameter B was proven in [65, Theorem 2.2], though [65] works only in the case of Dirichlet boundary condition, that is in the case $A = +\infty$. Hence, one needs to combine the arguments from [66] and [65] to deduce this result.

Using the stationary structure from Sect. 4.2 together with Claim 4.6, we obtain the following. Let $Z(x, t)$ be as in Claim 4.6 and assume that $A + B + 1 = 0$. Then, for any time $t > 0$, $Z(x, t)/Z(0, t)$ is the exponential of a Brownian motion with drift $-B - 1/2$.

We may also compute the expectation of $h(x, t) = \log Z(x, t)$ in the stationary case when $A + B + 1 = 0$. Using (4.7) and plugging there the scalings of Claim 4.6, we obtain that

$$ \mathbb{E}[h(x, t)] = -\frac{t}{12} + \left(B + \frac{1}{2}\right)^2 t - \left(B + \frac{1}{2}\right)x, \quad A + B + 1 = 0. \quad (4.18) $$

In particular, when $A = B = -1/2$, we have that $\mathbb{E}[h(0, t)] = -t/12$.

## 4.4 Moments of the Half-Space KPZ Equation with Brownian Initial Data

Using the moment formula from Proposition 4.2 and the convergence result from Claim 4.6, we obtain the following moment formula for the half-space stochastic heat equation $Z(x, t)$. Note that the formula is valid for any $x \ge 0$.

**Claim 4.7** Let $Z(x, t)$ be the solution to the half-space stochastic heat equation (Definition 2.1) with Brownian initial data with drift $-1/2 - B$ and boundary parameter A. Assume that $B > 0$, and $A + B > k - 1$. Then, we have

$$ \begin{aligned} \mathbb{E}[Z(x, t)^k] = & 2^k \frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - k)} \int_{r_1 + i\mathbb{R}} \frac{dz_1}{2i\pi} \cdots \int_{r_k + i\mathbb{R}} \frac{dz_k}{2i\pi} \prod_{1 \le a < b \le k} \frac{z_a - z_b}{z_a - z_b - 1} \frac{z_a + z_b}{z_a + z_b - 1} \\ & \times \prod_{i=1}^k \frac{z_i}{z_i + A} \frac{1}{B^2 - z_i^2} e^{tz_i^2 - xz_i}, \end{aligned} \quad (4.19) $$

where the contours are chosen so that $B > r_1 > r_2 + 1 > \dots, > r_k + k - 1 > k - 1 - A$.

**Remark 4.8** One may also compute mixed moments of $Z(x, t)$, that is $\mathbb{E}[Z(x_1, t) \dots Z(x_k, t)]$.

**Proof** We start from the moment formula given in Proposition 4.2. Under the scalings considered in Claim 4.6, the second line of (4.1) becomes

$$ \prod_{i=1}^{k} \frac{2z_i}{z_i - \alpha_o + 1/2} \frac{1}{(\alpha_1 - 1/2)^2 - z_i^2} \left( \frac{1}{\sqrt{n}/2 - z_i} \right)^{\frac{nt+\sqrt{nx}}{4}-1} \left( \frac{1}{\sqrt{n}/2 + z_i} \right)^{\frac{nt-\sqrt{nx}}{4}-1}. $$

Using dominated convergence, one readily obtains that

$$ \lim_{n \to \infty} \mathbb{E} \left( \frac{Z \left( \frac{nt + \sqrt{nx}}{4}, \frac{nt - \sqrt{nx}}{4} \right)}{\left( \frac{2}{\sqrt{n}} \right)^{nt/2 - 2}} \right)^k = \int_{r_1 + i\mathbb{R}} \frac{dz_1}{2i\pi} \cdots \int_{r_k + i\mathbb{R}} \frac{dz_k}{2i\pi} $$
$$ \prod_{1 \le a < b \le k} \frac{z_a - z_b}{z_a - z_b - 1} \frac{z_a + z_b}{1 + z_a + z_b} \prod_{i=1}^{k} \frac{2z_i}{z_i - A} \frac{1}{B^2 - z_i^2} e^{t z_i^2 + x z_i}. \quad (4.20) $$

Using the convergence in distribution from Claim 4.6 the left hand side in (4.20) converges to $m_k \mathbb{E}[Z(x, t)^k]$, where $m_k$ is the k-th moment of an inverse Gamma random variable with parameter $\alpha_o$ and $\alpha_1$ (the convergence in distribution implies the convergence of moments modulo some tail bounds, which, using Markov inequality, can be proven from our explicit moment formulae. Details of this argument are provided in a similar case in [100, Sect. 3]). It is well known that $m_k = \Gamma(\alpha_o + \alpha_1 - k)/\Gamma(\alpha_o + \alpha_1)$, so that

$$ \mathbb{E}[Z(x, t)^k] = \frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - k)} \times \text{R.H.S. of (4.20)} \quad (4.21) $$

Finally, we have used the change of variables $\tilde{z}_i = -z_{k-i+1}$ to obtain the statement of the Proposition. □

## 4.5 Symmetry Between Drift and Boundary Parameters

We now exploit one of the symmetries of the log-gamma polymer, arising from more general symmetries of so-called half-space Macdonald processes [51], which, in the KPZ limit, extends a result of Parekh [65].

**Claim 4.9** Let us denote here by $Z_A^B(x, t)$ (with explicit dependence in the parameters A, B) the solution to the multiplicative noise SHE with boundary parameter A and initial drift $-1/2 - B$ (Definition 2.1). For any fixed $t > 0$ and A, B $\in \mathbb{R}$, we have the equality in distribution

$$ Z_A^B(0, t) = Z_B^A(0, t). \quad (4.22) $$

*Proof* This is a consequence of [51, Proposition 8.1] which states that the law of the partition function $\mathcal{Z}(n, n)$ of the half-space log-gamma polymer is invariant under exchanging parameters $\alpha_o$ and $\alpha_1$ (recall Definition 4.1), along with the convergence result from Claim 4.6.

However, [51, Proposition 8.1] assumes that $\alpha_o + \alpha_1 > 0$ as in Definition 4.1, which would require the condition $A + B + 1 > 0$. Let us explain why we do not need to assume this condition. Recall Definition 4.1 and let us denote the partition function by $Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1}(n, n)$, where we indicate explicitly the dependence on parameters $\tilde{\alpha}_o, \tilde{\alpha}_1$. The result from [51, Proposition 8.1] implies that for $\tilde{\alpha}_o, \tilde{\alpha}_1$ such that $\tilde{\alpha}_o + \tilde{\alpha}_1 > 0$ we have the equality in distribution

$$ Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1}(n, n) = Z_{\tilde{\alpha}_1}^{\tilde{\alpha}_o}(n, n). $$

Notice that by Definition 4.1, we have

$$ Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1}(n, n) = w_{1,1} Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1, \text{stat}}, \quad Z_{\tilde{\alpha}_1}^{\tilde{\alpha}_o}(n, n) = w_{1,1} Z_{\tilde{\alpha}_1}^{\tilde{\alpha}_o, \text{stat}}, $$

where $Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1, \text{stat}}$ is defined as in Sect. 4.2 and $w_{1,1}$ has the same distribution in both cases. This implies that we have also the equality in distribution

$$ Z_{\tilde{\alpha}_o}^{\tilde{\alpha}_1, \text{stat}}(n, n) = Z_{\tilde{\alpha}_1}^{\tilde{\alpha}_o, \text{stat}}(n, n). \quad (4.23) $$

The distribution of the random variable in (4.23) depends on parameters $\tilde{\alpha}_o, \tilde{\alpha}_1, \alpha_2, \dots, \alpha_n$. The equality in distribution can be analytically extended to all parameters such that $\alpha_i + \alpha_j > 0$ for any $2 \le i < j \le n$, $\tilde{\alpha}_1 + \alpha_i > 0$ for any $2 \le i \le n$, and $\tilde{\alpha}_o + \alpha_i > 0$ for any $2 \le i \le n$. In particular, we do not require anymore that $\tilde{\alpha}_o + \tilde{\alpha}_1 > 0$. Passing to the limit in (4.23) using Claim 4.6, we obtain the desired result. □
## 4.6 Another Conjectural Identity in Law

The symmetry between parameters A and B stated in Claim 4.9 relies on a similar property for the log-gamma polymer (symmetry between $\alpha_0$ and $\alpha_1$) based on the theory of half-space Macdonald processes as stated in [51, Proposition 8.1]. This result was stated in [51] for the log-gamma polymer model where $\alpha_2 = \dots = \alpha_n$. However, the same property actually holds for general half-space Macdonald process and for any choice of parameters, see [51, Proposition 2.6], as long as we restrict to the partition function on the boundary. Furthermore, the law of $Z(n, n)$ is symmetric with respect to permutation of the parameters $\alpha_i$. Thus, we claim that one can also exchange the roles of the parameter $\alpha_0$ and the parameter $\alpha_2$.

When scaling parameters to the (multiplicative noise) stochastic heat equation in Claim 4.6, we set $\alpha_0 = A - \frac{1}{2}$, $\alpha_1 = B - \frac{1}{2}$ and $\alpha_i = \frac{\sqrt{n}}{2}$. If we exchange parameters $\alpha_0$ and $\alpha_2$, we expect that we will obtain the stochastic heat equation with boundary parameter equal to $+\infty$ (i.e. with Dirichlet boundary condition $Z(0, t) = 0$) and initial condition given by

$$Z(x, 0) = \int_{0}^{x} \exp (B_{1}(y) + B_{2}(x) - B_{2}(y)) dy, \quad (4.24)$$

where $B_1$ and $B_2$ are independent Brownian motions with respective drifts $-(B + 1/2)$ and $-(A + 1/2)$. Let us call $Z_{\text{Dir}}^{A,B}$ the solution to the heat equation with Dirichlet boundary condition and the initial condition given above in (4.24). We refer to [64] regarding the exact meaning of the Dirichlet boundary condition in this context. Following similar arguments as in the proof of [64, Theorem 1.1], we conjecture that for any $t > 0$, we have the identity in distribution

$$Z_{A}^{B}(0, t) = \lim_{x \to 0} \frac{Z_{\text{Dir}}^{A, B}(x, t)}{x}. \quad (4.25)$$

The identity in law (4.25) allows to predict the law of large numbers for $h(0, t) = \log Z_A^B(0, t)$ in the bound phase. Recall that we expect that $h(0, t)$ follows the asymptotics

$$h(0, t) \approx v_{\infty}^{A,B} t + t^{\beta} \chi \quad (4.26)$$

where $\chi$ is an $O(1)$ random variable, and $\beta$ the growth fluctuation exponent. Using (4.25), $\log \frac{Z_{\text{Dir}}^{A,B}(x,t)}{x}$ should also follow the same asymptotics. When $A < -1/2$ or $B < -1/2$, we see that $Z_{\text{Dir}}^{A,B}(x, 0)$ grows as $e^{\max\{|A+\frac{1}{2}|, |B+\frac{1}{2}|\}x}$. Hence, the polymer partition function will be dominated by paths from $(x, 0)$ to $(0, t)$ where $x$ is of order $t$. More precisely, since the point to point free energy from $(x, 0)$ to $(0, t)$ behaves asymptotically as $-\frac{t}{12} - \frac{x^2}{4t}$, the partition function will be dominated by paths leaving from $x = x_t^*$, where $x_t^* = \arg\max_{x>0}(-\frac{x^2}{4t} + \max\{A + \frac{1}{2}, B + \frac{1}{2}\}x) = 2 \max\{|A + \frac{1}{2}|, |B + \frac{1}{2}|\}t$. Thus, we have that for general A, B,

$$v_{\infty}^{A,B} = -\frac{1}{12} + \left( \min \{A + \frac{1}{2}, B + \frac{1}{2}, 0\} \right)^2.$$

We can even predict the fluctuation exponent $\beta$ and the nature of fluctuations. When $A < -1/2$ or $B < -1/2$ with $A \neq B$ (say $A < B$ for simplicity), the initial condition for $Z_{\text{Dir}}^{A,B}(x, t)$ in (4.24) will essentially be the Brownian motion, i.e. $B_2(x)$ with drift $-(A+1/2)$. The fluctuations of the initial condition at the optimal point $x_t^*$ are thus $\approx B_2(x_t^*)$ i.e. Gaussian on the scale $t^{1/2}$, and they will dominate the fluctuations in the partition function, hence we find that $\beta = 1/2$ and $\chi$ is Gaussian. The situation is completely similar when $B < A$.
When $A = B < -1/2$, the situation is a bit more delicate since we cannot approximate the initial condition (4.24) by a Brownian motion. Instead, we notice that (4.24) can be interpreted as the partition function in the O'Connell-Yor directed polymer model [101]. For large values of $x$, it behaves as

$$ \log Z_{\text{Dir}}^{A,B}(x, 0) = \log \int_{0}^{x} \exp (\mathcal{B}_{1}(y) + \mathcal{B}_{2}(x) - \mathcal{B}_{2}(y)) \, dy \quad (4.27) $$

$$ \approx \left| B + \frac{1}{2} \right| x + \sqrt{x} \max_{0 \le y \le 1} \{ \mathcal{B}_1(y) + \mathcal{B}_2(1) - \mathcal{B}_2(y) \} . \quad (4.28) $$

It was proved [102,103] that the latter quantity behaves asymptotically as the largest eigenvalue of a $2 \times 2$ GUE matrix in the scale $x^{1/2}$. We must now evaluate this quantity at the optimal point $x = x_t^* = |2B + 1|t$, hence we find that $\beta = 1/2$ and that $\chi$ has the same distribution as the top eigenvalue of a $2 \times 2$ GUE matrix, discussed for instance in [101, Sect. 6].

## 4.7 Residue Expansion

In this section, we restrict to $x = 0$ and denote $Z(0, t) = Z(t)$ as in the previous sections. The moment formula (4.19) is not convenient for asymptotic analysis because the contours are different for each variable, so that the complexity of the formula significantly increases with $k$. To overcome this issue, one has to deform the contours to all lie on a fixed vertical line and take into account the residues encountered during this contour deformation. This procedure was implemented in [55], but the computation of residues is very involved and the result relies on a conjectural combinatorial simplification. Applying this result [55, Conjecture 5.2] to the moment formula (4.19), we conjecture that for $A > k - 1$ and $B > k - 1$, we have

$$ \begin{aligned} \mathbb{E}[Z(t)^k] = {}& 2^k \frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - k)} \sum_{\substack{\lambda \vdash k \\ \lambda = 1^{m_1} 2^{m_2} \dots}} \frac{(-1)^{\ell(\lambda)}}{m_1! m_2! \dots} \int_{i\mathbb{R}} \frac{dw_1}{2i\pi} \dots \int_{i\mathbb{R}} \frac{dw_{\ell(\lambda)}}{2i\pi} \\ & \times \prod_{j=1}^{\ell(\lambda)} \frac{(w_j + 1/2)_{\lambda_j - 1}}{4(w_j)_{\lambda_j}} \mathrm{Pf} \left[ \frac{u_i - u_j}{u_i + u_j} \right]_{i, j=1}^{2\ell(\lambda)} \\ & \times E(w_1, w_1 + 1, \dots, w_1 + \lambda_1 - 1, \\ & \qquad w_2, \dots, w_2 + \lambda_2 - 1, \dots, w_{\ell(\lambda)}, \dots, w_{\ell(\lambda)} + \lambda_{\ell(\lambda)} - 1), \end{aligned} \quad (4.29) $$

where we use the Pochhammer notation for rising factorials $(w)_\lambda = w(w+1)\dots(w+\lambda-1)$, we define variables $u_i$ for $1 \le i \le 2\ell(\lambda)$ as

$$ \begin{aligned} & (u_1, \dots, u_{2\ell(\lambda)}) \\ &= (-w_1 + \frac{1}{2}, w_1 - \frac{1}{2} + \lambda_1, -w_2 + \frac{1}{2}, w_2 - \frac{1}{2} + \lambda_2, \dots, -w_{\ell(\lambda)} + \frac{1}{2}, w_{\ell(\lambda)} - \frac{1}{2} + \lambda_{\ell(\lambda)}), \end{aligned} \quad (4.30) $$

and where

$$ E(z_1, \dots, z_k) = \prod_{i=1}^{k} \frac{e^{zt_i^2}}{B^2 - z_i^2} \sum_{\sigma \in B C_k} \sigma \left( \prod_{1 \le j < i \le k} \frac{z_i - z_j - 1}{z_i - z_j} \frac{z_i + z_j - 1}{z_i + z_j} \prod_{i=1}^{k} \frac{z_i}{z_i + A} \right). \quad (4.31) $$
It turns out that the symmetrization can be performed using [55, Equation (54)], relying on the theory of BC-symmetric Hall-Littlewood polynomials [104], and one finds

$$E(z_1, \dots, z_k) = 2^k k! \prod_{i=1}^{k} \frac{e^{t z_i^2}}{B^2 - z_i^2 z_i^2 - A^2} \quad (4.32)$$

**Remark 4.10** It is now apparent that the moment formulae are invariant with respect to the transformation $(A, B) \mapsto (B, A)$.

Performing explicitly the evaluation into strings, we obtain

$$\mathbb{E}[Z(t)^k] = 4^k k! \frac{\Gamma(A+B+1)}{\Gamma(A+B+1-k)} \sum_{\substack{\lambda \vdash k \\ \lambda=1^{m_1} 2^{m_2} \dots}} \frac{(-1)^{\ell(\lambda)}}{m_1! m_2! \dots} \int_{i \mathbb{R}} \frac{dw_1}{2i\pi} \dots \int_{i \mathbb{R}} \frac{dw_{\ell(\lambda)}}{2i\pi} \\ \times \mathrm{Pf} \left[ \frac{u_i - u_j}{u_i + u_j} \right]_{i,j=1}^{2\ell(\lambda)} \prod_{j=1}^{\ell(\lambda)} \frac{e^{\mathrm{Tr}(w_i + \lambda_i)}}{e^{\mathrm{Tr}(w_i)}} \frac{(w_j + 1/2)_{\lambda_j - 1}}{4(w_j)_{\lambda_j}} \left( \frac{\Gamma(w_i + \lambda_i)}{\Gamma(w_i)} \right)^2 \\ \times \frac{\Gamma(B - w_j - \lambda_j + 1) \Gamma(B + w_j)}{\Gamma(B - w_j + 1) \Gamma(B + w_j + \lambda_j)} \frac{\Gamma(w_j - A) \Gamma(w_j + A)}{\Gamma(w_j + \lambda_j - A) \Gamma(w_j + \lambda_j + A)}, \quad (4.33)$$

where

$$\mathbb{G}(w) = \frac{w^3}{3} - \frac{w^2}{2} + \frac{w}{6},$$

so that

$$\mathbb{G}(w + \ell) - \mathbb{G}(w) = w^2 + (w + 1)^2 + \dots + (w + \ell - 1)^2. \quad (4.34)$$

Note that $B^2 - z^2 = (\pm B - z)(\pm B + z)$ and $z^2 - A^2 = -(\pm A - z)(\pm A + z)$, so that many choices are possible for the evaluation of $E(z_1, \dots, z_k)$ into strings. The most convenient choice seems to be the following formula.

**Claim 4.11** (based on [55, Conjecture 5.2]) For $A > k - 1$ and $B > k - 1$, we have

$$\mathbb{E}[Z(t)^k] = 4^k k! \frac{\Gamma(A+B+1)}{\Gamma(A+B+1-k)} \sum_{\substack{\lambda \vdash k \\ \lambda=1^{m_1} 2^{m_2} \dots}} \frac{(-1)^{\ell(\lambda)}}{m_1! m_2! \dots} \int_{i \mathbb{R}} \frac{dw_1}{2i\pi} \dots \int_{i \mathbb{R}} \frac{dw_{\ell(\lambda)}}{2i\pi} \\ \times \mathrm{Pf} \left[ \frac{u_i - u_j}{u_i + u_j} \right]_{i,j=1}^{2\ell(\lambda)} \prod_{j=1}^{\ell(\lambda)} \frac{e^{\mathrm{Tr}(w_i + \lambda_i)}}{e^{\mathrm{Tr}(w_i)}} \frac{(w_j + 1/2)_{\lambda_j - 1}}{4(w_j)_{\lambda_j}} \frac{\Gamma(-w_j + 1) \Gamma(w_j + \lambda_j)}{\Gamma(-w_j - \lambda_j + 1) \Gamma(w_j)} \\ \times \frac{\Gamma(B - w_j - \lambda_j + 1) \Gamma(B + w_j)}{\Gamma(B - w_j + 1) \Gamma(B + w_j + \lambda_j)} \frac{\Gamma(A - w_j - \lambda_j + 1) \Gamma(A + w_j)}{\Gamma(A - w_j + 1) \Gamma(A + w_j + \lambda_j)}. \quad (4.35)$$

Comparing with the formula (3.27) obtained from the replica Bethe ansatz, we see that (4.35) and (3.27) agree after the substitutions

$$w_j \rightarrow i k_j + \frac{1 - m_j}{2}, \quad \ell(\lambda) \rightarrow n_s, \quad \lambda_j \rightarrow m_j. \quad (4.36)$$

Indeed, under this change of variables, we have

$$(u_1, u_2, \dots, u_{2\ell(\lambda)-1}, u_{2\ell(\lambda)}) = (X_2, X_1, \dots, X_{2n_s}, X_{2n_s-1}) \quad (4.37)$$
where $X_j = m_p + 2ik_p$ for $j = 2p-1$, $X_j = m_p - 2ik_p$ for $j = 2p$. Thus we have that

$$(-1)^{\ell(\lambda)} \mathrm{Pf} \left[ \frac{u_i - u_j}{u_i + u_j} \right]_{i,j=1}^{2\ell(\lambda)} = \mathrm{Pf} \left[ \frac{X_i - X_j}{X_i + X_j} \right]. \quad (4.38)$$

One may also check that $e^{\mathcal{G}(w_i+\lambda_i)-\mathcal{G}(w_j)} = e^{(m_j^3-m_j)\frac{\ell}{12}-m_j k_j^2\ell}$. Using the reflection formula for the Gamma function, we may write the hyperbolic sine in the definition of $B_{k,m}$ in (3.20) as a product of Gamma functions so that we have (dropping indices of variables $k_j, m_j, \lambda_j, w_j$)

$$\frac{B_{k,m}}{4ik} = \frac{1}{2} \frac{\Gamma(2w - 2\lambda - 1)\Gamma(1 - 2w)\Gamma(B - w - \lambda + 1)\Gamma(B + w)\Gamma(A - w - \lambda + 1)\Gamma(A + w)}{\Gamma(2w + \lambda)\Gamma(1 - 2w - \lambda)\Gamma(B - w + 1)\Gamma(B + w + \lambda)\Gamma(A - w + 1)\Gamma(A + w + \lambda)}. \quad (4.39)$$

Using the duplication formula for the Gamma function and the reflection formula a few times, we arrive at

$$\frac{B_{k,m}}{4ik} = 2^{2\lambda} \frac{(w + 1/2)_{\lambda-1}}{4(w)_{\lambda}} \frac{\Gamma(-w + 1)\Gamma(w + \lambda)}{\Gamma(-w - \lambda + 1)\Gamma(w)} \times \frac{\Gamma(B - w - \lambda + 1)\Gamma(B + w)\Gamma(A - w - \lambda + 1)\Gamma(A + w)}{\Gamma(B - w + 1)\Gamma(B + w + \lambda)\Gamma(A - w + 1)\Gamma(A + w + \lambda)}. \quad (4.40)$$

Thus, we have shown that (4.35) and (3.27) agree as claimed.

## 5 Generating Function in Terms of a Fredholm Pfaffian

We will now write the moment generating function of $Z(t)$. We define, for $\zeta > 0$,

$$g(\zeta) = \mathbb{E} \left[ \exp(-\zeta e^{\frac{\ell}{12}} WZ(t)) \right]. \quad (5.1)$$

Ignoring the fact that the summation over $n$ cannot be exchanged with the expectation due to the divergence of moments, we will consider the following formal power series

$$1 + \sum_{n=1}^{\infty} \frac{(-\zeta e^{\frac{\ell}{12}})^n}{n!} \mathbb{E}[W^n Z(t)^n],$$

that we will again denote by $g(\zeta)$.

**Remark 5.1** Following [97], we may informally write

$$\mathbb{E}[W^n]\zeta^n = \frac{\Gamma(A + B - n + 1)}{\Gamma(A + B + 1)}\zeta^n = \frac{\Gamma(A + B - \zeta\partial_\zeta + 1)}{\Gamma(A + B + 1)}\zeta^n \quad (5.2)$$

so that

$$\mathbb{E} \left[ \exp(-\zeta e^{\frac{\ell}{12}} Z(t)) \right] = \frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - \zeta \partial_{\zeta})} \mathbb{E} \left[ \exp(-\zeta e^{\frac{\ell}{12}} WZ(t)) \right], \quad (5.3)$$

where the operator $\frac{\Gamma(A+B+1)}{\Gamma(A+B+1-\zeta\partial_{\zeta})}$ should be understood in the following sense: If

$$\frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - z)} = \sum_{n=0}^{+\infty} a_n z^n,$$
in a neighborhood of $z = 0$, then

$$ \frac{\Gamma(A + B + 1)}{\Gamma(A + B + 1 - \zeta \partial_{\zeta})} = \sum_{n=0}^{+\infty} a_n (\zeta \partial_{\zeta})^n. $$

Equation (5.3) is a positive temperature analogue of [105, Eq. (4.5)] which we will use in the next section. This type of arguments can be traced back to the work of Baik and Rains [58].

**Remark 5.2** Alternatively, one may use the density of the inverse Gamma distribution and write

$$ g(\zeta) = \int_{0}^{+\infty} \frac{du}{\Gamma(A+B+1)} u^{-A-B-2} e^{-1/u} \mathbb{E} \left[ \exp(-\zeta u e^{t/12} Z(t)) \right]. \quad (5.4) $$

We may rewrite this equation using the following representation of the Bessel function

$$ \int_{0}^{+\infty} du \, u^{-\nu} e^{1/u - x^2 u} = 2x^{\nu} K_{\nu}(2x), \quad (5.5) $$

for $x > 0$, where $K_{\nu}$ denotes the modified Bessel $K$ function. Exchanging the expectation with the integral over $u$ in (5.4), and using this integral representation, we arrive at

$$ \mathbb{E} \left[ 2 \left( \zeta Z(t) e^{\frac{t}{12}} \right)^{\frac{1+A+B}{2}} K_{1+A+B} \left( 2 \sqrt{\zeta Z(t)} e^{\frac{t}{12}} \right) \right] = g(\zeta) \Gamma(1+A+B), \quad (5.6) $$

where $K_{\nu}(z)$ denotes the modified Bessel $K$ function. Note that similar integral transforms involving the modified Bessel function $K$ appear in [31, Theorem 2.9]. It is plausible that this expression can be inverted to compute the distribution of $Z(t)$. For $A + B + 1 = 0$, an inversion formula is provided in [31, Appendix E]. We will see in Sect. 6 that we will actually not need to perform this inversion.

The constraint $\sum_{i=1}^{n_s} m_i = n$ in (3.27) can then be relaxed by reorganizing the series according to the number of strings:

$$ g(\zeta) = 1 + \sum_{n_s=1}^{\infty} \frac{1}{n_s!} Z(n_s, \zeta) \quad (5.7) $$

where $Z(n_s, \zeta)$ is the partition sum at fixed number of strings $n_s$, calculated below. We now show that one can write the generating function as a Fredholm Pfaffian. It will be possible thanks to the Schur Pfaffian identity, (3.24), given above. The partition sum at fixed number of strings, expressed in terms of the reduced variables $X_{2p-1} = m_p + 2ik_p$ and $X_{2p} = m_p - 2ik_p$ for $p \in [1, n_s]$, reads

$$ Z(n_s, \zeta) = \prod_{p=1}^{n_s} \sum_{m_p \ge 1} \int_{\mathbb{R}} \frac{dk_p}{2\pi} (-\zeta)^{m_p} \frac{B_{k_p, m_p}}{4ik_p} e^{-tm_p k_p^2 + \frac{t}{12} m_p^3} \mathrm{Pf}_{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \quad (5.8) $$

where $B_{k,m}$ was given in (3.20). The summation over the variables $m_p$ can be done using the Mellin-Barnes summation trick similarly to Refs. [30,97]. The barrier $A > (n-1)/2$ is overcome exactly as in Ref. [97] (see Lemma. 6 and the discussion therein) from an analytic continuation of Gamma functions included in the $B_{k,m}$ factor, the introduction of a particular contour $C_0$ and a final requirement for the drift $A + 1/2 > 0$. Indeed, define the contour $C_0 = a + i\mathbb{R}$ with $a \in (0, \min\{2B+1, 2A+1, 1\})$, then for any holomorphic function $f$
having sufficient decay at infinity and in particular denoting the summand of Eq. (5.8) by the function $f(m_p)$, we have

$$ \sum_{m \ge 1} (-\zeta)^m f(m) = - \int_{C_0} \frac{dw}{2i\pi} \zeta^w \frac{\pi}{\sin \pi w} f(w) = - \int_{\mathbb{R}} dr \frac{\zeta}{\zeta + e^{-r}} \int_{C_0} \frac{dw}{2i\pi} f(w) e^{-wr}. \quad (5.9) $$

For each $m_p$ we therefore introduce two variables $r_p$ and $w_p$ and we redefine the reduced variables $X_{2p}$ and $X_{2p-1}$ under the minimal replacement $m_p \to w_p$ imposed by the Mellin-Barnes formula, which we will apply despite the presence of poles on the right of the contour $C_0$. This is an a priori illegal step, but it will exactly turn the diverging moment generating series into a well-defined and converging series equal to the Laplace transform. We refer to [55, Sect. 6] where this procedure and its degree of rigor is discussed in great details. This leads to the following rewriting of the coefficient $Z(n_s, \zeta)$ as (see section 5 in [54] for similar manipulations)

$$ \begin{aligned} Z(n_s, \zeta) = &(-1)^{n_s} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-r_p}} \iint_{C_0^2} \frac{dX_{2p-1}}{4i\pi} \frac{dX_{2p}}{4i\pi} \frac{\sin(\frac{\pi}{2}(X_{2p} - X_{2p-1}))}{2\pi} \\ & \times \frac{\Gamma(A + \frac{1}{2} - \frac{X_{2p}}{2}) \Gamma(A + \frac{1}{2} - \frac{X_{2p-1}}{2}) \Gamma(B + \frac{1}{2} - \frac{X_{2p}}{2}) \Gamma(B + \frac{1}{2} - \frac{X_{2p-1}}{2})}{\Gamma(A + \frac{1}{2} + \frac{X_{2p}}{2}) \Gamma(A + \frac{1}{2} + \frac{X_{2p-1}}{2}) \Gamma(B + \frac{1}{2} + \frac{X_{2p}}{2}) \Gamma(B + \frac{1}{2} + \frac{X_{2p-1}}{2})} \\ & \times \Gamma(X_{2p-1}) \Gamma(X_{2p}) e^{-(X_{2p-1} + X_{2p})\frac{r_p}{2}} + t(\frac{x_{2p-1}^3}{24} + \frac{x_{2p}^3}{24}) \frac{\text{Pf}}{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \end{aligned} \quad (5.10) $$

**Remark 5.3** Note that the contour $C_0$ passes to the left of the poles of the Gamma function at $X = 2A + 1, 2B + 1$.

We observe that the integrals are almost separable in $X_{2p}$ and $X_{2p-1}$ except for the sine function which couples them. By anti-symmetrization and similarly to [54, Sect. 5], we can proceed to the replacement<fnref n="2" />

$$ \sin\left(\frac{\pi}{2}(X_{2p} - X_{2p-1})\right) \to 2 \sin\left(\frac{\pi}{2}X_{2p}\right) \cos\left(\frac{\pi}{2}X_{2p-1}\right). \quad (5.11) $$

The last manipulations consist in rescaling all variables $X$ by a factor 2 and replacing the contours of integration by $C = \frac{a}{2} + i\mathbb{R}$. Hence we have

$$ \begin{aligned} Z(n_s, \zeta) = &(-1)^{n_s} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-r_p}} \iint_{C_2} \frac{dX_{2p-1}}{2i\pi} \frac{dX_{2p}}{2i\pi} \frac{\sin(\pi X_{2p}) \cos(\pi X_{2p-1})}{\pi} \\ & \times \frac{\Gamma(A + \frac{1}{2} - X_{2p}) \Gamma(A + \frac{1}{2} - X_{2p-1}) \Gamma(B + \frac{1}{2} - X_{2p}) \Gamma(B + \frac{1}{2} - X_{2p-1})}{\Gamma(A + \frac{1}{2} + X_{2p}) \Gamma(A + \frac{1}{2} + X_{2p-1}) \Gamma(B + \frac{1}{2} + X_{2p}) \Gamma(B + \frac{1}{2} + X_{2p-1})} \\ & \times \Gamma(2X_{2p-1}) \Gamma(2X_{2p}) e^{-(X_{2p-1} + X_{2p})r_p + t(\frac{x_{2p-1}^3}{3} + \frac{x_{2p}^3}{3})} \frac{\text{Pf}}{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \end{aligned} \quad (5.12) $$

<footnote n="2">² This replacement is done using the trigonometric identity $\sin(x - y) = \sin(x) \cos(y) - \sin(y) \cos(x)$ and the antisymmetry of the Schur Pfaffian upon the relabelling $X_{2p} \leftrightarrow X_{2p-1}$. Note that similarly to [54], we could have the more general decomposition $\sin(\frac{\pi}{2}(X_{2p} - X_{2p-1})) \to 2 \sin(\frac{\pi}{2} X_{2p} + \theta) \cos(\frac{\pi}{2} X_{2p-1} + \theta)$. valid for all $\theta \in \mathbb{R}$ but we do not make use of it here.</footnote>
There are a few last steps before we introduce the Fredholm Pfaffian. First define the functions

$$ \phi_{2p}(X) = \frac{\sin(\pi X)}{\pi} \Gamma(2X) \frac{\Gamma(A + \frac{1}{2} - X) \Gamma(B + \frac{1}{2} - X)}{\Gamma(A + \frac{1}{2} + X) \Gamma(B + \frac{1}{2} + X)} e^{-r_p X + it \frac{X^3}{3}} \quad (5.13) $$

$$ \phi_{2p-1}(X) = \cos(\pi X) \Gamma(2X) \frac{\Gamma(A + \frac{1}{2} - X) \Gamma(B + \frac{1}{2} - X)}{\Gamma(A + \frac{1}{2} + X) \Gamma(B + \frac{1}{2} + X)} e^{-r_p X + it \frac{X^3}{3}} $$

Using a known property of Pfaffians (see De Bruijn [106]), we can rewrite the partition sum at fixed number of strings itself as a Pfaffian, i.e. we use that

$$ \prod_{\ell=1}^{2n_s} \int_C \frac{dX_\ell}{2i\pi} \Phi_\ell(X_\ell) \mathrm{Pf}_{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] = \mathrm{Pf}_{2n_s \times 2n_s} \left[ \iint_C^2 \frac{dw}{2i\pi} \frac{dz}{2i\pi} \Phi_i(w) \Phi_j(z) \frac{w-z}{w+z} \right] \quad (5.14) $$

This leads to the definition of a $2n_s \times 2n_s$ matrix $M$ such that

$$ M_{ij} = \iint_C^2 \frac{dw}{2i\pi} \frac{dz}{2i\pi} \Phi_i(w) \Phi_j(z) \frac{w-z}{w+z} \quad (5.15) $$

Since a variable $r_p$ will be shared between four elements of this matrix, it is more convenient to view $M$ as composed of $2 \times 2$ blocks which we denote $K$, whose elements are presented in Eqs. (2.12). Finally, the string-replicated partition function is given by an infinite series of Pfaffians

$$ g(\zeta) = 1 + \sum_{n_s=1}^{\infty} \frac{(-1)^{n_s}}{n_s!} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-r_p}} \mathrm{Pf}_{n_s \times n_s} (K(r_k, r_\ell)) \quad (5.16) $$

This series is a Fredholm Pfaffian,

$$ g(\zeta) = \mathbb{E}[\exp(-\zeta W e^{H(t)})] = \mathrm{Pf}(J - \sigma_{\zeta} K)_{\mathbb{L}^2(\mathbb{R})}, \quad (5.17) $$

where $K$ is given in (2.12), the function $\sigma_\zeta$ is given by $\sigma_\zeta(r) = \frac{\zeta}{\zeta+e^{-r}}$ and the $2 \times 2$ kernel $J$ is given by $J(r, r') = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} 1_{r=r'}$. For the precise definition and properties of Fredholm Pfaffians see Section 8 in [107], as well as e.g. Section 2.2. in [44], Appendix B in [108] and Appendix G in [24,25].

# 6 Large Time Limit of the Fredholm Pfaffian and the Distribution of the KPZ Height: Crossover Kernel

We will now study the large time limit of our kernel. To understand the scaling required at large time, let us recall the expression of the partition sum at fixed number of strings

$$ Z(n_s, \zeta) = (-1)^{n_s} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-r_p}} \iint_C^2 \frac{dX_{2p-1}}{2i\pi} \frac{dX_{2p}}{2i\pi} \frac{\sin(\pi X_{2p}) \cos(\pi X_{2p-1})}{\pi} \\ \times \frac{\Gamma(A + \frac{1}{2} - X_{2p}) \Gamma(A + \frac{1}{2} - X_{2p-1}) \Gamma(B + \frac{1}{2} - X_{2p}) \Gamma(B + \frac{1}{2} - X_{2p-1})}{\Gamma(A + \frac{1}{2} + X_{2p}) \Gamma(A + \frac{1}{2} + X_{2p-1}) \Gamma(B + \frac{1}{2} + X_{2p}) \Gamma(B + \frac{1}{2} + X_{2p-1})} \\ \times \Gamma(2X_{2p-1}) \Gamma(2X_{2p}) e^{-(X_{2p-1}+X_{2p})r_p + \frac{\zeta}{3}(X_{2p-1}^3+X_{2p}^3)} \mathrm{Pf}_{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \quad (6.1) $$
At large time, we want to eliminate the time factor in the exponential, hence we perform the change of variables

$$X = t^{-1/3} \tilde{X}, \quad r = t^{1/3} \tilde{r}, \quad A + \frac{1}{2} = at^{-1/3}, \quad B + \frac{1}{2} = bt^{-1/3}. \quad (6.2)$$

In the large time limit, the Gamma, cosine and sine functions simplify using that for small positive argument

$$\Gamma(x) \approx \frac{1}{x}, \quad \cos(x) \approx 1, \quad \sin(x) \approx x. \quad (6.3)$$

Under these simplifications, the partition sum at fixed number of strings reads in the limit $t \to +\infty$ (dropping all tildes)

$$Z(n_s, \zeta) = (-1)^{n_s} \prod_{p=1}^{n_s} \int_{\mathbb{R}} dr_p \frac{\zeta}{\zeta + e^{-t^{1/3}r_p}} \iint_{C^2} \frac{dX_{2p-1}}{2i\pi} \frac{dX_{2p}}{2i\pi} \frac{1}{4X_{2p-1}} \frac{a + X_{2p}}{a - X_{2p}} \frac{a + X_{2p-1}}{a - X_{2p-1}} \\ \times \frac{b + X_{2p}}{b - X_{2p}} \frac{b + X_{2p-1}}{b - X_{2p-1}} e^{-(X_{2p-1} + X_{2p})r_p + \frac{x_{2p-1}^3}{3} + \frac{x_{2p}^3}{3}} \text{Pf}_{2n_s \times 2n_s} \left[ \frac{X_i - X_j}{X_i + X_j} \right] \quad (6.4)$$

The contours $C$ have now to be understood as $C = \tilde{a} + i\mathbb{R}$, where $\tilde{a} \in (0, \min\{a, b\})$. We emphasize that the contours all lie at the left of the poles at $X = a$ and $X = b$. We may now use (5.14) and (5.16) to write the Laplace transform $g(\zeta)$ under the scalings in (6.2) as the Pfaffian $\text{Pf}[J - \sigma_\zeta K^{(a,b)}]_{L^2(\mathbb{R})}$. The matrix valued kernel $K^{(a,b)}$ reads in this limit

$$K_{11}^{(a,b)}(r, r') = \frac{1}{4} \iint_{C^2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} \frac{1}{wz} \frac{a+w}{a-w} \frac{a+z}{a-z} \frac{b+w}{b-w} \frac{b+z}{b-z} e^{-rw-r'z+\frac{w^3+z^3}{3}} \\ K_{22}^{(a,b)}(r, r') = \frac{1}{4} \iint_{C^2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} \frac{a+w}{a-w} \frac{a+z}{a-z} \frac{b+w}{b-w} \frac{b+z}{b-z} e^{-rw-r'z+\frac{w^3+z^3}{3}} \\ K_{12}^{(a,b)}(r, r') = \frac{1}{4} \iint_{C^2} \frac{dw}{2i\pi} \frac{dz}{2i\pi} \frac{w-z}{w+z} \frac{1}{w} \frac{a+w}{a-w} \frac{a+z}{a-z} \frac{b+w}{b-w} \frac{b+z}{b-z} e^{-rw-r'z+\frac{w^3+z^3}{3}} \quad (6.5)$$

**Remark 6.1** The kernel $K^{(a,b)}$ has a particular structure, indeed its elements are related through derivative identities: $K_{22}^{(a,b)}(r, r') = \partial_r \partial_{r'} K_{11}^{(a,b)}(r, r')$, $K_{12}^{(a,b)}(r, r') = -\partial_{r'} K_{11}^{(a,b)}(r, r')$ and $K_{22}^{(a,b)}(r, r') = -\partial_r K_{12}^{(a,b)}(r, r')$.

**Remark 6.2** The kernel $K^{(a,b)}$ can be obtained equivalently from the kernel (2.12) by rescaling, as in [56].

Finally, choosing the variable $\zeta$ as $\zeta = e^{-st^{1/3}}$, at large time we have $\lim_{t \to +\infty} \sigma_\zeta(r t^{1/3}) = \Theta(r-s)$, where $\Theta$ is the Theta Heaviside function. The Fredholm Pfaffian formula for the generating function then becomes in the limit

$$\lim_{t \to +\infty} g(\zeta = e^{-st^{1/3}}) = \text{Pf}(J - K^{(a,b)})_{L^2(s, +\infty)}. \quad (6.6)$$

On the other hand, at large time, the Laplace transform of the distribution of the exponential of the KPZ height converges towards the cumulative probability of the height (see [78, Lemma 4.1.39]), i.e.
$$ 
\begin{aligned}
g(\varsigma = e^{-st^{1/3}}) &= \mathbb{E}[\exp(-We^{H(t)-st^{1/3}})] \\
&\approx_{t \to +\infty} \mathbb{E}[\Theta(st^{1/3} - H(t) - \log W)] \\
&\approx_{t \to +\infty} \mathbb{P}\left(\frac{H(t) + \log W}{t^{1/3}} \le s\right)
\end{aligned} 
\quad (6.7) $$

where $\Theta$ is the Theta Heaviside function.

Note that since $W$ is an inverse Gamma random variable with parameter $A + B + 1 = t^{-1/3}(a+b)$, the random variable $\log(W)/t^{1/3}$ weakly converges to an exponential random variable with parameter $a+b$.

At this point, we obtain that for any $a, b > 0$,

$$ \lim_{t \to +\infty} \mathbb{P}\left(\frac{H(t)}{t^{1/3}} \le s - \varepsilon\right) = \mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)} \quad (6.8) $$

where $\mathcal{E}$ is an exponential random variable with parameter $a+b$ independent from $H(t)$, and the matrix kernel $K^{(a,b)}$ is given in (6.5). Using the density of the exponential distribution, and denoting $F^{(a,b)}(s) = \lim_{t \to \infty} \mathbb{P}\left(\frac{H(t)}{t^{1/3}} \le s\right)$, we may rewrite (6.8) as

$$ \int_0^{+\infty} dx (a+b)e^{-x(a+b)} F^{(a,b)}(s-x) = \mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)}. \quad (6.9) $$

Following [105, Eq. (4.3)] (see also Remark 5.1), we differentiate in $s$ in (6.9) and use integration by parts in the left hand side. We obtain

$$ 
\begin{aligned}
\partial_s \mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)} &= (a+b)F^{(a,b)}(s) - (a+b)^2 \int_0^{+\infty} dx e^{-x(a+b)} F^{(a,b)}(s-x) \\
&= (a+b)F^{(a,b)}(s) - (a+b)\mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)}.
\end{aligned} 
$$

Finally, we can write

$$ \lim_{t \to +\infty} \mathbb{P}\left(\frac{H(t)}{t^{1/3}} \le s\right) = \left(1 + \frac{\partial_s}{a+b}\right) \mathrm{Pf}(J - K^{(a,b)})_{\mathbb{L}^2(s, +\infty)} := F^{(a,b)}(s). \quad (6.10) $$

In the next sections, we show that the distribution $F^{(a,b)}$ interpolates between various known distribution as $b$ or $a$ goes to infinity. The most interesting case, corresponding to stationary growth, is when $b, a$ go to zero and will be studied in details in Sect. 7.

**Remark 6.3** Performing the large time limit at any fixed $A > -1/2$ and $B > -1/2$ corresponds – modulo an exchange of limits – to the scaling considered above with $a, b = +\infty$. One can show that in the limit $a, b \to +\infty$, the kernel (6.5) converges to the GSE matrix kernel, by the same manipulations as in [56, Sect. 4.1]. Indeed, as the contours of kernel $K^{(a,b)}$ are parallel to the imaginary axis and cross the real axis between 0 and $\min\{a, b\}$, we can take the limit $a, b \to \infty$ in the integrand without affecting the contours. All rational functions involving the parameter $a, b$ in the large time limit of the kernel $K^{(a,b)}$ in (6.5) converge to the value $-1$. Hence in this limit we obtain the kernel $K^\infty$ given in [56, Eq. (113)], which is precisely the kernel associated to the Gaussian Symplectic Ensemble (GSE) of random matrices as also given in Lemma 2.7. of [44]. Hence this shows that the distribution of the height at $x = 0$ converges at large time for boundary conditions such that $a, b \to \infty$ (e.g. for any fixed $A, B > -1/2$) to the GSE Tracy–Widom distribution, as we will also show below.

**Remark 6.4** In the limit $b \to +\infty$, the kernel $K^{(a,b)}$ in (6.5) converges to the kernel $K^\epsilon$ with $\epsilon = a$ obtained in [56, Eq. (64)] in the study of the droplet initial condition. The
Fredholm Pfaffian $F^{(a)}(s) := Pf(J - K^{(a,+\infty)})_{\mathbb{L}^2(s,+\infty)}$ interpolates between the CDF of the GSE Tracy–Widom distribution, $F_4(s)$, (at $a \to +\infty$) and CDF of the GOE Tracy–Widom distribution function, $F_1(s)$, (at $a = 0$). Note that the GOE kernel obtained in [56, Eq. (83)] was found to provide a new representation of the GOE Tracy–Widom distribution. If instead we take $a \to +\infty$ for fixed $b$, we obtain again the same kernel $K^b$ because of the symmetry $a \leftrightarrow b$. Physically, it describes the CDF of the distribution of the rescaled height of the KPZ equation with Brownian initial data and Dirichlet boundary condition.

# 7 From a Matrix Valued Kernel to a Scalar Kernel

## 7.1 Solution for the KPZ Generating Function at All Times for Generic A, B in Terms of a Scalar Kernel

The general kernel we have obtained in (2.12) has a particular structure in the form of a Schur Pfaffian. With this structure, the kernel verifies the hypothesis of Proposition B.2 of [54] recalled in Lemma B.2 in Appendix B. This proposition states that we can transform the Fredholm Pfaffian of (5.17) which involves a matrix valued kernel, into a Fredholm determinant of a scalar kernel. To proceed, let us first define the functions

$$ 
\begin{aligned}
G(w) &= \frac{\Gamma(A + \frac{1}{2} - w) \Gamma(B + \frac{1}{2} - w)}{\Gamma(A + \frac{1}{2} + w) \Gamma(B + \frac{1}{2} + w)} \Gamma(2w) \\
f_{\text{odd}}(r) &= \int_C \frac{dw}{2i\pi} G(w) \cos(\pi w) e^{-rw + t \frac{w^3}{3}} \\
f_{\text{even}}(r) &= \int_C \frac{dz}{2i\pi} G(z) \frac{\sin(\pi z)}{\pi} e^{-rz + t \frac{z^3}{3}}
\end{aligned} 
\quad (7.1) $$

and the kernel $\bar{K}_{t, \zeta}$ such that for all $(x, y) \in \mathbb{R}_+^2$

$$ 
\begin{aligned}
\bar{K}_{t, \zeta}(x, y) &= 2\partial_x \int_{\mathbb{R}} dr \frac{S}{S + e^{-r}} [f_{\text{even}}(r+x) f_{\text{odd}}(r+y) - f_{\text{odd}}(r+x) f_{\text{even}}(r+y)] \\
&= 2\partial_x \int_{\mathbb{R}} dr \frac{S}{S + e^{-r}} \\
&\quad \times \iint_{C^2} \frac{dw dz}{(2i\pi)^2} G(z) G(w) \frac{\sin(\pi(z-w))}{\pi} e^{-xz-yw-rw-rz+t \frac{w^3+z^3}{3}} \\
&= 2\partial_x \iint_{C^2} \frac{dw dz}{(2i\pi)^2} G(z) G(w) \frac{\sin(\pi(z-w))}{\sin(\pi(z+w))} e^{w+z} e^{-xz-yw+t \frac{w^3+z^3}{3}}
\end{aligned} 
\quad (7.2) $$

Then, the Laplace transform of the one-point distribution of the exponential of the KPZ height admits the following representation:

$$ g(\zeta) = \mathbb{E}[\exp(-\zeta We^{H(t)})] = Pf(J - \sigma_\zeta K)_{\mathbb{L}^2(\mathbb{R})} = \sqrt{\text{Det}(I - \bar{K}_{t, \zeta})_{\mathbb{L}^2(\mathbb{R}_+)}} \quad (7.3) $$

## 7.2 Large Time Limit of the Scalar Kernel

To deduce the large time asymptotics of $H(t)$ from the determinantal formula (7.3), one performs the same rescaling as in Sec. 6, namely one chooses $\zeta = e^{-t^{1/3}s}$ and one rescales $(w, z) \to t^{-1/3}(w, z)$, $r \to t^{1/3}r$. The kernel $\bar{K}_{t, \zeta}$ becomes
$$ \bar{K}^{(a,b)}(x, y) = \frac{1}{2} \iint_{C^2} \frac{dw dz}{(2i\pi)^2} \frac{a+w b+w a+z b+z w-z}{a-w b-w a-z b-z w+z w} \frac{1}{w} e^{-xz-yw+\frac{w^3+z^3}{3}}, \quad (7.4) $$

where the contour $C$ is an upwardly oriented vertical line with real part between 0 and $\min\{a, b\}$ as previously. Then, using the definition of $F^{(a,b)}$ from (6.10), the determinantal formula (7.3) implies the following: for any $a, b > 0$,

$$ F^{(a,b)}(s) = \lim_{t \to +\infty} \mathbb{P} \left( \frac{H(t)}{t^{1/3}} \le s \right) = \left( 1 + \frac{\partial_s}{a+b} \right) \sqrt{\mathrm{Det}(I - \bar{K}^{(a,b)})_{\mathbb{L}^2(s, +\infty)}} $$

where $\bar{K}^{(a,b)}$ is defined in (7.4). Using $\frac{1}{2} \frac{w-z}{(w+z)w} = \frac{1}{w+z} - \frac{1}{2w}$, we obtain that

$$ \bar{K}^{(a,b)}(x, y) = \int_{0}^{+\infty} d\lambda A^{(a,b)}(x + \lambda) A^{(a,b)}(y + \lambda) - \frac{1}{2} A^{(a,b)}(x) \int_{0}^{+\infty} A^{(a,b)}(y + \lambda) d\lambda, \quad (7.5) $$

where the function $A^{(a,b)}(x)$ is defined by

$$ A^{(a,b)}(x) = \int \frac{dz}{2i\pi} \frac{a+z}{a-z} \frac{b+z}{b-z} e^{-xz+z^3/3}, \quad (7.6) $$

where the contour is a vertical line with real part between 0 and $\min\{a, b\}$. Note that the function $A^{(a,b)}$ has exponential decay at $+\infty$, that is for any $c \in (0, \min\{a, b\})$, there exist $C \in \mathbb{R}$ such that $|A^{(a,b)}(x)| \le Ce^{-cx}$. Let us introduce an operator $\hat{A}_s$ acting on $\mathbb{L}^2(0, +\infty)$ with kernel

$$ \hat{A}_s(x, y) = A^{(a,b)}(x + y + s), \quad (7.7) $$

and an operator $\bar{K}_s^{(a,b)}$ acting on $\mathbb{L}^2(0, +\infty)$ with kernel $\bar{K}_s^{(a,b)}(x, y) := \bar{K}^{(a,b)}(x+s, y+s)$.

**Claim 7.1** For any $s \in \mathbb{R}$, and $a, b > 0$,

$$ \sqrt{\mathrm{Det}(I - \bar{K}_s^{(a,b)})} = \frac{1}{2} \left( \mathrm{Det}(I - \hat{A}_s) + \mathrm{Det}(I + \hat{A}_s) \right), \quad (7.8) $$

where all operators act on $\mathbb{L}^2(0, +\infty)$.

**Proof** Equation 7.5 implies that, as operators acting on $\mathbb{L}^2(0, +\infty)$, we have

$$ \bar{K}_s^{(a,b)} = \hat{A}_s^2 - \frac{1}{2} |\hat{A}_s \delta\rangle \langle 1 \hat{A}_s| $$

At this point, we recognize that the operator $\bar{K}_s^{(a,b)}$ has the same structure as in [87] and we may apply the same steps as Equations (32)–(35) therein. More precisely, we may use the matrix determinant Lemma to obtain that

$$ \mathrm{Det}(I - \bar{K}_s) = \mathrm{Det}(I - \hat{A}_s^2) \left( 1 + \frac{1}{2} \langle 1 | \frac{\hat{A}_s^2}{I - \hat{A}_s^2} | \delta \rangle \right). \quad (7.9) $$

Then, we use the decomposition

$$ \frac{\hat{A}_s^2}{I - \hat{A}_s^2} = -I + \frac{1}{2} \left( \frac{I}{I - \hat{A}_s} + \frac{I}{I + \hat{A}_s} \right), $$

and recall that $\langle 1 | I | \delta \rangle = 1$. Using (the proof of) [109, Proposition 1] we have

$$ \langle 1 | \frac{I}{I \pm \hat{A}_s} | \delta \rangle = \frac{\mathrm{Det}(I \mp \hat{A}_s)}{\mathrm{Det}(I \pm \hat{A}_s)}. \quad (7.10) $$
Plugging (7.10) into (7.9) yields the statement of the Lemma. □

**Remark 7.2** The identity (7.10) is true for any kernel of the type $B_s(x, y) = B(x + y + s)$ such that $B$ has sufficient decay at $+\infty$ so that $\text{Det}(I \pm B_s)$ and $\langle 1|\frac{I}{I+B_s}|\delta \rangle$ converges to 1 as $s$ goes to $+\infty$ (see the proof of [109, Proposition 1] for details). In our case, this condition is satisfied due to the exponential decay of the function $A^{(a,b)}$ for fixed $a, b > 0$).

Hence one has

$$F^{(a,b)}(s) = \frac{1}{2} \left( 1 + \frac{\partial_s}{a+b} \right) \left( \text{Det}(I - \hat{A}_s) + \text{Det}(I + \hat{A}_s) \right) \quad (7.11)$$

## 7.3 Limit $a, b \to 0$: The Critical Stationary Case

### 7.3.1 CDF in Terms of a Fredholm Determinant

Moving the contour to the right in the definition of $A^{(a,b)}$ in (7.6), we obtain

$$A^{(a,b)}(x) = \tilde{A}^{(a,b)}(x) + 2 \frac{a+b}{a-b} (h_b(x) - h_a(x)), \quad (7.12)$$

with $h_b(x) = be^{-xb+b^3/3}$ and

$$\tilde{A}^{(a,b)}(x) = \int \frac{dz}{2i\pi} \frac{a+z}{a-z} \frac{b+z}{b-z} e^{-xz+z^3/3} = \text{Ai}(x) + 2(a+b) \int_x^{+\infty} d\lambda \, \text{Ai}(\lambda) + \mathcal{O}(b^2, a^2, ba), \quad (7.13)$$

where in the integral over $z$, the contour passes to the right of $b, a$. We introduce the operator $\tilde{A}_s$ acting on $\mathbb{L}^2(0, +\infty)$ with kernel

$$\tilde{A}_s(x, y) = \tilde{A}^{(a,b)}(s + x + y).$$

We thus write

$$\hat{A}_s(x, y) = \tilde{A}_s(x, y) + 2 \frac{a+b}{a-b} (|f_b(x)\rangle\langle f_b(y)| - |f_a(x)\rangle\langle f_a(y)|), \quad (7.14)$$

with $f_b(x) = \sqrt{b}e^{b^3/6-bs/2-bx}$. Using the matrix determinant lemma, we have

$$\text{Det}(I \mp \hat{A}_s) = \text{Det}(I \mp \tilde{A}_s) \left( \left( 1 \mp 2 \frac{a+b}{a-b} I_{b,b} \right) \left( 1 \pm 2 \frac{a+b}{a-b} I_{a,a} \right) + 4 \left( \frac{a+b}{a-b} \right)^2 I_{b,a} I_{a,b} \right) \quad (7.15)$$

where

$$I_{\alpha, \beta} = \langle f_{\alpha} | f_{\beta} \rangle \pm \frac{\tilde{A}_s}{I \mp \tilde{A}_s} |f_{\beta}\rangle. \quad (7.16)$$

We now consider the limit $b, a \to 0$ with a fixed arbitrary ratio $r = a/b$. We use the exact expressions for the scalar products

$$\langle f_{\alpha} | f_{\beta} \rangle = \frac{\sqrt{\alpha \beta}}{\alpha + \beta} e^{\frac{\alpha^3}{6} + \frac{\beta^3}{6} - \frac{\alpha + \beta}{2} s} \quad (7.17)$$
as well as

$$ \langle f_{\alpha} | \frac{\tilde{A}_s}{I \mp \tilde{A}_s} | f_{\beta} \rangle = \sqrt{\alpha \beta} e^{\frac{\alpha^3}{6} + \frac{\beta^3}{6}} e^{-\frac{\alpha+\beta}{2}s} \langle e^{-\alpha x} | \frac{\tilde{A}_s}{I \mp \tilde{A}_s} | e^{-\beta x} \rangle \quad (7.18) $$

One has

$$ I_{b,b} = e^{-bs} \left( \frac{1}{2} \pm b R_s^\mp + \mathcal{O}(b^2) \right) , \quad R_s^\mp = \left\langle 1 \left| \frac{\tilde{A}_s}{I \mp \tilde{A}_s} \right| 1 \right\rangle \quad (7.19) $$

$$ I_{a,a} = e^{-as} \left( \frac{1}{2} \pm a R_s^\mp + \mathcal{O}(a^2) \right) \quad (7.20) $$

$$ I_{b,a} = I_{a,b} = \sqrt{ab} e^{-(b+a)s/2} \left( \frac{1}{a+b} \pm R_s^\mp + \mathcal{O}(a, b) \right) \quad (7.21) $$

Plugging these asymptotics in (7.15) we obtain

$$ \mathrm{Det}(I - \hat{A}_s) = \mathcal{O}(b^2) \quad (7.22) $$

$$ \mathrm{Det}(I + \hat{A}_s) = \mathrm{Det}(I + \tilde{A}_s) \times 2b(1+r)(2R_s^+ + s) + \mathcal{O}(b^2). \quad (7.23) $$

Furthermore, if we keep only the first order in $b$, we may replace $\tilde{A}_s(x, y)$ by $\mathrm{Ai}(x + y + s)$ as $a, b \to 0$. Thus, we have found that

$$ F(s) := F^{0,0}(s) = \partial_s [\mathrm{Det}(I + \mathrm{Ai}_s)(2R_s^{+,0} + s)] \quad (7.24) $$

where

$$ R_s^{+,0} = \left\langle 1 \left| \frac{\mathrm{Ai}_s}{I + \mathrm{Ai}_s} \right| 1 \right\rangle \quad (7.25) $$

and $\mathrm{Ai}_s$ is the operator with kernel $\mathrm{Ai}(x + y + s)$. Using the Sherman-Morrison formula, we have that

$$ \left\langle 1 \left| \frac{\mathrm{Ai}_s}{I + \mathrm{Ai}_s} \right| 1 \right\rangle = \left\langle 1 \left| \frac{I}{I + \mathrm{Ai}_s} \right| \mathrm{Ai}_s 1 \right\rangle = \frac{\mathrm{Det}(I + \mathrm{Ai}_s + |\mathrm{Ai}_s 1\rangle\langle 1|)}{\mathrm{Det}(I + \mathrm{Ai}_s)} - 1. \quad (7.26) $$

which yields the following alternative formula in terms of Fredholm determinants

$$ F(s) = \partial_s [2\mathrm{Det}(I + \mathrm{Ai}_s + |\mathrm{Ai}_s 1\rangle\langle 1|) + (s-2)\mathrm{Det}(I + \mathrm{Ai}_s)] . \quad (7.27) $$

which is given in the Introduction in (2.26).

### 7.3.2 CDF in Terms of the Solution of the Painlevé II Equation

In this Section, we show that the distribution $F(s)$ can also be written in terms of the Tracy-Widom distributions $F_1(s)$ and $F_2(s)$ only. We also provide formulae in terms of the Hastings-McLeod solution of Painlevé II equation (see Appendix C). Define $q(s)$ to be the solution of the Painlevé II equation for $s \in \mathbb{R}$,

$$ q''(s) = sq(s) + 2q(s)^3 \quad (7.28) $$

which satisfies the asymptotic condition $q(s) \sim_{s \to +\infty} \mathrm{Ai}(s)$. This solution is called the Hastings-McLeod solution of the Painlevé II equation [110].

Let us first recall some results from [111]. Combining the formula 4.18 and the one just below 4.52 in [111] we can write

$$ \left\langle 1 \left| \frac{I}{I + \mathrm{Ai}_s} \right| \delta \right\rangle = e^{-\int_s^{+\infty} dt q(t)}, \quad (7.29) $$
**Remark 7.3** One has, see [111, Eq. 4.22], $q(s) = \langle \delta | \frac{\text{Ai}_s}{1-K\text{Ai}_s} | \delta \rangle$.

Another result is obtained combining Remark 7.2 and (7.29) as

$$ \langle 1 | \frac{I}{I \pm \text{Ai}_s} | \delta \rangle = \frac{\text{Det}(I \mp \text{Ai}_s)}{\text{Det}(I \pm \text{Ai}_s)} = e^{\pm \int_s^{+\infty} dt \, q(t)} = \left( \frac{F_1(s)^2}{F_2(s)} \right)^{\pm 1}. \quad (7.30) $$

where the last identity is from [112].

The next result allows to rewrite the resolvant $R_s^{+,0}$ appearing in (7.24) and (7.25) in terms of $q$.

**Proposition 7.4** We have

$$ \langle 1 | \frac{\text{Ai}_s}{I + \text{Ai}_s} | 1 \rangle = \frac{1}{2} \left( \int_{-\infty}^{s} dr \, e^{-2 \int_{r}^{+\infty} dt \, q(t)} - s \right). \quad (7.31) $$

**Proof** Before proving this proposition, we need the following known lemma.

**Lemma 7.5** ([109, Lemma 3]) We have the following relation

$$ \partial_s \frac{\text{Ai}_s}{I + \text{Ai}_s} = - \frac{\text{Ai}_s}{I - K_{\text{Ai},s}} D - \frac{\text{Ai}_s}{I - K_{\text{Ai},s}} |\delta\rangle \langle\delta| \frac{I}{I + \text{Ai}_s} \quad (7.32) $$

where $D$ is the derivative operator.

Since $D|1\rangle = 0$ we have

$$ \begin{aligned} \partial_s \langle 1 | \frac{\text{Ai}_s}{I + \text{Ai}_s} | 1 \rangle &= - \langle 1 | \frac{\text{Ai}_s}{I - K_{\text{Ai},s}} | \delta \rangle \langle \delta | \frac{I}{I + \text{Ai}_s} | 1 \rangle \\ &= -\frac{1}{2} \langle 1 | \left( \frac{I}{I - \text{Ai}_s} - \frac{I}{I + \text{Ai}_s} \right) | \delta \rangle \langle \delta | \frac{I}{I + \text{Ai}_s} | 1 \rangle \\ &= -\frac{1}{2} \left( \frac{\text{Det}(I + \text{Ai}_s)}{\text{Det}(I - \text{Ai}_s)} - \frac{\text{Det}(I - \text{Ai}_s)}{\text{Det}(I + \text{Ai}_s)} \right) \times \frac{\text{Det}(I - \text{Ai}_s)}{\text{Det}(I + \text{Ai}_s)} \\ &= -\frac{1}{2} \left( 1 - \frac{F_1(s)^4}{F_2(s)^2} \right) \end{aligned} \quad (7.33) $$

Using the expressions of $F_1$ and $F_2$ in terms of the Painlevé II equation, we have

$$ \partial_s \langle 1 | \frac{\text{Ai}_s}{I + \text{Ai}_s} | 1 \rangle = -\frac{1}{2} \left(1 - e^{-2 \int_s^{+\infty} dt \, q(t)}\right) \quad (7.34) $$

Recalling the notation $R_s^{+,0} = \langle 1 | \frac{\text{Ai}_s}{1+\text{Ai}_s} | 1 \rangle$, we have

$$ \partial_s (2R_s^{+,0} + s) = e^{-2 \int_s^{+\infty} dt \, q(t)} = \frac{F_1(s)^4}{F_2(s)^2} \quad (7.35) $$

We integrate the last quantity between $-\infty$ and $s$ using that $q(s) \to +\infty$ for $s \to -\infty$ as in (C.4) and obtain

$$ 2R_s^{+,0} + s = \int_{-\infty}^{s} dr \, e^{-2 \int_{r}^{+\infty} dt \, q(t)} + \kappa. \quad (7.36) $$

Since $R_s^{+,0} \to 0$ in the limit $s \to +\infty$, and using (C.2) and the asymptotics (C.5) and (C.6), we obtain that $\kappa = 0$, which concludes the proof. □
Using Proposition 7.4 we arrive at the following equivalent expressions for $F(s)$ (defined in (7.24)). The first one reads

$$F(s) = \partial_s \left[ \frac{F_2(s)}{F_1(s)} \int_{-\infty}^{s} dt \frac{F_1(t)^4}{F_2(t)^2} \right]. \qquad (7.37)$$

where the first line was given in (2.28). The second formula is expressed in terms of the Hastings–McLeod solution of the Painlevé II equation (see Appendix C) and reads

$$ \begin{aligned} F(s) &= \partial_s \left[ e^{-\frac{1}{2} \int_s^{+\infty} dr [(r-s)q(r)^2 - q(r)]} \int_{-\infty}^{s} dr e^{-2 \int_r^{+\infty} dt q(t)} \right] \\ &= \partial_s \left[ e^{-\frac{1}{2} \int_s^{+\infty} dr [(r-s)q(r)^2 + 3q(r)]} (s - 2q'(s) + 2q(s)^2) \right]. \end{aligned} \qquad (7.38) $$

which was given in (2.29). A third formula, useful for the asymptotics, is obtained applying the derivative in front of the bracket and using Eqs. (C.2) and (C.3)

$$ F(s) = e^{-\frac{1}{2} \int_s^{+\infty} dr [(r-s)q(r)^2 + 3q(r)]} \times \left( 1 + \frac{1}{2}(q'(s))^2 - sq(s)^2 - q(s)^4 - q(s))(s - 2q'(s) + 2q(s)^2) \right) \qquad (7.39) $$

### 7.3.3 Properties of $F(s)$: First Moments

We check in Appendix C using the formulae (7.39) that the function $F(s)$ has the behaviour at $s \to \pm\infty$ that is required for a CDF, i.e. its limit at $s = -\infty$ is 0 and its limit at $s = +\infty$ is 1. The detailed asymptotics for $s \to \pm\infty$ is performed in Appendix D. The CDF takes the form $F(s) = \partial_s[F(s)]$. Provided $F(s) \to 0$ sufficiently fast for $s \to -\infty$ and $F(s) - s \to 0$ sufficiently fast for $s \to +\infty$, conditions which can be checked from Appendix D, integration by parts give the following formula for the k-th positive integer moment $M_k$ of the distribution $F(s)$

$$ M_k = k(k-1) \int_{\mathbb{R}} (F(s) - \max(s, 0)) s^{k-2} ds \qquad (7.40) $$

This formula can be used to obtain the moments and the cumulants through a numerical evaluation of $F(s)$. One notices that the mean vanishes, $M_1 = 0$, which indeed must be the case since $\mathbb{E}[H(t)] = 0$ in the critical stationary case (see Sect. 4.3 where we have computed the expectation of $h(x, t)$ using the stationary structure of the log-gamma polymer). We used two numerical methods to evaluate $F(s)$. The first one uses Eq. (7.27) where the Fredholm determinants are calculated using the method described in Ref. [76,113]. The second method uses the formula (7.37) and uses the Mathematica routines for $F_{1,2}(s)$, and is in agreement with the first one. The CDF $F(s)$ and its derivative $F'(s)$ are plotted in Fig. 3. The mean, variance, skewness and excess kurtosis are given in Table 1.

### 7.4 Limit $a, b \to +\infty$: Convergence to the GSE

As we already discussed in Remark 6.3 the limit $a, b \to +\infty$ can be performed on the Fredholm pfaffian formula and leads to GSE Tracy–Widom fluctuations. This limit can also be performed on the formula (7.11) involving the scalar kernel $\hat{A}_s$ defined in (7.7) in terms of the function $A^{(a,b)}(x)$ defined in (7.6). It is clear from the definition of $A^{(a,b)}(x)$ that if $a, b \to +\infty$ simultaneously, then $A^{(a,b)}(x)$ converges to the standard Airy function. The
CDF $F^{(a,b)}(s)$ in (7.11), for $a, b \to +\infty$ then takes the form of the GSE Tracy–Widom distribution found in [53, Eq (35)]. This result thus matches smoothly with the result (2.16) valid for any fixed $A, B > -1/2$ in the large time limit.

## 7.5 Limit (a, b) → (0, +∞): Convergence to the GOE

Another interesting limit, that we call $F^{(a)}(s) = \lim_{b \to +\infty} F^{(a,b)}(s)$, is the limit $b \to +\infty$ at fixed $a$ and by the $A \leftrightarrow B$ symmetry, the case $A \to +\infty$ at fixed $b$ is similar. In particular we consider now the limit $a \to 0$.

The manipulations follow closely the ones of Section. 7.3. We start by moving the contour to the right in the definition of $A^{(a,\infty)}$ in (7.6) already taking into account the $b \to +\infty$ limit. We obtain $A^{(a,\infty)}(x) = \tilde{A}^{(a,\infty)}(x) + 2h_a(y)$, with $h_a(x) = ae^{-xa+a^3/3}$ and

$$ \tilde{A}^{(a, \infty)}(x) = \int \frac{dz}{2i\pi} \frac{a+z}{a-z} e^{-xz+z^3/3} = -\mathrm{Ai}(x) - 2a \int_x^{+\infty} d\lambda \, \mathrm{Ai}(\lambda) + \mathcal{O}(a^2), \quad (7.41) $$

where in the integral over $z$, the contour passes to the right of $a$. We introduce the operator $\tilde{A}_s$ acting on $\mathbb{L}^2(0, +\infty)$ with kernel $\tilde{A}_s(x, y) = \tilde{A}^{(a,\infty)}(s+x+y)$. We thus write

$$ \hat{A}_s(x, y) = \tilde{A}_s(x, y) + 2|f_a(x)\rangle\langle f_a(y)|, \quad (7.42) $$

with $f_a(x) = \sqrt{a}e^{a^3/6-\frac{a^2}{2}-ax}$. Using the matrix determinant lemma, we have

$$ \mathrm{Det}(I \pm \hat{A}_s) = \mathrm{Det}(I \pm \tilde{A}_s) \left( 1 \pm 2\langle f_a | f_a \rangle - 2 \left\langle f_a \left| \frac{\tilde{A}_s}{I \pm \tilde{A}_s} \right| f_a \right\rangle \right) \quad (7.43) $$

We now consider the limit $a \to 0$ and we use the exact expressions for the scalar product $\langle f_a | f_a \rangle = \frac{1}{2}e^{\frac{a^3}{3}-as}$ as well as

$$ \left\langle f_a \left| \frac{\tilde{A}_s}{I \pm \tilde{A}_s} \right| f_a \right\rangle = a e^{\frac{a^3}{3}-as} \left\langle e^{-ax} \left| \frac{\tilde{A}_s}{I \pm \tilde{A}_s} \right| e^{-ax} \right\rangle \quad (7.44) $$

One has

$$ \mathrm{Det}(I \pm \hat{A}_s) = \mathrm{Det}(I \pm \tilde{A}_s) \left( 1 \pm e^{-as} - 2e^{-as}a \left\langle 1 \left| \frac{\tilde{A}_s}{I \pm \tilde{A}_s} \right| 1 \right\rangle + \mathcal{O}(a^2) \right) \quad (7.45) $$

Looking at formula (7.11) in the limit $b \to +\infty$ we see that we need only the following estimates from (7.45) up to $\mathcal{O}(a)$ as $a \to 0$

$$ \mathrm{Det}(I - \hat{A}_s) = \mathcal{O}(a) \quad (7.46) $$

$$ \mathrm{Det}(I + \hat{A}_s) = \mathrm{Det}(I - \mathrm{Ai}_s)(2 + \mathcal{O}(a)). \quad (7.47) $$

which, inserted in (7.11), lead to

$$ F^{(0)}(s) := F^{(0,\infty)}(s) = F^{(\infty,0)}(s) = \mathrm{Det}(I - \mathrm{Ai}_s) = F_1(s) \quad (7.48) $$

This coincides with the determinantal representation of the GOE Tracy–Widom CDF. This formula matches smoothly with the result (2.17) which states that the CDF of the one-point KPZ height field is given by the GOE Tracy–Widom distribution for $A = -1/2$ and any $B > -1/2$.
**Acknowledgements** Open access funding provided by Scuola Internazionale Superiore di Studi Avanzati - SISSA within the CRUI-CARE Agreement. AK and PLD acknowledge support from ANR grant ANR-17-CE30-0027-01 RaMaTraF. AK acknowledges support from ERC under Consolidator Grant Number 771536 (NEMO).

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

# Appendix

## A Overlap of the Half-Line Bethe States with the Brownian Initial Condition

Here we give some details on the calculation of the overlap $\langle\Psi_{\mu}|\Phi_0\rangle$ between the half-line Bethe states and the Brownian initial condition. We recall that $\Phi_0$, given in (3.5), is a fully symmetric function of its arguments, and that in the sector $0 \le x_1 \le \dots \le x_n$ it equals

$$ \Phi_0(x_1, \dots, x_n) = \exp \left( \frac{1}{2} \sum_{j=1}^{n} (2n - 2j + 1)x_j - (1/2 + B)x_j \right), \quad (A.1) $$

where b is the drift of the Brownian. Since we will find that the overlap is real, we will instead calculate its complex conjugate and use that $\langle\Psi_{\mu}|\Phi_0\rangle^* = \langle\Phi_0|\Psi_{\mu}\rangle = \langle\Psi_{\mu}|\Phi_0\rangle$. Since $\Psi_{\mu}$ is also a symmetric function of its arguments, by definition the overlap can thus be written as

$$ \langle\Phi_0|\Psi_\mu\rangle = n! \int_{0<y_1<y_2<\dots<y_n} dy_1 \dots dy_n \Psi_\mu(y_1, \dots, y_n) e^{\sum_{j=1}^n \left(\frac{1}{2}(2n+1-2j)y_j - (1/2+B)y_j\right)} \quad (A.2) $$

Inserting the explicit form of the Bethe eigenstate (3.7) as a superposition of plane waves we obtain

$$ \begin{aligned} \langle\Phi_0|\Psi_\mu\rangle = \frac{n!}{(2\mathbf{i})^n} & \sum_{P \in S_n} \prod_{p=1}^{n} \left( \sum_{\epsilon_p = \pm 1} \epsilon_p \right) \prod_{\ell=1}^{n} \left( 1 + \mathbf{i} \frac{\epsilon_\ell \lambda_P(\ell)}{A} \right) \\ & \times \prod_{k<\ell} \left( 1 + \frac{\mathbf{i}}{\epsilon_\ell \lambda_P(\ell) - \epsilon_k \lambda_P(k)} \right) \left( 1 + \frac{\mathbf{i}}{\epsilon_\ell \lambda_P(\ell) + \epsilon_k \lambda_P(k)} \right) \\ & G_{n,w}(\epsilon_1 \lambda_P(1), \dots, \epsilon_n \lambda_P(n)) \end{aligned} \quad (A.3) $$

where we have defined the integrals

$$ G_{n,w}(\lambda_1, \dots, \lambda_n) = \int_{0<y_1<y_2<\dots<y_n} dy_1 \dots dy_p e^{\sum_{j=1}^n \left(-B - \frac{1}{2} + i\lambda_j\right)y_j + \frac{1}{2}(2n+1-2j)y_j} \quad (A.4) $$
These integrals can be explicitly evaluated

$$G_{n,w}(\lambda_1, \dots, \lambda_n) = \prod_{j=1}^{n} \frac{-1}{-j(B + 1/2) + i\lambda_n + \dots + i\lambda_{n+1-j} + j^2/2} \quad (A.5)$$

Now in (A.3) for each permutation $P$ we can relabel all the $\varepsilon_p \to \varepsilon_{P(p)}$ and denoting by $\sum_{\varepsilon=\{\pm 1\}^n}$ the operation of summation over all the variables $\varepsilon_i$ (an operation independent of their labeling) we can rewrite (A.3) as

$$\langle \Phi_0 | \Psi_\mu \rangle = \frac{n!}{(2i)^n} \sum_{\varepsilon=\{\pm 1\}^n} \prod_{\ell=1}^{n} \varepsilon_\ell \left(1 + i \frac{\varepsilon_\ell \lambda_\ell}{A}\right) \prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_k \lambda_k + \varepsilon_\ell \lambda_\ell}\right) \\ \sum_{P \in S_n} \prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_{P(\ell)} \lambda_{P(\ell)} - \varepsilon_{P(k)} \lambda_{P(k)}}\right) G_{n,w}(\varepsilon_{P(1)} \lambda_{P(1)}, \dots, \varepsilon_{P(n)} \lambda_{P(n)}) \quad (A.6)$$

where we have used the fact that the products

$$\prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_{P(\ell)} \lambda_{P(\ell)} + \varepsilon_{P(k)} \lambda_{P(k)}}\right) = \prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_k \lambda_k + \varepsilon_\ell \lambda_\ell}\right) \quad (A.7)$$

and

$$\prod_{\ell=1}^{n} \left( 1 + i \frac{\varepsilon_{P(\ell)} \lambda_{P(\ell)}}{A} \right) = \prod_{\ell=1}^{n} \left( 1 + i \frac{\varepsilon_{\ell} \lambda_{\ell}}{A} \right) \quad (A.8)$$

are independent of the permutation $P$. Now we use the following symmetrization identity, given in [114] (this is a limit of [115, Eq. (9)] which was a slight generalization of [116, Eq. (1.6)]),

$$\sum_{P \in S_n} \prod_{k<\ell} \left(1 + \frac{i}{\lambda_{P(\ell)} - \lambda_{P(n)}}\right) G_{n,w}(\lambda_{P(1)}, \dots, \lambda_{P(n)}) = \prod_{j=1}^{n} \frac{1}{B - i\lambda_j} \quad (A.9)$$

Applying it to the set $\{\varepsilon_k \lambda_k\}$ we obtain

$$\langle \Phi_0 | \Psi_\mu \rangle = \frac{n!}{(2i)^n} \sum_{\varepsilon=\{\pm 1\}^n} \prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_k \lambda_k + \varepsilon_\ell \lambda_\ell}\right) \prod_{j=1}^{n} \frac{\varepsilon_j + i \frac{\lambda_j}{A}}{B - i \varepsilon_j \lambda_j} \quad (A.10)$$

It remains to perform the symmetrization over $\varepsilon$.

**Lemma A.1** For any set of complex numbers $\{\lambda_j\}_{1 \le j \le n}$,

$$\sum_{\varepsilon=\{\pm 1\}^n} \prod_{k<\ell} \left(1 + \frac{i}{\varepsilon_k \lambda_k + \varepsilon_\ell \lambda_\ell}\right) \prod_{j=1}^{n} \frac{\varepsilon_j + i \frac{\lambda_j}{A}}{B - i \varepsilon_j \lambda_j} = \left(\frac{2i}{A}\right)^n \frac{\Gamma(A + B + 1)}{\Gamma(A + B - n + 1)} \prod_{j=1}^{n} \frac{\lambda_j}{B^2 + \lambda_j^2} \quad (A.11)$$

**Remark A.2** In the limit $A \to +\infty$, this formula yields back the overlap for Dirichlet boundary condition, see [54].

**Proof** This identity was essentially known in the context of Hall-Littlewood polynomials of type BC. In particular, a generalization of it was proved in [104, Theorem 2.6]. We refer to [55, Eq. (54)] for more details about how to degenerate Venkateswaran's symmetrization
identity to the one we need (trigonometric to rational limit). For any parameters A, B, C ∈ C and complex variables $(z_i)_{1 \le i \le n}$, we have

$$ \sum_{\epsilon \in \{\pm 1\}^n} \sum_{P \in S_n} \prod_{k<l} \frac{\epsilon_{P(k)} z_{P(k)} + \epsilon_{P(l)} z_{P(l)} - C}{\epsilon_{P(k)} z_{P(k)} + \epsilon_{P(l)} z_{P(l)}} \frac{\epsilon_{P(k)} z_{P(k)} - \epsilon_{P(l)} z_{P(l)} - C}{\epsilon_{P(k)} z_{P(k)} - \epsilon_{P(l)} z_{P(l)}} \\ \times \prod_{j=1}^{n} \frac{(\epsilon_{P(j)} z_{P(j)} + A)(\epsilon_{P(j)} z_{P(j)} + B)}{\epsilon_{P(j)} z_{P(j)}} = 2^n n! \prod_{j=0}^{n-1} (A + B - j). \quad (A.12) $$

Let us perform first the symmetrization over $P \in S_n$. Since the first product in the left hand side of (A.12) is $S_n$-invariant, and using the symmetrization identity ([117, Chap. III, Eq. (1.4)])

$$ \sum_{P \in S_n} \frac{\epsilon_{P(k)} z_{P(k)} - \epsilon_{P(l)} z_{P(l)} - C}{\epsilon_{P(k)} z_{P(k)} - \epsilon_{P(l)} z_{P(l)}} = n!, $$

we obtain that

$$ \sum_{\epsilon \in \{\pm 1\}^n} \prod_{k<l} \frac{\epsilon_k z_k + \epsilon_l z_l - C}{\epsilon_k z_k + \epsilon_l z_l} \prod_{j=1}^{n} \frac{(\epsilon_j z_j + A)(\epsilon_j z_j + B)}{\epsilon_j z_j} = 2^n \prod_{j=0}^{n-1} (A + B - j). \quad (A.13) $$

We obtain (A.11) using the substitutions $C \to 1$, $z_j \to i\lambda_j$ in (A.13). □

Applying Lemma A.1 in (A.10), we obtain the formula for the overlap (3.17) given in the text.

We note that the overlap formula it is a priori valid before the insertion of the solution of the Bethe equations, i.e. it is valid for any set of complex $\lambda_j$ such that the overlap integral converges. The condition for that to be true can be read from (A.9) as

$$ \mathrm{Re}(i\lambda_j) < B \quad (A.14) $$

for all $j \in [1, n]$. Inserting a string state, labeled by $\{k_j, m_j\}_{j=1,...,n_s}$, the condition becomes $\max_j \frac{1}{2}(2m_j - 1) \le \frac{1}{2}(2n - 1) < B$, which leads to the condition $\frac{n}{2} < B + \frac{1}{2}$.

## B From Two-Dimensional Kernels to Scalar Kernels

We present in this section an equivalent representation of a class of Fredholm Pfaffians with 2 × 2 block kernels in terms of a Fredholm determinant with a scalar valued kernel. Consider a measure $d\mu$ on a contour C in the complex plane and another measure $d\nu_\zeta$ on the real line $\mathbb{R}$, depending on a real parameter $\zeta$. Consider the quantity $g(\zeta)$ defined by the series

$$ g(\zeta) = 1 + \sum_{n_s=1}^{\infty} \frac{(-1)^{n_s}}{n_s!} Z(n_s, \zeta) \quad (B.1) $$

and

$$ Z(n_s, \zeta) = \prod_{p=1}^{n_s} \int_{\mathbb{R}} d\nu_{\zeta}(r_p) \iint_{C^2} d\mu(X_{2p-1})d\mu(X_{2p}) \\ \phi_{\mathrm{odd}}(X_{2p-1})\phi_{\mathrm{even}}(X_{2p})e^{-r_p[X_{2p-1}+X_{2p}]} \mathrm{Pf} \left[ \frac{X_i - X_j}{X_i + X_j} \right]_{i,j=1}^{2n_s} \quad (B.2) $$
Then, following Ref. [54], we have the following equivalent representations for g(ζ).

**Lemma B.1** (Fredholm Schur Pfaffian) g(ζ) is equal to a Fredholm Pfaffian with a 2 × 2 matrix valued skew-symmetric kernel

$$g(\zeta) = \mathrm{Pf}(J - K)_{\mathbb{L}^2(\mathbb{R}, v_\zeta)} \quad (B.3)$$

For $(r, r') \in \mathbb{R}^2$ the matrix kernel K is given by

$$ 
\begin{align}
K_{11}(r, r') &= \iint_{C^2} d\mu(v) d\mu(w) \frac{v-w}{v+w} \phi_{\mathrm{odd}}(v) \phi_{\mathrm{odd}}(w) e^{-rv-r'w} \nonumber \\
K_{22}(r, r') &= \iint_{C^2} d\mu(v) d\mu(w) \frac{v-w}{v+w} \phi_{\mathrm{even}}(v) \phi_{\mathrm{even}}(w) e^{-rv-r'w} \nonumber \\
K_{12}(r, r') &= \iint_{C^2} d\mu(v) d\mu(w) \frac{v-w}{v+w} \phi_{\mathrm{odd}}(v) \phi_{\mathrm{even}}(w) e^{-rv-r'w} \nonumber \\
K_{21}(r, r') &= \iint_{C^2} d\mu(v) d\mu(w) \frac{v-w}{v+w} \phi_{\mathrm{even}}(v) \phi_{\mathrm{odd}}(w) e^{-rv-r'w} \tag{B.4}
\end{align}
$$

and the matrix kernel J is defined by $J(r, r') = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} 1_{r=r'}$.

**Lemma B.2** (Scalar Fredholm determinant - Proposition B.2 of Ref. [54]) g(ζ) is equal to the square root of a Fredholm determinant with scalar valued kernel

$$g(\zeta) = \sqrt{\mathrm{Det}(I - \bar{K})_{\mathbb{L}^2(\mathbb{R}_+)}} \quad (B.5)$$

where $\mathbb{L}^2(\mathbb{R}_+)$ is considered with the Lebesgue measure on $\mathbb{R}_+$. Introducing the functions

$$f_{\mathrm{odd}}(r) = \int_C d\mu(v) \phi_{\mathrm{odd}}(v) e^{-rv}, \quad f_{\mathrm{even}}(r) = \int_C d\mu(v) \phi_{\mathrm{even}}(v) e^{-rv} \quad (B.6)$$

which are assumed to be in $\mathbb{L}^2(\mathbb{R}_+)$, the scalar kernel $\bar{K}$ is given, for $(x, y) \in \mathbb{R}_+^2$, by

$$\bar{K}(x, y) = 2\partial_x \int_{\mathbb{R}} dv_\zeta(r) [f_{\mathrm{even}}(x+r)f_{\mathrm{odd}}(r+y) - f_{\mathrm{odd}}(x+r)f_{\mathrm{even}}(r+y)] \quad (B.7)$$

and the scalar kernel I is the identity kernel $I(x, y) = 1_{x=y}$.

# C Painlevé II Equation and Tracy–Widom Distributions

## C.1 The Hastings–McLeod Solution of the Painlevé II Equation

Define $q(s)$ to be the solution of the Painlevé II equation for $s \in \mathbb{R}$,

$$q''(s) = sq(s) + 2q(s)^3 \quad (C.1)$$

which satisfies the asymptotic condition $q(s) \sim_{s \to +\infty} \mathrm{Ai}(s)$. The unique smooth solution of the Painlevé II equation with that asymptotic condition is called the Hastings–McLeod solution [110]. Let us indicate the two following formula which we use in the manuscript. The first one is given in Ref. [118, Eq. (2.18)]

$$s - 2q'(s) + 2q(s)^2 = e^2 \int_{s}^{\infty} dr \, q(r) \int_{-\infty}^{s} dr \, e^{-2 \int_{r}^{\infty} dr \, q(r)} \quad (C.2)$$
The second one is given in Ref. [118, Eqs. (2.5), (2.6)]

$$ \int_s^{+\infty} dr \, q(r)^2 = q'(s)^2 - sq(s)^2 - q(s)^4 \quad (C.3) $$

### C.1.1 Left Asymptotics

Useful asymptotics are given [112] and

$$ q(s) = s \to -\infty \sqrt{-\frac{s}{2}} \left( 1 + \frac{1}{8s^3} - \frac{73}{128s^6} + \frac{10657}{1024s^9} + o\left(\frac{1}{s^9}\right) \right) \quad (C.4) $$

### C.1.2 Right Asymptotics

We also needed the following right asymptotics in the proof of Proposition 7.4, see [119, Sect. 1.3].

$$ \int_s^{+\infty} dr (r-s) q(r)^2 = s \to +\infty \frac{e^{-4/3s^{3/2}}}{15\pi s^{3/2}} \left( 1 - \frac{35}{24s^{3/2}} + \mathcal{O}\left(\frac{1}{s^3}\right) \right) \quad (C.5) $$

and

$$ \int_s^{+\infty} dr \, q(r) = s \to +\infty \frac{e^{-2/3s^{3/2}}}{2\sqrt{\pi}s^{3/4}} \left( 1 - \frac{41}{48s^{3/2}} + \mathcal{O}\left(\frac{1}{s^3}\right) \right) \quad (C.6) $$

## C.2 Relations Between the Tracy-Widom, Baik-Rains Distributions and the Hastings-McLeod Solution of the Painlevé II Equation

The Tracy-Widom distributions for $\beta = 2, 1, 4$ and the Baik-Rains distribution (denoted BR) are given by [109,118]

- For $\beta = 2$

$$ F_2(s) = \exp\left(-\int_s^{+\infty} dr (r-s) q^2(r)\right) = \text{Det}(I - \text{Ai}_s^2)_{\mathbb{L}^2(\mathbb{R}_+)} \quad (C.7) $$

- For $\beta = 1$

$$ F_1(s) = \exp\left(-\frac{1}{2} \int_s^{+\infty} dr [(r-s)q^2(r) + q(r)]\right) = \text{Det}(I - \text{Ai}_s)_{\mathbb{L}^2(\mathbb{R}_+)} \quad (C.8) $$

- For $\beta = 4$ (everywhere in the paper we use the conventions of [112])

$$ \begin{aligned} F_4(s) &= \exp\left(-\frac{1}{2} \int_s^{+\infty} dr (r-s) q^2(r)\right) \cosh\left(\frac{1}{2} \int_s^{+\infty} dr q(r)\right) \\ &= \frac{1}{2} \left( \text{Det}(I - \text{Ai}_s)_{\mathbb{L}^2(\mathbb{R}_+)} + \text{Det}(I + \text{Ai}_s)_{\mathbb{L}^2(\mathbb{R}_+)} \right) \end{aligned} \quad (C.9) $$

- For the Baik-Rains distribution

$$ \begin{aligned} F_{\text{BR}}(s) &= (1 + (s - 2q'(s) + 2q^2(s))(\int_s^{+\infty} dr q(r)^2)) \exp\left(-\int_s^{+\infty} dr [(r-s)q^2(r) + 2q(r)]\right) \\ &= \partial_s \left[ \int_{-\infty}^s dt \exp\left(-2 \int_t^{+\infty} dr q(r)\right) \exp\left(-\int_s^{+\infty} dr (r-s)q^2(r)\right) \right] \end{aligned} \quad (C.10) $$
It allows to obtain the following relation between the distributions:

$$F_4(s) = \frac{1}{2} \left( F_1(s) + \frac{F_2(s)}{F_1(s)} \right) \qquad (C.11)$$

and

$$F_{\text{BR}}(s) = \partial_s \left[ F_2(s) \int_{-\infty}^{s} dt \frac{F_1(t)^4}{F_2(t)^2} \right]. \qquad (C.12)$$

The next sections provide tail asymptotics for $F_1$ and $F_2$.

### C.2.1 Left Asymptotics of $F_1$ and $F_2$

The left asymptotics of the Tracy–Widom $\beta = 1, 2$ are given by (see e.g. [120] and references therein)

$$F_1(s) = 2^{-11/48} e^{\zeta'(-1)/2} \exp \left[ -\frac{|s|^3}{24} - \frac{|s|^{3/2}}{3\sqrt{2}} - \frac{\log|s|}{16} - \frac{|s|^{-3/2}}{24\sqrt{2}} + \frac{3|s|^{-3}}{128} - \frac{73|s|^{-9/2}}{1152\sqrt{2}} + \frac{63|s|^{-6}}{512} + \mathcal{O}(|s|^{-15/2}) \right] \qquad (C.13)$$

$$F_2(s) = 2^{1/24} e^{\zeta'(-1)} \exp \left[ -\frac{|s|^3}{12} - \frac{\log|s|}{8} + \frac{3|s|^{-3}}{64} + \frac{63|s|^{-6}}{256} + \mathcal{O}(|s|^{-9}) \right]$$

### C.2.2 Right Asymptotics of $F_1$

The right asymptotics of the Tracy–Widom $\beta = 1, 2$ are given by (see e.g. [119] and references therein)

$$1 - F_1(s) = \frac{e^{-\frac{2s^{3/2}}{3}}}{4\sqrt{\pi} s^{3/4}} \left[ 1 - \frac{41}{2^4 \cdot 3} s^{-3/2} + \frac{9241}{2^9 \cdot 3^2} s^{-3} - \frac{5075225}{2^{13} \cdot 3^4} s^{-9/2} + \frac{5153008945}{2^{19} \cdot 3^5} s^{-6} - \frac{1674966309205}{2^{23} \cdot 3^6} s^{-15/2} + \frac{3985569631633205}{2^{28} \cdot 3^8} s^{-9} + \mathcal{O}(s^{-21/2}) \right] \qquad (C.14)$$

and

$$1 - F_2(s) = \frac{e^{-\frac{4s^{3/2}}{3}}}{16\pi s^{3/2}} \left[ 1 - \frac{35}{2^3 \cdot 3} s^{-3/2} + \frac{3745}{2^7 \cdot 3^2} s^{-3} - \frac{805805}{2^{10} \cdot 3^4} s^{-9/2} + \frac{289554265}{2^{15} \cdot 3^5} s^{-6} - \frac{31241084875}{2^{18} \cdot 3^6} s^{-15/2} + \frac{23604769513325}{2^{22} \cdot 3^8} s^{-9} + \mathcal{O}(s^{-21/2}) \right] \qquad (C.15)$$

## D Asymptotics of $F(s)$

We compute the asymptotics of $F(s)$ for large positive and negative values of $s$ and plot in Fig. 6 the overlap between the complete PDF $F'(s)$ and the asymptotics obtained.
## D.1 Right Tail Using the Determinantal Formula

Let us perform the trace expansion for $s \to +\infty$ on the form

$$F(s) = \partial_s [2\mathrm{Det}(I + \mathrm{Ai}_s + |\mathrm{Ai}_s 1\rangle\langle 1|) + (s-2)\mathrm{Det}(I + \mathrm{Ai}_s)] . \quad (\mathrm{D}.1)$$

Here we can perform an expansion in powers of Airy functions that is in $e^{-\frac{2k}{3}s^{3/2}}$, $k = 1, 2, \dots$ using the first two orders of the trace expansion of the Fredholm determinant

$$\mathrm{Det}(I + M) = 1 + \mathrm{Tr}M + \frac{1}{2}((\mathrm{Tr}M)^2 - \mathrm{Tr}M^2) + \mathcal{O}(M^3) \quad (\mathrm{D}.2)$$

We have, where $\mathcal{O}(K^\ell)$ indicate the higher order traces in the expansion. Up to the second order in powers of Airy functions, we obtain

$$\begin{aligned} 2\mathrm{Det}(I + \mathrm{Ai}_s + |\mathrm{Ai}_s 1\rangle\langle 1|) + (s-2)\mathrm{Det}(I + \mathrm{Ai}_s) &= s + s\mathrm{Tr}\mathrm{Ai}_s + 2\mathrm{Tr}|\mathrm{Ai}_s 1\rangle\langle 1| + \mathcal{O}(\mathrm{Ai}^2) \\ &= s + s \int_0^{+\infty} dx \mathrm{Ai}(2x + s) + 2 \int_0^{+\infty} \int_0^{+\infty} dx dy \mathrm{Ai}(s + x + y) + \mathcal{O}(\mathrm{Ai}^2) \end{aligned} \quad (\mathrm{D}.3)$$

Differentiating allows to find back the CDF $F(s)$

$$\begin{aligned} F(s) &= 1 + \int_0^{+\infty} dx \mathrm{Ai}(2x + s) - \frac{s}{2}\mathrm{Ai}(s) - 2 \int_0^{+\infty} dx \mathrm{Ai}(x + s) + \mathcal{O}(\mathrm{Ai}^2) \\ &= 1 - \frac{s}{2}\mathrm{Ai}(s) - \frac{3}{2} \int_0^{+\infty} dx \mathrm{Ai}(x + s) + \mathcal{O}(\mathrm{Ai}^2) \end{aligned} \quad (\mathrm{D}.4)$$

## D.2 Right Tail Using the Asymptotics of $q(s)$ and the Tracy-Widom Distributions

We can now compare with the formula (7.39)

$$F(s) = e^{-\frac{1}{2} \int_{s}^{+\infty} dr [(r-s)q(r)^2 + 3q(r)]} \times \left( 1 + \frac{1}{2}(q'(s)^2 - sq(s)^2 - q(s)^4 - q(s))(s - 2q'(s) + 2q(s)^2) \right) \quad (\mathrm{D}.5)$$

which we recast into

$$F(s) = \frac{F_1(s)^3}{F_2(s)} \left( 1 + \frac{1}{2}(q'(s)^2 - sq(s)^2 - q(s)^4 - q(s))(s - 2q'(s) + 2q(s)^2) \right) \quad (\mathrm{D}.6)$$

Since $q(s)$ behaves asymptotically for large positive $s$ as the Airy function, we obtain at first order in $q(s)$ the expansion

$$F(s) = \frac{F_1(s)^3}{F_2(s)} \left(1 - \frac{s}{2}q(s)\right) + \mathcal{O}(q(s)^2) \quad (\mathrm{D}.7)$$

Using the asymptotics of the Tracy-Widom distributions (C.14) and (C.15) and the asymptotics of the Airy function, we obtain

$$1 - F(s) = \frac{s^{3/4} e^{-\frac{2s^{3/2}}{3}}}{4\sqrt{\pi}} \left[ 1 + \frac{139s^{-3/2}}{48} - \frac{11423s^{-3}}{4608} + \frac{3907027s^{-9/2}}{663552} - \frac{2886147455s^{-6}}{127401984} + o(s^{-6}) \right] \quad (\mathrm{D}.8)$$
<figure><img src="image_8.png" /><figcaption>Fig. 6 Overlap of the left and right tails of the PDF of the critical stationary case (derivative of Eqs. (D.10) and (D.4)) with the complete PDF (derivative of Eq. (7.37)). **Top.** True scale. **Bottom.** Logarithmic scale on the vertical axis</figcaption></figure>

### D.3 Left Tail Using the Asymptotics of $q(s)$ and the Tracy–Widom Distributions

Now we investigate the behaviour of $F(s)$ when $s \to -\infty$. Using once again that

$$F(s) = \frac{F_1(s)^3}{F_2(s)} \left( 1 + \frac{1}{2}(q'(s)^2 - sq(s)^2 - q(s)^4 - q(s))(s - 2q'(s) + 2q(s)^2) \right) \quad (D.9)$$

and reading the asymptotics of the Painlevé transcendent (C.4) and of the Tracy–Widom distributions (C.13) we obtain the left tail of $F(s)$ as

$$F(s) = 2^{-203/48} e^{\zeta'(-1)/2} \exp \left[ -\frac{|s|^3}{24} - \frac{|s|^{3/2}}{\sqrt{2}} + \frac{23}{16} \log |s| + \frac{91}{8\sqrt{2}|s|^{3/2}} - \frac{3957}{128|s|^3} + \frac{28717}{128\sqrt{2}|s|^{9/2}} - \frac{469683}{512|s|^6} + o(s^{-6}) \right] \quad (D.10)$$
E Extended Kernel and the Kadomtsev-Petviashvili Equation

Recently it was shown that the height CDF at the (large time) KPZ fixed point for the full space problem is related to scale-invariant solutions of the Kadomtsev-Petviashvili (KP) equation [121]. A related observation was made for the periodic KPZ fixed point [122]. This connection to the KP equation extends to the generating function at arbitrary time for the KPZ equation in full space, for some particular initial conditions, droplet, half-Brownian [121] and Brownian [123]. A similar relation was also obtained for a class of linear statistics associated to the Airy process [123]. However at present no such relation is known for the half-space problem.

Here we provide an extended version of our kernel which can be related to the KP equation. It involves however an additional variable which plays the role of a “fictitious space”. Although at this stage we could not see a physical interpretation for this variable, we believe this fact is curious enough to be reported.

We recall that for the half-space problem at all times, the generating function of the exponential of the KPZ height at the origin, $x = 0$, reads (7.3)

$$ \mathbb{E} \left[ \exp(-\zeta W e^{H(t)}) \right] = \sqrt{\mathrm{Det}(I - \bar{K}_{t,\zeta})_{\mathbb{L}^2(\mathbb{R}_+)}}, \quad (E.1) $$

Denoting $\zeta = e^{-r}$, one rewrites the kernel $\bar{K}$ as

$$ \bar{K}_{t, e^{-r}}(x, y) = 2\partial_x \iint_{C^2} \frac{dw dz}{(2i\pi)^2} G(z)G(w) \frac{\sin(\pi(z-w))}{\sin(\pi(z+w))} e^{-(x+r)z-(y+r)w+t \frac{w^3+z^3}{3}}, \quad (E.2) $$

where the function $G$ reads

$$ G(z) = \frac{\Gamma(A + \frac{1}{2} - z)}{\Gamma(A + \frac{1}{2} + z)} \frac{\Gamma(B + \frac{1}{2} - z)}{\Gamma(B + \frac{1}{2} + z)} \Gamma(2z). \quad (E.3) $$

The kernel $\bar{K}_{t, e^{-r}}$ verifies two simple identities

$$ \left\{ \begin{aligned} \partial_t \bar{K}_{t, e^{-r}}(x, y) &= -\frac{1}{3} [\partial_x^3 + \partial_y^3] \bar{K}_{t, e^{-r}}(x, y) \\ \partial_r \bar{K}_{t, e^{-r}}(x, y) &= [\partial_x + \partial_y] \bar{K}_{t, e^{-r}}(x, y) \end{aligned} \right. \quad (E.4) $$

We can now extend the kernel by introducing a fictitious variable $u$ in the following way

$$ \mathcal{K}_{t, e^{-r}, u}(x, y) = 2\partial_x \iint_{C^2} \frac{dw dz}{(2i\pi)^2} G(z)G(w) \frac{\sin(\pi(z-w))}{\sin(\pi(z+w))} e^{-(x+r)z-(y+r)w+t \frac{w^3-z^3}{3}} \quad (E.5) $$

so that

$$ \mathcal{K}_{t, e^{-r}, u=0}(x, y) = \bar{K}_{t, e^{-r}}(x, y) \quad (E.6) $$

The new kernel $\mathcal{K}_{t, e^{-r}, u}$ verifies the same set of differential equations (E.4) and in addition verifies a third one

$$ \partial_u \mathcal{K}_{t, e^{-r}, u}(x, y) = \frac{1}{2} [\partial_y^2 - \partial_x^2] \mathcal{K}_{t, e^{-r}, u}(x, y) \quad (E.7) $$

It was shown in [121], and known earlier in [124] (see discussion in [123, Appendix E]), that the three conditions (E.4) and (E.7) imply that the Fredholm determinant associated to $\mathcal{K}$ is a $\tau$-function of the KP equation. Thus, denoting $\mathcal{F}(r, t, u) = \sqrt{\mathrm{Det}(I - \mathcal{K}_{t, e^{-r}, u})_{\mathbb{L}^2(\mathbb{R}_+)}},$ the
function $\phi(r, t, u) := 2\partial_r^2 \log(\mathcal{F}(r, t, u))$ solves the KP equation for $(r, t, u)$

$$ \partial_t \phi + \phi \partial_r \phi + \frac{1}{12} \partial_r^3 \phi + \partial_r^{-1} \partial_u^2 \phi = 0 \qquad (E.8) $$

Note that the knowledge of $\phi(r, t, u = 0)$ is equivalent to the knowledge of the half-space generating function for the KPZ equation. For the full-space KPZ equation, the variable $u$ was interpreted as a spatial variable whereas in our case, it is a fictitious variable with no obvious direct interpretation.

## References

1. Kardar, M., Parisi, G., Zhang, Y.-C.: Dynamic scaling of growing interfaces. Phys. Rev. Lett. 56, 889 (1986)

2. Baik, J., Rains, E.M.: Symmetrized random permutations, and in Random matrix models and their applications, vol. 40 of Math. Sci. Res. Inst. Publ., pp. 1–19. Cambridge University Press (2001) arXiv:math/9910019

3. Prahofer, M., Spohn, H.: Current fluctuations for the totally asymmetric simple exclusion process, progress in Probability, vol. 51, edited by V. Sidoravicius (Birkhauser, Boston, 2002) 185, arXiv:cond-mat/0101200

4. Spohn, H.: Exact solutions for KPZ-type growth processes, random matrices, and equilibrium shapes of crystals. Physica A 369(1), 71 (2006)

5. Ferrari, P.L.: From interacting particle systems to random matrices. J. Stat. Mech. P10016 (2010)

6. Quastel, J.: Introduction to KPZ. J. Quastel (2012)

7. Corwin, I.: The Kardar-Parisi-Zhang equation and universality class. Random Matrices 1, 1130001 (2012)

8. Corwin, I.: Macdonald processes, quantum integrable systems and the Kardar-Parisi-Zhang universality class. In: Proceedings of the ICM. arXiv:1403.6877

9. Borodin, A., Gorin, V.: Lectures on integrable probability. Lecture notes (2012). arXiv:1212.3351

10. Borodin, A., Petrov, L.: Integrable probability: from representation theory to Macdonald processes. Probab. Surveys 11, 1 (2014)

11. Quastel, J., Spohn, H.: The one-dimensional KPZ equation and its universality class. J. Stat. Phys. 160, 965–984 (2015)

12. Takeuchi, K.A.: An appetizer to modern developments on the Kardar-Parisi-Zhang universality class. Physica A 504, 77–105 (2016)

13. Borodin, A., Petrov, L.: Lectures on integrable probability: stochastic vertex models and symmetric functions (2016). arXiv:1605.01349

14. Quastel, J., Matetski, K.: From the totally asymmetric simple exclusion process to the KPZ fixed point (2017). arXiv:1710.02635

15. Calabrese, P., Le Doussal, P., Rosso, A.: Free-energy distribution of the directed polymer at high temperature. Europhys. Lett. 90, 20002 (2010)

16. Dotsenko, V.: Bethe ansatz derivation of the Tracy-Widom distribution for one-dimensional directed polymers. EPL 90, 20003 (2010)

17. Dotsenko, V.: Replica Bethe ansatz derivation of the Tracy–Widom distribution of the free energy fluctuations in one-dimensional directed polymers. J. Stat. Mech. P07010 (2010)

18. Dotsenko, V., Klumov, B.: Bethe ansatz solution for one-dimensional directed polymers in random media. J. Stat. Mech. P03022 (2010)

19. Sasamoto, T., Spohn, H.: One-dimensional Kardar-Parisi-Zhang equation: an exact solution and its universality. Phys. Rev. Lett. 104, 230602 (2010)

20. Sasamoto, T., Spohn, H.: Exact height distributions for the KPZ equation with narrow wedge initial condition. Nucl. Phys. B 834, 523 (2010)

21. Sasamoto, T., Spohn, H.: The crossover regime for the weakly asymmetric simple exclusion process. J. Stat. Phys. 140, 209 (2010)

22. Amir, G., Corwin, I., Quastel, J.: Probability distribution of the free energy of the continuum directed random polymer in 1 + 1 dimensions. Commun. Pure Appl. Math. 64, 466 (2011)

23. Calabrese, P., Kormos, M., Le Doussal, P.: From the sine-Gordon field theory to the Kardar-Parisi-Zhang growth equation. EPL 107, 10011 (2014)
24. Calabrese, P., Le Doussal, P.: Exact solution for the Kardar-Parisi-Zhang equation with flat initial conditions. Phys. Rev. Lett. **106**, 250603 (2011)

25. Le Doussal, P., Calabrese, P.: The KPZ equation with flat initial condition and the directed polymer with one free end. J. Stat. Mech. P06001 (2012)

26. Le Doussal, P.: Crossover from droplet to flat initial conditions in the KPZ equation from the replica Bethe ansatz. J. Stat. Mech. P04018 (2014)

27. De Luca, A., Le Doussal, P.: Crossing probability for directed polymers in random media: exact tail of the distribution. Phys. Rev. E **93**, 032118 (2016)

28. Dotsenko, V.: Replica Bethe ansatz derivation of the GOE Tracy-Widom distribution in one-dimensional directed polymers with free boundary conditions. J. Stat. Mech. P11014 (2012)

29. Ortmann, J., Quastel, J., Remenik, D.: Exact formulas for random growth with half-flat initial data Ann. Appl. Probab. **26**, 507 (2016)

30. Imamura, T., Sasamoto, T.: Exact solution for the stationary Kardar-Parisi-Zhang equation. Phys. Rev. Lett. **108**, 190603 (2012)

31. Borodin, A., Corwin, I., Ferrari, P.L., Veto, B.: Height fluctuations for the stationary KPZ equation. Math. Phys. Anal. Geom. **18**, 20 (2015)

32. Corwin, I., Quastel, J., Remenik, D.: Renormalization fixed point of the KPZ universality class. J. Stat. Phys. **160**, 815 (2015)

33. Tracy, C.A., Widom, H.: Level-spacing distributions and the Airy kernel. Commun. Math. Phys. **159**(1), 151–174 (1994)

34. Tracy, C.A., Widom, H.: On orthogonal and symplectic matrix ensembles. Commun. Math. Phys. **177**(3), 727–754 (1996)

35. Takeuchi, K.A., Sano, M.: Universal fluctuations of growing interfaces: evidence in turbulent liquid crystals. Phys. Rev. Lett. **104**, 230601 (2010)

36. Takeuchi, K.A., Sano, M., Sasamoto, T., Spohn, H.: Growing interfaces uncover universal fluctuations behind scale invariance. Sci. Rep. **1**, 34 (2011)

37. Takeuchi, K.A., Sano, M.: Evidence for geometry-dependent universal fluctuations of the Kardar-Parisi-Zhang interfaces in liquid-crystal turbulence. J. Stat. Phys. **147**, 853–890 (2012)

38. Takeuchi, K.A.: Crossover from growing to stationary interfaces in the Kardar-Parisi-Zhang class. Phys. Rev. Lett. **110**, 210604 (2013)

39. Halpin-Healy, T., Takeuchi, K.A.: A KPZ cocktail-shaken, not stirred: toasting 30 years of kinetically roughened surfaces. J. Stat. Phys. **160**, 794 (2015)

40. De Nardis, J., Le Doussal, P., Takeuchi, K.A.: Memory and universality in interface growth. Phys. Rev. Lett. **118**, 125701 (2017)

41. Ito, I., Takeuchi, K.A.: When fast and slow interfaces grow together: connection to the half-space problem of the Kardar-Parisi-Zhang class. Phys. Rev. E **97**, 040103 (2018)

42. Kardar, M.: Depinning by quenched randomness. Phys. Rev. Lett. **55**, 2235 (1985)

43. Sasamoto, T., Imamura, T.: Fluctuations of a one-dimensional polynuclear growth model in a half space. J. Stat. Phys. **115**, 749 (2004)

44. Baik, J., Barraquand, G., Corwin, I., Suidan, T.: Pfaffian Schur processes and last passage percolation in a half-quadrant. Ann. Probab. **46**(6), 3015 (2018)

45. Baik, J., Barraquand, G., Corwin, I., Suidan, T.: Facilitated Exclusion Process, Computation and Combinatorics in Dynamics, Stochastics and Control Abelsymposium 2016, vol. 13. Springer, New York (2017)

46. Borodin, A., Rains, E.M.: Eynard-Mehta theorem, Schur process, and their Pfaffian analogs. J. Stat. Phys. **121**(3–4), 291–317 (2005)

47. Tracy, C., Widom, H.: The Bose gas and asymmetric simple exclusion process on the half-line. J. Stat. Phys. **150**, 1 (2013)

48. Tracy, C.A., Widom, H.: The asymmetric simple exclusion process with an open boundary. J. Math. Phys. **54**(10), 103301 (2013)

49. Barraquand, G., Borodin, A., Corwin, I., Wheeler, M.: Stochastic six-vertex model in a half-quadrant and half-line open ASEP. Duke Math. J. **167**(13), 2457–2529 (2018)

50. O’Connell, N., Seppäläinen, T., Zygouras, N.: Geometric RSK correspondence, Whittaker functions and symmetrized random polymers. Invent. math. **197**, 361–416 (2014)

51. Barraquand, G., Borodin, A., Corwin, I.: Half-space Macdonald processes. Forum Math. (2018). arXiv:1802.08210

52. Bisi, E., O’Connell, N., Zygouras, N.: The geometric Burge correspondence and the partition function of polymer replicas (2020). arXiv:2001.09145

53. Gueudré, T., Le Doussal, P.: Directed polymer near a hard wall and KPZ equation in the half-space. Europhys. Lett. **100**, 26006 (2012)
54. Krajenbrink, A., Le Doussal, P.: Large fluctuations of the KPZ equation in a half-space. SciPost Phys. **5**, 032 (2018)

55. Borodin, A., Bufetov, A., Corwin, I.: Directed random polymers via nested contour integrals. Ann. Phys. **368**, 191–247 (2016)

56. Krajenbrink, A., Le Doussal, P.: Replica Bethe Ansatz solution to the Kardar-Paris–Zhang equation on the half-line. SciPost Phys. **8**, 035 (2020)

57. de Nardis, J., Krajenbrink, A., Le Doussal, P., Thiery, T.: Delta-Bose gas on a half-line and the KPZ equation: boundary bound states and unbinding transitions (2019). arXiv:1911.06133

58. Baik, J., Rains, E.M.: Limiting distributions for a polynuclear growth model with external sources. J. Stat. Phys. **100**, 523–542 (2000)

59. Baik, J., Ferrari, P.L., Péché, S.: Limit process of stationary TASEP near the characteristic line. Commun. Pure Appl. Math. **63**, 1017–1070 (2010)

60. Aggarwal, A.: Current fluctuations of the stationary ASEP and six-vertex model. Duke Math. J. **167**, 269–384 (2018)

61. Seppäläinen, T.: Scaling for a one-dimensional directed polymer with boundary conditions. Ann. Probab. **40**(1), 19 (2012)

62. Betea, D., Ferrari, P.L., Occelli, A.: Stationary half-space last passage percolation. Commun. Math. Phys. **377**, 421–467 (2020)

63. Corwin, I., Shen, H.: Open ASEP in the weakly asymmetric regime. Commun. Pure Appl. Math. **71**(10), 2065 (2018)

64. Parekh, S.: The KPZ limit of ASEP with boundary. Commun. Math. Phys. **365**(2), 569–649 (2019)

65. Parekh, S.: Positive random walks and an identity for half-space SPDEs (2019). arXiv:1901.09449

66. Wu, X.: Intermediate disorder regime for half-space directed polymers (2018). arXiv:1804.09815

67. Liggett, T.M.: Ergodic theorems for the asymmetric simple exclusion process. Trans. Am. Math. Soc. **213**, 237 (1975)

68. Derrida, B., Evans, M.R., Hakim, V., Pasquier, V.: Exact solution of a 1D asymmetric exclusion model using a matrix formulation. J. Phys. A **26**(7), 1493 (1993)

69. Grosskinsky, S.: Phase transitions in nonequilibrium stochastic particle systems with local conservation laws PhD thesis, TU Munich (2004)

70. Derrida, B., Enaud, C., Lebowitz, J.: The asymmetric exclusion process and Brownian excursions. J. Stat. Phys. **115**(1–2), 365–382 (2004)

71. Bryce, W., Wang, Y.: Limit fluctuations for density of asymmetric simple exclusion processes with open boundaries. Ann. l'Inst. Henri Poincaré Probab. Stat. **55**, 2169–2194 (2019)

72. Tracy, C., Widom, H.: Matrix kernels for the Gaussian orthogonal and symplectic ensembles. Ann. l'inst. Fourier **55**(6), 2197–2207 (2005)

73. Forrester, PJ.: Painleve transcendent evaluation of the scaled distribution of the smallest eigenvalue in the Laguerre orthogonal and symplectic ensembles (2000). arXiv:nlin/0005064

74. Baik, J., Rains, E.M.: The asymptotics of monotone subsequences of involutions. Duke Math. J. **109**(2), 205 (2001)

75. Wang, D.: The largest sample eigenvalue distribution in the rank 1 quaternionic spiked model of Wishart ensemble. Ann. Probab. 1273–1328 (2009)

76. Bornemann, F.: On the numerical evaluation of distributions in random matrix theory: a review. Markov Process. Relat. Fields **16**(2010), 803–866 (2009)

77. Halpin-Healy, T., Lin, Y.: Universal aspects of curved, flat, and stationary-state Kardar-Paris–Zhang statistics. Phys. Rev. E **89**(1), 010103 (2014)

78. Borodin, A., Corwin, I.: Macdonald processes. Probab. Theor. Rel. Fields **158**(1–2), 225–400 (2014)

79. Borodin, A., Corwin, I., Ferrari, P.: Free energy fluctuations for directed polymers in random media in 1+1 dimension. Commun. Pure Appl. Math. **67**(7), 1129 (2014)

80. Borodin, A., Corwin, I., Petrov, L., Sasamoto, T.: Spectral theory for interacting particle systems solvable by coordinate Bethe ansatz. Commun. Math. Phys. **339**(3), 1167 (2015)

81. Borodin, A., Corwin, I., Petrov, L., Sasamoto, T.: Spectral theory for the q-Boson particle system. Compos. Math. **151**, 1 (2015)

82. Kardar, M.: Replica Bethe ansatz studies of two-dimensional interfaces with quenched random impurities. Nucl. Phys. B **290**, 582 (1987)

83. Lieb, E.H., Liniger, W.: Exact analysis of an interacting Bose gas. I. The general solution and the ground state. Phys. Rev. **130**, 1605 (1963)

84. Gaudin, M.: Boundary energy of a Bose gas in one dimension. Phys. Rev. A **4**, 386 (1971)

85. Oelkers, N., Batchelor, M.T., Bortz, M., Guan, X.W.: Bethe Ansatz study of one-dimensional Bose and Fermi gases with periodic and hard wall boundary conditions. J. Phys. A **39**, 1073 (2006)

86. Gaudin, M., Caux, J.-S.: The Bethe wavefunction. Cambridge University Press, Cambridge (2014)
87. Gueudré, T., Le Doussal, P.: Directed polymer near a hard wall and KPZ equation in the half-space. *Europhys. Lett.* **100**, 26006 (2012)

88. Castillo, I.P., Dupic, T.: Reunion probabilities of *N* one-dimensional random walkers with mixed boundary conditions. *J. Stat. Phys.* **3**(156), 606–616 (2014)

89. van Diejen, J.F., Emsiz, E.: Orthogonality of Bethe Ansatz eigenfunctions for the Laplacian on a hypercubical Weyl alcove. *Commun. Math. Phys.* **350**, no. 3, 1017 (2017)

90. van Diejen, J.F., Emsiz, E., Zurrian, I.N.: Completeness of the Bethe Ansatz for an open q-boson system with integrable boundary interactions. *Ann. Henri Poincaré*, **19**, 1349 (2018)

91. Emsiz, E.: Completeness of the Bethe ansatz on Weyl alcoves. *Lett. Math. Phys.* **91**, 61–70 (2010)

92. Gutkin, E., Sutherland, B.: Completely integrable systems and groups generated by reflections. *PNAS* **76**, 6057 (1979)

93. Heckman, G.J., Opdam, E.M.: Yang's system of particles and Hecke algebras. *Ann. Math.* **145**, 139–173 (1997)

94. P. Le Doussal (2016) *Unpublished notes*

95. Hao, Y., Zhang, Y., Liang, J.Q., Chen, S.: Ground-state properties of one-dimensional ultracold Bose gases in a hard-wall trap. *Phys. Rev. A* **73**, 063617 (2006)

96. McGuire, J.B.: Study of exactly soluble one-dimensional N-body problems. *J. Math. Phys.* **5**, 622 (1964)

97. Imamura, T., Sasamoto, T.: Stationary correlations for the 1D KPZ equation. *J. Stat. Phys.* **150**, 908–939 (2013)

98. Knuth, D.E.: Overlapping Pfaffians. *Electron. J. Comb.* (1995)

99. Alberts, T., Khanin, K., Quastel, J.: The intermediate disorder regime for directed polymers in dimension 1 + 1. *Ann. Probab.* **42**(3), 1212 (2014)

100. Ghosal, P.: Moments of the SHE under delta initial measure (2018). arXiv:1808.04353

101. O'Connell, N., Yor, M.: Brownian analogues of Burke's theorem. *Stoch. Process. Appl.* **96**(2), 285 (2001)

102. Baryshnikov, Y.: GUEs and queues. *Probab. Theory Rel. Fields* **119**(2), 256 (2001)

103. Gravner, J., Tracy, C.A., Widom, H.: Limit theorems for height fluctuations in a class of discrete space and time growth models. *J. Stat. Phys.* **102**(5–6), 1085 (2001)

104. Venkateswaran, V.: Symmetric and nonsymmetric Hall-Littlewood polynomials of type BC. *J. Algebr. Comb.* **42**, 331 (2015)

105. Ferrari, P.L., Spohn, H.: Scaling limit for the space-time covariance of the stationary totally asymmetric simple exclusion process. *Commun. Math. Phys.* **265**, 1–44 (2006)

106. De Bruijn, N.: On some multiple integrals involving determinants. *J. Indian Math. Soc.* **19**, 133–151 (1955)

107. Rains, E.M.: Correlation functions for symmetrized increasing subsequences (2000). arXiv:math/0006097

108. Ortmann, J., Quastel, J., Remenik, D.: A Pfaffian representation for flat ASEP. *Commun. Pure Appl. Math.* **70**(1), 3 (2015)

109. Ferrari, P.L., Spohn, H.: A determinantal formula for the GOE Tracy-Widom distribution. *J. Phys. A* **38**(33), L557 (2005)

110. Hastings, S.P., Mcleod, J.B.: A boundary value problem associated with the second Painlevé transcendent and the Korteweg-De Vries equation. *Arch. Ration. Mech. Anal.* **73**(1), 31–51 (1980)

111. Imamura, T., Sasamoto, T.: Fluctuations of the one-dimensional polynuclear growth model with external sources. *Nucl. Phys. B* **699**(3), 503–544 (2004)

112. Baik, J., Buckingham, R., DiFranco, J.: Asymptotics of Tracy-Widom distributions and the total integral of a Painlevé II function. *Commun. Math. Phys.* **280**(2), 463–497 (2008)

113. Bornemann, F.: On the numerical evaluation of Fredholm determinants. *Math. Comput.* **79**, 871 (2010)

114. Imamura, T., Sasamoto, T.: Replica approach to the KPZ equation with half Brownian motion initial condition. *J. Phys. A* **44**, 385001 (2011)

115. Tracy, C.A., Widom, H.: On ASEP with step Bernoulli initial condition. *J. Stat. Phys.* **137**(5–6), 825 (2009)

116. Tracy, C.A., Widom, H.: Integral formulas for the asymmetric simple exclusion process. *Commun. Math. Phys.* **279**(3), 815–844 (2008)

117. Macdonald, I.G.: Symmetric Functions and Hall Polynomials, vol. 354. Clarendon Press, Oxford (1995)

118. Baik, J., Rains, E.M.: Limiting distributions for a polynuclear growth model with external sources. *J. Stat. Phys.* **100**, 523–541 (2000)

119. Borot, G., Nadal, C.: Right tail asymptotic expansion of Tracy-Widom beta laws. *Random Matrices* **1**(03), 1250006 (2012)

120. Borot, G., Eynard, B., Majumdar, S.N., Nadal, C.: Large deviations of the maximal eigenvalue of random matrices. *J. Stat. Mech. P11024* (2011)
121. Quastel, J., Remenik, D.: KP governs random growth off a one dimensional substrate (2019). arXiv:1908.10353

122. Prolhac, S.: Riemann surfaces for KPZ with periodic boundaries. SciPost Phys. 8(1) (2020)

123. Le Doussal, P.: Large deviations for the KPZ equation from the KP equation (2019). arXiv:1910.03671

124. Poppe, C.: General determinants and the tau function for the Kadomtsev-Petviashvili hierarchy. Inverse Probl. 5(4), 613 (1989)

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.