Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in recent years , the investigations of non - intersecting brownian motion and random walk paths focused on the description of a critical process called the tacnode process . this process appears when two groups of trajectories are asymptotically supported in two ellipses in the time - space plane such that the ellipses touch each other creating a tacnode ( self - touching point ) of the macroscopic boundary . the aim is to describe the behaviour of the paths near the touching point . a series of recent papers by different groups of authors studied the tacnode process in parallel using various methods . the first result in this direction is due to adler , ferrari and van moerbeke : in @xcite , a model of non - intersecting random walk paths is considered such that the paths form a symmetric tacnode . in @xcite , delvaux , kuijlaars and zhang studied the non - symmetric case of non - intersecting brownian trajectories and they expressed the critical correlation kernel in terms of the solution of a @xmath2 riemann  hilbert problem . a different approach to the brownian case is due to johansson who gave a formula for the tacnode kernel in terms of the resolvent of the airy kernel in @xcite in the symmetric case . the latter approach was extended by ferrari and vet in @xcite to the general non - symmetric case . the tacnode process was also obtained in the tiling problem of the double aztec diamond by adler , johansson and van moerbeke in @xcite . it was not a priori clear that the various formulas for the tacnode kernel give rise to the same limit process , since the results in @xcite , @xcite and @xcite contain airy resolvent type formulas whereas the kernel in @xcite is expressed with the solution of a riemann  hilbert problem . it was shown in @xcite that the formulation of @xcite and @xcite are equivalent . a more recent result @xcite gives the equivalence of the riemann  hilbert formulas in @xcite and the airy resolvent formulas in @xcite and @xcite .    apart from systems of non - intersecting brownian motions and random walk paths , the study of non - intersecting squared bessel paths is also natural due to their representation as eigenvalues of the laguerre process , see  @xcite . this process is a positive definite matrix valued brownian bridges . since the transition kernel of the squared bessel paths can be given explicitly via modified bessel functions , see , a karlin  mcgregor type formula like can be applied . non - intersecting squared bessel paths were also studied by katori and tanemura in @xcite , but the first description of the hard edge tacnode process is due to delvaux @xcite . the formation of the tacnode in this case is slightly different : instead of two touching groups of trajectories in the time - space plane , only one group of paths is considered and its boundary macroscopically touches the hard edge , i.e.  the @xmath3 line . this configuration is referred to as the hard edge tacnode and the critical process in the neighbourhood of the touching point is called the hard edge tacnode process . in @xcite , the correlation kernel of the hard edge tacnode process is expressed in terms of the solution of a @xmath2 riemann  hilbert problem which is different from the one that appears in the case of the brownian trajectories . see the left - hand side of figure  [ fig : simulation ] for non - intersecting squared bessel paths with the tacnode scaling of the starting and endpoints . for further figures , see also  @xcite .    in theorem  [ thm : hardedgetacnode ] which is the main result of the present paper we provide an explicit airy type formula for the multi - time correlation kernel of the hard edge tacnode process , under the assumption that the parameter @xmath4 of the bessel process is a ( non - negative ) integer . the kernel is expressed as a double complex integral , see . the formula is completely new in the literature . it is reminiscent to the usual tacnode kernel , but the airy resolvent operator on @xmath0 is replaced by an operator which acts on the direct sum space of @xmath0 and a finite dimensional space . a different phenomenon appears when the starting points of the non - intersecting squared bessel paths are taken to be @xmath1 and the endpoints are scaled linearly with the number of paths . in this case , there is a critical time such that for any earlier time the lowest path stays close to @xmath1 whereas for any later time the distance of the lowest path from @xmath1 is macroscopic . after rescaling around the critical time , one obtains the hard edge pearcey process , see the right - hand side of figure  [ fig : simulation ] . the single time hard edge pearcey kernel was already described by desrosiers and forrester @xcite and later in a different formulation by kuijlaars , martnez - finkelshtein and wielonsky @xcite . we give the multi - time correlation kernel of the hard edge pearcey process in different formulations in theorems  [ thm : hardedgepearcey ] and  [ thm : pearcey : alternative ] which also shows the equivalence of the formulas in @xcite and @xcite . in corollary [ cor : bk ] , we obtain the kernel by borodin and kuan @xcite in the special case of one - dimensional squared bessel paths , i.e.  of absolute values of one - dimensional brownian motions . the main steps how the convergence of the kernel of @xmath5 non - intersecting squared bessel paths to that of the hard edge tacnode is proved are the following . first we write the kernel as a double complex contour integral using the representation of the modified bessel function . the integrand can be transformed into the ratio of two toeplitz determinants of sizes @xmath6 and @xmath5 respectively . the toeplitz determinants have symbols with non - zero winding numbers around @xmath1 , hence we apply a generalized borodin  okounkov formula due to bttcher and widom in @xcite to obtain theorem  [ theorem : toeplitz ] if the parameter @xmath7 of the squared bessel paths is a ( non - negative ) integer . then we obtain the ratio of two fredholm type determinants where the two operators are rank one perturbations of each other . this yields a resolvent type formula in theorem  [ theorem : kernel : n ] which is suitable for asymptotic analysis . the functions that appear in the finite @xmath5 kernel can be rewritten in terms of bessel functions , hence the asymptotic analysis relies on the convergence of bessel functions which is proved separately in the appendix . now we introduce the model of non - intersecting squared bessel paths that we consider . the squared bessel process depends on a parameter @xmath8 . the transition probability of the squared bessel process for any time @xmath9 is defined by @xmath10 for @xmath11 and @xmath12 where @xmath13 is the modified bessel function which can be given by the series @xmath14 the transition probability @xmath15 is obtained by taking the limit @xmath16 . if @xmath17 is an integer , then the squared bessel process can be obtained as the squared absolute value of a @xmath18-dimensional brownian motion . in this case , we call @xmath18 the dimension of the squared bessel process . in the present paper , we consider non - intersecting squared bessel paths which start from one fixed point and end at another fixed point . in order to construct this system of paths , we first take @xmath5 non - intersecting squared bessel paths with fixed different starting points @xmath19 at time @xmath20 and fixed different ending points @xmath21 at time @xmath22 . the paths are conditioned to be non - intersecting in the time interval @xmath23 . it is well known that this defines an extended determinantal point process . that is , the joint probability at a sequence of times @xmath24 can be expressed via the determinant of an extended correlation kernel @xmath25 . the kernel @xmath26 is defined in terms of the transition probability @xmath27 in . namely as given also in ( * ? ? ? * eq .  ( 1.1 ) ) , we have @xmath28 for any positions @xmath29 and times @xmath30 with @xmath31 denoting the characteristic function of @xmath32 and with @xmath33 defined as the @xmath34 matrix @xmath35    next we take the confluent limit of the starting and ending points @xmath36 and @xmath37 . if @xmath38 and @xmath39 are suitably scaled with the number @xmath5 of the paths , we can create a picture in the time - space plane with a cusp or a tacnode at the hard edge . in this paper , we obtain the limiting extended correlation kernel of non - intersecting squared bessel paths near the tacnode and the cusp which we call the hard edge tacnode process and the hard edge pearcey process . the paper is organized as follows . we first state our main results in section  [ s : mainres ] . the correlation kernel @xmath26 is expressed as a ratio of two toeplitz determinants and in a pre - asymptotic form for finite @xmath5 in subsection  [ ss : finiten ] . then two different scalings and the corresponding limit processes are considered : the hard edge tacnode process is introduced and discussed along with our results on the convergence in subsections  [ ss : hetacnode ] and [ ss : alt_tacnode ] . the hard edge pearcey process with different formulations of its correlation kernel and our results are given in subsection  [ ss : hepearcey ] . the finite @xmath5 formulas for the kernel @xmath26 are proved in section  [ s : derivation_finite ] . the asymptotic analysis for the hard edge tacnode process is performed in section  [ s : analysis_tacnode ] , the one for the hard edge pearcey process is in section  [ s : analysis_pearcey ] . the alternative hard edge pearcey formulas are proved in section  [ s : alternative_pearcey ] . the proof of proposition  [ prop : limschur ] which is needed for the existence of the given formulation of the hard edge tacnode kernel is postponed to section  [ s : asymptinv ] . section  [ s : bessel_limit ] of the appendix contains a statement about the convergence of the derivatives of bessel functions to those of the airy function and useful tail bounds . we first report our formulas for @xmath5 non - intersecting squared bessel paths for @xmath5 finite . next we describe the limits under the tacnode and the pearcey scaling .      in what follows , we will assume that the parameter @xmath7 of the squared bessel process is a ( non - negative ) integer . see however remark  [ remark : nonint ] where we state an extension of one of our results for finite @xmath5 to the case of non - integer @xmath7 . define the weight function @xmath40 on the unit circle @xmath41 with parameters @xmath42 . we abbreviate @xmath43 where @xmath44 denote the starting and ending point of the non - intersecting squared bessel paths . we also define the weight function @xmath45 where @xmath46 are free parameters for the moment , but they will depend on the integration variables in as given by below . first , we express the kernel of @xmath5 non - intersecting squared bessel paths as a double integral where the integrand is a ratio of two toeplitz determinants . the proof of the theorem is given in section  [ subsection : proof : theorem : toeplitz ] . [ theorem : toeplitz ] let @xmath7 be an integer and consider @xmath5 non - intersecting squared bessel paths of parameter @xmath7 with starting point @xmath47 and ending point @xmath39 . then the extended correlation kernel can be written as @xmath48 where @xmath49 denotes the unit circle in the complex plane , oriented counterclockwise and @xmath50 is a ratio of two toeplitz determinants with weight functions @xmath51 and @xmath52 defined in  with @xmath53    [ remark : nonint ] theorem  [ theorem : toeplitz ] can be extended to arbitrary real values of the parameter @xmath7 ( @xmath8 ) , not necessarily integer . to that end , it suffices to replace both integration contours @xmath49 in by a contour @xmath54 where @xmath54 is a contour encircling the origin in counterclockwise direction , beginning at and returning to @xmath55 , and never intersecting the negative real line except at @xmath55 . the contour @xmath54 is the negative of the standard ` hankel contour ' . we then assume that all powers in the formulas have a branch cut along the negative real axis . note that in the special case where @xmath7 is integer , the integration over the contour @xmath54 can be replaced by an integration over the unit circle @xmath49 and then we retrieve the formula stated in theorem  [ theorem : toeplitz ] above . the proof of the extension to non - integer @xmath7 is similar to the one of theorem  [ theorem : toeplitz ] given in section  [ subsection : proof : theorem : toeplitz ] , by noting that holds for non - integer values of @xmath4 provided that the integration contour @xmath49 is replaced by the contour @xmath54 described above . we omit the details . unfortunately , we were unable to apply an asymptotic analysis to the toeplitz determinants for this generalized setting with @xmath4 non - integer , therefore we will always assume below that @xmath4 is integer .    in order to calculate the asymptotics of the toeplitz determinants in , we apply a borodin  okounkov type formula . note that the symbols of our toeplitz matrices have winding number @xmath7 around the origin , which is in general non - zero . a borodin  okounkov type formula in this setting is described by bttcher  widom @xcite . the formulas in @xcite allow to write each toeplitz determinant as the fredholm determinant of an operator on @xmath56 times a determinant of size @xmath57 . alternatively , we will see that the toeplitz determinant can be expressed via the fredholm determinant of a _ single _ operator acting on the direct sum space @xmath58    fix @xmath59 . we define a block matrix operator of size @xmath60 @xmath61 which we view as the matrix representation of a linear operator acting on the space @xmath62 with respect to the natural basis of @xmath62 . the block matrix decomposition is compatible with the direct sum decomposition of @xmath62 . let us describe the four blocks of . we will denote by @xmath63 the circle with radius @xmath64 . * @xmath33 is a semi - infinite matrix with @xmath65th entry @xmath66 for @xmath67 . * @xmath68 is a matrix of size @xmath69 with @xmath65th entry given by @xmath70 for @xmath71 . * @xmath54 is a toeplitz matrix of size @xmath72 with @xmath65th entry @xmath73 * @xmath74 is a toeplitz matrix of size @xmath75 with @xmath65th entry given by @xmath76    next we define two vectors in the space @xmath62 . the vectors will depend on parameters @xmath77 which we consider for the moment to be fixed numbers such that @xmath78 and @xmath79 . denote again by @xmath63 the circle of radius @xmath80 which we now take such that @xmath81 . we define @xmath82 to be the column vectors of length @xmath83 and @xmath4 respectively with @xmath84th entry @xmath85 let @xmath86 be the row vector of length @xmath83 with @xmath87th entry @xmath88 and we define @xmath89 as the row vector of length @xmath4 given by @xmath90 with these notations , let @xmath91 where the matrix is the same as the one given in and the entries of the two vectors are defined in  and  respectively . the block matrix and block vector notation is again compatible with the direct sum decomposition of the space @xmath62 given in .    the inverse of the block matrix in the above formula can again be partitioned as a block matrix of size @xmath60 : @xmath92 where @xmath93 is the schur complement of @xmath33 . hence the invertibility of the block matrix in follows from the invertibility of @xmath33 and @xmath94 . it is not clear a priori that the block matrix in the definition of @xmath95 in is invertible , but as it will be seen later in proposition  [ prop : limschur ] for the tacnode scaling and in for the pearcey scaling , we get invertible matrices in these two limits . therefore , @xmath95 is certainly well - defined if the parameters are close enough to any of these limits . [ theorem : inner : product ] the ratio of toeplitz determinants in can be written as @xmath96 where @xmath95 is given above by . the proof of the theorem is given in section  [ subsection : proof : theorem : bo ] . [ def : contour : xi : eta ] we denote by @xmath97 and @xmath98 counterclockwise oriented closed contours as follows . @xmath97 is a loop encircling @xmath1 but not the points @xmath99 and @xmath100 . @xmath98 consists of two pieces : a loop encircling @xmath1 lying at the inside of @xmath97 , and a small loop surrounding the point at @xmath101 lying at the outside of @xmath97 . a particular choice is shown later in figure  [ fig : xietacontours ] . [ theorem : kernel : n ] under the same assumptions as in theorem  [ theorem : toeplitz ] , the correlation kernel @xmath26 for @xmath5 non - intersecting squared bessel paths can be written as @xmath102 where @xmath95 is defined by and the contours @xmath98 and @xmath97 are given in definition  [ def : contour : xi : eta ] . it turns out that it is more convenient for our purposes to work with @xmath103 when we consider the hard edge tacnode process . the free parameter @xmath104 will parametrize the location of the hard edge tacnode process , however @xmath105 will disappear from the limit process . [ def : tacscale ] let the _ tacnode scaling _ be the scaling when time and space are scaled according to @xmath106 and the starting and ending points are rescaled by @xmath107 where @xmath108 is the _ temperature parameter_.    for @xmath108 , we define the shifted airy kernel by @xmath109 to state the result about the limiting kernel , we define the function @xmath110 and the two - variate functions ( with @xmath111 denoting the @xmath84th derivative of the airy function ) @xmath112 the matrix component @xmath113 always appears with a subscript throughout this paper in order to avoid confusions with the set of complex numbers . let @xmath114 and for any function @xmath115 or @xmath116 , we define the column vector on the space @xmath117 @xmath118=\begin{pmatrix } ( f(x))_{x\in{\mathbb{r}}_+ } & f(x)|_{x=0 } & \frac{{\partial}}{{\partial}x } f(x)|_{x=0 } & \ldots & \frac{{\partial}^{{\alpha}-1}}{{\partial}x^{{\alpha}-1 } } f(x)|_{x=0 } \end{pmatrix}^t.\ ] ] finally define @xmath119.\ ] ]    [ def : contours ] the contours depend on real parameters @xmath120 but we will not indicate this dependence in the notation . let @xmath121 be a jordan arc in the complex plane which comes from @xmath122 , goes to @xmath123 and it crosses the real axis on the right of @xmath124 and @xmath125 . let @xmath126 consist of two pieces : a jordan arc in the complex plane coming from @xmath127 and going to @xmath128 , lying to the right of @xmath121 , and a clockwise oriented small loop around @xmath124 , lying to the left of @xmath121 . a possible choice of these contours can be seen in figure  [ fig : hxicontours ] . let @xmath129 be the hard edge tacnode kernel where the integration contours are given in definition  [ def : contours ] and the function @xmath130 is defined by . the next theorem will be proved in section  [ s : analysis_tacnode ] . [ thm : hardedgetacnode ] let @xmath7 be a non - negative integer .    1 .   under the tacnode scaling given by definition  [ def : tacscale ] , the sequence of rescaled kernels converges @xmath131 uniformly as @xmath132 and @xmath133 are in a compact subset of @xmath134 where @xmath135 and the kernel on the right - hand side of is defined by . 2 . as a consequence , the hard edge tacnode process @xmath136 exists as the limit of @xmath5 non - intersecting squared bessel processes under the tacnode scaling . it is characterized by the following gap probabilities . for any fixed @xmath84 and @xmath137 and for any compact set @xmath138 , @xmath139    in the definition of the kernel @xmath140 in  , we implicitly assume that the limiting block matrix is invertible . equivalently , the limiting schur complement @xmath141 should be invertible . this follows from proposition  [ prop : limschur ] below for @xmath142 large enough . let @xmath143 be the integral operator on @xmath144 with kernel @xmath145 note that @xmath146 and the matrices in  can also be easily expressed in terms of the operator @xmath143 as @xmath147    we will use the following general notation in the sequel : if @xmath148 is a smooth function on @xmath149 , we denote the kernel of the integral operator @xmath150 on @xmath144 by @xmath151 . in particular , we have @xmath152 in below . straightforward substitution of  into yields the following . [ lemma : schur : repr ]    1 . the schur complement @xmath153 in can be represented as @xmath154 2 . the schur complement @xmath153 is a symmetric matrix .    in particular , for @xmath155 the schur complement takes the scalar form @xmath156 it is well - known @xcite that this is equal to @xmath157 where @xmath105 is the hastings  mcleod solution to the painlev  ii equation . the invertibility of the schur complement then amounts to the fact that @xmath158 for all @xmath159 . we conjecture that it holds more generally . [ conj : schur : inv ] the schur complement is invertible for all @xmath159 .    let us prove conjecture [ conj : schur : inv ] under the assumption that @xmath160 is large enough . note that for @xmath161 the operator @xmath162 tends exponentially fast to the identity operator . hence the @xmath65th entry of the matrix is approximated by @xmath163 we will prove the following proposition in section  [ s : asymptinv ] .    [ prop : limschur ]    1 . let @xmath5 be fixed . we have the asymptotics @xmath164 for @xmath165 where @xmath166 means that the ratio of the two sides goes to @xmath167 . the schur complement @xmath153 is invertible for all @xmath160 large enough . we take the derivative of the kernel   with respect to the temperature parameter @xmath142 . it turns out that this derivative has a convenient form : it is a rank 1 kernel . the proof of this fact is a rather lengthy calculation which we omit here , hence we give the next theorem without proof . [ theorem : temp : der ] the derivative of the hard edge tacnode kernel with respect to @xmath142 is a rank 1 kernel , that is @xmath168 if @xmath169 , then @xmath170 and they both are given by the contour integral @xmath171 if @xmath172 , then the function @xmath173 is given by @xmath174\ ] ] and @xmath175 can be represented as @xmath176\right].\end{gathered}\ ] ]    the functions @xmath173 and @xmath175 are the same . the conjecture holds for @xmath169 trivially , the @xmath177 case can still be checked by rewriting the functions in schur complement forms . for general @xmath7 , it is more complicated . the equality of the two functions @xmath178 and @xmath179 follows from the hidden symmetry of the kernel . integration of with respect to the temperature yields the following . the hard edge tacnode kernel can be written as @xmath180    finally , we show how the block matrix notations in the above formulas can be avoided . we show this first for the function @xmath130 . let us use the notation @xmath181 for the exponential function for any @xmath182 and let @xmath183 be the usual scalar product of functions @xmath148 and @xmath184 in @xmath0 . observe that with this notation @xmath185    [ lemma : mschur ] we can rewrite the function @xmath130 in with the notation of and as @xmath186 here the two expressions between large parentheses @xmath187 denote a row and a column vector respectively , both of length @xmath4 . in a very similar way , one can express the functions @xmath178 and @xmath179 without the block matrix notation . [ rem : cafasso ] in the @xmath169 case after the change of variables @xmath188 and @xmath189 , the kernel in restricted to single time reads @xmath190 the contour @xmath191 is a clockwise curve that has a cusp at the origin such that it leaves the origin in the @xmath192 direction , it crosses the positive real axes and it returns to the origin from the @xmath193 direction . the other contour @xmath194 consists of two parts : the first one is a clockwise loop inside @xmath191 and it has a cusp at the origin such that it leaves the origin in the @xmath195 direction and it returns to the origin from the @xmath196 direction ; the other part is a counterclockwise circle around @xmath191 . the kernel in with a further change of variables is a special case of the single time kernels of the form @xmath197 where @xmath198 is a rational function in @xmath199 which depends on the parameter @xmath132 linearly , i.e.  @xmath200 . in this case , @xmath201 is a third order rational function . in the hard edge pearcey case , it is a second order function , see corollary  [ cor : kmw ] below and also @xcite . as discussed in @xcite , to the type of kernels with some further restrictions on @xmath201 , one can associate a riemann  hilbert problem which was used in @xcite to express the gap probabilities in terms of the tau function of the riemann  hilbert problem . it is natural to ask about the existence of a corresponding riemann  hilbert problem in the present case . here we consider the non - intersecting squared bessel paths which all start from @xmath1 . more precisely , let @xmath202 be the confluent limit of the kernel of @xmath5 non - intersecting squared bessel processes as the starting points tend to @xmath1 . this means that @xmath203 is the correlation kernel for @xmath5 non - intersecting squared bessel processes conditioned to start at @xmath1 at time @xmath1 and end at positions @xmath204 at time @xmath167 . this is similar to the situation considered by kuijlaars , martnez - finkelshtein and wielonsky ( * ? ? ? [ prop : confluent : a0 ] for any real @xmath8 , the kernel @xmath205 in has the double integral representation @xmath206 where the @xmath207-contour @xmath208 is an counterclockwise loop surrounding the points @xmath209 and being disjoint from the @xmath210-contour @xmath211 $ ] . the branch for @xmath212 is defined by @xmath213 and along the contour @xmath208 , we use the extension of @xmath214 from the positive real axis .    this proposition is closely related to a result by katori  tanemura ( * ? ? ? * theorem 2.1 ) . it can also be considered as the hard edge analogue of tracy  widom ( * ? ? ? * eq .  ( 2.11 ) ) . for completeness we will provide a proof of the proposition in section  [ s : analysis_pearcey ] the confluent limit of the endpoints @xmath215 for @xmath216 is simply obtained by substituting it in . if we choose @xmath217 where @xmath104 is a free parameter , then the region filled by the trajectories of the bessel processes behaves as follows . for any time @xmath218 , the lowest path stays close to @xmath1 whereas for @xmath219 , its distance from @xmath1 is macroscopic . our main focus is the neighbourhood of @xmath220 under the following scaling where the hard edge pearcey process is observed . [ def : peascaling ] let the _ pearcey scaling _ be the following scaling of the parameters . we consider @xmath5 non - intersecting squared bessel paths , take the confluent limit as all the starting points are @xmath1 and the endpoints are equal to @xmath217 with a free parameter @xmath104 and we rescale time and space variables @xmath221 and @xmath133 as @xmath222 respectively . ( we make this rescaling for all the considered time and space variables . ) we define the hard edge pearcey kernel by @xmath223 where @xmath54 is the contour which consists of two rays : one from @xmath224 to @xmath1 and one from @xmath1 to @xmath225 . in the single time case @xmath226 , the above kernel essentially reduces to the one of desrosiers  forrester  ( * ? ? ? 5 ) which was also discussed in the last paragraph of ( * ? ? ? * sec .  1 ) . applying an asymptotic analysis to the double integral formula in proposition  [ prop : confluent : a0 ] similarly to ( * ? ? ? iii ) , we obtain the multi - time extension of the hard edge pearcey process . [ thm : hardedgepearcey ] let @xmath8 . the hard edge pearcey process @xmath227 is obtained as the limit of @xmath5 non - intersecting squared bessel processes starting at @xmath1 at time @xmath1 and ending at @xmath217 at time @xmath167 under the time - space scaling . it is defined by the following gap probabilities . for any fixed @xmath84 and @xmath137 and for any compact set @xmath138 , @xmath228 where @xmath229 is the hard edge pearcey kernel given by . [ rem : pearceysigmadep ] it is not hard to see from the proof of theorem  [ thm : hardedgepearcey ] that if we set @xmath230 for some fixed @xmath159 , then the correlation kernel of the hard edge pearcey process becomes @xmath231 . similarly as before , one can think of @xmath142 as a _ temperature _ parameter .    for non - negative integer values of @xmath7 , we can use our methods involving toeplitz determinants to obtain an alternative formulation of the kernel @xmath232 . the following result is proved in section  [ s : alternative_pearcey ] . [ thm : pearcey : alternative ] if @xmath7 is a non - negative integer , then we can write @xmath229 in as @xmath233 where @xmath234 is a clockwise oriented circle surrounding the singularity at @xmath124 , and @xmath235 is chosen such that the contour @xmath236 passes to the right of the singularity at @xmath125 and to the right of the contour @xmath234 . [ rem : pearcy : nonint ] although we derive it only if @xmath4 is a non - negative integer , formula   makes sense for any real @xmath237 with a few minor modifications . in that case , we replace @xmath234 by a clockwise oriented closed loop which intersects the real line at a point to the right of @xmath124 , and also at @xmath124 itself where it has a cusp at angle @xmath238 . we then take the principal branches of the powers @xmath239 , @xmath240 , i.e. , with a branch cut along the negative half - line . see also remark  [ remark : kmw : noninteger ] .    in the single time case @xmath241 , the kernel @xmath232 in can be further reduced to the kernel of kuijlaars , martnez  finkelshtein and wielonsky , as stated in  ( * ? ? ? * eq .  ( 1.19 ) ) : @xmath242 where @xmath243 is a clockwise circle in the left half - plane touching zero along the imaginary axis , and @xmath244 is a counterclockwise loop surrounding @xmath243 . [ cor : kmw ] in the single time case @xmath241 , the kernel @xmath232 in reduces to @xmath245    as a consequence of theorem  [ thm : hardedgepearcey ] in the special case of @xmath246 , we recover the process described by borodin and kuan @xcite via the kernel @xmath247 as follows . [ cor : bk ] for @xmath246 , let us choose @xmath248 then the hard edge pearcey process with the above change of the space and time variables is the same as the borodin  kuan process , that is , their kernels satisfy the relation @xmath249    we assume above that the non - intersecting paths start at @xmath250 and end at @xmath217 . instead we could let them start at @xmath251 and end at @xmath252 . this leads to an alternate hard - edge pearcey process that we could denote by the superscript @xmath253 . the corresponding kernels are related by @xmath254 note that the prefactor @xmath255 is just a conjugation of the kernel , hence the two kernels give rise to essentially the same correlation functions . note that this prefactor also appears in the right - hand side of . the equivalence of the two formulas and for the kernel @xmath229 can be seen directly as follows . in both formulations , it can be shown that the derivative of the kernel with respect to the parameter @xmath142 factorizes as @xmath256 where the functions @xmath257 solve the third - order ode s @xmath258 where the prime denotes the derivative with respect to the space variable @xmath132 or @xmath133 . by using this fact and also some extra information to determine the relevant solutions to the ode s , one can prove that the functions @xmath257 in are the same for both formulations of the kernel @xmath232 , and hence the kernels and themselves are equal as well . we leave the details to the interested reader . we note that if @xmath259 the functions @xmath257 in also satisfy the pde s @xmath260 in this section , we prove the different formulations of the kernel @xmath26 stated in theorems  [ theorem : toeplitz ] , [ theorem : inner : product ] and [ theorem : kernel : n ] . the proof of theorem  [ theorem : toeplitz ] consists of two steps . first we express the quantity @xmath261 as a christoffel  darboux kernel , then we rewrite it as the ratio of two toeplitz determinants . [ [ step-1 . ] ] step 1 . + + + + + + +    let us take the confluent limit of the starting and ending points @xmath262 and @xmath37 . the expression reduces to @xmath263 where @xmath264    since we assumed that @xmath7 is an integer , the integral representation for the modified bessel function @xmath265 holds where @xmath49 denotes the unit circle . therefore , the transition probability can be written as @xmath266 with the weight function on the unit circle @xmath267 . we can then evaluate the quantities in  : @xmath268 and @xmath269 inserting these expressions in the kernel and using the linearity of matrix multiplication , we find @xmath270 with @xmath271 where @xmath272 and @xmath273 are arbitrary polynomials of degree @xmath274 and @xmath84 respectively . note that the factors @xmath275 in the denominators in  cancel each other . so we obtain the formula , but with @xmath276 there replaced by the christoffel  darboux kernel @xmath277 with again the notation and . by taking @xmath278 and @xmath279 , this reduces to @xmath280 where @xmath281 note that @xmath33 is a toeplitz matrix . [ [ step-2 . ] ] step 2 . + + + + + + +    next we express the christoffel  darboux kernel @xmath276 in as a ratio of two toeplitz determinants . we discuss this in a broader context . let @xmath51 be an arbitrary weight function on the unit circle @xmath49 . define the moments @xmath282 for @xmath283 . obviously the moments @xmath284 depend only on the difference @xmath285 . so they satisfy the toeplitz property @xmath286 but we will not use this for the moment . let @xmath287 be the moment matrix , i.e. , the toeplitz matrix .    with the setting of the last paragraph , consider the christoffel  darboux kernel . it can be written as a ratio of two toeplitz determinants @xmath288 with the moments @xmath284 given in and @xmath289 which also satisfy the toeplitz property @xmath290 .    the definition of the christoffel  darboux kernel is equivalent with @xmath291 which in turn yields the determinantal formula @xmath292 the denominator of is a toeplitz determinant . we want to achieve the same for the determinant in the numerator . to this end we apply a series of elementary row and column operations . from the @xmath274th row we subtract @xmath293 times the @xmath294st row for @xmath295 , and from the @xmath84th column we subtract @xmath296 times the @xmath297st column for @xmath298 . these operations kill all the entries in the last row and column except for the entries @xmath299 and @xmath300 . hence the determinant factorizes as @xmath301 times the determinant obtained by skipping the last two rows and columns . the lemma follows . note that @xmath302 are expressed as moments @xmath303 where the weight function @xmath52 is given in . this proves , which ends the proof of theorem  [ theorem : toeplitz ] . we prove theorem  [ theorem : inner : product ] in three steps . first we apply a borodin  okounkov type formula for the toeplitz determinants , then we rewrite the ratio of determinants , finally we use the rank one structure to obtain the assertion . [ [ step-1.-1 ] ] step 1 . + + + + + + +    we express both toeplitz determinants in via a borodin  okounkov type formula , in the formulation of bttcher  widom @xcite . let @xmath5 be a fixed integer . recall the semi - infinite matrix @xmath33 in , viewed as an operator on @xmath56 , and similarly define @xmath304 by . both operators @xmath33 and @xmath304 are a trace class perturbation of the identity operator on @xmath56 so their fredholm determinants @xmath305 , @xmath306 are well - defined . then we obtain the following factorizations . [ thm : borodin : okounkov ] with the above notations , the following is true .    1 .   the toeplitz determinant in the numerator of can be factorized as @xmath307 where @xmath308 and @xmath309 will be defined in . 2 . the toeplitz determinant in the denominator of can be factorized as @xmath310 where @xmath311 and @xmath312 will be defined in the proof . the second part follows from the first by taking the limit @xmath313 , @xmath314 and replacing @xmath5 by @xmath315 . so it suffices to prove the first part of the lemma . it will be convenient for us to replace the weight function @xmath52 in the left hand side of by @xmath316 . note that this is equivalent to a conjugation of the toeplitz matrix by the diagonal matrix @xmath317 , which does not affect the value of the determinant . moreover , we multiply the weight function @xmath316 by the prefactor @xmath318 ; this prefactor multiplies the toeplitz determinant in by the factor @xmath319 , but this factor will cancel out later . we use the setting in bttcher  widom @xcite . with the notations in that paper , the weight function of the toeplitz determinant is @xmath320 on the unit circle @xmath41 . note the changes of signs compared to  by the replacement of @xmath52 with @xmath316 , and note that @xmath207 turned into @xmath321 due to the conventions in @xcite . we have the wiener  hopf factorization @xmath322 with @xmath323 where @xmath324 are analytic inside and outside the unit circle respectively and we normalized them as in bttcher  widom ( * ? ? ? * eq .  ( 1 ) ) . following ( * ? ? ? * page  2 ) , set @xmath325 hence @xmath326 let @xmath327 , @xmath328 be the semi - infinite hankel matrices associated to the symbols @xmath329 and @xmath330 respectively . the row and column indices of these matrices are labelled by @xmath331 .    from ( * ? ? ? * theorems  1.1 and 1.3 ) ( with @xmath5 replaced by @xmath6 and with @xmath332 ) , we obtain the required formula for appropriately defined @xmath333 ( see below ) , but with the final factor @xmath334 in replaced by the fredholm determinant @xmath335 where @xmath336 denotes the projection matrix @xmath337 with @xmath338 being the sequence of @xmath1s of length @xmath84 . ( note that the prefactor @xmath339 in ( * ? ? ? * theorem  1.3 ) cancels out with the one at the end of the first paragraph of the present proof , see above . ) we calculate the product of hankel matrices @xmath340 in . its @xmath65th entry , @xmath341 , is given by @xmath342 times @xmath343 hence the @xmath65th entry of @xmath344 is given by @xmath345 which corresponds with our definition of @xmath304 in up to a shift of the indices @xmath346 by @xmath347 . therefore the fredholm determinant equals @xmath348 as desired . next , the expressions for @xmath349 and @xmath350 are easily calculated from bttcher  widom @xcite ( keeping in mind formula above ) : @xmath351 and @xmath352 as desired . finally , we discuss the quantity @xmath309 . from ( * ? ? ? * eq .  ( 4 ) ) , we have @xmath353 where @xmath354 where again @xmath355 denotes the projection matrix @xmath356 with @xmath357 being the sequence of @xmath1s of length @xmath5 . the quantity @xmath358 can be defined similarly . it is also obtained as the limit of @xmath359 as @xmath360 and @xmath314 .    from lemma  [ thm : borodin : okounkov ] we obtain    [ theorem : bo ] the ratio of toeplitz determinants in can be written as @xmath361    [ [ step-2.-1 ] ] step 2 . + + + + + + +    we express the last factor in corollary  [ theorem : bo ] via the determinants of certain block matrices , i.e.  determinants on the space @xmath62 . fix @xmath59 . we define a block matrix operator @xmath362 of size @xmath60 . let us describe the four blocks . denote again by @xmath63 the circle with radius @xmath80 which we choose such that @xmath81 . * @xmath304 is a semi - infinite matrix with @xmath65th entry @xmath363 for @xmath67 . * @xmath364 is a matrix of size @xmath69 with @xmath65th entry given by , but with the @xmath365 term removed and the factor @xmath366 replaced by @xmath367 for @xmath368 and @xmath369 . * @xmath370 is a toeplitz matrix of size @xmath72 with @xmath65th entry @xmath371 for @xmath372 and @xmath373 . * @xmath374 is a toeplitz matrix of size @xmath75 with @xmath65th entry given by , but with the factor @xmath375 replaced by @xmath376 for @xmath377 .    by taking the limit @xmath378 and replacing @xmath5 by @xmath315 in the above formulas for @xmath379 , we retrieve the formulas for the matrices @xmath380 . [ prop : block : matrix ] with the above notations , the last factor in corollary [ theorem : bo ] can be written as @xmath381 where the constant @xmath382 is the same in both formulas . we only prove , the other identity being similar . we recall the general matrix identity @xmath383 for arbitrary matrices @xmath384 with compatible dimensions and @xmath385 invertible . we will apply this identity with the matrices @xmath386 in and @xmath387 . then the left - hand side of reduces to @xmath388 by virtue of . on the other hand , the right - hand side of can be reduced to@xmath389 with @xmath390 . [ [ step-3 . ] ] step 3 . + + + + + + +    we relate the block matrices via a rank 1 formula . let @xmath391 be the upper triangular toeplitz matrix @xmath392    [ prop : rank1:alpha ] the block matrices in proposition  [ prop : block : matrix ] satisfy the identity @xmath393 where the right - hand side of is a rank @xmath167 matrix and we used the notations . recall the matrices @xmath370 and @xmath391 in and respectively . the @xmath65th entry of @xmath394 is given by @xmath395\notag\\ & = \frac{1}{2\pi { \mathbf{i}}}\int_{s_1}{\,\mathrm{d}}w\,\frac{w^{-(k+n+{\alpha}+1)}}{w+\eta}\exp\left(\frac{bw - aw^{-1}}{2}\right)\notag\\ & \qquad\qquad\times\left[(w+\xi)w^l-(\xi-\eta)(w^{l}-(-\xi)^l)\right]\label{proof : rank1:1}\end{aligned}\ ] ] for @xmath372 and @xmath373 . on the other hand , the @xmath65th entry of @xmath54 is given by @xmath396.\ ] ] subtracting from and cancelling terms , we find @xmath397.\ ] ] this proves the matrix equality @xmath398 on account of and . in a similar way , one shows that @xmath399 which yields the second block column of . finally , from the identity @xmath400 one obtains the matrix equation @xmath401 on account of . this yields the first block column of . [ lemma : rankone : doublesum ] the last factor in corollary [ theorem : bo ] can be written as @xmath402    this follows from propositions  [ prop : block : matrix ] and [ prop : rank1:alpha ] . theorem  [ theorem : inner : product ] is now an easy consequence of corollary  [ theorem : bo ] and corollary  [ lemma : rankone : doublesum ] . this ends the proof of theorem  [ theorem : inner : product ] .      in what follows , we show that @xmath403 where @xmath95 is defined in where in the definitions of @xmath404 in  we now take the circle radius @xmath80 such that @xmath405 and where the contours @xmath406 and @xmath407 are as follows . let @xmath406 and @xmath407 denote counterclockwise oriented closed contours such that @xmath78 and @xmath79 holds for all @xmath408 and @xmath409 and such that @xmath406 encircles @xmath1 but it does not encircle the points @xmath99 and @xmath100 and @xmath407 does surround the points @xmath99 and @xmath100 . see figure  [ fig : primecontours ] for illustration .    as a first step for proving , we combine theorems  [ theorem : toeplitz ] and [ theorem : inner : product ] along with and to get @xmath410 where we recall that @xmath46 are defined in terms of @xmath411 by . next we show that the circle radius @xmath80 in the definitions of @xmath404 can be taken as in . recall that we had @xmath81 in the original definitions of @xmath404 in . if we enlarge @xmath80 as in , then a residue term will be picked up in each of these definitions . these residue terms ( ignoring the signs ) are given by @xmath412 we claim that the contributions from these residue terms all drop out when inserting them in the double integral via . note that the dependence on @xmath293 and @xmath413 ( that is , on @xmath414 and @xmath415 ) factorizes in the second term of if we substitute it in . hence in the second term , we can separate the @xmath414 and the @xmath415 integrals . first , we show that the contribution coming from the residue @xmath416 vanishes . note that there are essential cancellations when we multiply by the @xmath415 ( and @xmath413 ) dependent factors in . the resulting contribution written as an @xmath413 integral is @xmath417 where the integration contour is the image of the unit circle under , that is the circle of radius @xmath418 around @xmath99 . note that the integrand is analytic outside this circle ( also at @xmath83 ) , hence we can blow up the contour and conclude that the integral is @xmath1 by cauchy s theorem . the argument for @xmath419 is similar , for @xmath420 after the change from @xmath414 to @xmath293 according to , we can shrink the contour to @xmath1 . next we change the integration variables @xmath411 in the double integral in to @xmath46 as defined in . the integration contours also transform under , i.e.  the contour for @xmath293 becomes the circle of radius @xmath421 around @xmath422 and the contour for @xmath413 will be the circle of radius @xmath418 around @xmath423 . note that for any @xmath424 , the point @xmath100 is outside the @xmath293 contour and for any @xmath425 , the point @xmath99 is inside the @xmath413 contour . since the @xmath293 integral has no singularity at @xmath99 and the @xmath413 integral has no singularity at @xmath100 , the contours can be deformed to @xmath406 and @xmath407 as shown in figure  [ fig : primecontours ] by cauchy s theorem and reduces to . now we are ready to derive theorem  [ theorem : kernel : n ] from . note that the order of the @xmath293 and @xmath413 integration in the double integral can be interchanged . this is due to fubini s theorem and the fact that the contours @xmath406 and @xmath407 are disjoint . we deform the integration contours in formula from @xmath406 and @xmath407 to the new contours @xmath97 and @xmath98 in the following steps . first we enlarge @xmath407 to @xmath426 if necessary and we move the contour @xmath406 to @xmath97 inside @xmath426 . then , for any fixed @xmath427 , we pick up a residue contribution when deforming @xmath426 to @xmath98 due to the @xmath428 term in @xmath95 . the residue term can be written as @xmath429 the integrand above is analytic away from @xmath99 and @xmath100 , in particular inside @xmath97 , hence it vanishes thanks to cauchy s theorem . in this section , we prove theorem  [ thm : hardedgetacnode ] . the proof of the first part follows the lines of the classical method of steep descent analysis . the leading terms in the exponent of the kernel @xmath26 can be expressed by the function @xmath430 we define the following contours with counterclockwise orientation . if @xmath431 and close enough to @xmath167 , then let @xmath432 where @xmath433 and @xmath434 . if @xmath435 and close to @xmath167 , then let @xmath436 where @xmath437 and @xmath438 . the next lemma is about the steep descent properties of these paths . [ lemma : steep ]    1 . there is an @xmath439 close enough to @xmath167 for which @xmath440 defines a steep descent contour for the function @xmath441 with respect to the critical point at @xmath442 . 2 .   there is an @xmath443 close enough to @xmath167 for which @xmath444 defines a steep descent contour for the function @xmath445 with respect to the critical point at @xmath442 .    to show the steep descent property , we consider the following derivatives . in the neighborhood of @xmath442 , we have @xmath446 if @xmath447 . furthermore , @xmath448 if @xmath449 . on the circular part of the contour , @xmath450 if @xmath431 , is negative for @xmath451 $ ] and positive for @xmath452 $ ] , whereas for @xmath435 , the sign of is opposite . this proves the required property of the paths . we use the method of steep descent analysis . if we substitute  and the rescaled space and time variables into the double integral in , one sees that the leading term in the exponent is @xmath453 . we also substitute for @xmath95 . next we specify the initial contours @xmath98 and @xmath97 to be the steep descent contours @xmath454 and @xmath455 . @xmath455 coincides with @xmath456 except for an @xmath457 neighbourhood of @xmath458 and there it crosses the real axis to the right of @xmath99 and @xmath100 . @xmath454 consists of two parts : the first one coincides with @xmath459 except for an @xmath457 neighbourhood of @xmath458 and it remains inside @xmath455 . the other part is a small ( order @xmath457 ) circle around the singularity at @xmath101 which is also in the @xmath457 neighbourhood of @xmath458 due to the scaling . see figure  [ fig : xietacontours ] . now we are ready for the steep descent analysis . consider the kernel with the integration contours replaced by @xmath454 and @xmath455 . let us choose a small but fixed @xmath235 and denote @xmath460 and @xmath461 . for @xmath462 small enough , these contours do not contain the circular part of the original steep descent contour coming from the last part in and , but they completely contain the order @xmath457 modifications if @xmath463 is large enough . by neglecting the integral over @xmath464 and @xmath465 , we make an error of order @xmath466 .    in the @xmath462 neighborhood of @xmath458 , we can replace @xmath467 by its taylor series . now we can do the change of variables given by @xmath468 the reader should not be confused by the fact that @xmath413 and @xmath293 were initially defined by @xmath414 and @xmath415 in where @xmath414 and @xmath415 played different role than in the rescaling . for the factor giving the main exponential term in , we have @xmath469 using also . the other factors after substituting the scaling on the left - hand side of become @xmath470 and @xmath471 finally , the change of variables gives a factor @xmath472 and @xmath473 goes to @xmath167 . together with proposition [ prop : detdetlimit ] , this shows that the double integral in multiplied by the prefactor @xmath474 is equal to @xmath475 up to an error of order @xmath476 in the exponent coming from the taylor approximation and the additive @xmath466 error . by the identity @xmath477 , the error in the exponent can be removed by an additive @xmath476 error . finally , we extend the integral over @xmath478 to infinity in the @xmath479 directions and the integral over @xmath480 to infinity in the @xmath481 directions on the cost of an @xmath466 error again . this proves the convergence towards the double integral on the right - hand side of . the term @xmath482 appears by simple substitution . if @xmath132 and @xmath133 are confined in a compact subset of @xmath134 , then all the above error estimates can be given uniformly in @xmath132 and @xmath133 yielding the stated uniform convergence . this completes the proof of the first part of the theorem . the second part follows directly from the first one since the uniform convergence of the kernel on bounded intervals implies the convergence of fredholm determinants on compact sets by the dominated convergence theorem . [ prop : detdetlimit ] under the scaling , we have @xmath483 as @xmath484 , where @xmath95 is given by and @xmath130 is defined in . furthermore , @xmath485 . first we will show that @xmath486 is independent of the parameter @xmath105 . to see this , define the diagonal matrices @xmath487 in equation , the matrix product on the right - hand side can be conjugated by the invertible matrices @xmath488 . this yields @xmath489 let us show that the right - hand side of is independent of @xmath105 . the idea is to perform a change of variables @xmath490 and @xmath491 in the defining integrals of @xmath492 , @xmath493 , @xmath494 , @xmath495 , @xmath496 , @xmath497 and @xmath498 , see  and . it is not hard to see that after this change of variables , @xmath499 becomes @xmath492 with @xmath413 replaced by @xmath500 and with @xmath501 and @xmath502 defined for @xmath503 , i.e.  @xmath504 in . @xmath505 will be @xmath506 with @xmath413 replaced @xmath500 and with @xmath501 and @xmath502 defined for @xmath503 in . similarly , @xmath507 is the same as @xmath494 with @xmath293 replaced by @xmath508 and with @xmath501 and @xmath502 defined for @xmath503 in and @xmath509 is equal to @xmath510 with @xmath293 replaced by @xmath508 in . by the same change of variables , one sees that the block matrix on the right - hand side of is equal to @xmath511 with @xmath501 and @xmath502 defined for @xmath503 in . note that by , replacing @xmath413 and @xmath293 by @xmath500 and @xmath508 means the same scaling as taking @xmath413 and @xmath293 for @xmath503 . hence we have shown that @xmath486 is independent of @xmath105 . therefore we can and do assume that @xmath503 in what follows . denote by @xmath512 ( resp . @xmath62 ) the lower triangular ( resp . upper triangular ) pascal triangle type matrix @xmath513 in equation , the matrix product on the right - hand side can be further conjugated by the invertible matrices @xmath514 . this yields ( recall @xmath503 ) @xmath515 now follows from propositions  [ prop : functionlimits ] , [ prop : operatorlimits ] , [ prop : functopbounds ] below and the dominated convergence theorem . note that the right - hand side of can be rewritten by using the block matrix inversion formula . the matrix inverse in the limit is well - defined since the operator @xmath516 is invertible for any @xmath108 , and for @xmath142 large enough the limiting schur complement @xmath517 is also invertible by proposition  [ prop : limschur ] . finally , the second assertion follows from the estimates given in proposition  [ prop : functopbounds ] . [ prop : functionlimits ] assume that @xmath503 . let @xmath518 be fixed and choose @xmath519 . we have the following limits for the functions in : @xmath520 as @xmath484 pointwise . [ prop : operatorlimits ] assume that @xmath503 . fix the integers @xmath518 and let @xmath519 . then the operators in converge as follows : @xmath521 as @xmath484 pointwise .    we have the following decay properties of the functions and operators . [ prop : functopbounds ] assume that @xmath503 . for any @xmath518 and any constant @xmath522 , there are constants @xmath523 such that @xmath524 uniformly for @xmath525 and as @xmath526 for some thresholds @xmath527 and @xmath528 . the strategy for proving propositions  [ prop : functionlimits ] and [ prop : operatorlimits ] is the following . to keep the computations simpler , we first prove for @xmath529 , and we point out the changes for general @xmath530 . we first rewrite our functions and operators using the bessel function @xmath531 defined by @xmath532 it turns out that some of the operators and functions can be written even more conveniently using the function @xmath533 which is an approximation of the @xmath84th derivative of the airy function as seen from proposition  [ prop : jconv ] . we remark that one can also rewrite @xmath534 by its definition and the binomial theorem as @xmath535 then the limit statement of proposition  [ prop : jconv ] and the bound of proposition  [ prop : jbound ] imply the claimed convergence results in propositions  [ prop : functionlimits ] and [ prop : operatorlimits ] by dominated convergence . the bounds in proposition  [ prop : functopbounds ] follow from the bounds given for the dominated convergence argument in the proofs of propositions  [ prop : functionlimits ] and [ prop : operatorlimits ] . suppose first that @xmath529 . in the definition of @xmath492 we write @xmath536 and use  ( with @xmath503 and @xmath529 ) to obtain @xmath537 on account of . if we substitute the scaling of @xmath413 from ( with @xmath503 ) and @xmath538 and replace the summation by an integral for @xmath539 , then the above expression up to @xmath457 error is @xmath540 the factor @xmath541 converges to @xmath542 and it is bounded by @xmath543 uniformly in @xmath463 . it follows from propositions  [ prop : jconv ] and [ prop : jbound ] for @xmath544 that the factor @xmath545 converges to @xmath546 pointwise and that its decay in @xmath547 ( and in @xmath132 ) is faster than any exponential , in particular , faster than @xmath548 . this shows that the dominated convergence theorem applies and it yields for @xmath529 .    for general @xmath530 , we have @xmath549 instead of @xmath550 on the right - hand side of . then we use proposition  [ prop : jconv ] as @xmath551 after the change of variable @xmath552 in the integral obtained as the limit of @xmath553 , one gets . suppose again @xmath529 . in the definition of @xmath493 , we write @xmath554 as a geometric series as in to get @xmath555 by the definition of the matrix @xmath512 and by that of the function @xmath556 , we have @xmath557 by the same dominated convergence argument as above using propositions  [ prop : jconv ] and [ prop : jbound ] for @xmath558 , one can conclude the proof of for @xmath529 . the general case follows as before . the proof of is similar . we write @xmath559 as a geometric series . next we substitute @xmath560 . this allows for expressing @xmath494 with bessel functions as @xmath561 the dominated convergence argument and the extension to general @xmath530 is the same as for @xmath492 . finally , the proof of is straightforward . since the steps of this proof are rather similar to those of proposition  [ prop : functionlimits ] , we only give a sketch here , from which the reader can complete the details . we also restrict ourselves to the @xmath529 case .    by rewriting @xmath562 as a geometric series , substituting @xmath560 and using  ( with @xmath503 ) , we get @xmath563 for @xmath538 and @xmath564 , an application of propositions  [ prop : jconv ] and [ prop : jbound ] for @xmath544 and the dominated convergence theorem gives . similarly , @xmath565 as in by applying the matrix @xmath512 , one gets @xmath566 a dominated convergence argument as before yields . it is easy to see that @xmath567 if we apply the matrix @xmath62 from the right , the convergence follows from the equation @xmath568    similarly , @xmath569 after applying the matrix @xmath512 from the left and the matrix @xmath62 from the right , one can rewrite as @xmath570 then follows . this completes the proof of proposition  [ prop : operatorlimits ] . the assertions easily follow from the bounds given in the proofs of propositions  [ prop : functionlimits ] and [ prop : operatorlimits ] for the dominated convergence argument , which were obtained from proposition  [ prop : jbound ] and from obvious estimates . in this section , we first prove proposition  [ prop : confluent : a0 ] which is the starting formula for the asymptotic analysis for the hard edge pearcey process . then , we turn to the actual asymptotic analysis , i.e.  we prove theorem  [ thm : hardedgepearcey ] . we use the following lemma for the proof of proposition  [ prop : confluent : a0 ] . [ lemma : kernel : joh ] the kernel @xmath26 for @xmath5 non - intersecting squared bessel paths in can be written as @xmath571 for any @xmath29 and @xmath30 where we use the notation @xmath572    we follow johansson ( * ? ? ? * ( 2.7)(2.10 ) ) . define the column vector @xmath573 applying cramer s rule to , we find @xmath574 where @xmath575 denotes the matrix @xmath33 in with its @xmath84th column replaced by the vector @xmath576 . the following markov type formula holds for the transition probabilities of the squared bessel paths given in for negative time spans : @xmath577 for any @xmath29 and @xmath578 . here we define the branches of the powers in  by @xmath579 and @xmath580 when @xmath207 is negative . the markov type formula is the analogue of ( * ? ? ? * ( 2.5)(2.6 ) ) . from , we get @xmath581 with the notation . by inserting into and interchanging the sum and the integral , we obtain . this proves the lemma . consider the ratio of determinants in the right hand side of . from the definition of the transition probability and using the linearity of the determinant , we can write it as @xmath582 where we also used . the factors @xmath583 in the above determinants depend only on the row index @xmath584 , therefore they can be extracted from the determinants cancelling each other completely . but it will be convenient to keep these factors inside the determinants . when taking the confluent limit @xmath585 of , certain derivatives appear . this is discussed by tracy  widom @xcite in a similar context . if we apply this idea to , we obtain @xmath586_{a=0}\right)_{i , j } } { \det\left(\frac{\partial^{i-1}}{\partial a^{i-1}}\left [ a^{-\alpha/2 } i_{\alpha}\left(\sqrt{a b_j}\right)\right]_{a=0}\right)_{i , j}}\ ] ] where again @xmath587 run from @xmath167 to @xmath5 .    we evaluate the derivatives . from the series expansion of the modified bessel function , we find @xmath588_{a=0}=\frac{1}{\gamma(i+\alpha)}\left(\frac{b}{4}\right)^{i-1+\alpha/2}.\ ] ] using this for both matrices on the right - hand side of and cancelling again the common factors , we obtain @xmath589 where the second line follows from the evaluation of the vandermonde determinants and cancellation of the common factors recalling . , by inserting into and using the residue theorem , we obtain the double integral formula . this proves the proposition . next we prove theorem  [ thm : hardedgepearcey ] which immediately follows from the following convergence statement for the kernels . let @xmath590 be fixed numbers . it holds that @xmath591 and the convergence is uniform for @xmath132 and @xmath133 from compact subsets of @xmath134 . we use again the method of steep descent analysis . first we substitute the scaled parameter values of the pearcey scaling to and use the definition of the transition kernel to get @xmath592 we perform the change of variables @xmath593 and @xmath594 in the double integral above . since @xmath595 , the following exponential integral representation holds : @xmath596 the advantage of this formula is that it separates the dependence of the integrand of on the variables @xmath210 and @xmath207 , so the @xmath597 limit of the @xmath210-integral and the @xmath207-integral can be taken separately , as we will proceed . at the end , we will redo the operation . in fact , one could perform the asymptotic analysis without using the identity , but then one would have even longer formulas in the proof . using in , we arrive at @xmath598 we start with the asymptotic analysis for the @xmath210-integral . we are to take the limit of @xmath599 where @xmath600 means an error term which is bounded in absolute value by constant times @xmath601 for @xmath5 large enough . the constant may depend on other parameters : @xmath133 and @xmath221 here . the term @xmath602 has a similar meaning in . we use the method of steep descent analysis to determine the limit of as explained also in @xcite and @xcite . note that by the exponential asymptotics @xmath603 as @xmath604 , see e.g.  9.7.1 in @xcite , the bessel function in corresponds to a term of order @xmath605 in the exponent . the first step of the analysis is to consider the function that is multiplied by the highest power of @xmath5 in the exponent @xmath606 we show that @xmath607 is a steep descent path for the function @xmath608 which means that it reaches its maximum at @xmath609 and that @xmath608 is monotone along @xmath607 . it is easily seen by taking the derivative @xmath610 for all @xmath611 .    for a @xmath235 small enough which will be given later , restricting the integral in to @xmath612 gives an error of order @xmath613 as @xmath597 because of the steep descent property of the path . we substitute @xmath467 by its taylor expansion given in , and now we choose @xmath462 small enough such that the error of the taylor expansion is much smaller than the main terms in the exponent . we do a change of variables @xmath614 . as @xmath597 , the error term of the taylor expansion is @xmath602 . in order to compare the integral with and without this error term in the exponent , we use the inequality @xmath477 and we get that the difference of the two integrals is @xmath602 . in the integral without the error term in the exponent , we extend the steepest descent path to @xmath607 again where an error of @xmath613 is made . this procedure of steep descent analysis proves that converges to @xmath615 as @xmath597 . it is also seen from the above argument that the convergence is uniform if @xmath133 is in a compact interval .    with a similar method for the @xmath207-integral , one can show that @xmath616 as @xmath597 uniformly as @xmath132 is in a finite interval where @xmath617 is @xmath618 with opposite orientation . the error terms here have similar meaning as in . in this argument , one has to first specify the contour @xmath619 . before taking the @xmath597 limit , let @xmath619 to be the counterclockwise oriented sector contour @xmath620\}\cup\{e^{{\mathbf{i}}\theta}r,\theta\in[-\pi/2+\varepsilon,\pi/2-\varepsilon]\}\cup\{e^{{\mathbf{i}}(\pi/2-\varepsilon)}t , t\in[r,0]\}$ ] where @xmath621 $ ] refers to opposite orientation . here @xmath622 is a small value and @xmath623 is sufficiently large such that the contour surrounds the point @xmath624 . as @xmath625 , the integrand on the left - hand side of along @xmath626\}$ ] becomes exponentially small in @xmath623 for large enough but fixed @xmath5 , hence we can replace @xmath619 by the contour @xmath627 consisting of two rays one coming from @xmath628 to @xmath1 and one from @xmath1 to @xmath629 . the reason why we use the contour @xmath627 for some small @xmath622 for the analysis instead of @xmath630 is that by deforming the contour on the left - hand side of , we do not have enough control on the decay of the integrand along @xmath631\}$ ] . the path @xmath627 is a steep descent contour for the function @xmath632 , since @xmath633 for all @xmath611 if @xmath622 is small enough . one can perform a similar asymptotic analysis as for the @xmath210-integral using the steep descent contour @xmath627 for the @xmath207-integral . one gets the convergence with @xmath617 replaced by @xmath627 on the right - hand side , but these two integrals are the same due to the decay coming from the quadratic term in the exponent .    to show that the second term in the right - hand side of without the minus sign converges to @xmath634 we use a dominated convergence argument for the @xmath547 integral of . following the lines of the proof of proposition  [ prop : jbound ] , but actually easier as in lemma 4.1 of @xcite , one can give an exponential bound in @xmath547 on the expression in and on the left - hand side of uniformly in @xmath5 . we omit the proof here . this justifies that the second term on the right - hand side of converges to . what is left to show is that can be rewritten as the double integral in . to this end , it is convenient to split the integration contour @xmath618 in into the half - lines @xmath635 and @xmath636 . let us first consider the contribution from @xmath636 . by substituting @xmath207 by @xmath637 and @xmath210 by @xmath638 and taking into account the branches of the powers , we get the contribution @xmath639 where we also used the identity @xmath640 for any @xmath641 for the bessel functions ( see e.g.  9.6.3 in @xcite ) . similarly , the contribution from the half - line @xmath635 in can be written in exactly the same form as , but with @xmath635 replaced by @xmath636 . to get this , we also use the identity @xmath642 for @xmath643 and @xmath644 which is a consequence of the series expansion of the bessel function ( see e.g.  9.1.10 in @xcite ) .    combining these two contributions , we conclude that equals @xmath645 making the change of variables @xmath646 and @xmath647 , we rewrite this as @xmath648 where we recall the definition of the contour @xmath54 . hence we obtain the double integral on the right - hand side of . in this section , we prove proposition  [ thm : pearcey : alternative ] and corollaries  [ cor : kmw ] and [ cor : bk ] . the strategy for proving the proposition is the following . if @xmath7 is an integer , we can give another formulation of the kernel @xmath229 in proposition  [ prop : pearcey : alternative2 ] below via the approach used to derive the hard edge tacnode kernel . then we show in proposition  [ prop : pearcey : rank1 ] that the temperature derivative of the kernel has a rank one structure . finally , in lemma  [ lemma : rank1:eq ] , we see that the functions which appear in the rank one derivative can be written in a simpler form given in . the proof of proposition  [ thm : pearcey : alternative ] is easy after this . let @xmath649 let @xmath650 be the matrix of size @xmath75 with @xmath65th entry @xmath651 where the superscript denotes the @xmath652th derivative . let the row vector @xmath653 of size @xmath4 be defined as before by , and let @xmath654 be the column vector of size @xmath4 with @xmath84th entry @xmath655 then we put @xmath656 note that we use the same symbols @xmath654 and @xmath130 in this section for different functions as in the sections about the hard edge tacnode process , but this should not lead to any confusion . the invertibility of the matrix @xmath650 follows from the identity @xmath657 which is easily proved by induction on @xmath4 , see below . now we are ready to state    [ prop : pearcey : alternative2 ] if @xmath7 is a non - negative integer , then we can write @xmath229 in as @xmath658 where @xmath234 is a clockwise oriented circle surrounding the singularity at @xmath124 and @xmath235 is chosen such that the contour @xmath236 passes to the right of the singularity at @xmath125 and to the right of the contour @xmath234 .    [ prop : pearcey : rank1 ] the derivative of the kernel with respect to @xmath142 is the rank 1 kernel @xmath659 where @xmath660,\\ \label{g2:pea : schur } g_2(t , y ) & = \frac1{2\pi{\mathbf{i}}}\int_{\delta+{\mathbf{i}}{\mathbb{r}}}{\mathrm d}{u}\ , e^{\frac{{u}^2}2-{\sigma}{u}-\frac y{{u}+t}}({u}+t)^{{\alpha}-1 } \left[{u}^{{\alpha}}-{\boldsymbol { \bf u}}{\mathcal d_{\alpha}}^{-1}\left ( b^{(k+{\alpha})}(-{\sigma } ) \right)_{k=0}^{{\alpha}-1}\right].\end{aligned}\ ] ] if @xmath259 then the expressions between square brackets above are understood to be @xmath167 . the proofs of propositions  [ prop : pearcey : alternative2 ] and [ prop : pearcey : rank1 ] are very similar to the ones for the analogous statements for the hard edge tacnode kernel . we leave this to the interested reader . contrary to the tacnode case , it turns out that the above formulas can be considerably simplified . [ lemma : rank1:eq ] the functions given in  can be written in the simplified form @xmath661    an important ingredient for the proof of lemma  [ lemma : rank1:eq ] is the following elementary analysis lemma which we give with a draft of the proof . [ lemma : int : parts : repeated ] for any @xmath662 and any function @xmath148 with sufficient decay at infinity in the direction of the contour @xmath126 , we have @xmath663    this is trivial for @xmath259 . for @xmath155 , it follows by applying integration by parts to the left - hand side of . for higher values of @xmath4 , it follows from a repeated ( @xmath4-fold ) use of integration by parts writing each time @xmath664 and then integrating by parts . this allows to express the left - hand side of as @xmath665 times a sum of several terms . then one notes that these terms cancel out each other completely except for one surviving term which is precisely the right - hand side of . this is obvious for @xmath259 . for general @xmath4 , the first step in the proof is to use lemma  [ lemma : int : parts : repeated ] . then it remains to prove two identities , one for @xmath666 and one for @xmath667 . we start with the identity for @xmath666 which reads @xmath668 to prove this identity , we first derive the formula @xmath669 as follows . note that the derivatives of the function @xmath68 satisfy the recursion relation @xmath670 for all @xmath671 which is closely related to the recursion for the hermite polynomials . by taking the @xmath478 derivative of and using , one gets @xmath672 integration by parts yields the equation @xmath673 which substituted into exactly gives . if we apply one derivation on the left - hand side of to the vector of entries @xmath674 , then by , we get @xmath675 if @xmath177 , then @xmath676 is @xmath1 and follows easily . otherwise ( if @xmath677 ) after substituting into the left - hand side of and cancelling @xmath650 with its inverse , the scalar product of the two vectors is @xmath1 , hence the left - hand side of can be written as @xmath678    at this point we apply the following elementary column operations to the matrix @xmath679 . for @xmath680 we update column @xmath87 by subtracting @xmath142 times column @xmath681 and adding @xmath681 times column @xmath682 . by with @xmath683 , the resulting matrix has zeros in its zeroth row , except for the @xmath684 entry , which is @xmath685 . we should apply the same operations to the row vector @xmath686 , but this vector does not change under the operations . then reduces to @xmath687 where @xmath688 denotes an unimportant column vector . taking into account the sparsity pattern of the involved matrices and vectors and using the block inversion formula , this in turn can be reduced to @xmath689^{-1 } \frac{{\partial}^{{\alpha}-1}}{{\partial}{v}^{{\alpha}-1}}\left(e^{-{v}^2/2+{\sigma}{v}}{l h_{l-1}}\right)_{l=1}^{\alpha-1}.\ ] ] finally , we note that the factors @xmath297 in the matrix rows and @xmath87 in the entries of the column vector cancel each other . so we get @xmath690 but this is just what we get from the identity for @xmath691 in , so the proof of ends by induction yielding the equivalence of @xmath666 in and in . now we turn to the identity for @xmath667 . using lemma  [ lemma : int : parts : repeated ] , we will prove the identity @xmath692 by induction on @xmath4 . for @xmath259 the identity is trivial since both sides are @xmath167 . for @xmath693 , denote the left - hand side of by @xmath694 and the right - hand side by @xmath695 . note that @xmath694 is basically a hermite polynomial and it satisfies the recursion @xmath696 where the prime denotes the @xmath480-derivative . we will show that @xmath695 satisfies the same recursion . the expression @xmath695 has a schur complement form and therefore it can be written as a ratio of two nested determinants : @xmath697 note that the last column of @xmath698 nicely follows the pattern of the previous columns . we note that the entries in the topmost rows of @xmath698 are basically the hermite polynomials and they satisfy the recursion with @xmath683 . we use the following column operations for @xmath698 : for @xmath680 we update column @xmath87 by subtracting @xmath142 times column @xmath681 and adding @xmath681 times column @xmath682 . the resulting matrix has zeros in its zeroth row , except for the @xmath684 entry , which is @xmath685 , and so we can reduce the size of the determinant . this gives us @xmath699 we take out the factor @xmath700 from the @xmath84th row yielding a total contribution @xmath701 next we split the determinant by linearity with respect to its last row , recognizing the smaller matrix @xmath702 in one of the terms , and its @xmath480-derivative @xmath703 in the other term . this yields @xmath704.\ ] ] similar ( but actually simpler ) manipulations yield @xmath705 using these relations in , we find that @xmath695 satisfies the required recursion relation . this ends the proof of the identity , hence the two definitions and of @xmath667 coincide .    by proposition  [ prop : pearcey : rank1 ] , the derivative of the kernel @xmath229 with respect to @xmath530 can be written in a rank one form with @xmath666 and @xmath667 given by . note that the derivative of the right - hand side of is also rank one but with @xmath666 and @xmath667 given by . hence lemma  [ lemma : rank1:eq ] shows that the two derivatives are the same . moreover , both double integrals converge to zero in the @xmath706 limit . this ends the proof . now we turn to the proofs of the two corollaries .    similarly to remark  [ rem : cafasso ] for the tacnode kernel , by introducing the new integration variables @xmath707 and @xmath708 in , the kernel can be written alternatively as @xmath709 where @xmath243 is a clockwise circle in the left half - plane touching zero along the imaginary axis , and @xmath244 is a counterclockwise loop surrounding @xmath243 . it is now straightforward to get . [ remark : kmw : noninteger ] similarly to remark  [ rem : pearcy : nonint ] , also extends to any real @xmath237 . in this case , we take @xmath243 as above but we replace @xmath244 by a counterclockwise contour coming from @xmath710 from above , encircling @xmath243 , and returning to @xmath710 from below . we then take the powers @xmath711 , @xmath712 with a branch cut along the positive half - line , i.e. , with @xmath713 , see also ( * ? ? ? * eq .  ( 1.19 ) ) .    for @xmath246 , we have @xmath714 if @xmath715 or @xmath716 . substituting this , we get the following expression for the kernel @xmath717 to conclude the proof , one has to do the change of variables @xmath718 and @xmath719 and substitute the new variables from along with the derivative @xmath720 , the rest is straightforward . in this section , we prove proposition  [ prop : limschur ] . it relies on lemmas  [ lemma : detairewrite ] , [ lemma : fklasympt ] and [ lemma : detan ] . in order to state them , we introduce the following notation . then we give the proof of the proposition , finally , we prove the lemmas . inspired by the asymptotics of the airy function @xmath721 as @xmath165 , see e.g.  10.4.59 in @xcite , let us denote @xmath722 for @xmath723 . [ lemma : detairewrite ] for any positive integer @xmath5 , one has the equality @xmath724 where @xmath725 is the @xmath87th derivative of @xmath726 .    by the definition @xmath727 of the airy function , all higher derivatives of @xmath728 can be expressed in terms @xmath729 and @xmath730 . by and the asymptotics of @xmath730 , see 10.4.61 in @xcite , it follows that @xmath731 as @xmath165 with @xmath732 where the empty product is defined to be @xmath167 . note that already yields the following    [ lemma : fklasympt ] we have the asymptotics @xmath733 as @xmath165 . the matrix in the lemma has the form @xmath734    [ lemma : detan ] the determinant of the @xmath34 matrix defined by is @xmath735    1 . the asymptotics readily follows if one puts together , and . 2 .   the assertion is a consequence of the first part and the discussion in the paragraph before the statement of the proposition .    by definition , one immediately gets with @xmath736 replaced by @xmath737 in the determinant on the right - hand side . hence one only needs to show that @xmath738    it follows from that the recursion relation @xmath739 holds for any @xmath740 integer . in , @xmath741 is expressed in terms of @xmath726 and its derivative where the coefficient of the latter is @xmath167 . by applying and its derivative on the right - hand side of again , one can express @xmath742 in terms of @xmath726 and its first and second derivative where the coefficient of the second derivative is @xmath167 . in general , one gets by induction that @xmath743 for all @xmath346 non - negative integers where @xmath744 are coefficients that depend on @xmath132 but not on @xmath84 . the equality of the determinants now follows by simple column transformations . we first represent the entries on the left - hand side of as given in . then we perform the following step for @xmath745 : from the @xmath87th column of the matrix , we subtract @xmath746 times the @xmath1 column , @xmath747 times the first column , up to @xmath748 times the @xmath681st column . after these transformations , we recover the right - hand side of , which completes the proof . we apply elementary column operation on the matrix @xmath749 in order to get @xmath1 everywhere in the first row except for the first position as follows . starting from the last column and proceeding from right to left until the second column , we add an appropriate multiple of the previous column to each of the columns such that we get @xmath1 in the top position of the column . with this procedure , we get that @xmath750 where we expanded the determinant along the first row in the last step . by dividing the @xmath84th row on the right - hand side of by @xmath751 , we get that the determinant of the matrix @xmath749 is @xmath752 times a similar @xmath753 determinant as that of @xmath749 . the statement of the lemma now follows by induction . [ prop : jconv ] fix a non - negative integer @xmath84 . then for any @xmath754 , @xmath755 pointwise as @xmath756 .    in the integral representation of @xmath556 in , the coefficient of the leading term in the exponent as @xmath756 is the function @xmath757 for which we introduce the following steep descent path . it follows from the more general derivative calculation in that the circle @xmath758 is a steep descent path for @xmath759 for any @xmath622 . let us choose a small @xmath760 such that @xmath758 replaced in a small neighborhood of @xmath167 by a segment of the path @xmath761 such that they form a closed loop around the origin .    for the contour given above , we perform the usual steep descent analysis , so we first neglect the integral on the right - hand side of except for a small neighborhood of @xmath167 . after the change of variable @xmath762 , we get that @xmath763 to show that the right - hand side of and that of are the same , we apply a dominated convergence argument . with the notation @xmath764 we can write @xmath765 where we used the dominated convergence theorem in the second step above . the integrand can be dominated uniformly in @xmath766 because of the fast decay of @xmath767 along the @xmath768 contour . [ prop : jbound ] fix a non - negative integer @xmath84 . for any @xmath769 , there is a finite constant @xmath54 and thresholds @xmath770 , @xmath771 such that @xmath772 holds uniformly for @xmath773 if @xmath774 . we adapt the method of proof of proposition 5.3 in  @xcite or that of lemma 4.1 in  @xcite to our setting . in order to investigate the large values of @xmath775 , we rescale it as @xmath776 and we define the function @xmath777 which gives the leading contribution in the exponent of the integrand on the right - hand side of if @xmath775 is of order @xmath778 .    for small values of @xmath779 , this function has two critical points at @xmath780 at first order , and we will pass through the larger one . hence define @xmath781 for some small @xmath622 to be chosen later .    on the right - hand side of , we can change the integration contour by the cauchy theorem to a circle of radius @xmath7 . using the definition , this yields @xmath782 the path @xmath783 is steep descent for the function @xmath784 , since @xmath785 which is negative for @xmath786 and @xmath787 . we show that the integral in can be bounded by the value of the integrand at @xmath7 . to this end , define @xmath788 let @xmath789 for some small @xmath235 . by the steep descent property of @xmath783 , the contribution of the integral over @xmath790 in is bounded by @xmath791 where @xmath792 does not depend on @xmath221 .    to bound the integral over @xmath793 , we first observe by series expansion that @xmath794 we also use that for any @xmath795 @xmath796 holds for some @xmath797 and @xmath512 constants by if @xmath775 is larger than some @xmath771 . if @xmath798 is such that @xmath799 , then the stronger bound @xmath800 applies . putting the estimates , and together , we get that the integral on the right - hand side of taken only over @xmath793 can be bounded by @xmath801 } { \mathrm d}\theta e^{-\frac{\alpha-\alpha^{-1}}2 \theta^2 t(1+{\mathcal o}(\theta))}k + \int_{\substack{[-\delta,\delta]\setminus\\ [ -t^{-\frac13},t^{-\frac13 } ] } } { \mathrm d}\theta e^{-\frac{\alpha-\alpha^{-1}}2 \theta^2 t(1+{\mathcal o}(\theta))}kt^{k/3}\bigg]\end{gathered}\ ] ] after the change of variable @xmath798 . for @xmath221 large enough , the error terms in the exponents on the right - hand side of and factor @xmath802 in the second integral can be removed by replacing @xmath803 by @xmath804 . thus we get that @xmath805 by bounding the gaussian integral . the estimate above is the largest if @xmath7 is close to @xmath167 . note that , by and , @xmath806 which is large if @xmath775 is large enough , therefore is at most constant times @xmath807 . the authors thank patrik ferrari for many illuminating discussions at different stages of the project and mattia cafasso for the discussion which led to remark  [ rem : cafasso ] . b.  v.  thanks the invitation of lun zhang and the hospitality of the ku leuven where the present research was initiated . he is grateful for the generous support of the humboldt research fellowship for postdoctoral researchers during his stay at the university of bonn and for the postdoctoral fellowship of the hungarian academy of sciences . his work was partially supported by otka ( hungarian national research fund ) grant k100473 . s.  d.  is a postdoctoral fellow of the fund for scientific research  flanders ( belgium ) .
Abstract:

Answer: a system of non - intersecting squared bessel processes is considered which all start from one point and they all return to another point . under the scaling of the starting and ending points when the macroscopic boundary of the paths touches the hard edge , a limiting critical process is described in the neighbourhood of the touching point which we call the hard edge tacnode process . we derive its correlation kernel in an explicit new form which involves airy type functions and operators that act on the direct sum of @xmath0 and a finite dimensional space . as the starting points of the squared bessel paths are set to @xmath1 , a cusp in the boundary appears . the limiting process is described near the cusp and it is called the hard edge pearcey process . we compute its multi - time correlation kernel which extends the existing formulas for the single - time kernel . our pre - asymptotic correlation kernel involves the ratio of two toeplitz determinants which are rewritten using a borodin  okounkov type formula .