Resurgence of one-point functions in a matrix model for 2D type IIA superstrings

In the previous papers, the authors pointed out correspondence between a supersymmetric double-well matrix model and two-dimensional type IIA superstring theory on a Ramond-Ramond background. This was confirmed by agreement between planar correlation functions in the matrix model and tree-level amplitudes in the superstring theory. Furthermore, in the matrix model we computed one-point functions of single-trace operators to all orders of genus expansion in its double scaling limit, and found that the large-order behavior of this expansion is stringy and not Borel summable. In this paper, we discuss resurgence structure of these one-point functions and see cancellations of ambiguities in their trans-series. More precisely, we compute both series of ambiguities arising in a zero-instanton sector and in a one-instanton sector, and confirm how they cancel each other. In case that the original integration contour is a finite interval not passing through a saddle point, we have to choose an appropriate integration path in order for resurgence to work.


Introduction
Resurgence [1][2][3][4][5][6][7][8] has attracted lots of attention by its intriguing property to make intimate connection between perturbative and nonperturbative quantities. From data of higher-order perturbative expansion, resurgence enables us to extract nonperturbative aspects. It has been investigated how resurgence works in each of various quantum theories, 1 whereas we still do not know much about a unified picture or classification of their resurgence structure. 2 In addition, even in a specific model or theory, we have not clarified how resurgence structure changes depending on its physical observables. Namely, some observables may be strongly affected by nonperturbative effects and resurgence plays an important role in extracting nonperturbative information, while other observables may not and their perturbative series may behave too well to get insight into nonperturbative aspects by resurgent analysis. Furthermore, we have not known much about relation between resurgence structure and physics. In general, nonperturbative dynamics is important to explain how interesting physical phenomena like quark confinement or symmetry breaking take place. When resurgence extracts some information on nonperturbative effects from perturbation theory, we 1 Resurgence structure has been studied in various models and theories based on several motivations: see e.g. in quantum mechanics , string theories [36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52] as well as quantum field theories . 2 Recent progress in this direction have been made, e.g. in [34].

JHEP05(2019)138
expect that some insights into such interesting physics are available. Thus, it is desirable to find various examples where resurgence helps us understand nonperturbative aspects of physics. Based on these motivations, in this paper we study resurgence structure in a supersymmetric double-well matrix model with the action which lead to the nilpotency: Q 2 =Q 2 = {Q,Q} = 0. One of the most interesting features of this model is that the supersymmetry is spontaneously broken by nonperturbative effects in a certain large-N limit called as double scaling limit (defined in (2.2)) [90,91]. We also have proposed in [92] that this model under the double scaling limit gives nonperturbative formulation of type IIA superstring theory in two dimensions on a Ramond-Ramond background. From these, we can regard this model as an invaluable example of spontaneously broken target-space supersymmetry in string theory. In this paper, we concentrate on one-point functions of powers of matrix φ: 1 N tr φ n . In the previous work [93], it is shown that the operators with n even (n ∈ 2N ) are essentially supersymmetric, and 1/N or genus expansions of their one-point functions are polynomials in the parameter µ 2 , terminating at some genus, which do not lead to any nonanalytic behavior in the double scaling limit. In [90], taking account of effects nonperturbative in 1/N , we have calculated the one-point function 1 N tr (φ 2 − µ 2 ) (or equivalently 1 N tr B ) as an order parameter of the spontaneous supersymmetry breaking. On the other hand, the odd-power operators (n ∈ 2N − 1) are not supersymmetric, and genus expansions of their one-point functions exhibit stringy growth of the expansion coefficients as (2h)! as genus h grows. In this paper, we consider the one-point functions of the odd-power operators and study their resurgence structure. Since instantons in the matrix model (1.1) contribute to the order parameter 1 N tr (φ 2 − µ 2 ) and trigger spontaneous supersymmetry breaking [90], our main interest is to clarify how resurgence structure for the odd-power operators is related to such nonperturbative physics.
Another advantage of considering our model (1.1) is that the existence of the Nicolai mapping [90]. Although resurgence requires data of large-order perturbation series, it is not easy to obtain such data in general. In case that a theory is supersymmetric, we may compute its perturbative expansion to all orders, but in turn it may be Borel summable and have trivial resurgence structure. One of nice approaches to overcome this issue is to introduce a parameter explicitly breaking supersymmetry [27][28][29][30][31]. In [96] and this paper, we propose another way to obtain perturbative expansion to all orders for JHEP05(2019)138 resurgence: we consider nonsupersymmetric quantities in a supersymmetric model. In fact, even in calculation of nonsupersymmetric quantities in (1.1), the Nicolai mapping is available reflecting the existence of supersymmetry in the action. This kind of idea will be useful, in particular, in supersymmetric field theories as in [82][83][84][85][86][87][88].
The organization of this paper is as follows. In the next section, we give a brief review of the supersymmetric double-well matrix model. Correlation functions are expressed in terms of eigenvalues of the matrix φ and are defined in each instanton sector. In section 3, we explain how to compute the one-point functions of odd powers of φ by utilizing the Nicolai mapping. In section 4, we consider contribution from the zero-instanton sector to the one-point functions, and find that there exists a series of ambiguities after applying the Borel resummation technique. Then, in section 5 we see that contribution from the one-instanton sector also has another series of ambiguities, and confirm that these series exactly cancel each other at the leading and next-to-leading orders. The last section is devoted to conclusions and discussions.

