Replica Bethe Ansatz solution to the Kardar-Parisi-Zhang equation on the half-line

We consider the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height h(x, t) on the positive half line with boundary condition ∂xh(x, t)|x=0 = A. It is equivalent to a continuum directed polymer (DP) in a random potential in half-space with a wall at x = 0 either repulsive A > 0, or attractive A < 0. We provide an exact solution, using replica Bethe ansatz methods, to two problems which were recently proved to be equivalent [Parekh, arXiv:1901.09449]: the droplet initial condition for arbitrary A > −1/2, and the Brownian initial condition with a drift for A = +∞ (infinite hard wall). We study the height at x = 0 and obtain (i) at all time the Laplace transform of the distribution of its exponential (ii) at infinite time, its exact probability distribution function (PDF). These are expressed in two equivalent forms, either as a Fredholm Pfaffian with a matrix valued kernel, or as a Fredholm determinant with a scalar kernel. For droplet initial conditions and A > − 2 the large time PDF is the GSE Tracy-Widom distribution. For A = 1 2 , the critical point at which the DP binds to the wall, we obtain the GOE Tracy-Widom distribution. In the critical region, A+ 1 2 = ǫt −1/3 → 0 with fixed ǫ = O(1), we obtain a transition kernel continuously depending on ǫ. Our work extends the results obtained previously for A = +∞, A = 0 and A = −12 .


Introduction
The continuum Kardar-Parisi-Zhang (KPZ) equation in one dimension [1][2][3][4][5][6][7] describes the stochastic growth of an interface parameterized by a height field h(x, t) at point x ∈ and time t. In appropriate units it reads where ξ(x, t) is a centered Gaussian white noise with ξ(x, t)ξ(x , t ) = δ(x − x )δ(t − t ). An important question is to determine the probability distribution function (PDF) for the height at one point, h(0, t), given an initial condition h(x, t = 0). Exact solutions, valid for all times t > 0, have been obtained for several initial conditions, notably droplet, flat and Brownian (including stationary) [8][9][10][11][12][13][14][15][16]. Remarkably they are expressed using Fredholm determinants or Pfaffians. From them it was shown that, in the large time limit, the PDF's of the shifted height fluctuations, H(t) = h(0, t) + t 12 ∼ t 1/3 , are described by the Tracy Widom distributions [17,18], i.e. the distributions of the largest eigenvalues of standard Gaussian random matrix ensembles [19,20]. Some of these results were obtained using the replica Bethe ansatz (RBA) method. This method, pioneered by Kardar [21], consists of several steps. First one maps, via the Cole-Hopf transformation Z(x, t) = e h(x,t) , the KPZ equation for h(x, t) to the stochastic heat equation (SHE) for Z(x, t), which can then be interpreted as the partition function of a continuum directed polymer of length t in a random potential with final point at x. Next one expresses the n-th integer moment of Z(x, t) using (i) a mapping to the attractive delta Bose gas with n bosons (the replica) in one dimension, and (ii) its exact solution from the Bethe ansatz. From the moments one then obtains the generating function (i.e. the Laplace transform of the PDF of Z(x, t)), and at large time, the PDF of h(x, t).
Here, we consider the KPZ equation on the half-line, where Eq. (1) is considered for x ∈ + along with the Neumann boundary condition (b.c.) where A is a real parameter which describes the interaction with the boundary (a wall at x = 0). This problem was considered in a pioneering paper by Kardar [22] in the equivalent representation in terms of a directed polymer in a half space bounded by a wall. The wall is repulsive for A > 0 and attractive for A < 0. The case A = +∞ imposes Z(x = 0, t) = 0, an infinitely repulsive wall. It can also be seen as absorbing wall, while A = 0 can be seen as a reflecting wall. A binding transition to the wall was predicted for A = −1/2 as A is decreased, from heuristic considerations on the ground state of the associated delta Bose gas in presence of a wall. It was later observed in numerical simulations of a discretized model [23]. More recently, exact results for three specific values of A, i.e. A = +∞, 0, −1/2, were obtained [24][25][26] for the droplet IC defined as What was studied is the height at the origin, H(t), i.e. H(t) = h(κ = 0, t) for finite A, and H(t) = h(κ, t) − log κ with κ → 0 + since a regularization is needed for A = +∞ (see below). In all three cases the solution can be expressed in terms of a Fredholm Pfaffian involving a matrix kernel [24][25][26][27]. For A = +∞, the infinitely repulsive wall, it was found [24] that the PDF of the scaled height, H(t)/t 1/3 , converges at large t to the Tracy-Widom distribution associated to the Gaussian Symplectic Ensemble (GSE) of random matrices [18][19][20]. Note that the original RBA solution of [24] involves a scalar kernel, but recently an equivalent solution was obtained by us in terms of a matrix kernel [27]. For A = 0, although the finite time matrix kernel differs from the case A = +∞, it was found that the large time limit of the PDF also corresponds to the TW-GSE distribution [25]. Both cases used the mapping to the delta Bose gas, with use of respectively the RBA for A = +∞ and nested contour integral representations of the moments for A = 0. The case A = − 1 2 , i.e. the critical case, was solved instead using a continuum limit from the ASEP model with an open boundary [26]. It was found that at large time the PDF converges to the Tracy Widom distribution associated to the Gaussian Orthogonal Ensemble (GOE) which thus describes the large time behavior at the transition.
Until now the KPZ equation for general A has resisted the analysis. The RBA approach of [24] has been extended using that the associated Bose gas is integrable for general A. This led to an explicit formula for the n-th integer moment of Z [28], which will be fully confirmed here, but valid however only in some restricted range of n and A, namely A + 1/2 > n/2. The origin of this restriction is that for any finite A, a complicated structure of multi particle bound states to the wall arises in the spectrum of the associated delta Bose gas. This problem did not impact the study in [24] since for A = +∞ there are no boundary bound states. Extracting the PDF from the moments in the restricted range has eluded previous attempts, but here we will provide the solution to this puzzle. Another approach was followed recently: the full structure of these bound states was obtained [29] leading to an improved moment formula a priori valid for all A, n. A (quite formal) Pfaffian formula was obtained from these moments. Among the partial results obtained, the large time limit as a GSE-TW distribution for all A > −1/2, and the convergence to a Gaussian PDF in the bound phase A < −1/2 [29]. This study however failed until now to produce the full finite time answer, and also to obtain the large time behavior in the critical region.
The half-space problem has also been addressed for other solvable models in the KPZ universality class, mostly in the mathematics literature. In a pioneering paper, Baik and Rains [30] studied the longest increasing sub-sequences (LIS) of symmetrized random permutations. The problem maps to a discrete zero temperature model of a directed polymer in a half-space, with a tunable parameter α which makes the boundary more attractive as α increases. They found, in the limit of large polymer length t, a transition when α reaches the critical value α c = 1. For α < α c the PDF of the fluctuations of the DP energy (analog to the height in a growth model) is given by the Tracy-Widom GSE distribution [18] on the characteristic KPZ scale t 1/3 . For α > α c the PDF is Gaussian on the scale t 1/2 , as the DP paths are bound to the diagonal line. At the critical point, α = α c the PDF is given by the GOE Tracy-Widom distribution on the t 1/3 scale. A similar transition was found for the height distribution in the discrete PNG growth model on a half-line, with a source at the origin, as the the nucleation rate at the origin is increased above a threshold [31,32]. Results were also obtained for the TASEP in a half-space [33]. Finally, formula were obtained for the (finite temperature) log-gamma DP with symmetric weights [34,41] and half-quadrant geometries [35] with, however, not yet rigorous asymptotic results. In the case of the ASEP model in a half space a Bethe ansatz study was carried in Ref. [36] without however asymptotic results. In the case A = −1/2, the ASEP and its KPZ limit were studied using half-space Macdonald processes [41] and TW GOE asymptotics were proved, see Refs. [26,58].
It is natural to conjecture that the transition for the KPZ equation at A = −1/2 is in the same universality class, in the large time limit, as the one discovered by Baik and Rains in [30] and that this universality is common to the full KPZ class, see Ref. [58]. Baik and Rains performed a detailed analysis on a scale α − 1 = wt −1/3 around the transition. They found that the PDF depends continously on w and interpolates between the GSE/GOE/Gaussian distributions as w is increased. This transition PDF was obtained as a solution of Painlevé type system of equations. Further results were obtained recently using Pfaffian-Schur processes, for variants of TASEP models and last passage percolation in a half-quadrant [37,38]. Not only the onepoint, but also the multi-point height distributions where studied (the extended process) and for arbitrary positions with respect to the wall. A Fredholm Pfaffian formula was obtained with an explicit expression for the (extended) matrix kernel around the GSE/GOE/Gaussian transition. One may conjecture that it is compatible with the Painlevé system of [30] but this has not yet been verified. Finally numerical studies have addressed the half-line problem. In [39] the convergence to the GSE was explored in a half-space geometry aimed to open the way for an experimental confirmation. In [23] the Baik Rains transitional PDF was verified and connections to conductance fluctuations in Anderson localization were explored.
The aim of this paper is to provide an exact solution for the KPZ equation on a half-line for generic values of A using the replica Bethe ansatz. We study the height at x = 0 and obtain (i) at all time the generating function (i.e. the Laplace transform of the PDF of Z = e h ) (ii) at infinite time, its exact probability distribution function (PDF). These are expressed in two equivalent forms, either as a Fredholm Pfaffian with a matrix valued kernel, or as a Fredholm determinant with a scalar kernel. For droplet initial conditions and A > − 1 2 the large time PDF is the GSE Tracy-Widom distribution. We obtain an explicit formula for the transition kernel in the critical region, A + 1 2 = εt −1/3 → 0 with fixed ε = O(1), continuously depending on ε. Although it is different from the expression obtained in [37,38], we conjecture that the two expressions are in fact equivalent, and check it by a systematic expansion in traces to third order. In the limit ε → 0 we obtain that the PDF converges to the GOE TW distribution. Our results thus show universality and that the KPZ equation in the half-line is in the KPZ class. Note however that for a technical reason, we are able to present the result only for the case A −1/2, i.e. ε 0.
Our strategy is the following. It is inspired by a theorem shown recently: Remark 1.2. This theorem comes as a limit of an identity proved in Ref. [41], Proposition 8.1, on half-space Macdonald processes.
We thus choose to study, rather than our original problem with parameter A, the KPZ problem on the half-line with Dirichlet boundary conditions, but with Brownian initial conditions. To be able to apply the theorem we choose the drift of the Brownian to be A + 1 2 . Since the boundary conditions are Dirichlet (i.e. A = +∞) we can apply the same RBA method as in [24] (there are no boundary bound states for A = +∞). The only technical difficulty is the calculation of the "overlap" of the Brownian initial condition with the eigenstates of the delta Bose gas, which we are able to perform. We then obtain a formula for the n-th integer moment of Z Br (x, 0), which then leads us to a formula for the moments of the droplet initial condition for generic A using the theorem. Although obtained via a completely different calculation the formula is, in the end, identical to the one obtained in [28] (which, in a sense, confirms the theorem). It is also valid only in the same restricted range A + 1 2 > n/2 (as the overlap diverges beyond). This is a well known difficulty, which arises already on the full line, but was surmounted in Refs [14,15] to obtain the solution of KPZ with Brownian initial conditions on the full line. We thus follow the same strategy, and it leads us to the abovementioned results.
Let us close by mentioning that these exact formulae for the KPZ equation at all times are very useful to obtain exact results for the large deviations of the PDF of the KPZ height field both at large time [46][47][48][49] and short time [50][51][52][53] and in particular in the half-space [27] with excellent agreeement with numerics [54,55]. These exact solutions have also been used in the mathematics community to prove exact bounds on the tails of the PDF of the KPZ height, see Refs. [42][43][44][45]. Although we will study here only typical fluctuations and not large deviations, the new formulae obtained in this work should also allow for such results for generic A in the future.