Review of the supersymmetric matrix model
In this section, we give a brief review of the supersymmetric double-well matrix model defined by the action (1.1), which has been proposed as a nonperturbative formulation of type IIA superstring theory in two dimensions.

Supersymmetry and large-N limit
After integrating out the auxiliary variable B in (1.1), the scalar potential of φ reads (2.1) In the planar limit (N → ∞ with µ 2 fixed) of the matrix model, there are infinitely degenerate supersymmetric vacua parametrized by filling fractions (ν + , ν − ) for µ 2 ≥ 2. The filling fractions represent configurations that ν ± N of the eigenvalues of φ are around the minimum ±|µ| of the double-well potential (2.1) [94,95]. On the other hand, for µ 2 < 2 we have a unique vacuum which breaks the supersymmetry. The boundary µ 2 = 2 is a critical point at which the third-order phase transition occurs. In the planar limit, it is explicitly seen in [92,93] that the result of several types of correlation functions in the matrix model reproduces the tree amplitudes in two-dimensional type IIA superstring theory on a nontrivial Ramond-Ramond background. In addition, we have considered the following double scaling limit [90] that approaches the critical point from the inside of the supersymmetric phase: This limit of the matrix model is expected to provide a nonperturbative formulation of the superstring theory with string coupling constant g s proportional to s − 3 2 . From this viewpoint the planar limit mentioned above is regarded as g s → 0 limit. In fact, in [96] one-point functions for the single-trace operators of powers of φ are explicitly calculated at JHEP05(2019)138 arbitrary genus and found to be finite at each genus under the double scaling limit (2.2). In [90,91], contribution from matrix-model instantons (isolated eigenvalues of φ located around the top of the effective potential) to the free energy is found to be also finite and to have a factor exp (−C/g s ) with a constant C of O(1). This form is typical of solitonic objects in string theory (D-branes). The correspondence between isolated eigenvalues and solitons is also observed in well-established bosonic noncritical string theories [101][102][103][104][105][106][107][108]. The most interesting feature of the model ever found is that these instantons cause spontaneous supersymmetry breaking in the matrix model, which implies violation of target-space supersymmetry induced by nonperturbative effects in the corresponding superstring theory. The aim of this paper is to investigate connection between contribution from higher genus in the absence of the instanton (in the zero-instanton sector) and that from the one-instanton sector through the one-point functions.