Presentation of the main results
In this paper we obtain in a single unified calculation the statistics of (i) Z(t) = e h(t) and H(t) = h(0, t) + t 12 where h(0, t) is the height at x = 0 with droplet initial condition in presence of a wall of parameter A ; (ii)Ẑ Br (t) = e h Br (t) = lim κ→0 + e h Br (κ,t) and H Br (t) = h Br (κ, t) + t 12 where h Br (κ, t) − log κ is the height at x = κ with Brownian initial condition with a drift A+ 1 2 in presence of a hard-wall. From the theorem of Parekh in [40] one has the equality of the two generating functions, for any ς > 0 KPZ,Brownian exp(−ςe H Br (t) ) = KPZ exp(−ςe H(t) ) .
The expected value of the left hand side is taken over the white noise of the KPZ equation and the Brownian initial condition while the expected value of the right hand side is taken over the white noise of the KPZ equation.

Finite time
Our main result valid for all time t 0 and all A > − 1 2 is that the generating function for ς > 0 can be written as a Fredholm Pfaffian where kernel K is matrix valued and represented by a 2 × 2 block matrix with elements In this formula the contours C are parallel to the imaginary axis and cross the real axis between 0 and A + 1 2 . For definition and more details about Fredholm Pfaffians see the end of Section 3.3. For a visual illustration of the block structure in the Pfaffians of matrix valued kernel, appearing in (6), see e.g. (146).
Remark 2.1. This formula reproduces the known cases. For A → +∞, the ratio of Γ functions containing A behaves as ∼ e −2(w+z) log A . A change of variable r =r − 2 log A then recovers the formula that we obtained in Eq. (5) and (6) of [27] with the correspondence H 1 = H(t) + 2 log A (up to a trivial rescaling involving time explained below Eq. (62) in [27]). There it was shown that this formula is equivalent to the result obtained in [24] in terms of a scalar kernel. Our result for A = 0, using some elementary manipulations 1 identifies with the one in [25]. The limit A = − 1 2 of our formula is more involved, and we discuss it below. 1 We first rewrite (7) using sin πz , Γ (2z) = Γ (z)Γ (z + 1 2 )2 2z−1 / π for z and for w. We can now compare with formula (5) and below in [25] (arXiv version). One finds K 22 (r, r ) =K 11 (−r, −r )/π (using the change (w 1 , w 2 ) → (−w, −z)) K 11 (r, r ) =K 22 (−r, −r )π (using the change (s 1 , s 2 ) → (w, z)) and The series in Eq. (6) can also be interpreted as a Fredholm Pfaffian, see Eq. (54), implying a duality between the generating function of the exponential KPZ height and an average of a "Fermi factor" over a Pfaffian point process (see Ref. [27] Eq. 7 for instance for further details on this type of duality) where the set {a i } forms a Pfaffian point process with kernel K.