Correlation functions in fixed filling fraction
In this subsection, we define correlation functions of our model (1.1) in a fixed filling fraction. First, the partition function is expressed in terms of the eigenvalues of φ as follows [90,93]: where the normalization of the integration measure is fixed as [95], and (x) stands for the Vandermonde determinant for eigenvalues x i (i = 1, · · · , N ): . By dividing the integration region of each λ i according to the filling fraction, the total partition function can be expressed as a sum of each partition function with a fixed filling fraction: (2.5) By changing the integration variables λ j → −λ j (j = ν + N + 1, · · · , N ), it is easy to find that Z (ν + ,ν − ) = (−1) ν − N Z (1,0) and the total partition function vanishes.

JHEP05(2019)138
Next, we define the correlation function of K single-trace operators 1 N tr O a (φ) (a = 1, · · · , K) in the filling fraction (ν + , ν − ) as and express its connected part in the 1/N -expansion: denotes the connected correlation function on a handle-h random surface with the N -dependence factored out; i.e., the quantity of O(N 0 ). Let us consider the case where O a (φ) are polynomials of φ.
. For these operators, multipoint functions at the planar level (h = 0) and higher-genus one-point functions do not exhibit any nonanalytic behavior as s → 0 [93,96], which is characteristic of protection by supersymmetry. On the other hand, operators of odd powers: 1 N tr φ 2k+1 (k ∈ N ∪ {0}) are not invariant under either of Q orQ, and their correlation functions have nontrivial dependence on s [96] as we will mention in the next section. For simplicity, in the following we focus on the one-point function of the odd-power operators: (2.7) with K = 1, O 1 (φ) = φ 2k+1 in the filling fraction 3 (ν + , ν − ) = (1, 0).

Correlation functions in fixed instanton sector
In this subsection, we divide correlation functions in the (1, 0) sector into contributions from definite instanton numbers as done in [90]. In (2.5), the partition function Z (1,0) with the filling fraction (1, 0) is expressed as the integrations of N eigenvalues along the positive real axis. The eigenvalue distribution in the planar limit is given as [93,94] with a = µ 2 − 2 and b = µ 2 + 2, which means that all the eigenvalues are confined in the interval [a, b]. Dividing the integration region of each eigenvalue R + = [0, ∞) into the inside and outside of the support:

9)
3 It is shown in [93] that at least at the planar level (h = 0) and up to the three-point functions (1 ≤ K ≤ 3), it is easy to recover filling fraction dependence of correlation functions from those in (ν+, ν−) = (1, 0).

JHEP05(2019)138
we decompose the partition function as (2.10) Each contribution with fixed p is regarded as the partition function in the p-instanton sector. In fact, an instanton in our model corresponds to a saddle point of effective potential V eff (λ i ) with respect to a single eigenvalue λ i , which is obtained by integrating out all the eigenvalues other than λ i in (2.3). Its saddle point turns out to be the origin λ i = 0 [90]. For large s (small g s ) under the double scaling limit (2.2), the main contribution from the outside of the support R + \ [a, b] is provided by such an instanton located at the origin. Then as mentioned in Introduction, leading contribution from the instanton takes the form of exp (−C/g s ). Correlation functions in the p-instanton sector can also be defined in a similar manner: denotes an expectation value of O within the p-instanton configurations normalized by Z (1,0) p-inst. . According to [90,91], the partition functions behave as in the double scaling limit with s finite but large, whereas O (1,0) p-inst. has no such exponential suppression. Hence (2.11) is a trans-series expanded by the instanton weight  find that for odd n, the Borel resummation of this expansion has a series of ambiguities, but explicitly show that it is indeed canceled by another series of ambiguities arising in the one-instanton sector given in the second line in (2.13) up to the next-to-leading order. This provides a strong support that there is no ambiguity in the trans-series form up to the instanton number one.
In the following, we use the notation: . (2.14)

One-point functions via Nicolai mapping
In this section, following the derivation in [96], we express the one-point function 1 N tr φ 2k+1 (1,0) in terms of quantities in the Gaussian matrix model. Let us first consider the φ 2 -resolvent In terms of the eigenvalues, R 2 (z 2 ) becomes and 1/ (z − λ i ) (1/ (z + λ i )) has poles only on the positive (negative) real axis for the filling fraction (1, 0). This leads to where C 0 is a contour which encloses only the poles at z = λ i for ∀ i counterclockwise. In particular, in the zero-instanton sector, λ i 's are all confined in the interval [a, b], and therefore, where C denotes a contour encircling the interval [a, b] counterclockwise as depicted in figure 1. C does not contain z = 0 inside, and hence contribution from the instanton at the origin is not included in (3.4). Note that the φ 2 -resolvent is mapped to the resolvent in the Gaussian matrix model. In fact, the Nicolai mapping x i = µ 2 − λ 2 i (i = 1, · · · , N ) recasts the partition function Z (1,0) and the one-point function (2.6) with K = 1 in the filling fraction (1, 0) as JHEP05(2019)138 Figure 1. Integration contour C on the complex z-plane.
respectively. Differently from the standard Gaussian matrix model, the integrals of the eigenvalues x i are not over the entire real axis, but are bounded from the above by µ 2 . The superscript (G') indicates a quantity in such a Gaussian matrix model. By introducing an N × N Hermitian matrix M whose eigenvalues are x i (i = 1, · · · , N ), (3.6) can be written as and where R M (x) ≡ 1 N tr 1 x−M and the expectation value · (G') is taken in the Gaussian matrix model (3.5).
It is also useful to express the one-point function by introducing the orthogonal polynomials P n (x) (n = 0, 1, · · · ) associated with the Gaussian matrix model (3.5) [90]: with p (i) n coefficients, which satisfies the orthogonality  12) and Likewise eigenvalue distribution of the Gaussian matrix model From (3.13) and (3.15), we find a formula of the one-point function as an integral of the eigenvalue distribution In [90], the orthogonal polynomials P n (x) are expressed in terms of the orthogonal polynomials P n by taking account of the boundary effect of (3.10) in an iterative manner: (1) n (x) +P (2) n (x) + · · · , S n =S (1) n +S (2) n + · · · , R n =R (1) n +R (2) n + · · · , h n =h (1) n +h (2) n + · · · , (3.18)

JHEP05(2019)138
where the numbers (1), (2), · · · denote the steps of the iteration and turn out to correspond to the instanton numbers of the contributions [90]. Namely, quantities with the number (p) are suppressed by the factor exp −4ps 3/2 /3 as s grows as in (2.12).
Applying this expansion to the eigenvalue distribution (3.15) with we obtainρ where ρ

(G)
M (x) is the eigenvalue distribution of the standard Gaussian matrix model, and the ellipsis represents contribution from higher instanton numbers (p ≥ 2) [90]. Then the one-point function (3.16) is decomposed as From straightforward calculation similar to what is done in section 5 in [90], the second term of the r.h.s. turns out to be a quantity with higher instanton numbers (p ≥ 2), and can be neglected as far as cancellation between the zero-and one-instanton sectors is concerned: By the same reason, we can replace the resolvent in (3.8) by that of the standard Gaussian matrix model (with the superscript (G)):

Ambiguities in the zero-instanton sector
In [96], the all-order result of genus expansion of the one-point functions 1 N tr φ 2k+1 (1,0) is obtained at the zero-instanton sector in the double scaling limit (2.2). In this section, we apply the Borel resummation technique to the result and find that ambiguities arise.

Genus expansion to all orders
For the resolvent of the standard Gaussian matrix model R M (z) (G) , the expression of genus expansion is obtained at arbitrary genus in the literature e.g. [97]. In [96], utilizing