Large time limit
At large time, we have obtained the following limit behavior of the solution of the Kardar-Parisi-Zhang equation: • For any A > − 1 2 , the one-point KPZ height fluctuations follow the Tracy-Widom GSE distribution • For A = − 1 2 , the one-point KPZ height fluctuations follow the Tracy-Widom GOE distribution.
is the cumulative distribution function (CDF) of the GSE-TW distribution, as defined in [56]. Another convention, which we denoteF 4 with F 4 (s) =F 4 ( s 2 ), is given in [18,57]. Near the transition point A = −1/2, there is critical regime for A + 1 2 → 0 and t → +∞ simultaneously, with the crossover parameter ε = (A+ 1 2 )t 1/3 being kept fixed and O(1). From the above formula (7) taking ς = e −st 1/3 one obtains the large time limit of the matrix kernel. This limit kernel, which depends continuously on ε, is called the transition matrix kernel. Its expression, K ε , is given in two equivalent forms in Eqs. (64) and in Eqs. (71) for ε > 0. In the limit ε → +∞ this transition kernel becomes equal to the standard matrix kernel of the GSE (the GSE limit is also obtained by taking the large time limit for any fixed A > −1/2 directly from (7)). The transition at A = −1/2 for the KPZ equation is believed to be in the same universality class than the one for last passage percolation in a half-quadrant. For the latter an explicit Fredholm Pfaffian formula was obtained in Ref. [37,38]. Although the kernel there is different from ours, we have checked by expansion in traces to third order that for ε > 0 their Fredholm Pfaffian coincides (see Appendix C).
In Ref. [27] we presented a general procedure to construct a scalar kernel from a matrix valued kernel with a Schur Pfaffian structure. Here we use this method to obtain the CDF of the one-point KPZ height in terms of a Fredholm determinant involving the following scalar kernel K 12 (r , r) =K 12 (−r, −r )π (using the change (s, w) → (w, −z)), whereK denotes the kernel displayed in formula (5) and below in [25]. One has also performed a shift r, r → r, r + 2 log 2, which together with the identification ζ/4 ≡ −ς shows that Z (0) in [25] is identical to e H here. The extra (−1) n s = (−1) L factor amounts to the permutation of the column and lines of the 2 × 2 block Pfaffian. Finally, t/2 there is t here.
where the transition kernel for ε > 0 is given bȳ Here, the contours C are parallel to the imaginary axis and cross the real axis between 0 and ε.
For ε → +∞, we check that this kernel converges towards the scalar kernel associated to the GSE, obtained in [24] for the droplet initial condition with A = +∞. The limit ε → 0 + is more delicate to handle, but via a careful analysis we are able to show that it converges towards the known scalar kernel of the GOE-TW distribution. We have introduced the notation F (ε) (s) for further purpose. Hence we show here that F (0) (s) = F 1 (s).
In Section 5.4 we extend our calculation and obtain the scalar kernel at any finite time. Finally, although we do not address the case A < −1/2 or ε < 0 the conjectured equivalence with the kernel of Ref. [37,38] (which holds for ε > 0) suggests that for any A < − 1 2 , the onepoint KPZ height has Gaussian fluctuations.
which is also the conclusion of the RBA analysis of [29], arising there from the contribution of the boundary bound states. Physically it is expected indeed in the phase where the DP is bound to the wall.

Bethe ansatz calculation for Brownian initial condition with Dirichlet boundary condition
Let us recall that via the Cole-Hopf mapping the partition sum Z(x, t) = e h(x,t) , where h(x, t) is solution of the KPZ equation (1), satisfies the stochastic heat equation (SHE) understood here with the Ito prescription. Here we consider the problem on the half-line x 0. Let us denote Z Br (x, t) the solution of the SHE (14) with the Brownian initial condition in presence of a drift A + 1 where B(x) is the standard Brownian (i.e. with B(0) = 0), and with Dirichlet boundary condition Z Br (x = 0, t) = 0. We will first calculate the integer moments of Z Br (x, t), from which we will obtain the moments of the limitẐ Br (t) = lim κ→0 + Z Br (κ, t)/κ. Then we will use that KPZ,B Ẑ Br (t) n = KPZ Z(t) n (16) to obtain the moments of Z(t) = Z(x = 0, t), which denotes here the solution of the SHE for the droplet initial condition and a wall of parameter A, our problem of main interest. Note that on the left hand side the average is over the Brownian B(x) and the KPZ noise while on the r.h.s. it is over the KPZ noise only.
The general equal time moments of the solution of the SHE, Z Br (x, t), over the KPZ noise can be expressed [21] as quantum mechanical expectation values of the evolution operator in imaginary time of the Lieb Liniger model [59] Here H n is the Hamiltonian of the Lieb Liniger model [59] for n quantum particles with attractive delta function interactions of strength c = −c < 0 with here an below, in our unitsc = 1. The initial state is such that 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 |Ψ(t = 0)〉 as A simple calculation shows that Φ 0 (x 1 , . . . , x n ) is the fully symmetric function which in the sector 0 x 1 < · · · x n takes the form