JHEP05(2019)138
the result there to (3.23), we have arrived at the expression in the double scaling limit (2.2), where [x] denotes the greatest integer less than or equal to x. We can see that the infinite series in the bracket on the r.h.s. gives the genus expansion where the power of g 2 s ∝ s −3 counts the number of handles. The suffix "univ." on the l.h.s. means that the most dominant nonanalytic term at s = 0 is taken in the limit (2.2) (the universal part). The overall factor N − 2 3 (k+2) can be absorbed in the "wave function renormalization" of the operator 1 N tr φ 2k+1 .

Borel resummation
The second line on the r.h.s. in (4.1) is a series exhibiting factorial growth as Γ(3h−k−2) h! ∼ (2h)!, which is a characteristic feature of string perturbation series and gives further support that the matrix model describes a string theory in the double scaling limit [98].
The factorial growth means that (4.1) is a divergent series with convergence radius zero. In order to try to make the series well-defined, let us apply the Borel resummation technique to (4.1). It amounts to inserting 1 = 1 Γ (2h + 1) for large h. In addition, a binomial coefficient α h with α / ∈ Z asymptotically behaves

= 1 4π
s k+2 Here the ellipsis stands for integrals containing higher powers of 1 − z 2 z 2 0 and lower genus contributions up to h = 1 3 (k + 2) . The former contributes to a series of ambiguities as well as the first two terms, whereas the latter is a finite sum providing nothing ambiguous. The integrals in (4.6) yield ambiguities due to the cut of the integrands [z 0 , +∞) on the integration contour. As depicted in figure 2, there are two ways to avoid the cut in the integrals of z, but the result will change depending on which choice we take. This means that the divergent series cannot be made well-defined (non-Borel summable) for the zeroinstanton sector alone. However, by taking account of contributions from nonzero instanton numbers, there is a possibility that ambiguity arising there cancels the ambiguity from the zero-instanton sector. Therefore, the theory is free from ambiguity and well-defined as a whole. This is the idea of resurgence. In the following, we will actually see that a series of ambiguities from (4.6) cancels that from one-instanton contribution (up to the next-to-leading order).
In order to identify the precise form of ambiguities, let us give a tiny imaginary part to s: s → s ± i with > 0. For s + i (s − i ), the integration contour of z in (4.6) shifts slightly below (above) the positive real axis to C + (C − ) as in figure 2. Then the ambiguities are given as difference between contribution from s + i (C + ) and that from 4 For m = 0, 1, 2, · · · , k+  s − i (C − ), i.e., integrated discontinuity of the integrands across the cut: (Ambiguities of (4.6)) ≡ (4.6)| s→s+i − (4.6)| s→s−i where the ellipsis in the last line stands for terms of modified Bessel function K k+3+ν (z 0 ) with ν ≥ 2. They are accompanied with z −ν 0 as in the second term and hence suppressed by s −3ν/2 compared to the first term. By using the asymptotic form we finally find (Ambiguities of (4.6)) = i(−1) k Γ k +  Notice that the exponential factor e − 4 3 s 3 2 is characteristic of the one-instanton contribution and its exponent comes from the value of the branch point z 0 . This opens a profound connection between perturbative ambiguities and nonperturbative effects [99]. In addition, ambiguities from higher powers of z 2 z 2 0 − 1 seem to be related to higher corrections by holes and handles which are created by D-branes and closed strings, respectively.

JHEP05(2019)138 5 Ambiguities in the one-instanton sector
In the previous section, we have seen that the one-point functions of the odd-power operators provide divergent string perturbation series in the zero-instanton sector, and explicitly computed a series of ambiguities of its Borel resumed series at the leading and next-toleading orders. Here we find that another series of ambiguities appears in the one-instanton sector, and that these two series with different origins indeed cancel each other.

One-point functions in the one-instanton sector
By using the fact that the eigenvalue distribution of the standard Gaussian matrix model becomes the (diagonal) Airy kernel in the double scaling limit (2.2) [91,100]: we obtain the expression of (3.22) in the double scaling limit as From the relation of the Nicolai mapping 2 + N − 2 3 ξ = x = µ 2 − λ 2 , we see that the region [a, b] for λ in (2.10) is mapped to (−∞, 0] for ξ, while the region [0, a) is mapped to (0, s] in the double scaling limit. (ξ = s corresponds to the location of the instanton λ = 0.) Thus the latter gives contribution from the one-instanton sector:

Saddle point method
We consider contribution to the integral (5.3) from ξ ∼ s 1, which is expected to be comparable with the ambiguities (4.10). From asymptotic behavior of the Airy function, we find that the Airy kernel behaves as where we take account of the expansion up to the next-to-leading order in order to compare the result in the zero-instanton sector (4.10). This leads us to consider the following integral

JHEP05(2019)138
For a while, we do not specify the lower bound of the integral region in (5.5), since our concern is contribution around ξ = s. Let us evaluate the integral (5.5) by a saddle point method. The saddle point around ξ = s is found to be where v(s) represents corrections by subleading contributions. We see that the saddle point corresponding to the instanton ξ * = s slightly shifts due to the presence of the operator depending on k.
In the standard saddle point method, we first find a saddle point and rotate an original integration contour so that it will go along the steepest descent path passing through the saddle point. In our case, the original integration contour is not an infinite line, but a finite interval [0, s], which is major difference from the standard case. The region [0, s] is along the steepest descent path, but it terminates at the branch point ξ = s. The saddle point (5.7) is not in the interval [0, s], but on the branch cut [s, ∞). This can be treated by the shift s → s ± i as in section 4. Here we rotate the contour by the angle π (or −π) around ξ = s. The rotated contour goes in the opposite direction decreasing ξ on the real axis, passes through ξ * and ends at ξ = s. We should choose the π rotation or the −π rotation of the contour of ξ in accordance with s → s + i or s → s − i , respectively.
Here it is worth noticing that what resurgence implies in the present setting. We first note that (3.22) itself is a well-defined real quantity without ambiguity, and we are interested in it after taking the double scaling limit. Because of technical difficulty, without explicit computation of the integral (3.22), we are trying to deduce its expression in the double scaling limit in the form of trans-series. Quantity in each instanton sector would have ambiguities, but resurgence ensures that reflecting the well-definedness of the original expression, such ambiguities are expected to cancel among instanton sectors. Even if the integration region does not contain a saddle point, we should develop perturbation theory around it in order to construct trans-series. As we will see later, information of the original contour can be included in how to rotate the contour and the end point of the rotated contour.
Let us go back to the computation. In calculating f (ξ * ), we should use the shift s → s ± i to the term of ln(s − ξ * ):

JHEP05(2019)138
while ξ-derivatives have no ambiguity: The Taylor expansion of f (ξ) is summed up as Looking at the factor of the Gaussian integral f (ξ * ) in (5.11), we can regard x as a quantity at most O(s − 1 2 ). From this, we find that all the terms in the Taylor expansion are of the same order and should be kept.
By changing the integration variable to We end up with (Ambiguities of (5.3)) ≡ (5.14)| s→s+i − (5.14)| s→s−i 5 We can check that the integrand does not depend on the contour rotation by π or −π. Setting s − ξ = re iθ , θ is supposed to rotate from 0 to ±π in accordance with s → s ± i as mentioned before. Then