Bethe ansatz formula for the moments
We can now rewrite (17) at coinciding points using the decomposition of the evolution operator e −t H n in terms of the eigenstates of H n as Here the un-normalized eigenfunctions of H n are denoted Ψ µ (of norm denoted ||µ||) with eigenenergies E µ . 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 µ denotes a sum over all bosonic eigenstates of the Lieb-Liniger model, also called delta Bose gas, and 〈Ψ µ |Φ 0 〉 denotes the overlap, i.e. the hermitian scalar product of the initial state (21) with the eigenstate Ψ µ .
We should remember now that H n is defined on the half-line x 0. The boundary condition (2) 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 = +∞, i.e. for Dirichlet boundary condition as needed here [36,60,61] (see also section 5.1 of [62]). Note that it can also be solved for arbitrary A, [25,[62][63][64][65][66][67][68] which led to the formula in [28] and [29], but here we circumvent this, using instead the Brownian with Dirichlet boundary condition.
From the Bethe ansatz the eigenstates Ψ µ are thus Bethe states, i.e. superpositions of plane waves over all permutations P of the n rapidities λ j for j ∈ [1, n] with an additional summation over opposite pairs ±λ j due to the infinite hard wall. The bosonic (fully symmetric) eigenstates can be obtained everywhere from their expression in the sector x 1 < · · · < x n , which reads This wavefunction automatically satisfies both 1. The matching condition arising from the δ(

The hardwall boundary condition
The allowed values for the rapidities λ 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 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 {λ j }, depend on the choice of boundary condition at x = L. However, in the L → +∞ limit, these details do not matter. For simplicity we choose another hardwall at x = L. The Bethe equations then read Remark 3.1. This set of equations is invariant by λ j → −λ j for any j.
In the case of the infinite hardwall, these equations are also given in Ref. [61] and their solutions in the large L limit were studied in Ref. [69]. 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 1 particles with n = n s j=1 m j . Each bound state is a perfect string [70] , i.e. a set of rapidities where a = 1, . . . , m j labels the rapidities within the string. Such eigenstates have momentum and energy 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 λ j → −λ j , i.e. k j → −k j . From now on, we will denote the wavefunctions of the string states as Ψ {k ,m } . The inverse of the squared norm of an arbitrary string state was obtained in Ref.
[24] as 2 with S k,1 = 1. Note that we have only kept the leading term in L as L → +∞. Inserting the norm formula Eq. (28) into Eq. (22), we obtain the starting formula for the integer moments of the partition sum with Brownian weight on the endpoint in the limit L → +∞ Here the Kronecker delta enforces the constraint n s j=1 m j = n with m j 1 and in the summation over states we used k j → m j L dk 2π 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. [9,12,13,24,29]. Since we are interested in we can simplify the factor Ψ {k ,m } (x, · · · , x) in (29). For the general Bethe state (23) (before insertion of the string solution), the small x limit reads Inserting the string solution we see that we can replace in (29) at small x To obtain the n-th moment of (30) from (29) we still need to calculate the overlap 〈Ψ {k ,m } |Φ 0 〉 where Φ 0 is given in (21). In general it involves sums over permutations and leads to complicated expressions unless there is some kind of "miracle", known to happen in full space only for a few special initial conditions (droplet, half-flat, Brownian). Here, as we find in the Appendix A, the final result in the half-space for Brownian initial conditions is quite simple Inserting the rapidities λ j of the string state one see that the denominator of this formula reads while the numerator was already calculated in (32). We can thus define C k, j = A 2 k, j E k, j and putting all together we obtain the starting expression for the integer moments: where the expressions for C, S, D are given above and where we recall the constraint n s j=1 m j = n. Let us usec = 1 from now on. Denoting .
The starting formula for the moments of our two equivalent problems (droplet initial conditions with any A, and Brownian initial condition with drift A + 1/2 and Dirichlet) reads where B k,m is given in (37) and D k i ,m i ,k j ,m j is given in (28) and where we recall the constraint n s j=1 m j = n. This formula is identical to the one obtained by a completely different calculation in [28] for the problem of the droplet initial condition with any A, consistent with the theorem of Ref. [40]. For A → +∞ one recovers the expression (11) in [24], i.e the same without the ratio of Gamma functions involving A, using that Z there is lim A→+∞ A 2 Z(t) here.
It is important to note that, strictly, the formula (34) for the overlap of the string states with the Brownian initial condition is valid only when the multiple integrals in the scalar product are convergent. This requires (see discussion in the Appendix A) the condition n/2 < A + 1 2 . Hence the above formula for the integer moments KPZ,B Ẑ Br (t) n is valid only when the drift is large enough, for each value of n. This requirement for the drift is identical to the one obtained in the full space, for the Brownian initial condition in [14,15] and for the half-Brownian initial condition in [71]. There, it was shown how to nevertheless use these restricted moment formula to construct the full solution for the generating function of the moments at finite time (see definition below). Here we follow the same strategy as in these works and we refer the reader to Ref. [15] Section 4.3 for further details. The main idea behind the continuation of the moment formula is that it is possible to avoid all divergences in the complex plane by finely tuning the contour integrals of the Mellin-Barnes summation formula, which is possible as long as A + 1/2 0.
The resulting solution obtained in [14,15,71] for the full space holds for any positive drifts. The limit of zero drift is quite delicate, but can be performed, and leads to the Baik Rains distribution for stationary KPZ. The analog in the half-space is that we obtain below a formula valid for positive drift, A > −1/2. With some efforts, we are able to extend it to A = −1/2 the critical case, leading to the GOE-TW distribution. We have not attempted to obtain the solution for A < −1/2, i.e. negative drift.
Since the two problems are in correspondence, due to the theorem of Parekh in [40], the above formula for the moments KPZ Z(t) n for the droplet initial condition and wall parameter A is also valid only for n/2 < A + 1 2 . In that case, this restriction arises from an a priori different origin. As shown in [29], the string states discussed above do not form a complete basis for generic A: new boundary bound states arise. To obtain a formula for the moments valid for any n, A is possible but requires to include all boundary states, as done in [29]. This is illustrated here in the Appendix B for the first moment of Z(t). We will not follow that route here, although connections between the two approaches must exist.

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 ∈ [1, n s ], the norm reads where we recall that the Pfaffian of an anti-symmetric matrix A of size N × N is defined by and that the Schur Pfaffian is given by (see Ref. [78]) Hence the starting formula for the moments now becomes:

Generating function in terms of a Fredholm Pfaffian
We will now write the generating function for the moments of Z(t), i.e. focusing on the droplet initial condition with generic A (the result being identical for the Brownian with Dirichlet boundary conditions). It is defined, for ς > 0, as The constraint n s i=1 m i = n in Eq. (42) can then be relaxed by reorganizing the series according to the number of strings: where Z(n s , ς) is the partition sum at fixed number of strings, 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, Eq. (39), given above. The partition sum at fixed number of strings, expressed in terms of the reduced variables where B k,m was given in (37). The summation over the variables m p can be done using the Mellin-Barnes summation trick similarly to Refs. [14,15]. The barrier A > (n−1)/2 is overcome exactly as in Ref.
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 → w p imposed by the Mellin-Barnes formula. This leads to the following rewriting of the partition sum at fixed number of strings (see section 5 in [27] for similar manipulations) Remark 3.2. Note that the contour C 0 passes to the left of the pole of the Gamma function at X = 2A + 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 Ref. [27] section 5, we can proceed to the replacement The last manipulations consist in rescaling all variables X by a factor 2 and replacing the contours of integration by C = a 2 + i . Hence we have There are a few last steps before we introduce the Fredholm Pfaffian. First define the functions Using a known property of Pfaffians (see De Bruijn [72]), we can rewrite the partition sum at fixed number of strings itself as a Pfaffian, i.e. we use that This leads to the definition of a 2n s × 2n s matrix M such that 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 × 2 blocks which we denote K, whose elements are presented in Eqs. (7) (for visualization see formula (146)). Finally, the string-replicated partition function is given by an infinite series of Pfaffians This series admits a closed form in terms of a Fredholm Pfaffian, which is our main result for the generating function at finite time together with the explicit expression of the kernel given in Eqs.
(56) For the rest of the paper we will call this the α-decomposition. To study the limit A → +∞, we will use α = 0 and for the limit A → − 1 2 , we will use α = A + 1 2 . In particular, from the same symmetrization argument than above, the partition function for n s strings reads, using the α-decomposition: This leads to an α-extension of our functions φ, which in the special case α = A + 1 2 read We see that this choice α = A+ 1 2 can be interesting since it suppresses the pole at X = A+ 1 2 in the function φ 2p , a property that we will use below. At the end we obtain a the same Fredholm Pfaffian expression for the generating function, with the minimal replacement in the kernel (7) sin(πX ) → sin(π(X − (A + 1 2 ))), We emphasize that all the kernels parametrized by α yield the same Fredholm Pfaffian by construction.

Large time limit of the Fredholm Pfaffian and the distribution of the KPZ height
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 At large time, we want to eliminate the time factor in the exponential, hence we perform the change of variables In the large time limit, the Gamma, cosine and sine functions simplify using that for small positive argument Under these simplifications, the partition sum at fixed number of strings reads in the limit t → +∞ (dropping all tildes) The contours C have now to be understood as C =ã + i , whereã ∈]0, ε[. We emphasize that the contours all lie at the left of the poles at X = ε. Now that the rescaling at large time is well understood, we can reconstruct the matrix valued kernel K ε which reads in this limit Remark 4.1. The kernel K ε has a particular structure, indeed its elements are related through derivative identities: K ε 22 (r, r ) = ∂ r ∂ r K ε 11 (r, r ), K ε 12 (r, r ) = −∂ r K ε 11 (r, r ) and K ε 22 (r, r ) = −∂ r K ε 12 (r, r ).
Remark 4.2. The kernel K ε can be obtained equivalently from the kernel (7) by the same rescaling of r p = t 1/3r p and the changes (w, z) → t −1/3 (w, z) in the integrals.
Finally, choosing the variable ς as ς = e −st 1/3 , at large time we have where Θ is the Theta Heaviside function. The Fredholm Pfaffian formula for the generating function then becomes in the limit 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, i.e.

KPZ exp(−ςe H(t) ) = KPZ exp(−e H(t)−st
where Θ is the Theta Heaviside function. From this, we obtain the CDF of the height distribution in the large time limit as a Fredholm Pfaffian in terms of the matrix kernel given in (64), also called the transition kernel, as it describes the critical region around the transition. Anticipating a bit, it should describe the crossover from GSE/GOE/Gaussian discussed in the introduction, expected from universality arguments. Although we will indeed show most of it, our formula for the kernel is, however, limited at this stage to ε > 0.

Large ε behavior of the matrix valued kernel and convergence to the GSE-TW distribution
Performing the large time limit at any fixed A > −1/2 corresponds to the previous calculation with ε = +∞. Hence we now investigate the large ε behavior of the kernel at large time. As the contours of kernel K ε are parallel to the imaginary axis and cross the real axis between 0 and ε, we can push the limit ε → ∞ without any ambiguity. All rational functions involving the parameter ε in the large time limit of the kernel K ε in Eq. (64) converge to the value −1.
Hence in this limit we obtain a kernel which we denote K ∞ , and it reads This is precisely the kernel associated to the Gaussian Symplectic Ensemble (GSE) of random matrices as given in Lemma 2.7. of Ref. [37]. Hence this shows that the distribution of the height at x = 0 converges at large time for boundary conditions such that ε → ∞ (e.g. for any fixed A > −1/2) to the GSE Tracy-Widom distribution

Alternative expression for the large-time transition matrix-valued kernel
In this section we will present an equivalent form of the large time transition kernel for the KPZ equation, better suited to study general ε. In particular, we will compare it with the transition kernel obtained in in Ref. [37] from the solution of discrete models, i.e. last passage percolation and facilitated TASEP.
To study the large time behavior for general ε, i.e. the critical region, it is useful to return to the α-decomposition of Section 3.4. Let us choose, as discussed there, α = A + 1 2 . Inserting the replacement (59) into the finite time kernel (7) and repeating the same steps one obtains a large time kernel equivalent to (64) (which, for simplicity, we will also denote by K ε ). In that limit the replacement only amounts to multiply by (w − ε)(z − ε)/(wz) the integrand in K ε 22 , and by (z − ε)/z the integrand in K ε 12 , leading to As advertized above, this choice of α = A + 1 2 has decreased the number of poles for w or z equal to ε, which will be useful below.
We now slightly rewrite this kernel to make more apparent its asymptotic behavior for large arguments r, r → +∞. For this purpose it is useful to move the contours of integration. The contour C crosses the real axis between 0 and ε. We now choose to move all the contours to the right of ε. This will pick some residues that we will evaluate. Doing so, and introducing a contourĈ parallel to the imaginary axis and crossing the real axis at the right of ε, we find the decomposition This form is useful to evaluate the transitional CDF F (ε) (s) at large s, i.e. the Fredholm Pfaffian in a trace expansion at large s. The Fredholm Pfaffian involves integrals of the variables r, r on [s, +∞[. Its behavior for large s is thus controled by the region of large r, r > s. In that region, each term in (72) can be evaluated by saddle point: indeed the new contourĈ can be deformed to include the saddle point. Hence the asymptotics is of the standard Airy type, i.e. each cubic exponential in the integrand leads to a decay ∼ exp(− 2 3 s 3/2 ) and the total decay of each term is thus of the form ∼ exp(− 2 3 ps 3/2 ), where the integer p = 1, 2 is the number of cubic exponentials (akin to the number of Airy functions). To be systematic one thus labels the terms in (72) in terms of their homogeneity in cubic exponentials (not counting the terms in ε 3 ). We introduce the notation K (n) when the element K comprises n cubic exponentials and hence we write the kernel as 12 . (73) is a rank 1 operator.
An important property of our result is that our transition kernel K ε in (72) is equivalent to the cross-over kernel 3 K cross in Ref.
[37] Definition 2.9, i.e. we conjecture that for ε > 0 Let us first recall that from Ref.
[37] Definition 2.9, the cross-over kernel K cross is the sum of two kernels I and R, the latter having a single non-zero component R 22 . In our notations and contour conventions, these kernels read where we recall that the contour C is parallel to the imaginary axis and crosses the real axis between 0 and ε and that the contourĈ is parallel to the imaginary axis and crosses the real axis at the right of ε.
Remark 4.4. With these definitions the above formula holds for ε > 0. The kernel of Ref. [37] is also valid for ε < 0, with the following changes. For ε < 0, the contours of I 11 are unchanged, for I 12 the contour of z crosses the real axis at the left of ε, the contour of w crosses the real axis at the right of 0 such that w + z > 0. For I 22 both contours of w and z cross the real axis at the right of 0. For R 22 the contours cross the real axis between the left of −ε and the right of ε.
It is important to note that for ε > 0, the first two terms of R 22 can be absorbed in I 22 by changing the contourĈ → C in I 22 . Indeed the first two terms of R 22 are nothing more than the residues of the poles of I 22 at w = ε and z = ε. Hence, in later uses, we will redefine for convenience R 22 to be only its last term and include the first two terms in I 22 . Since only the sum I 22 + R 22 matters this is immaterial.
Since it is already proved in Ref. [37] that in the large positive ε → +∞ limit, their kernel is the GSE kernel, we see that in this limit it is equivalent to ours. For general ε > 0 we have expanded the Fredholm Pfaffians on both sides of Eq. (74) in series of their traces, using the counting in cubic exponentials explained above. It is equivalent to expanding both predictions for the CDF for the height in the transition regimes at large s. We found perfect matching up to, and including, third order, i.e. exp(−2s 3/2 ). Thus we have a strong case for the conjecture of the full equivalence of the Fredholm Pfaffians. A complete proof of this identity seems for the moment out of reach.

Small ε behavior of the matrix-valued kernel
In this Section, we will take the small ε limit of the two-dimensional kernel. Since the kernel K cross converges to the Tracy-Widom GOE kernel as ε → 0, see Ref. [37], the limit of K ε for ε → 0 will provide a new kernel for the Tracy-Widom GOE distribution in accordance to our conjecture in Eq. (74). The starting point is the expression of K ε in Eq. (72) and we will show that the ε → 0 limit can be taken without any ambiguity within the different elements of the kernel.
Whenever there is an Airy term of the type e −rz−rε−z 3 /3 , the ε → 0 limit is clear and we can set ε = 0 directly. In the case of the different projectors, where there is solely a term of the type e −rε in Eq. (72), the limit has to be investigated properly because of possible convergence issue when calculating the Fredholm Pfaffian which involves products of the elements of the matrix-valued kernel, see Eq. (146). Hence, with this remark and expressing the projectors in terms of Airy functions, we rewrite the kernel K ε as If one defines the following vectors then the kernel K ε→0 admits the reduced representation where K ∞ is the GSE kernel of Eq. (69). In order to treat the ε dependence in the term |u ε 〉 〈v ε | − |v ε 〉 〈u ε |, we treat it as a rank one perturbation to a Fredholm Pfaffian.
By a generalization of the matrix determinant lemma and the Sherman Morrison formula for Fredholm determinants 4 , we have The remaining goal is to show that the ε → 0 limit can be taken in the inner product in Eq. (81). It is possible to expand the resolvant (I + J K ∞ ) −1 , we will show on the lowest order term of the expansion that the ε → 0 limit can be taken easily and we will conjecture the same conclusion for higher order terms.
From the first to the second line, we proceeded to the change of variable (r, λ) → (u = r + λ ∈ + , v = r − λ ∈ [−u, u]) and relabeled u → r.
Our conclusion is that we can simply take the ε → 0 limit from the beginning in the vectors |u ε 〉 and |v ε 〉 so that the final resulting infinite-time kernel for A = −1/2 reads Remark 4.5. Note that in the kernel K 0 , the derivative relation between its elements is conserved: K 0 22 (r, r ) = ∂ r ∂ r K 0 11 (r, r ), K 0 12 (r, r ) = −∂ r K 0 11 (r, r ) and K 0 22 (r, r ) = −∂ r K 0 12 (r, r ).
According to Ref. [37] and our conjecture in Eq. (74), the Fredholm Pfaffian of the kernel K 0 describes the Tracy-Widom GOE distribution. To the best of our knowledge, this kernel K 0 has not appeared before in the literature and therefore provides an alternative expression for the Tracy-Widom GOE distribution. Since it enjoys the same simple derivative structure as the GSE kernel, see Remark 4.5 above, without involving any sign function, contrary to K cross , we hope this kernel will find further applications.
Having understood the ε → 0 limit on the infinite time matrix-valued kernel, we extend the calculation to study the critical case A = −1/2 at finite time. To this aim, we first rewrite the kernel Eq. (7) using the α-prescription with α = A + 1 2 .
We recall that the contours C are parallel to the imaginary axis and cross the real axis between 0 and A + 1 2 . Let us now push this contour to the right of A + 1 2 by picking the related pole and call the new contourĈ. We will directly write the result at A = −1/2 without further details to conclude this Section on the matrix-valued kernel.
This final kernel is the one that describes the critical case A = −1/2 at all times.
It was shown in Ref. [27] that Fredholm Pfaffians with matrix-valued kernel of Schur type can also be expressed as Fredholm determinants with a scalar kernel. We now dedicate the rest of this work to the understanding of the structure of the scalar-valued kernel.

Solution for the KPZ generating function at all times for generic A in terms of a scalar kernel
The general kernel we have obtained in Eq. (7) has a particular structure in the form of a Schur Pfaffian. With this structure, the kernel verifies the hypothesis of Proposition B.2 of [27]. This proposition states that we can transform the Fredholm Pfaffian of Eq. (54) which involves a matrix valued kernel, into a Fredholm determinant of a scalar kernel. To proceed, let us first define the functions and the kernelK t,ς such that for all (x, y) ∈ 2 We recall that in this formula the contours C are parallel to the imaginary axis and cross the real axis between 0 and A + 1 2 .
Then, the Laplace transform of the one-point distribution of the exponential of the KPZ height admits the following representation: Remark 5.1. In the case of the scalar kernel, there is no α-decomposition as any α prescription in the functions f odd and f even will eventually lead to the same kernel (87).
Large time limit. To obtain the large time limit of (88) one performs the same rescaling as in Sec. 4, namely one chooses ς = e −t 1/3 s and one rescales (w, z) → t −1/3 (w, z), r → t 1/3 r. The Fermi factor σ ς produces a Heaviside function Θ(r − s) and then (88) becomes (11), in terms of the large time scalar kernelK given in Eq. (12).

Large ε behavior of the scalar kernel and convergence to the GSE
We recall that at large time the critical cross-over parameter is given by ε = (A + 1 2 )t 1/3 . For large ε, we use the exact same procedure as in Section 4.1 and we obtain from Eq. (12) the kernelK which was proved to be the one-dimensional GSE kernel in Refs. [24,27]. Performing the integrals one can equivalently write it as where K Ai is the standard Airy kernel. This is the form in which it was obtained in [24] for the case A = +∞. Here we see that the GSE-TW holds for any fixed A > −1/2 at large time, consistent with the conclusion obtained above from the matrix valued kernel.

Small ε limit of the scalar kernel and convergence to the GOE
To investigate the small ε behavior of the scalar kernel, one proceeds as in the matrix case by moving the contours C of the scalar kernel of Eq. (12) toĈ which is to the right of ε, and collecting all residues along the way. Doing so, we rewrite the scalar kernel as The new contourĈ is now well suited for the small ε limit in the first term. However one cannot simply set ε = 0 in the last two terms. Indeed the operator product of the last term with the second one leads to a divergent integral. We now study this delicate limit in details.
To further study the small ε limit, one first conjugates the kernel on the left by the multiplication of e −εx and on the right by the multiplication of e ε y . This operation leaves the associated Fredholm determinant unchanged. The new kernel obtained by this manipulation is the following We split this kernel into three components and we will analyze them separately. (93) • The small ε limit of K 0 is taken easily as ε does not play any particular role in the convergence of the integrals, it reads • Similarly, the small ε limit of K 2 is taken straightforwardly and reads lim ε→0 Indeed, the Airy function will always make any integral over x convergent.
• The convergence of K 1 is more complicated and as the devil hides in the details, we have to analyze this part of the kernel very carefully. One could be tempted to write that This is wrong and the whole subtlety of the small ε limit. The convergence of the kernel with respect to the variable y is controlled by the Airy function, hence we are allowed to take the limit e ε y+ ε 3 3 → 1. We are now going to show the distributional identity When we multiply the kernel K 1 on the left by an arbitrary element K, we have to calculate the following product Proceeding to the change of variablex = 2εx, we have In the limit of small ε, this integral should be equal to Naively one would claim that this limit is always zero due to some exponential decay of the kernels. However, the kernel K 2 does not decay lim X →+∞ hence products of K 2 and K 1 give a non zero result for ε → 0.
To summarize this discussion, we have shown that the kernel for ε = 0 is equal to The last effort of this Section now consists in proving that the Fredholm determinant of K is indeed identical to the one of the GOE-TW distribution. As it is common for such manipulations, we introduce the brackets notations where we have implicitly promoted the Airy function to an operator of kernel Ai(x, y) = Ai(x+ y).
Remark 5.2. With this notation, we have K GLD = K Ai − 1 2 |Ai〉 〈1| Ai. Under these conventions, our kernel for ε = 0 is written as The kernel K can be interpreted as the K GLD kernel with a rank 2 perturbation. By the matrix determinant lemma coupled to the Sherman-Morrison formula, we can evaluate the Fredholm determinant of K and extract the two rank 1 perturbations as Since K GLD decreases fast enough at infinity in both arguments, we have for any integer n 1 the equality (K GLD ) n |δ ∞ 〉 = 0.
On the contrary, the identity keeps |δ ∞ 〉 invariant, i.e. I |δ ∞ 〉 = |δ ∞ 〉 and hence by a regular series expansion we obtain the action of δ ∞ on the resolvant of K GLD .
Similarly to the above considerations, by the decay property of the Airy function, we have that 〈1| Ai|δ ∞ 〉 = 0 and we also have readily the inner product 〈1|δ ∞ 〉 = 1. Implementing these results in the Fredholm determinant of Eq. (105) leads to In the first line, we simplified Eq. (105) using the action of δ ∞ , from the first to the second line we used the Matrix determinant lemma and from the second to the third line we introduced the expression of K GLD with brackets notations. The last Fredholm determinant was shown by Forrester in Ref. [77] to be related to the GOE Tracy-Widom cumulative distribution function Hence, we finally conclude that for ε → 0 + , the cumulative probability of the one-point KPZ height field is given by the GOE Tracy-Widom distribution.

Solution for finite time at A = −1/2
Now that we have understood how to handle the limit ε → 0 + on the infinite time kernel we can extend the same method to study A = −1/2 at finite time. We first rewrite the finite time scalar kernel (87) for any A > −1/2 as follows where we have performed the derivative ∂ x , rearranged the Gamma functions and used the alternative form of the Mellin Barnes formula (46). We consider now the case 0 < A + 1 2 < 1 (since we will perform the limit A → −1/2 below this is sufficient for our purpose). Then we can move the integration contour of both variables from C toĈ which, at finite time iŝ C = a + i with a ∈]A + 1 2 , 1[. This allows to take the A = −1/2 limit, following similar steps as in the previous Section.
To study the A → − 1 2 limit, we will first denoteÃ = A + 1 2 and similarly to the infinite time limit, we conjugate the kernelK t,ς on the left by the multiplication of e −Ãx and on the right by the multiplication of eÃ y and we denote K t,ς the new kernel. This operation leaves again the associated Fredholm determinant unchanged. Evaluating the two residues of the kernel splits it in three parts K t,ς = KÃ 0 + KÃ 1 + KÃ 2 . The first part is the same expression as (111) with C replaced byĈ. We will write it here with the limit A → −1/2 already taken We now determine the residues at w, z =Ã obtained by pushing the contour C intoĈ.
• The residue at z =Ã gives the projector with We only need theÃ = 0 limit of ψÃ( y) given as where in that limitĈ crosses the real axis in the interval ]0, 1[. As in the infinite time limit, the kernel KÃ 1 converges to Remark 5.3. Introducing the Fourier transform of the Airy function and the Barnes integral e t w 3 3 = dr Ai(r)e t 1/3 wr , it is possible to express ψ 0 ( y) as This is reminiscent of the result obtained in Refs. [26,58] about the solution at all times for A = −1/2, see below.
• The residue at w =Ã gives another projector In the limit A → − 1 2 , this kernel converges to To summarize, the A → − 1 2 limit of the kernel K t,ς is given by The structure of this kernel is similar to the one for infinite time in (104). Hence we can follow the same steps and treat the two projectors and the δ ∞ term using the Sherman-Morrison formula and the matrix determinant lemma which allows us to conclude that which translates for the generating function of the exponential of the KPZ height as Our final result for the generating function at all time for A = −1/2 can thus be written more explicitly as with where K 0 0 (x, y) is given in (112) and where in all formula the contourĈ runs parallel to the imaginary axis and crosses the real axis in the interval ]0, 1[. In the large time limit one sets ς = e −t 1/3 s and rescales the integration variables w → wt −1/3 , z → z t −1/3 and perform a similarity transformation on the kernel. The function φ 0 then leads to the Airy function, and one recovers the large time result of the previous Section.
where the set {a i } forms a Pfaffian GOE point process. The equivalence between the two formulae is under investigation.

Conclusion
In this paper we have extended the replica Bethe ansatz solution to the KPZ equation in a half-space for droplet initial condition near the wall, previously obtained for wall parameter A = +∞ and A = 0, to generic value of A. This is an important problem which is believed, via universality, to exhibit a rich GSE/GOE/Gaussian transition first discovered by Baik and Rains. Despite much progress in the topic of exact solutions for the KPZ equation, this has remained a challenge for several years. A recent theorem which maps this problem to the (easier) case of A = +∞ (Dirichlet) but with a Brownian initial condition was inspiring in the solution. Indeed we could apply the same type of summation as was developped in the solution of the Brownian initial condition in the full space. We used the known results for the Bethe ansatz eigenstates of the delta Bose gas in the half-line with Dirichlet boundary conditions. At present, for a technical reason, our solution is valid for all times but limited to A −1/2, i.e. the unbound phase. It allows to demonstrate the convergence to the GSE and GOE Tracy Widom distributions at large time in this phase. We have also obtained the crossover or transition kernel valid in the critical region A + 1 2 = εt −1/3 . The form of our transition kernel is novel, and we have shown that for ε 0 it agrees to a high order with the results of the kernel obtained by completely different methods for discrete models (facilitated TASEP and last passage percolation). This non-trivial check leads us to conjecture that these two kernel are fully equivalent. Additionally, let us note that the phase diagram we have presented for the solution to the KPZ equation echoes the one obtained in Ref. [41] Section 8 for the log-gamma polymer.
In a companion study [29] we have approached the problem differently by elucidating the structure of the boundary bound states which appear in the half-line Bose gas with generic parameter A. It quite easily yields results for the Gaussian phase A < −1/2 where these states dominate. It would be nice to find a bridge between these studies and to obtain a more complete solution of the KPZ equation in the half-space for all values of A. Let us close by mentioning that the physics content of the mapping between the Brownian initial condition and the droplet initial condition with different interactions with the wall remains to be understood deeper.

A Overlap of the half-line Bethe states with the Brownian initial condition
Here we give some details on the calculation of the overlap 〈Ψ µ |Φ 0 〉 between the half-line Bethe states and the Brownian initial condition. We recall that Φ 0 , given in (21), is a fully symmetric function of its arguments, and that in the sector 0 x 1 < · · · x n it equals where w = (A + 1 2 ) 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 〈Ψ µ |Φ 0 〉 * = 〈Φ 0 |Ψ µ 〉 = 〈Ψ µ |Φ 0 〉. Since Ψ µ is also a symmetric function of its arguments, by definition the overlap can thus be written as 〈Φ 0 |Ψ µ 〉 = n! 0< y 1 < y 2 <···< y n d y 1 . . . d y n Ψ µ ( y 1 , . . . , y n )e n j=1 1 2 (2n+1−2 j) y j −w y j .
Now in (130) for each permutation P we can relabel all the p → P(p) and denoting by { }=±1 the operation of summation over all the variables (an operation independent of their labeling) we can rewrite (130) as ) G n,w ( P(1) λ P(1) , . . . , P(n) λ P(n) ), where we have used the fact that the product is independent of the permutation P. Now we use the following "miracle" equality, proved in [71] P∈S n k< Applying it to the set { k λ k } we obtain Let us denote w = A + 1 2 . We now conjecture the following identity valid for any set of complex numbers λ j which we have checked to high order n = 5 using Mathematica. Inserting in (136) we obtain the formula for the overlap (34) 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 λ j , invariant by complex conjugation, such that the overlap integral converges. The condition for that to be true can be read from (135) as for all j ∈ [1, n]. Inserting a string state, labeled by {k j , m j } j=1,...,n s , and specifying w = A + 1 2 , the condition becomes 1 2 max j m j n 2 < A + 1 2 .

B Calculation of the first moment of Z(t)
Starting from the general formula (38) for the moments of Z(t) we can read the explicit expression for KPZ Z(t) , i.e. n = 1, for A > 0.
This result however is not correct for A < 0 as we now show. The Lieb-Liniger model for a single particle reduces to the standard heat equation on the half-line with boundary coefficient A, which is easily solved. The (un-normalized) eigenfunctions are the • For all A ∈ , the following plane wave solves the eigenvalue equation with eigenenergy Its norm on a line of length L is ||Ψ|| 2 = L 2 (1 + k 2 A 2 ). The real momentum k is quantized on the line [0, L]. In the L → ∞ limit, the summation over eigenstates becomes k → L dk 2π . The contribution of this plane wave to the first moment of Z(t) is • For all A < 0, an additional state, i.e. a boundary bound state (BS), with energy −A 2 appears Its norm on the infinite line is ||Ψ BS || 2 = − 1 2A and its contribution to the first moment of where Θ is the Heaviside function.
The total result valid for all A ∈ is thus This can also be obtained by solving the heat equation in time, leading to the same result. We observe that the second contribution Z BS is absent from our result (38). It confirms that our moment formula is valid only in a certain range of values, i.e. A + 1/2 > n/2 for the n-th moment of Z(t). The complete formula thus requires including the boundary bound states.

C Comparison of the Fredholm Pfaffians
As announced in the main text, we now test our conjecture (74), i.e. we compare the Fredholm Pfaffian based on our transition kernel, given in (72), and the one based on the crossover kernel of Ref.

C.2 Kernel with all contours at the right of ε
To be able to track systematically the degree of homogeneity of each integral in the expansion at large s as explained in the text, we move all integration contours toĈ which is at the right of ε. For our kernel this was done in (72), leading to the schematic form (73) where the upper label (p) , p = 1, 2 denotes the number of independent integrals with a cubic exponential, leading to the large positive s decay ∼ exp( 2 3 ps 3/2 ). We need to put the kernel K cross in (75) The important observation is that there is an anomalous term, of homogeneity zero p = 0, coming from R 22 .

C.3 Expansion of the Pfaffians in degree of homogeneity
Consider first the expansion of Pf(J − K ε ) 2 ([s,+∞[) using the formula (146) (and if needed (160) and (148)). We expand in the total degree of homogeneity (i.e. the degree of stretched exponential decay at large s) so we need to pick up terms in various orders in the expansion of the Pfaffian. Fortunately, to expand to third order in homogeneity this order is not too large.
Recalling that Below we will prove the derivatives w.r.t. s of these equalities , which is simpler in general and sufficient since all terms vanish at s = +∞. Once we have shown the first two equalities, we can use them to simplify the third one, so we will need to show 1 2 TrK (1)

C.3.1 Comparison of terms of homogeneity 1
One sees from (72) and (149) The first equality in (154) then holds.

C.3.2 Comparison of terms of homogeneity 2
From (154) we now have to compare We recall that K (2) 12 (r, r ) = and is easily checked to be true using mathematica. Hence the second equality in (154) holds.
Putting now all terms together as in (155) and symmetrizing the preexponential factor in (w, z) as we did for the terms of homogeneity 2, we find, after a tedious calculation (using mathematica) that the symmetrized form vanishes. We believe that the agreement of the first three orders is a rather non-trivial check, which comes in very strong support of our conjecture (74).