JHEP05(2019)138
which precisely cancels the series of ambiguities in the zero-instanton sector (4.10) regarding the leading and next-to-leading contributions. In the above derivation, we identify the origin of ambiguities as the saddle point value of the integrand, more precisely, of the operator 1 N tr φ 2k+1 . In fact, the imaginary ambiguities come from the logarithmic term (5.8) whose origin is the operator (s − ξ) k+ 1 2 after the Nicolai mapping. It is reasonable because the partition function or even-power operators 1 N tr φ 2k do not have any ambiguity and hence the existence of ambiguities must not depend on the eigenvalue distribution or the Airy kernel, but on the kind of operators.
In (5.12), we manage to sum up the Taylor series to obtain the logarithmic term which is a key to quickly derive the Gamma function Γ k + 3 2 in the ambiguities (5.14). If we perform the ordinary saddle point calculation, i.e., Gaussian integral over the whole real axis by bringing down the higher terms of n ≥ 3 in (5.12) from the exponent, the Gamma function will appear in the form of the asymptotic series as k + 1 2 grows. 6 A similar situation was found in [27] where difference between a Gamma function factor and its form of Stirling's formula is supplemented by terms whose order is higher than quadratic. There, it is necessary to go beyond the one-loop determinant since quasi zero-modes appear. In our case, there is no such zero-mode since f (ξ * ) is positive definite, but still all order terms needs to be taken into account in order to confirm resurgence. It would be interesting that in both cases we need all order terms in the saddle point method for different reasons.
In order to derive the trans-series around the instanton saddle for the finite interval, we took a prescription to rotate the integration contour so that it will pass the saddle point with taking care of the direction and the end point of the original contour. We have explicitly seen that this prescription realizes the cancellation not only at the leading order but also at the next-to-leading order. This result supports validity of our prescription. Some powerful technique will be necessary to check the cancellation to all orders or at the level of higher instanton numbers. It is also desirable to give more justification of this prescription from the viewpoint of general theory on resurgence applied to an interval.
Finally, concerning our motivation mentioned in Introduction, we make a comment on a relation to physics, in particular spontaneous supersymmetry breaking. As shown in [90] (eqs. (5.25) and (5.26) there 7 ), its order parameter is given by the Airy kernel (5.1) as where the ellipsis stands for contribution from higher instantons. As explicitly seen from (5.4), this expression can be interpreted as contribution of the instanton, namely an isolated eigenvalue at the top of the potential. It would be important to understand physical aspects of a connection between the ambiguity computed here and the order parameter (5.16).  In this paper, we have investigated resurgence structure in the one-point functions of nonsupersymmetric operators in the supersymmetric double-well matrix model which is proposed as nonperturbative formulation of two-dimensional type IIA superstring theory. In the zero-instanton sector, the Borel resummation technique is applied to a divergent string perturbation series, and a series of ambiguities arises depending on how to avoid the cut on the Borel plane. In the one-instanton sector, a special care is necessary for the integration contour of a finite interval which does not pass through a saddle point. Another series of ambiguities arise from the integrand itself evaluated at the saddle point. We have confirmed that these two series of ambiguities cancel each other both at the leading and next-to-leading order. Our prescription of the integration contour in the one-instanton sector is worth studying further and needs to be understood from the viewpoint of resurgence theory extended to cases of integrals over finite intervals.
Another interesting question is how resurgence structure changes for other correlation functions in the same model. In fact, we have developed a derivation of multi-point functions of odd-power operators in [96] by using the result in the Gaussian matrix model [97]. There, it should be possible to read off large-order behavior of genus expansion for the twopoint (or multi-point) functions of the odd-power operators [109]. It would be interesting to check that it again shows the stringy growth of the expansion coefficient as (2h)! and to find structure of singularities on the Borel plane. On the other hand, in the one-instanton sector, we need to treat off-diagonal components of the Airy kernel K Ai (ξ, η), and it is expected that we find richer resurgence structure with variety of saddle points and steepest descent paths (Lefschetz thimbles).
It is also anticipated that our analysis here can be extended to other models where the Nicolai mapping is available. If they are mapped to the Gaussian matrix model via the Nicolai mapping, correlation functions will be calculated from the result in the Gaussian matrix model as in (3.7). If one is interested in the soft edge scaling limit, studies on the Airy kernel and resurgence structure in this paper will be useful there as well.