in the recent study on quantum field theories and quantum mechanics , topologically neutral soliton molecules , which are locally composed of ( fractional ) instantons and anti - instantons , have been attracting a great deal of attention in relation to the ir - renormalon@xcite .
imaginary ambiguities arising in amplitudes of such topologically neutral configurations can cancel out those arising in non - borel - summable perturbative series ( ir - renormalon ) in quantum theories under certain conditions on the spacetime manifold @xcite . in field theories on compactified spacetime with a small compact dimension ,
these objects are termed as  bions " @xcite .
it is expected that full semi - classical expansion including perturbative and non - perturbative sectors as bions , which is called  resurgent " expansion @xcite , leads to unambiguous and self - consistent definition of field theories in the same manner as the conjecture in quantum mechanics @xcite .
the resurgence in theoretical physics was at first investigated in the matrix model and topological string theory @xcite .
then , the study on the topic has been extended to abjm theory @xcite , string and supersymmetric gauge theories @xcite , and general quantum systems @xcite . bions and resurgence in non - susy field theories , especially in the low - dimensional models , have been extensively investigated for the @xmath0 model @xcite , the grassmann sigma model @xcite , the principal chiral model @xcite , and the @xmath2 model @xcite . according to these studies , the leading - order renormalon ambiguity @xmath3 arising in non - borel - summable perturbative series , which corresponds to the singularity closest to the origin on the borel plane ,
is compensated by the amplitude of neutral bions . on the other hand ,
it is expected but not verified that the ambiguities corresponding to singularities further from the origin ( @xmath4 , @xmath5 , ... ) are cancelled by amplitudes of bion molecules with more than four instanton constituents .
in the case of quantum mechanics , not only the sector of zero instanton charge but also those of nonzero instanton charge contribute to physical observables such as the energy levels .
the authors in refs .
@xcite investigated quantum mechanics with several types of potential including the sine - gordon type .
they showed that the leading instanton contributions are consistent with the perturbative calculation , and conjectured the explicit equation connecting perturbative and instanton contributions , which they call the generalized quantization condition .
recently the authors in refs .
@xcite adopted the uniform - wkb method based on the boundary condition , which is equivalent to the quantization condition in refs .
@xcite , and pointed out the general relation between perturbative and non - perturbative contributions .
explicit calculations of multi - instanton amplitudes at each configuration level are expected to clarify the structure of resurgence and to verify the conjectured relation between perturbative and non - perturbative contributions @xcite .    in this paper , we focus on a quantum mechanical system with the sine - gordon potential , and we calculate the multi - instanton amplitude by explicitly integrating quasi moduli parameters corresponding to separations of instanton - constituents in a semi - classical limit , in comparison with the uniform wkb calculations @xcite .
we adopt an extension of bogomolnyi  zinn - justin prescription @xcite for multi - instanton configurations with an appropriate subtraction scheme for divergent parts .
we calculate contributions to the energy eigenvalue of the lowest band from each multi - instanton configuration in a semi - classical limit ( @xmath6 ) . for the configurations with only instantons ( anti - instantons ) such as
@xmath7 $ ] , @xmath8 $ ] and @xmath9 $ ] , we have unambiguous results without imaginary parts . here , we have denoted an instanton ( anti - instanton ) as @xmath10 ( @xmath11 ) . for configurations containing both instantons and anti - instantons such as
@xmath12 $ ] , @xmath13 $ ] and @xmath14 $ ] , the results contain ambiguous imaginary parts , which depend on the path of analytic continuation .
these imaginary parts correspond to the large - order behavior of perturbation series around the saddle point without the @xmath12 $ ] pair .
for instance , we show explicitly that the imaginary part of the multi - instanton amplitude @xmath13 $ ] cancels the imaginary part of the borel resummation of the large - order perturbation series around the nontrivial background with a single instanton @xmath15 $ ] . by investigating the uniform - wkb calculations in detail
, we verify that all of our results agree completely with those based on the uniform - wkb calculations up to a four - instanton order .
while the sine - gordon quantum mechanics is worth to study on its own , another strong motivation lies in its close relationship to small circumference limit of the two - dimensional @xmath0 model on @xmath1 ( circumference @xmath16 ) with the @xmath17-symmetric twisted boundary condition @xcite .
however , we observe that some of the field configurations of the @xmath0 model are not faithfully represented by means of the sine - gordon quantum mechanics .
one can derive the sine - gordon quantum mechanics from the two - dimensional @xmath0 model by applying the scherk - schwarz dimensional reduction , which requires a particular dependence of the phases of fields on the coordinate @xmath18 of compactified dimension ( @xmath19 ) .
it is important to realize that only parts of field configuration of @xmath0 model can be consistent with this @xmath18 dependence .
for instance , the bps solution of two fractional instantons is not consistent with the scherk - schwarz reduction , and hence its small circumference limit can not be described by the sine - gordon quantum mechanics . on the other hand , two adjacent instantons in the sine - gordon
quantum mechanics are mutually non - bps , although each individual instanton may be understood as a limit of bps fractional instanton ( with a different @xmath18 dependence ) . even in the instanton and anti - instanton configurations ,
@xmath0 model has a significant difference compared to the sine - gordon quantum mechanics : the phase moduli of the fractional instantons in the @xmath0 model are neglected in the sine - gordon quantum mechanics . for the configuration of a neutral bion composed of a fractional instanton and an anti - fractional instanton , we find that the interaction between them strongly depends on the relative phase of constituents .
we calculate the neutral bion contribution in the @xmath0 model , based on the interaction potential with the quasi moduli parameter corresponding to the relative phase between the fractional instanton and anti - instanton .
we find that this calculation gives a correction factor compared to the neutral bion amplitude obtained in the sine - gordon quantum mechanics @xcite .
this paper is organized as follows . in sec .  [
sec : sg ] , we review instantons and their interactions in the quantum mechanics with sine - gordon potential and the borel summation . in sec .  [
sec : ic ] we calculate amplitudes of multi - instanton configurations in sine - gordon quantum mechanics by integrating out the moduli parameters . in sec .
[ sec : wkb ] we discuss the results from the uniform wkb calculations , and show that they completely agree with the instanton moduli calculations . in sec .
[ sec : cpn ] we discuss the neutral bion contributions in the compactified @xmath0 model based on the interaction potential including the relative phase parameter .
section [ sec : sd ] is devoted to a summary and discussion . in appendix  [ app1 ]
we give some details of four - instanton calculations .
in this article , we focus on the sine - gordon quantum mechanics described by the schrdinger equation @xmath20 where we follow the notation in refs .
@xcite except @xmath21 is replaced by @xmath22 here .
the euclidian lagrangian for the sine - gordon quantum mechanics is given by , and their euclidean lagrangian is @xmath23 . ]
@xmath24 in the @xmath25 limit , it reduces to the schrdinger equation of the harmonic oscillator .
the energy eigenvalues of periodic potentials split into bands of states . within each band , they are labeled by the bloch angle @xmath26 $ ] defined by @xmath27 in this article , we are interested in the lowest band , although excited bands can be treated similarly .
the energy eigenvalue @xmath28 of the lowest band can be expressed in terms of the path - integral @xmath29 for weak coupling , the path - integral has contributions @xmath30 around the perturbative vacuum , as well as contributions @xmath31 from nonperturbative saddle points @xmath32    perturbation series in powers of coupling constant @xmath33 in quantum field theories or in quantum mechanics are extremely useful , but are usually factorially divergent @xmath34 it is useful to define the borel transform @xmath35 @xmath36 the borel resummation @xmath37 of the divergent series @xmath30 is defined as an integral of the borel transform along the positive real axis in the complex borel plane @xmath38 @xmath39    if the factorially divergent series is alternating , the borel transform has no singularities along the positive real @xmath38 axis , and the borel resummation becomes well - defined ( borel - summable ) . for the potential with degenerate minima such as the sine - gordon quantum mechanics , however , the perturbation series is non - alternating factorially divergent . in that case
, the borel transform is convergent with the finite radius of convergence , but the borel resummation is ill - defined because of singularities in the complex borel plane . since the series become alternating and the borel resummation @xmath37 is unambiguous for @xmath40 , we can analytically continue it from @xmath40 to the physical region @xmath41 to obtain a real analytic function @xmath37 .
if there is no complex singularities , we obtain a branch cut along the positive real axis of complex @xmath33 plane .
the imaginary part @xmath42 at @xmath41 is related to the large - order behavior ( @xmath43 ) of perturbation series @xmath30 in eq.([eq : perturbation - series ] ) through the dispersion relation @xcite @xmath44}{(g^2)^{k+1}}. \label{eq : dispersion_rel}\ ] ] this large - order behavior corresponds to the singularities of the borel transform @xmath45 in the complex borel plane @xmath38 .
of course this ambiguous ( path - dependent ) imaginary part is unacceptable , and should disappear , since the energy eigenvalue @xmath28 should be real , and ambiguity due to the choice of path is unphysical .
in fact , it has been found that the leading term of the imaginary ambiguities is cancelled by the contributions from non - perturbative saddle points associated with neutral objects composed of instantons @xcite .
this phenomenon is called the resurgence of the perturbation series .
let us now consider non - perturbative saddle points . by rescaling the variable @xmath46 the euclidean lagrangian in eq.([eq : lagrangian ] ) can be rewritten as @xmath47 and the instanton number as a topological charge may be defined by @xmath48 single instanton solution ( @xmath49 ) is given by . ]
@xmath50 whereas single anti - instanton solution ( @xmath51 ) is given by @xmath52 with the euclidean action @xmath53 the moduli parameter @xmath54 is a zero mode ( moduli ) associated to the breakdown of translation , representing the location of the ( anti-)instanton . for even @xmath55 , the solutions ( [ eq : instanton ] ) and ( [ eq : anti - instanton ] ) satisfy the following bps equation of
the ( anti-)bps equation ( [ eq : bps - bound ] ) and ( [ eq : anti - bps - bound ] ) is twice as large as the periodicity @xmath56 of the lagrangian in eq.([eq : lagrangian ] ) . ]
saturating the bps bound for @xmath57 @xmath58 for odd @xmath55 , they satisfy the anti - bps equation saturating the anti - bps bound for @xmath57 @xmath59    by integrating over the translational zero mode @xmath54 , one finds the contribution @xmath60 of single instanton @xmath15 $ ] to the energy as @xmath61=   - \left({e^{-s_{i}}\over{\sqrt{\pi g^{2}}}}\right)\ , e^{i\theta}\,.\ ] ]    suppose , for instance , we have a bps instanton in eq.([eq : instanton ] ) with @xmath62 and wish to place another instanton or anti - instanton to its right , we are forced to take either instanton with @xmath63 in eq.([eq : instanton ] ) or anti - instanton with @xmath63 in eq.([eq : anti - instanton ] ) , both of which are anti - bps configurations
. therefore two successive ( anti-)instantons are inevitably non - bps .
the energy of the non - bps configuration of two successive instantons should be more than the sum of individual instanton energies .
they are found to repel each other with the potential  @xcite for large separations @xmath64 @xmath65\ , .
\label{eq : interaction - i - i}\ ] ] the non - bps configuration of successive instanton and anti - instantons are found to attract each other with the potential for large separations @xmath64 @xmath66\ , .
\label{eq : interaction - i - antii}\ ] ]    for later convenience , we introduce the uniform - wkb ansatz by following ref .  @xcite . with the coordinate variable @xmath67 in eq.([eq : rescaled_variable ] ) , eq.([seq1 ] ) can be rewritten as @xmath68 we define the potential as @xmath69 . by using the parabolic cylinder function @xmath70 satisfying the differential equation @xmath71 we introduce an ansatz for the wave function@xcite @xmath72 where the parameter @xmath73 is the shift of energy eigenvalue @xmath28 from the ground state energy of the harmonic oscillator ( @xmath74 limit ) .
then the schrdinger equation ( [ seq2 ] ) becomes @xmath75 with @xmath76 . in the @xmath77 limit , eq.([eq : nonlinear_eq_u ] )
just reduces to @xmath78 , whose solution @xmath79 is @xmath80 which gives the zeroth - order argument of the parabolic cylinder function in eq.([dansatz ] ) and solves the schrdinger equation of the harmonic oscillator .
in this section we calculate multi - instanton amplitudes in sine - gordon quantum mechanics .
we need to integrate over the distances @xmath81 between instantons and ( anti-)instantons as quasi - moduli . since the interaction ( [ eq : interaction - i - i ] ) and ( [ eq : interaction - i - antii ] ) between instantons and ( anti-)instantons vanish at large distances , we need to regulate the integral by introducing a factor @xmath82 into the effective potential @xmath83 \,=\ , \pm { 2\over{g^{2}}}\exp(-r ) \,+\ , \epsilon r\ , , \label{pot2}\ ] ] where @xmath84 is for the instanton - instanton repulsive interaction and @xmath85 for the instanton  anti - instanton attractive interaction .
the regularization parameter @xmath82 can be identified as the number @xmath86 of fictitious fermions@xcite . after subtracting divergences
, we need to take the limit @xmath87 .    even after eliminating the divergence arising from large separations ( @xmath88 )
, we have another source of divergence for the case of the attractive instanton  anti - instanton interaction : the integrand @xmath89 $ ] becomes divergent as @xmath90 , contrary to the repulsive case .
therefore the moduli - integral gets divergent contributions from small @xmath81 regions , and is ill - defined in the semi - classical limit ( @xmath6 ) .
this is why we need to introduce the bogomolnyi  zinn - justin ( bzj ) prescription @xcite : we first regard @xmath91 as real positive ( @xmath92 ) to make the integral well - defined in the semi - classical limit , and then we analytically continue @xmath92 back to @xmath93 in the complex @xmath22 plane at the end of the calculation .
the energy eigenvalue of the lowest band has contributions @xmath94 from the amplitude @xmath95_{\rm all}$ ] of @xmath55-instanton and @xmath96-anti - instanton configuration as @xmath97_{\rm all}\,,\ ] ] @xmath98_{\rm all}\,=\ ,    \left({e^{-s_{\mathcal i}}\over{\sqrt{\pi g^{2}}}}\right)^{n+m}\ , e^{i(n - m)\theta}\ ,   \int dr_{1}dr_{2} ... dr_{n+m-1 } \ ,   e^{-v[r_{1}]-v[r_{2}]-\cdot\cdot\cdot v[r_{n+m-1 } ] } \ , , \ ] ] where @xmath99_{\rm all}$ ] stands for the sum of configurations which can be composed of @xmath55 instantons and @xmath96 anti - instantons in all possible orderings . as shown in ref .
@xcite , the contribution contains @xmath100 with @xmath101 being the instanton charge since the bloch angle @xmath102 shows up in a topological term @xmath103 in the euclidian action in eq.([eq : path - integral ] ) .
we perform the quasi - moduli integral taking only interactions between neighboring instantons among the ( @xmath104 ) instantons .
we should perform this multi - integral in the semi - classical region @xmath6 , and subtract the divergent parts appropriately at each level of the multi - integral .
we will evaluate them explicitly from the next subsection .
each horizontal line stands for the vacuum in the sine - gordon potential .
, scaledwidth=40.0% ]    the amplitude of two instantons shown in fig .
[ fig : ii ] is obtained as @xmath105e^{-2i\theta}\xi^{-2}\,=\ , \int_{0}^{\infty}dr
\exp\left(-{2\over{g^2}}e^{- r}-\epsilon r\right ) & \,=\ , \left({g^{2}\over{2}}\right)^{\epsilon } \int_{0}^{2/g^{2}}ds\ , e^{-s } s^{\epsilon-1 }   \nonumber\\[5pt ]   & \,\overset{|g^{2}|\ll 1}{\longrightarrow}\ , \left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon )   \nonumber\\[5pt ]   & \,=\ , -\left(\gamma\,+\,\log{2\over{g^{2}}}\right)\,+\ ,   o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon)\ , , \label{iia}\end{aligned}\ ] ] where @xmath106 is the euler constant and @xmath107 is an instanton factor defined by @xmath108 here we have neglected terms of order @xmath109 or higher . to simplify the formula , we divide the amplitude by @xmath110 and @xmath111 . precisely speaking
, the interaction energy between instantons at small separation @xmath112 may not be precisely represented by the potential in eq.([pot2 ] ) . however , our result is unchanged as long as @xmath6 is satisfied .
we need to subtract the divergent term @xmath113 while the @xmath114 term disappears in the @xmath115 limit .
the contribution from this amplitude to the energy eigenvalue of the lowest band is then given by @xmath116 where the superscript @xmath117 stands for two - instanton and zero  anti - instanton amplitude .
we note that the contribution from the two anti - instanton amplitude is obtained by replacing @xmath111 by @xmath118 .
the amplitude of one instanton and one anti - instanton amplitude is composed of two configurations @xmath12 $ ] and @xmath119 $ ] , as shown in fig .
[ fig : iibar ] . in these cases ,
the interaction between the two constituents is attractive , and the quasi moduli integral is ill - defined .
therefore we introduce the bogomolnyi  zinn - justin ( bzj ) prescription @xcite : we first evaluate the integral by taking @xmath92 , and then we analytically continue the result from @xmath92 back to @xmath93 in the complex @xmath22 plane .
this procedure provides the imaginary ambiguity depending on the path of the analytic continuation as @xmath120 .    , [ \bar{\mathcal{i}}\mathcal{i}]$ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential .
, title="fig:",scaledwidth=40.0% ] , [ \bar{\mathcal{i}}\mathcal{i}]$ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential .
, title="fig:",scaledwidth=40.0% ]    the amplitude of one - instanton and one anti - instanton configuration @xmath12 $ ] corresponding to the left of fig .
[ fig : iibar ] is obtained as @xmath121\xi^{-2}&\,=\ , \int_{0}^{\infty}dr\ , \exp\left(-{2\over{-g^2}}e^{- r}-\epsilon r\right ) \overset{|g^{2}|\ll 1}{\longrightarrow}\ , \left({-g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon )   \nonumber\\[5pt ]   & \,\overset{-g^{2}= e^{\mp i\pi}g^{2}}{\longrightarrow}\ , -\left(\gamma\,+\,\log{2\over{e^{\mp i\pi}g^{2}}}\right )    \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon ) \nonumber\\[5pt ]   & \,=\ , -\left(\gamma\,+\,\log{2\over{g^{2}}}\right ) \mp i\pi     \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon)\ , , \label{eq : i - bari-1}\end{aligned}\ ] ] where we perform the integral in the first line by considering @xmath92 , and in the second line analytically continue @xmath92 back to @xmath93 in the complex @xmath22 plane @xcite .
the third line shows a two - fold ambiguous expression of @xmath91 depending on the path of analytic continuation as @xmath120 . as with the two - instanton case
, we have subtracted the divergent part @xmath113 while the @xmath114 term disappears in the @xmath115 limit .
another amplitude of one - instanton and one anti - instanton configurations @xmath119 $ ] corresponding to the right of fig .
[ fig : iibar ] turns out to give identical contribution as that in eq.([eq : i - bari-1 ] ) .
the total contribution is given by the sum of them with dropping @xmath113 and @xmath114 terms , @xmath122\,+\,[\bar{\mathcal{i}}\mathcal{i } ] ) \xi^{-2}\,=\ , -2\left(\gamma\,+\,\log{2\over{g^{2}}}\right)\,\mp\ ,   2i\pi\,.\ ] ] its contribution to the energy eigenvalue of the lowest band is then given by @xmath123\ , .
\label{iibar}\ ] ]    if the resurgence idea is valid , this imaginary ambiguity should cancel the imaginary ambiguity of the non - borel summable divergent series of the perturbative contribution . @xmath124
+   { \rm im}[e_{\rm pert}]=0 . \label{imag_cancel}\ ] ] if we insert eqs.([iibar ] ) and ( [ imag_cancel ] ) into the dispersion relation in eq.([eq : dispersion_rel ] ) , we should be able to reproduce the large - order behavior of the perturbation series @xmath125\over{(g^2)^{k+1 } } } \,=\ , { 1\over{\pi}}\int_{0}^{\infty } \,d g^2 \ , { { \rm i m } [ \triangle e^{(1,1)}]\over{(g^2)^{k+1 } } } \nonumber\\ & \,=\ , -{1\over{\pi}}\int_{0}^{\infty } \,d(g^2 ) \ , { 2 e^{-1/g^2}\over{(g^2)^{k+2 } } } \,=\ , -{2\over{\pi } } \,k!\ , \quad\quad\quad\quad\quad\ , ( k\geq 2)\,,\end{aligned}\ ] ] in accordance with the leading large - order behavior of the perturbation series @xcite .
thus we find that the imaginary ambiguity of the instanton
anti - instanton amplitude correctly cancels the imaginary ambiguity of the borel resummation of the ( non - borel summable ) perturbation series .
each horizontal line stands for the vacuum in the sine - gordon potential.,scaledwidth=40.0% ]    for the three - instanton amplitude shown in fig .
[ fig : iii ] , we have two quasi modulus @xmath126 corresponding to the separations between adjacent instantons . for multiple moduli integral of each given configuration , we need to specify subtraction scheme explicitly , and propose the following :    1 .
enumerate possible ordering of quasi moduli integrations , such as @xmath127 and @xmath128 for the three instanton case .
2 .   subtract possible poles like @xmath129 for the first integration , and then perform the next integration successively , and retain the finite piece .
3 .   average the results of all possible orderings .
incorporating the repulsive potentials between adjacent instantons , we obtain the three instanton amplitude as ) , but our subtraction scheme specify the successive subtraction which gives a different result from the product of the result of single moduli integral in eq.([iia ] ) . ]
@xmath130e^{-3i\theta}\xi^{-3}&\,=\ , \int_{0}^{\infty}dr_{1}\int_{0}^{\infty}dr_{2 } \exp\left[-{2\over{g^2}}(e^{- r_{1}}+e^{- r_{2}})-\epsilon(r_{1}+r_{2})\right ] \nonumber\\[5pt ]   & \,\overset{|g^{2}|\ll 1}{\longrightarrow}\ , \left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left[\left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon)-{1\over{\epsilon}}\right ]   \nonumber\\[5pt ]   & \,=\ , { 3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 \,+ \ , { \pi^2\over{12 } } \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon)\,.\end{aligned}\ ] ] in the three instanton case , we have two possible orderings of the multi integral , each of which gives the identical contribution .
then the three instanton contribution to the energy eigenvalue of the lowest band is given by @xmath131\ , .
\label{iii}\ ] ] we note that the contribution from the three anti - instanton amplitude is obtained by replacing @xmath132 by @xmath133 .
the two - instanton and one  anti - instanton amplitudes consist of three types of configurations , as shown in fig .
[ fig : iiibar ] .    ,   [ \mathcal{i } \mathcal{i}\bar{\mathcal{i } } ] ,   [ \bar{\mathcal{i}}\mathcal{i } \mathcal{i}]$ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ] + ,   [ \mathcal{i } \mathcal{i}\bar{\mathcal{i } } ] ,   [ \bar{\mathcal{i}}\mathcal{i } \mathcal{i}]$ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ] ,   [ \mathcal{i } \mathcal{i}\bar{\mathcal{i } } ] ,   [ \bar{\mathcal{i}}\mathcal{i } \mathcal{i}]$ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ]    the first one is @xmath134 $ ] , where the anti - instanton is sandwiched between two instantons .
for this type of configuration , adjacent constituent ( anti-)instantons attract each other .
therefore we first take @xmath40 in order to apply the bzj prescription to the integral , and follow our subtraction prescription @xmath135e^{-i\theta}\xi^{-3}&\,=\ , \int_{0}^{\infty}dr_{1}dr_{2 } \exp\left[-{2\over{-g^2}}(e^{- r_{1}}+e^{- r_{2 } } ) -\epsilon(r_{1}+r_{2})\right ] \nonumber\\[5pt ]   & \,\overset{|g^{2}|\ll 1}{\longrightarrow}\ , \left({-g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left[\left({-g^{2}\over{2}}\right)^{\epsilon } \gamma(\epsilon)-{1\over{\epsilon}}\right ]   \nonumber\\[5pt ]   & \,\overset{-g^{2}=g^{2 } e^{\mp i\pi } } { \longrightarrow}\ , { 3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 \ , - \,{17\pi^2\over{12 } }   \pm 3i\pi\left(\gamma+\log{2\over{g^{2}}}\right )   \,+\,o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon)\,,\end{aligned}\ ] ] where we again subtract @xmath129 in the first integral , before integrating the second quasi - moduli .
the second type of configuration is @xmath136 $ ] , where one pair of constituents is repulsive and the other attractive . in order to apply the bzj prescription only to the attractive part of the interaction , we temporarily distinguish two coupling constants @xmath22 for the repulsive interaction and @xmath137 for the attractive interaction , and apply the bzj prescription to @xmath137 , but not to @xmath22 .
after analytic continuation , we identify @xmath137 as the original @xmath22 at the end of the calculation .
the moduli - integral of @xmath136 $ ] is given by @xmath138e^{-i\theta}\xi^{-3}\ , = \
, \int_{0}^{\infty}dr_{1}dr_{2}\ , \exp\left[-{2\over{g^2}}e^{- r_{1}}-{2\over{-\tilde{g}^2 } } e^{- r_{2}}-\epsilon(r_{1}+r_{2})\right]\,.\ ] ] for this configuration @xmath136 $ ] , we have two possible orderings of moduli integral .
the first ordering is to integrate over @xmath139 and to subtract the @xmath129 pole , before integrating over @xmath140 .
we call this ordering as @xmath141 @xmath142   \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi } } { \longrightarrow}\ , { 3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 \ , - \,{5\pi^2\over{12 } }   \pm 2i\pi\left(\gamma+\log{2\over{g^{2}}}\right ) \,+\ , o\left(\frac{1}{\epsilon}\right ) + o(\epsilon ) \,.\end{aligned}\ ] ] we call the result of another ordering as @xmath143 , which is given by integrating over @xmath140 first and then over @xmath139 @xmath144   \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , { 3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 \,+ \,{\pi^2\over{12 } }   \pm i\pi\left(\gamma+\log{2\over{g^{2}}}\right ) \,+\ , o\left(\frac{1}{\epsilon}\right ) \,+\ , o(\epsilon ) .\end{aligned}\ ] ] in the final expressions for @xmath145 , we implicitly put back @xmath137 to @xmath33 after analytic continuation .
we obtain the amplitude of @xmath146 $ ] as an average of @xmath147 and @xmath148 with dropping @xmath113 and @xmath114 terms @xmath149e^{-i\theta}\xi^{-3 } & \,=\ , ( f_{1}\,+\,f_{2})/2 \nonumber\\ & \,=\ ,   { 3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 \ , - \,{\pi^2\over{6 } }   \pm { 3\over{2}}i\pi\left(\gamma+\log{2\over{g^{2}}}\right)\ , .
\label{eq : iibari}\end{aligned}\ ] ] the third type of configuration is @xmath150 $ ] , which gives the identical result as the second type in eq.([eq : iibari ] ) : @xmath151= [ \mathcal{i}\mathcal{i } \bar{\mathcal{i}}]$ ] .    by taking the sum of all three types of configurations
, we end up with @xmath152\ , + \,[\mathcal{i } \mathcal{i } \bar{\mathcal{i } } ] \,+\,[\bar{\mathcal{i}}\mathcal{i }
\mathcal{i}])e^{-i\theta}\xi^{-3 } \,=\ , { 9\over{2}}\left[\left(\gamma+\log{2\over{g^{2}}}\right)^2 \ , - \,{7\pi^2\over{18 } }   \pm { 4\over{3}}i\pi\left(\gamma+\log{2\over{g^{2}}}\right)\right]\,.\end{aligned}\ ] ] its contribution to the energy eigenvalue of the lowest band is then given by @xmath153\ , .
\label{iiibar}\ ] ] we note that the two  anti - instanton and one - instanton contribution is obtained by replacing @xmath154 by @xmath155 .
resurgence implies that the two  anti - instanton and one - instanton amplitude should correspond to the large - order behavior of the perturbation series around the one - instanton saddle point . using the cancellation between imaginary ambiguities of borel resummed perturbation series around one - instanton saddle point and of the two  anti - instanton and one - instanton amplitude ,
the large - order behavior of the perturbation series around the one instanton saddle point can be estimated from the imaginary part of ( [ iiibar ] ) by means of the dispersion relation as @xmath156 \over{(g^2)^{k+1 } } } \,=\ , { 6\over{\pi}}\int_{0}^{\infty } \,dg^2   \ , { e^{-1/g^2}(\gamma+\log(2/g^2))\over{(g^2)^{k+2 } } } \nonumber\\ & \,=\ , { 6\over{\pi } } k!\,\left(\log 2\,+\,{s(k+1,2)\over{k ! } } \right )   \quad\quad\quad\quad\quad\ , ( k\geq 2)\ , , \label{ak3}\end{aligned}\ ] ] where @xmath157 is the stirling number of the first kind , which is the solution of the recurrence relation @xmath158 .
the first few numbers are given as @xmath159 for @xmath160 .
@xcite , the large - order perturbative series around one instanton is numerically calculated as @xmath161 which is consistent with eq .
( [ ak3 ] ) for large @xmath162 .
each horizontal line stands for the vacuum in the sine - gordon potential.,scaledwidth=40.0% ]    for the four - instanton amplitude , the interaction between each pair of instantons is repulsive , as shown in fig .
[ fig : iiii ] .
since all the possible orderings of multi - moduli integral has the identical contribution , we easily obtain the amplitude using the polygamma function defined as @xmath163 @xmath164e^{-4i\theta } \xi^{-4}&\,=\,\int_{0}^{\infty}dr_{1}dr_{2}dr_{3}\ , \exp\left[-{2\over{g^2}}(e^{- r_{1}}+e^{- r_{2}}+e^{-r_{3 } } ) -\epsilon(r_{1}+r_{2}+r_{3})\right ] \nonumber\\[5pt ]   & \,\overset{|g^{2}| \ll 1}{\longrightarrow}\ , \left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left[\left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left\{\left({g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) -{1\over{\epsilon}}\right\ } + { \gamma+\log(2/g^2)\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,=\ , -{8\over{3 } } \left [ \left(\gamma+\log{2\over{g^{2}}}\right)^3 \ , + \,{\pi^2\over{8}}\left ( \gamma+\log{2\over{g^{2}}}\right )   -{1\over{16}}\psi^{(2)}(1 ) \right ] \,+\,o\left(\frac{1}{\epsilon}\right)\,+\ , o(\epsilon)\,,\end{aligned}\ ] ] where we have first subtracted @xmath129 in the first integral @xmath165 , and subtracted @xmath166 in the next integral @xmath167 , and finally dropped the @xmath113 term in the integral @xmath168 .
its contribution to the eigenvalue of the lowest band is then given by @xmath169\ , .
\label{iiii}\ ] ] the four anti - instanton contribution is obtained by replacing @xmath170 by @xmath171 .
the three - instanton and one  anti - instanton amplitudes consist of four types of configurations @xmath172 $ ] , @xmath173 $ ] , @xmath174 $ ] and @xmath175 $ ] , as shown in fig .
[ fig : iiiibar ] .    $ ] , @xmath173 $ ] , @xmath174 $ ] and @xmath175 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ] $ ] , @xmath173 $ ] , @xmath174 $ ] and @xmath175 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ] $ ] , @xmath173 $ ] , @xmath174 $ ] and @xmath175 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ] $ ] , @xmath173 $ ] , @xmath174 $ ] and @xmath175 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=40.0% ]    in the first configuration @xmath172 $ ] , two pairs of adjacent instantons have repulsive interactions , and the other is attractive .
we again apply the bzj prescription only to the integral for the pair with the attractive interaction , where we denote the coupling as @xmath137 , evaluate the integral at @xmath176 , and analytically continue to @xmath177 , while keeping @xmath93 for the repulsive interactions .
the multi - moduli integral is given by @xmath178 \over{e^{2i\theta}\xi^{4}}}\,=\ , \int_{0}^{\infty}dr_{1}dr_{2}dr_{3 } \,\exp\left[-{2\over{g^2}}e^{- r_{1}}-{2\over{g^2 } } e^{- r_{2}}-{2\over{-\tilde{g}^2}}e^{-r_{3 } } -\epsilon(r_{1}+r_{2}+r_{3})\right]\ , .
\label{g123}\ ] ] here we have three orderings of the multi - integral , distinguished by the ordering of @xmath165 ( attractive interaction ) relative to @xmath179 .
we denote the results @xmath180 , @xmath181 , and @xmath182 as the first , the second , and the third integral , respectively . referring the calculations given in appendix  [ app1 ] ,
we just show the result @xmath183   g_{2}(g^2)&\,=\,-{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{1\over{6 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{5\over{2 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{12}}\big]\ , , \nonumber\\[5pt ]   g_{3}(g^2)&\,=\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,- \,{1\over{3 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{3\over{2 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{12}}\big]\,.\end{aligned}\ ] ]    we obtain the amplitude of @xmath184 $ ] as the average of @xmath185 , @xmath181 and @xmath182 , @xmath186\over{e^{2i\theta}\xi^{4 } } } & \,=\ , ( g_{1}(g^2)\,+\,g_{2}(g^2)\,+\ , g_{3}(g^2))/3 \nonumber\\ & \,=\ ,   -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \ , + \,{1\over{3}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{18}}\big]\ , .
\label{eq : iiibari}\end{aligned}\ ] ]    the amplitude of the second configuration @xmath173 $ ] is found to give identical result as the above configuration : @xmath187 = [ \mathcal{i}\mathcal{i}\mathcal{i }   \bar{\mathcal{i}}]$ ] .    in the third configuration @xmath174 $ ] , the two pairs of adjacent instantons have attractive interactions , while the other is repulsive .
we again apply the bzj prescription only to the attractive pairs , by denoting the couplings as @xmath137 .
the multi - integral is given as @xmath188\over { e^{2i\theta}\xi^{4 } } } \,=\,\int_{0}^{\infty}dr_{1}dr_{2}dr_{3}\ , \exp\left[-{2\over{g^2}}e^{- r_{1}}-{2\over{-\tilde{g}^2 } } e^{- r_{2}}-{2\over{-\tilde{g}^2}}e^{-r_{3 } } -\epsilon(r_{1}+r_{2}+r_{3})\right]\ , .
\label{g456}\ ] ]    we have three orderings of the multi - integral , distinguished by the ordering of @xmath168 ( repulsive interaction ) relative to @xmath189 .
we denote the results @xmath190 , @xmath191 , and @xmath192 as the first , the second , and the third integral , respectively . referring the calculations given in appendix  [ app1 ] ,
we just show the result @xmath193\ , , \nonumber\\[5pt ]   g_{5}&\,=\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \ , + \,{19\over{6}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{11\over{2 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2 -{5\pi^2\over{12}}\big]\ , , \nonumber\\[5pt ]   g_{6}&\,=\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{14\over{3 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \,\mp\,i\pi\big[{13\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2   -{11\pi^2\over{12}}\big]\,.\end{aligned}\ ] ]    we obtain the amplitude of @xmath194 $ ] as the average of @xmath195 , @xmath196 and @xmath197 , @xmath198\over{e^{2i\theta}\xi^{4 } } } & \,=\ , ( g_{4}(g^{2})\,+\,g_{5}(g^{2})\,+\ , g_{6}(g^{2}))/3 \nonumber\\ & \,=\ ,   -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,3\pi^{2 } \left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{16\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2   -{7\pi^2\over{18}}\big]\ , .
\label{eq : ibariii}\end{aligned}\ ] ] the amplitude of the fourth configuration @xmath175 $ ] turns out to be identical to the above third configuration : @xmath175 = [ \mathcal{i}\bar{\mathcal{i}}\mathcal{i}\mathcal{i } ] $ ] .
the sum of all the four configurations gives @xmath199 + [ \bar{\mathcal{i}}\mathcal{i}\mathcal{i}\mathcal{i}]+ [ \mathcal{i}\bar{\mathcal{i}}\mathcal{i}\mathcal{i}]+ [ \mathcal{i}\mathcal{i}\bar{\mathcal{i}}\mathcal{i}])e^{-2i\theta}\xi^{-4 } \nonumber\\ & \,=\ , -{32\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{20\over{3 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 2\over{3}}\psi^{(2)}(1)\,\mp\,i\pi\big[16\left(\gamma + \log{2\over{g^{2}}}\right)^2   -{2\pi^2\over{3}}\big]\,.\end{aligned}\ ] ] therefore the contribution of three instantons and one anti - instanton to the energy eigenvalue of the lowest band is given by @xmath200 \right]\ , .
\label{iiiibar}\end{aligned}\ ] ] we note that the three anti - instanton and one - instanton contribution is obtained by replacing @xmath111 by @xmath118 .
resurgence implies that the three - instanton and one  anti - instanton amplitude should correspond to the large - order behavior of the perturbation series around the two - instanton saddle point . using the cancellation between imaginary ambiguities of borel resummed perturbation series around two - instanton saddle point and of the three - instanton and one  anti - instanton amplitude ,
the large - order behavior of the perturbation series around the two instanton saddle point can be estimated from the imaginary part of ( [ iiiibar ] ) by means of the dispersion relation as @xmath201 \over{(g^2)^{k+1 } } } \nonumber\\ & \,=\ , -{16\over{\pi}}\int_{0}^{\infty } \,dg^2 \ , { e^{-1/g^2}\over{(g^2)^{k+2 } } } \left (     ( \gamma+\log 2)^{2}-{\pi^{2}\over{24 } } \ , - \ , 2(\gamma+\log 2)\log(g^2)\,+\ , \log^{2}(g^2 )   \right)\,.\end{aligned}\ ] ] this integral can be performed numerically , and the first few results are @xmath202 for @xmath203 .
one should be able to check this large - order behavior , by performing the perturbation around the two - instanton saddle point to high orders .
the two - instanton and two - anti - instanton amplitudes consist of six types of configurations @xmath204 $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] , as shown in fig .
[ fig : iiibaribar ] .    $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ] $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ] + $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ] $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ] + $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ] $ ] , @xmath205 $ ] , @xmath14 $ ] , @xmath206 $ ] , @xmath207 $ ] and @xmath208 $ ] ) .
each horizontal line stands for the vacuum in the sine - gordon potential.,title="fig:",scaledwidth=32.0% ]    the first one is @xmath204 $ ] , where all three adjacent pairs of constituents have attractive interactions .
thus we apply the simple bzj prescription to all the integration .
since all the orderings of the multi integral have the same contributions , the amplitude can be easily calculated as @xmath209\xi^{-4 } & \,=\ , \int_{0}^{\infty}dr_{1}dr_{2}dr_{3 } \exp\left[-{2\over{-g^2}}(e^{- r_{1}}+e^{- r_{2}}+e^{-r_{3 } } ) -\epsilon(r_{1}+r_{2}+r_{3})\right ] \nonumber\\[5pt ]   & \,\overset{|g^{2}|\ll 1}{\longrightarrow}\ , \left({-g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left[\left({-g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) \left\{\left({-g^{2}\over{2}}\right)^{\epsilon}\gamma(\epsilon ) -{1\over{\epsilon}}\right\ }    + { \gamma+\log(2/g^2)\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-g^{2}\to g^{2}e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{23\over{3 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ &
\qquad\qquad\qquad\mp\,i\pi\big[8\left(\gamma + \log{2\over{g^{2}}}\right)^2 -{7\pi^2\over{3}}\big ] \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,.\end{aligned}\ ] ]    the second configuration @xmath205 $ ] gives identical result as the first one : @xmath205= [ \mathcal{i}\bar{\mathcal{i}}\mathcal{i}\bar{\mathcal{i}}]$ ] .
for the third @xmath14 $ ] and fourth configurations @xmath206 $ ] , the two pairs of interactions are repulsive , but the other pair is attractive . since these moduli - integrals are the same as that of @xmath172 $ ] in eq.([eq : iiibari ] ) , we obtain @xmath210 \xi^{-4 } = [ \bar{\mathcal{i}}\bar{\mathcal{i}}\mathcal{i}\mathcal{i}]\xi^{-4 } \nonumber\\ & \,=\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{1\over{3 } } \pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{18}}\big]\,.\end{aligned}\ ] ]    for the fifth @xmath207 $ ] and the sixth configurations @xmath208 $ ] , the two pairs of interactions are attractive , but the other pair is repulsive .
since these moduli - integrals are the same as @xmath174 $ ] in eq.([eq : ibariii ] ) , we obtain @xmath211\xi^{-4 } = [ \bar{\mathcal{i}}\mathcal{i}\mathcal{i}\bar{\mathcal{i}}]\xi^{-4 } \nonumber\\ & \,=\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,3\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1)\,\mp\,i\pi\big[{16\over{3}}\left(\gamma+\log{2\over{g^{2}}}\right)^2   -{7\pi^2\over{18}}\big]\,.\end{aligned}\ ] ]    the sum of all the six configurations gives @xmath212 + [ \bar{\mathcal{i}}\mathcal{i}\bar{\mathcal{i}}\mathcal{i}]+ [ \mathcal{i}\mathcal{i}\bar{\mathcal{i}}\bar{\mathcal{i}}]+ [ \bar{\mathcal{i}}\bar{\mathcal{i}}\mathcal{i}\mathcal{i}]+ [ \mathcal{i}\bar{\mathcal{i}}\bar{\mathcal{i}}\mathcal{i}]+ [ \bar{\mathcal{i}}\mathcal{i}\mathcal{i}\bar{\mathcal{i}}])\xi^{-4 } \nonumber\\ \,=\ , & -16 \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,22\pi^{2 } \left ( \gamma+\log{2\over{g^{2}}}\right )   + \psi^{(2)}(1)\,\mp\,i\pi\big[32\left(\gamma+\log{2\over{g^{2}}}\right)^2   -{16\pi^2\over{3}}\big]\,.\end{aligned}\ ] ]    therefore the contribution of two instantons and two anti - instantons to the energy eigenvalue of the lowest band is given by    @xmath213 \right]\ , .
\label{iiibaribar}\end{aligned}\ ] ]    according to the resurgence , the imaginary part of the two - instanton and two anti - instanton amplitudes should cancel the imaginary part of the borel resummation of the divergent perturbation series . using the dispersion relation
, we can obtain the large - order behavior of the perturbation series corresponding to the two - instanton and two anti - instanton amplitude in eq .
( [ iiibaribar ] ) @xmath214 \over{(g^2)^{k+1 } } } \nonumber\\ & \,=\ , -{32\over{\pi^2}}\int_{0}^{\infty } \,dg^2   \ , { e^{-2/g^2}\over{(g^2)^{k+3 } } } \left (     ( \gamma+\log 2)^{2}-{\pi^{2}\over{6 } } \ , - \ , 2(\gamma+\log 2)\log(g^2)\,+\ , \log^{2}(g^2 )   \right)\ , .
\label{4iak}\end{aligned}\ ] ] this integral can be performed numerically , and the first few results are @xmath215 for @xmath203 . because of the factor @xmath216 associated with the four - instanton action @xmath217 , the large - order ( @xmath218 ) behavior of @xmath219 in eq .
( [ 4iak ] ) can be estimated with a constant @xmath220 as @xmath221 \,,\ ] ] exhibiting the @xmath222 factor besides the factorial growth @xmath223 .
this large - order behavior corresponds to the singularity of the borel transform @xmath45 at @xmath224 in the borel plane @xmath225 this is consistent with the known and expected results in quantum mechanics .
we comment that we can also calculate the large - order perturbative behavior around the one - instanton and one  anti - instanton vacuum as @xmath226/(g^{2})^{k+1}$ ] .      for general cases with @xmath55 instantons and @xmath96 anti - instantons
, we have @xmath227 configurations depending on the arrangement of constituents .
we can classify these @xmath227 configurations into classes based on how many attractive and repulsive interactions are among the @xmath104 nearest - neighbor interactions .
( we can classify them into @xmath228 classes for @xmath229 , and @xmath230 classes for @xmath231 with @xmath232 . )
we use the expression @xmath233 as the number of configurations including @xmath162 attractive and @xmath234 repulsive interactions satisfying @xmath235 , and show examples of classification for @xmath236 , @xmath237 , @xmath238 and @xmath239 as following , @xmath240 when calculating the amplitude of each in the @xmath233 configurations , we have @xmath241 possible orderings of the multi moduli integrals in relation to the subtraction scheme which is described in sec .
[ sec : general ] and [ sec:3instanton ] .
after calculating these @xmath241 integrals , we average the results and obtain a contribution from one of @xmath233 configurations . as such
, we calculate each contribution , and by summing up all the integral results , we end up with the semi - classical amplitude of the configuration with @xmath55 instantons and @xmath96 anti - instantons .
a systematic method called the uniform - wkb method has been applied extensively to study quantum mechanics for potentials with degenerate minima , such as the sine - gordon quantum mechanics  @xcite . in this method ,
one combines the ansatz in eq.([dansatz ] ) with a global boundary condition to take account of the other minima than the one to do perturbative computation .
this boundary condition enables one to go beyond the ordinary perturbative computation , and obtain all nonperturbative contributions leading to the resurgence . defining the even @xmath242 and odd @xmath243 functions , and the wronskian @xmath244 , @xmath245\ , , \\
f_{2}(y ) \,&=\ , { 1\over{\sqrt{u'(y)}}}\left [ d_{\nu}(u(y)/g)\ , -\ , d_{\nu}(-u(y)/g ) \right]\ , , \\ \mathcal{w}\,&=\ , f_{1}f'_{2}-f'_{1}f_{2}\ , = \,-\sqrt{8\pi\over{g^{2}}}{1\over{\gamma(-\nu)}}\,,\end{aligned}\ ] ] the following uniform - wkb boundary condition @xcite
is imposed at the midpoint @xmath246 of the adjacent minima of the sine - gordon potential @xmath247 where parameter @xmath102 is the bloch angle in eq.([eq : bloch - angle ] ) to label the state within the energy band .
this boundary condition can be rewritten as @xmath248 where @xmath107 is the instanton factor with the instanton action @xmath249 in eq.([eq : instanton_action ] ) @xmath250\,,\ ] ] the function @xmath251 describing the perturbative fluctuations around the instanton is defined in terms of hypergeometric functions in ref .
@xcite . we just show the necessary expression at the zeroth order of @xmath22 ( @xmath252 ) : @xmath253 this boundary condition of the uniform - wkb approximation is also derived as the quantization condition in ref .
@xcite .    since @xmath73 is the shift of energy eigenvalue @xmath28 from the harmonic oscillator ground state energy , the expansion of eq .
( [ uwkbc ] ) up to the @xmath254 order with respect to @xmath255 for the energy eigenvalue of the lowest band is @xmath256 \nonumber\\ & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\pm\,{i\pi\xi^2 \over{2 } } \left [   h_{0}^2 + \nu(2h_{0}h_{0}'+h_{0}^{2}q_{1}^{\pm } ) + \nu^{2}(2(h_{0}')^2 + 2h_{0}h_{0}''-h_{0}^{2}q_{2}^{\pm } ) \right ] \nonumber\\ & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,=\xi\cos\theta \left[h_{0 } + \nu h_{0 } ' + \nu^{2}h_{0}''+\nu^{3}h_{0 } ' '' \right]\ , , \label{expansion}\end{aligned}\ ] ] with @xmath257 with @xmath258 .
for @xmath259 , @xmath260 and @xmath261 , we just replace @xmath262 by @xmath263 .
we iteratively solve the expansion in eq .
( [ expansion ] ) in the semi - classical limit ( @xmath264 ) , then obtain the expression of @xmath255 as @xmath265 \nonumber\\ & \,+\,\xi^{3}h_{0}^{3}\left [ - ( 2q_{1}^2 + q_{2})\cos^{3}\theta \,+   \left ( { \pi^{2}\over{2 } } \mp i\pi{3\over{2}}q_{1 } \right ) \cos \theta \right ] \nonumber\\ & \,+\ , \xi^{4}h_{0}^{4}\big [ \left(5q_{1}^{3}+5q_{1}q_{2}+q_{3 } \right)\cos^{4}\theta \nonumber\\ & \,\,\,\,\,\,\,\,+\ , \left ( -2\pi^{2}q_{1 } \pm i\pi({9\over{2}}q_{1}^{2}+q_{2}-{\pi^{2}\over{4 } } ) \right ) \cos^{2}\theta   \,-\ , \left({\pi^{2}\over{2}}q_{1}\pm i\pi{\pi^{2}\over{4 } } \right ) \big]\ , .
\label{delnu}\end{aligned}\ ] ] for @xmath266 , the lowest energy level @xmath267 is expressed as @xmath268      from eq .
( [ delnu ] ) , we can calculate the contributions from each order of @xmath107 and @xmath154 .
the order of @xmath107 is interpreted as the number of instanton and anti - instanton constituents in the corresponding configurations while the order of @xmath154 is interpreted as the instanton number . in the followings ,
we show that the coefficients at each order of @xmath107 and @xmath154 are completely consistent with the contributions derived from the instanton calculations in the previous section .
we first start with the @xmath110 order , where we have two instanton constituents . in the @xmath110 order ,
the terms proportional to @xmath111 are @xmath269\,.\ ] ] it agrees with the contribution from the two - instanton amplitude in eq .  ( [ ii ] ) . on the other hand ,
the terms proportional to @xmath270 are given by @xmath271\,.\ ] ] these terms composed of the real and imaginary parts precisely agree with the result of the one - instanton and one  anti - instanton amplitude in eq .
( [ iibar ] ) .
we next look into the @xmath272 order . in the @xmath272 order , the terms with @xmath132 are @xmath273\,,\ ] ] which agree with the contributions from the three - instanton amplitude in eq .
( [ iii ] ) .
the terms proportional to @xmath154 are given by @xmath274\,.\ ] ] the real and imaginary parts in this result are identical to the results from the two - instanton and one  anti - instanton amplitudes in eq .
( [ iiibar ] ) .
we finally investigate the @xmath275 order . in this order ,
the terms proportional to @xmath170 are given by @xmath276\,,\ ] ] which agrees with the four - instanton contribution in eq .
( [ iiii ] ) .
next , the terms with @xmath111 are @xmath277 \right]\,.\ ] ] these terms are the same ones as the contribution obtained from the three - instanton and one  anti - instanton amplitudes in eq .
( [ iiiibar ] ) .
finally , the terms proportional to @xmath270 are given by @xmath278 \right]\,,\ ] ] which precisely agree with the contribution from two - instanton and two  anti - instanton amplitudes shown in eq .
( [ iiibaribar ] ) .
in this section , we discuss fractional instantons and neutral bions in the @xmath0 model on @xmath1 with the @xmath279 twisted boundary conditions in comparison to the sine - gordom quantum mechanics .
the @xmath0 model is described in terms of the @xmath280-component complex field @xmath281 satisfying the constraint @xmath282 with a constant @xmath283 , whose euclidean lagrangian is given by @xmath284 ,   \quad { \cal d}_\mu=\partial_\mu+i a_\mu ,   \quad   a_\mu = \frac{i}{v^2 }   { \partial}_\mu h h^{\dagger}.   \label{eq : cpn_lag}\ ] ] the parameter @xmath285 serves as a coupling constant manifold is denoted as @xmath286 , and is related to the notation of dunne and unsal notation @xcite as @xmath287 , @xmath288 .
] of the @xmath0 model , which is asymptotically free .
the bps solutions are obtained if a complex @xmath280-component vector @xmath289 called the moduli matrix is holomorphic ( independent of @xmath290 ) @xcite @xmath291 whereas the moduli matrix @xmath289 for anti - bps solutions should be anti - holomorphic ( @xmath292 independent of @xmath293 ) .
here we have defined the complex coordinate @xmath294 with @xmath295 .
the @xmath17 symmetric twisted boundary condition can be written as @xmath296 the moduli matrix for a fractional instanton as a bps solution is given by @xmath297 for which we find the topological charge to be fractional @xmath298 @xcite ( see also subsequent works @xcite ) .
this solution can be regarded as a kink connecting two neighboring vacua @xcite , where real parameters @xmath299 are their modulus , representing the relative phase @xmath300 of the neighboring vacua and the position of the kink at @xmath301 .
the action and topological charge densities of the single fractional instanton solution has no dependence on the coordinate @xmath18 of the compactified dimension , thus at the first glance the situation seems quite similar to that of the ( euclidean ) quantum mechanics with the periodic potential , namely the sine - gordon quantum mechanics .    to clarify the similarities and differences with the sine - gordon quantum mechanics
, we consider a small compactification length limit ( @xmath302 ) of the two - dimensional @xmath0 model .
let us first take the @xmath303 model for simplicity . by a stereographic projection , we can parametrize @xmath304 target space in terms of two fields @xmath305 and @xmath306 corresponding to the zenith and azimuth angles of @xmath307 @xmath308 one should note that only the ratio of the first and second components @xmath309 is needed to parametrize @xmath307 as the inhomogeneous coordinate of the @xmath303 field space @xmath310 choosing @xmath311 in eq.([eq : cpn_lag ] ) , the lagrangian of the two - dimensional @xmath303 model can be rewritten in terms of @xmath305 and @xmath306 @xmath312 .
\label{eq : cp1_lag}\ ] ] the scherk - schwarz dimensional reduction assumes the following ansatz of a particular @xmath18 dependence for the fields @xmath313 with constants @xmath314 .
one should note that we have restricted @xmath315 to @xmath18 independent field and ignored the fluctuation of field @xmath316 . by inserting the ansatz ( [ eq : ss_ansatz ] ) and integrating over @xmath18
, we obtain the euclidean action @xmath57 as an integral of the euclidean lagrangian @xmath16 over the euclidean time @xmath38 as @xmath317 = \int dt l. \label{eq : ssred_cp1_lag}\ ] ] we notice that the euclidean lagrangian @xmath16 is identical to that of the sine - gordon quantum mechanics in eq.([eq : y_lagrangian ] ) with the identification @xmath318 , @xmath319 , and the sine - gordon coupling @xmath21 as @xmath320 since the @xmath279-twisted boundary condition ( @xmath311 for @xmath303 model ) requires @xmath321 we finally find the coupling of the sine - gordon quantum mechanics @xmath21 in terms of the parameters of the @xmath0 model @xmath322 which is independent of the compactification scale @xmath16 . in fig .
[ fig : sg - cp1 ] we show how the sine - gordon model is embedded into the @xmath323 model .
the target space of the @xmath323 model is @xmath324 in which the @xmath325 as a great circle of @xmath307 is the target space of the sine - gordon model .
the two fixed points of the action of the twisted boundary conditions , the north and south poles denoted by n and s , respectively , corresponding to @xmath326 and @xmath327 ( modulo @xmath328 ) in the lagrangian in eq .
( [ eq : ssred_cp1_lag ] ) , are two vacua in the reduced sine - gordon model ( see fig .  [ fig : sg - cp1 ] ) : @xmath329     model .
the great circle @xmath325 denotes the target space of the sine - gordon model .
n and s are two vacua @xmath330 and @xmath331 modulo @xmath328 of the sine - gordon model . , scaledwidth=20.0% ]    to see the correspondence between instantons in the @xmath323 model and reduced sine - gordon model , let us first examine the bps one - instanton solution in eq.([eq:1instanton ] ) , whose inhomogeneous coordinate is given as @xmath332 @xmath333 we see that this solution is consistent with the assumption of the scherk - schwarz reduction in eq.([eq : ss_ansatz ] ) with eq.([eq : ssphase ] ) .
this solution gives the vertical path starting from n at @xmath334 and reaching s at @xmath335 with @xmath336 when @xmath337 . with varying from @xmath337 to @xmath338
, this vertical path revolves in @xmath316 to end up at @xmath338 with another vertical path connecting n ( at @xmath334 ) to s ( at @xmath335 ) with @xmath339 .
thus the paths sweep precisely half of @xmath307 , as illustrated in fig .
[ sg123 ] .
the solution @xmath315 in eq .
( [ eq : inhomo_1instanton ] ) is identical to the single instanton solution of the sine - gordon quantum mechanics in eq.([eq : instanton ] ) .
therefore the fractional instanton solution of the @xmath303 model is captured correctly by the sine - gordon quantum mechanics as the instanton solution .
the second homotopy group @xmath340 for the @xmath323 fractional instanton is 1/2 , while the first homotopy group @xmath341 for the reduced sine - gordon model is also 1/2 because it corresponds to a half orbit of the great circle .     in the reduced quantum mechanics
is depicted .
it corresponds to a single line from the north to the south pole which is rotated over the half of @xmath342 homogeneously .
the figure depicts the rotation of the line around the half sphere .
the red arrows denote paths depending on @xmath18 with a constant @xmath343 , while the blue arrows denote the @xmath343 dependence of such paths .
, scaledwidth=70.0% ]    in fig .
[ fig : fractional ] , we classify how all possible ( anti-)instanton configurations in the reduced sine - gordon model correspond to the fractional instantons in the original @xmath323 model . as shown in eqs .
( [ eq : bps - bound ] ) and ( [ eq : anti - bps - bound ] ) , a sine - gordon soliton connecting from a vacuum labeled by even @xmath55 to a vacuum labeled by odd @xmath55 is bps while that from odd @xmath55 to even @xmath55 is anti - bps .
when the sine - gordon quantum mechanics is embedded into the @xmath323 model , a soliton from n to s on @xmath307 is bps while that from s to n is anti - bps , consistent with bps and anti - bps kinks in the @xmath323 model .
in other words , configurations with positive values of the second homotopy class @xmath344 are bps while those with negative values are anti - bps both in the @xmath323 model and reduced sine - gordon model .
thus , ( a ) and ( d ) in fig .
[ fig : fractional ] are bps while ( b ) and ( c ) are anti - bps .    [ cols="^,^,^,^,^ " , ]     here we point out that the bps solution of two fractional - instantons in the @xmath303 model can not be described in the reduced sine - gordon model .
the bps solution of two fractional - instantons , which contains the ordinary one - instanton bps solution ( @xmath49 ) in the limit of small separation of two fractional instantons is given by @xmath345 this is a composite of ( a ) and ( d ) in fig.[fig : fractional ] .
the inhomogeneous coordinate of @xmath303 now reads @xmath346 which can not satisfy the assumption ( [ eq : ss_ansatz ] ) of the scherk - schwarz reduction .
this is because in the reduced sine - gordon model a configuration starting from n and ending at s [ ( a ) in fig .
[ fig : fractional ] ] can not be connected to another configuration starting from n and ending at s [ ( d ) in fig .
[ fig : fractional ] ] .
the former can be connected only to a configuration starting from s and ending at n. therefore the bps two fractional instanton solution can not be described by the sine - gordon quantum mechanics even in the limit of small @xmath16 .
more generally , all the bps multi - fractional - instanton solutions are inconsistent with the scherk - schwarz reduction and hence the sine - gordon quantum mechanics fails to capture them .
this is consistent with the fact that configurations containing @xmath55-instantons ( @xmath347 ) are always non - bps in the sine - gordon quantum mechanics .
of neutral bion configurations for @xmath348 and @xmath349 in the @xmath303 model on @xmath1 .
the same action density is depicted in two ways , as a function of @xmath350 ( left ) and @xmath343 ( right ) .
there is no @xmath18 dependence in the action density , with @xmath18 being a coordinate of the compactied dimension.,title="fig:",scaledwidth=55.0% ]   of neutral bion configurations for @xmath348 and @xmath349 in the @xmath303 model on @xmath1 .
the same action density is depicted in two ways , as a function of @xmath350 ( left ) and @xmath343 ( right ) .
there is no @xmath18 dependence in the action density , with @xmath18 being a coordinate of the compactied dimension.,title="fig:",scaledwidth=40.0% ]    next let us consider the non - bps configuration of neutral bion , which is a composite of a fractional instanton and a fractional anti - instanton as depicted in fig .
[ 025pi ] . we can write down an ansatz for the moduli matrix @xmath289 @xmath351 which is guaranteed to become an exact solution of the field equation in the limit of large separation ( @xmath352 ) between constituents .
the inhomogeneous coordinate of @xmath303 now reads @xmath353 which satisfy the assumption ( [ eq : ss_ansatz ] ) of the scherk - schwarz reduction if and only if @xmath354 . in that case
, we obtain the angular coordinate fields of @xmath307 as @xmath355 this configuration starts from n at @xmath334 . for the upper sign
, it goes through s at @xmath356 and reaches to n with @xmath357 at @xmath335 , namely it winds once around the great circle .
the configuration represents the double instanton configuration of the sine - gordon quantum mechanics as shown in fig.[fig : ii ] . for the lower sign
, the configuration returns back to n with @xmath330 at @xmath335 approaching but never reaching s at any point in @xmath358 .
this clearly represents the instanton and anti - instanton configuration @xmath12 $ ] of the sine - gordon quantum mechanics , as shown in the left panel of fig .
[ fig : iibar ] .
the sine - gordon quantum mechanics captures only field configurations that can cover the ( part of ) @xmath307 in the following specific fashion : when @xmath343 is varied with fixed @xmath18 , @xmath315 goes along the great circle ( namely fixed @xmath316 ) , whereas @xmath18 variation with fixed @xmath343 makes a rotation of @xmath316 with the constant velocity by an amount @xmath327 at fixed @xmath315 .
the first homotopy group @xmath341 for the sine - gordon model is one for the upper sign and zero for the lower sign , but the second homotopy group @xmath340 for the @xmath323 model is zero for the both cases . in fig .
[ sgii ] , we show the instanton
anti - instanton and instanton - instanton configurations in the sine - gordon quantum mechanics corresponding to @xmath354 in eq .
( [ eq : ratio_2instanton ] ) , and how the corresponding configuration of the @xmath303 model in eq.([eq : inhomo_ibari ] ) cover the sphere @xmath307 . here
, each of fractional instanton again corresponds to the line between the north and south poles sweeping around the half sphere .    on the other hand , generic configurations of the neutral bion of the @xmath303 model in eq .
( [ eq : ratio_2instanton ] ) exhibit complicated ways of covering ( part of ) @xmath307 in terms of @xmath359 . in terms of the @xmath360 field , however , they are quite similar to those special configurations with @xmath354 describable by the sine - gordon quantum mechanics except for an important new feature : they have the relative phase @xmath361 between fractional instanton and anti - instanton as an additional moduli of the neutral bion configuration . this relative phase moduli introduces a striking physical effect into interactions between fractional instanton and anti - instanton constituents even in the limit of @xmath302 , as we see immediately .
model are depicted.,scaledwidth=80.0% ]    the absence of the relative phase moduli between the instanton constituents in the bion configuration is a crucial drawback of the sine - gordon quantam mechanical description . to clarify how crucial it is ,
we discuss the neutral bion configuration in the @xmath0 model on @xmath1 with the @xmath279 twisted boundary conditions .
generalizing eq.([eq : bion_cp1 ] ) in the @xmath303 model , the neutral bion configuration with the relative phase moduli between the two constituents is given by @xcite @xmath362 where @xmath363 are real parameters corresponding to positions of constituent fractional instanton and fractional anti - instanton , and @xmath364 is the relative phase . here , the separation between the fractional instanton and anti - instanton is given by @xmath365 .
we study the interaction potential as a function of separation @xmath81 of constituents for the fixed relative phase @xmath364 . as we vary this relative phase @xmath364 from @xmath366 to @xmath328 , we find that the interaction between the instanton and anti - instanton constituents is attractive for @xmath367 and @xmath368 , whereas it is repulsive for @xmath369 .
we compare the energy density for various relative phases in fig .
[ thetadep ] .
the interaction potential between fractional instanton and fractional anti - instanton is derived by both numerically @xcite and analytically as @xmath83 \,=\ , -{4\kappa l v^2\cos \phi } \ , \exp(-\kappa r)\ , , \label{potnb}\ ] ] where we note @xmath370 and the coupling of @xmath0 model is denoted as @xmath285 .
we here follow the definition in our previous work @xcite : our coupling @xmath285 is related to the coupling @xmath371 adopted in @xcite as @xmath372 .
we note that the corresponding configuration in the sine - gordon quantum mechanics has no relative phase moduli and the interaction potential between the instanton and anti - instanton constituents is given by eq .
( [ pot2 ] ) with the minus sign ( attraction ) .
therefore the neutral bion contribution in the @xmath0 model can not be calculated in the reduced quantum mechanics correctly .     of neutral bion configurations for @xmath373
are shown for several values of relative phase moduli @xmath364 : @xmath374 ( left ) , @xmath375 ( center ) and @xmath327 ( right ) .
the total action increases as the relative phase @xmath364 is varied from @xmath366 to @xmath327 , while configurations in a far - separated situation are almost the same for the three cases .
the attractive force decreases in the region from @xmath376 to @xmath377 , then becomes repulsive for @xmath378.,title="fig:",scaledwidth=32.0% ]   of neutral bion configurations for @xmath373 are shown for several values of relative phase moduli @xmath364 : @xmath374 ( left ) , @xmath375 ( center ) and @xmath327 ( right ) .
the total action increases as the relative phase @xmath364 is varied from @xmath366 to @xmath327 , while configurations in a far - separated situation are almost the same for the three cases .
the attractive force decreases in the region from @xmath376 to @xmath377 , then becomes repulsive for @xmath378.,title="fig:",scaledwidth=32.0% ]    of neutral bion configurations for @xmath373 are shown for several values of relative phase moduli @xmath364 : @xmath374 ( left ) , @xmath375 ( center ) and @xmath327 ( right ) .
the total action increases as the relative phase @xmath364 is varied from @xmath366 to @xmath327 , while configurations in a far - separated situation are almost the same for the three cases .
the attractive force decreases in the region from @xmath376 to @xmath377 , then becomes repulsive for @xmath378.,title="fig:",scaledwidth=32.0% ]    the neutral bion amplitude has been computed previously by ignoring the relative phase moduli ( @xmath376 ) .
the result can be given in terms of the dimensionless separation variable @xmath379 @xcite @xmath380^{\phi=0 } & = & c e^{-2s_{i}/n}i(v^2 ) ,   \nonumber \\
i(v^2)&=&\int_0^\infty d\tau \exp\left(4\kappa l v^2e^{-\tau } -\epsilon \tau\right ) = \int_0^\infty d\tau \exp\left(\frac{8\pi v^2}{n}e^{-\tau } -\epsilon \tau\right ) \nonumber \\   & = & -\left(\gamma + \log \left(\frac{8\pi v^2}{n}\right)\right ) \mp i\pi ,   \label{eq : bion_no_phase}\end{aligned}\ ] ] where @xmath381 is the instanton action and @xmath220 denotes the numerical coefficient when the relative phase moduli is ignored ( @xmath376 ) .
based on the interaction potential we obtained in eq .
( [ potnb ] ) , we find that the corrected contributions of the neutral bion is given by taking account of the relative phase moduli as @xmath382 = c e^{-2s_{i}/n}\int_0^{2\pi } d\phi \,i(v^2 \cos\phi).\ ] ] since the moduli integrals for @xmath383 $ ] and @xmath384 $ ] are identical , we will double the result for @xmath383 $ ] .
we apply the bzj prescription to the integral in the case of the attractive interaction for @xmath385 $ ] , while we simply perform the integral for the repulsive case @xmath386 $ ] . using eqs.([iia ] ) and ( [ eq : i - bari-1 ] ) , we obtain @xmath387   \nonumber \\ & = -2\pi\left(\gamma + \log \left(\frac{4\pi v^2}{n}\right ) \right ) \,\mp\ , i\pi^{2}. \end{aligned}\ ] ] thus we obtain the neutral bion contribution as @xmath388 \,=\ ,   ce^{-\frac{2s_{i}}{n } } \left[-2\pi\left(\gamma + \log \left(\frac{4\pi v^2}{n}\right ) \right ) \,\mp\ , i\pi^{2 } \right ] \ , , \label{e0nb}\ ] ] in contrast to the result in eq.([eq : bion_no_phase ] ) where the relative phase is ignored . we have obtained quantitative corrections to that in ref .
@xcite , by taking account of the effects of the integral of the relative phase moduli . to establish the absolute magnitude of the instanton contribution definitely , it is desirable to examine the one - loop determinant around the fractional instanton and anti - instanton background for the two - dimensional @xmath0 model which has an explicit ( weak ) dependence on @xmath18 .
we consider these as future works .    before closing this section , we note an alternative possibility to relate the sine - gordon quantum mechanics and the @xmath0 model without compactification .
even if the compactification length is nonzero , the @xmath323 model can be deformed so as to produce the sine - gordon instanton solutions . to clarify this point
, we move to the @xmath389 non - linear sigma model equivalent to the @xmath303 model . in this model
, we have three real scalar fields @xmath390 subjected by the constraint @xmath391 .
we introduce the two potential terms @xmath392 and @xmath393 , with a mass hierarchy @xmath394 . for the parameter region @xmath395 ,
the potential @xmath396 admits two discrete vacua @xmath397 and a @xmath323 domain wall solution interpolating these two vacua @xcite with the width @xmath398 .
let us place it perpendicularly to the @xmath399 direction . with a small @xmath400 ,
the above vacua are shifted and the domain wall is deformed accordingly . in this case , the sine - gordon model is induced on the @xmath323 domain wall as the effective theory @xcite . then , an @xmath389 ( @xmath401 ) instanton is restricted to the @xmath323 domain wall and becomes a sine - gordon instanton with the width @xmath402 ( in the @xmath403 direction ) , in the domain wall effective theory @xcite .
this setting gives a precise relation between the @xmath323 model and the sine - gordon quantum mechanics . by sending @xmath404 to zero , we may be able to investigate how the contribution of the relative complex phase disappears in the reduction process from the @xmath323 model to the sine - gordon quantum mechanics . in the case of the @xmath0 model ,
multiple @xmath405 parallel domain walls in it @xcite play the role to connect the @xmath0 model instantons and instantons in a sine - gordon - like model .
in this paper we have calculated multi - instanton contributions in the quantum mechanics with the sine - gordon potential by integrating out separation moduli parameters between instantons and anti - instantons in the semi - classical limit ( @xmath6 ) . we have adopted an extended bogomolnyi  zinn - justin prescription for multi - instanton configurations and the step - by - step subtraction scheme for divergent parts .
we show that the imaginary parts of the multi - instanton amplitudes cancel those arising from the borel resummation of the large - order perturbation series .
we verify that our results completely agree with those based on the uniform - wkb calculations @xcite up to a four - instanton order .
we have also shown that the neutral bion amplitude in @xmath0 model based on the potential including the relative phase moduli parameter gives corrections to the results obtained in the sine - gordon quantum mechanics @xcite .
our main results , that the multi - instanton amplitudes in the sine - gordon quantum mechanics are consistent with the large - order behavior of the perturbative calculations , and are completely reproduced by the uniform - wkb boundary conditions , strongly indicate the following facts : the uniform - wkb boundary condition provides the correct link between the perturbation series around the perturbative vacuum and non - perturbative instanton effects , and furthermore the perturbative calculation _ knows _ non - perturbative aspects of the quantum system in the first place .
we can take the result as an evidence at least in the quantum mechanics , in favor of the resurgence conjecture , which gives an unambiguous and self - consistent definition of the quantum theory .
as for the neutral bion in the @xmath0 model on @xmath1 , we have calculated , for the first time , the neutral bion amplitude beyond the sine - gordon quantum mechanics , with the moduli integral of the relative phase parameter , which does not exist in the sine - gordon quantum mechanics @xcite . as a future study , we will perform perturbative calculations and multi - instanton calculations in the @xmath0 model with keeping the compactified - direction dependence , then we may be able to show that its large - order behavior of the perturbations originates in the imaginary part of the neutral - bion and bion - molecule amplitude _ at the field theory level_.    we here make a comment on the relation of the configuration eq .
( [ bion ] ) to the lefschetz thimble @xcite . in the recent attempt to perform the path integral in the extended theory @xcite , one complexifies the coupling constant ( @xmath406\not=0 $ ] ) and
the field variables , then the standard path integral can be replaced by the integral along the steepest descent curves in the complex configuration space , or lefschetz thimbles with the imaginary part of the action being constant @xmath407 = { \rm const.}$ ] . in terms of this method , the configuration eq .
( [ bion ] ) corresponds to the special thimble for @xmath406 = 0 $ ] , called the stokes line , where the two critical points , in which fractional instantons are infinitely - separated or completely compressed , are directly connected by the configuration .
although , for a general case @xmath406\not=0 $ ] , the two critical points should belong to two different lefschetz thimbles , our study on the stokes line @xmath406 = 0 $ ] could be a good starting point for investigating the general lefschetz thimbles .
part of future work will be devoted to this study .
we also note that the relation of lefschetz thimbles and the topologically neutral configurations has been investigated in terms of the hidden topological angles in ref .
@xcite .
we also comment that the cancellation of imaginary ambiguities in the sine - gordon quantum mechanics , or the resurgent structure in the system , can be understood as a special case of the generic resummation structure , called  median resummation " @xcite .
it is expected that this resummation is applicable to any quantum theories including quantum mechanics and field theory .
thus , we consider that it is intriguing to study the application of the median resummation to the sine - gordon quantum mechanics and verify its structure as a future work .
we are grateful to mithat nsal and gerald dunne for their interest and valuable comments and correspondences on their related work during the entire course of our study .
t.  m.  and n.  s.  thank cern theory institute 2014 ,  resurgence and transseries in quantum , gauge and string theories " for the fruitful discussion and useful correspondence .
this work is supported in part by the japan society for the promotion of science ( jsps ) grant - in - aid for scientific research ( kakenhi ) grant numbers ( 26800147 ( t.  m. ) , 25400268 ( m.  n. ) and 25400241 ( n.  s. ) ) .
in this appendix we show the details of calculation of the functions @xmath410 ( @xmath411 ) .
the functions @xmath412 are based on the three subtraction patterns of the integral in eq .
( [ g123 ] ) .
they are given by @xmath413    + { \gamma+\log(2/(-\tilde{g}^2))\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{7\over{6}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\,\mp\,4i\pi\left(\gamma+\log{2\over{g^{2}}}\right)^2 \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,,\end{aligned}\ ] ]    @xmath414    + { \gamma+\log(2/g^2)\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{1\over{6}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\,\mp\,i\pi\big[{5\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{12}}\big ] \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,,\end{aligned}\ ] ]    @xmath415    + { \gamma+\log(2/g^2)\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,- \,{1\over{3}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\mp\,i\pi\big[{3\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{12}}\big ] \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,.\end{aligned}\ ] ]      @xmath417    + { \gamma+\log(2/g^2)\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{7\over{6}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\,\mp\,i\pi\big[4\left(\gamma+\log{2\over{g^{2}}}\right)^2 + { \pi^2\over{6}}\big ] \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,,\end{aligned}\ ] ]    @xmath418    + { \gamma+\log(2/(-\tilde{g}^2))\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{19\over{6}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\,\mp\,i\pi\big[{11\over{2 } } \left(\gamma+\log{2\over{g^{2}}}\right)^2   -{5\pi^2\over{12}}\big]\,+\ , o\left(\frac{1}{\epsilon}\right)\ , + \,o(\epsilon)\,,\end{aligned}\ ] ]    @xmath419    + { \gamma+\log(2/(-\tilde{g}^2))\over{\epsilon}}\right ] \nonumber\\[5pt ]   & \,\overset{-\tilde{g}^{2}= \tilde{g}^{2 } e^{\mp i\pi}}{\longrightarrow}\ , -{8\over{3 } } \left(\gamma+\log{2\over{g^{2}}}\right)^3 \,+ \,{14\over{3}}\pi^{2}\left ( \gamma+\log{2\over{g^{2}}}\right )   + { 1\over{6}}\psi^{(2)}(1 ) \nonumber\\ & \qquad\qquad\qquad\,\mp\,i\pi\big[{13\over{2}}\left(\gamma+\log{2\over{g^{2}}}\right)^2 -{11\pi^2\over{12}}\big ] \,+\ , o\left(\frac{1}{\epsilon}\right)\,+\,o(\epsilon)\,.\end{aligned}\ ] ]        m.  nsal , `` abelian duality , confinement , and chiral symmetry breaking in qcd(adj ) , '' phys .
lett .   * 100 * , 032005 ( 2008 ) [ arxiv:0708.1772 [ hep - th ] ] .
m.  nsal , `` magnetic bion condensation : a new mechanism of confinement and mass gap in four dimensions , '' phys .
d * 80 * , 065001 ( 2009 ) [ arxiv:0709.3269 [ hep - th ] ] .
m.  shifman and m.  nsal , `` qcd - like theories on r(3 ) x s(1 ) : a smooth journey from small to large r(s(1 ) ) with double - trace deformations , '' phys .
d * 78 * , 065004 ( 2008 ) [ arxiv:0802.1232 [ hep - th ] ] .
e.  poppitz and m.  nsal , `` conformality or confinement : ( ir)relevance of topological excitations , '' jhep * 0909 * , 050 ( 2009 ) [ arxiv:0906.5156 [ hep - th ] ] . m.  m.  anber and e.  poppitz , `` microscopic structure of magnetic bions , '' jhep * 1106 * , 136 ( 2011 ) [ arxiv:1105.0940 [ hep - th ] ] .
e.  poppitz , t.  schaefer and m.  nsal , `` continuity , deconfinement , and ( super ) yang - mills theory , '' jhep * 1210 * , 115 ( 2012 ) [ arxiv:1205.0290 [ hep - th ] ] .
p.  argyres and m.  nsal , `` a semiclassical realization of infrared renormalons , '' phys .
lett .   * 109 * , 121601 ( 2012 ) [ arxiv:1204.1661 [ hep - th ] ] .
p.  c.  argyres and m.  nsal , `` the semi - classical expansion and resurgence in gauge theories : new perturbative , instanton , bion , and renormalon effects , '' jhep * 1208 * , 063 ( 2012 ) [ arxiv:1206.1890 [ hep - th ] ] . g.  v.  dunne and m.  nsal , `` resurgence and trans - series in quantum field theory : the cp(n-1 ) model , '' jhep * 1211 * , 170 ( 2012 ) [ arxiv:1210.2423 [ hep - th ] ] .
g.  v.  dunne and m.  nsal , `` continuity and resurgence : towards a continuum definition of the cp(n-1 ) model , '' phys .
d * 87 * , 025015 ( 2013 ) [ arxiv:1210.3646 [ hep - th ] ] .
r.  dabrowski and g.  v.  dunne , `` fractionalized non - self - dual solutions in the cp(n-1 ) model , '' phys .
d * 88 * , 025020 ( 2013 ) [ arxiv:1306.0921 [ hep - th ] ] .
g.  v.  dunne and m.  nsal , `` generating non - perturbative physics from perturbation theory , '' phys .
d * 89 * , 041701 ( 2014 ) [ arxiv:1306.4405 [ hep - th ] ] .
a.  cherman , d.  dorigoni , g.  v.  dunne and m.  nsal , `` resurgence in qft : unitons , fractons and renormalons in the principal chiral model , '' phys .
lett .   * 112 * , 021601 ( 2014 ) [ arxiv:1308.0127 [ hep - th ] ] .
g.  basar , g.  v.  dunne and m.  nsal , `` resurgence theory , ghost - instantons , and analytic continuation of path integrals , '' jhep * 1310 * , 041 ( 2013 ) [ arxiv:1308.1108 [ hep - th ] ]
. g.  v.  dunne and m.  nsal , `` uniform wkb , multi - instantons , and resurgent trans - series , '' phys .
d * 89 * , 105009 ( 2014 ) [ arxiv:1401.5202 [ hep - th ] ] .
a.  cherman , d.  dorigoni and m.  nsal , `` decoding perturbation theory using resurgence : stokes phenomena , new saddle points and lefschetz thimbles , '' arxiv:1403.1277 [ hep - th ] .
a.  behtash , t.  sulejmanpasic , t.  schaefer and m.  unsal , `` hidden topological angles and lefschetz thimbles , '' arxiv:1502.06624 [ hep - th ] .
s.  bolognesi and w.  zakrzewski , `` clustering and decomposition for non bps solutions of the @xmath420 models , '' phys .
d * 89 * , 065013 ( 2014 ) [ arxiv:1310.8247 [ hep - th ] ] .
t.  misumi , m.  nitta and n.  sakai , `` neutral bions in the @xmath0 model , '' jhep * 1406 * , 164 ( 2014 ) [ arxiv:1404.7225 [ hep - th ] ] ; t.  misumi , m.  nitta and n.  sakai , `` neutral bions in the @xmath421 model for resurgence , '' j.  phys .
ser .   * 597 * , no . 1 , 012060 ( 2015 ) [ arxiv:1412.0861 [ hep - th ] ] .
t.  misumi and t.  kanazawa , `` adjoint qcd on @xmath422 with twisted fermionic boundary conditions , '' jhep * 1406 * , 181 ( 2014 ) [ arxiv:1405.3113 [ hep - ph ] ] .
t.  misumi , m.  nitta and n.  sakai , `` classifying bions in grassmann sigma models and non - abelian gauge theories by d - branes , '' ptep * 2015 * , 033b02 ( 2015 ) [ arxiv:1409.3444 [ hep - th ] ] .
m.  nitta , `` fractional instantons and bions in the principal chiral model on @xmath423 with twisted boundary conditions , '' arxiv:1503.06336 [ hep - th ]
. m.  nitta , `` fractional instantons and bions in the o@xmath424 model with twisted boundary conditions , '' jhep * 1503 * , 108 ( 2015 ) [ arxiv:1412.7681 [ hep - th ] ] .
g.  v.  dunne and m.  unsal , `` resurgence and dynamics of o(n ) and grassmannian sigma models , '' arxiv:1505.07803 [ hep - th ] .
g.  t hooft , `` can we make sense out of quantum chromodynamics ? , '' subnucl
.  ser .   * 15 * , 943 ( 1979 ) .
v.  a.  fateev , v.  a.  kazakov and p.  b.  wiegmann , `` principal chiral field at large n , '' nucl .
b * 424 * , 505 ( 1994 ) [ hep - th/9403099 ] .
v.  a.  fateev , p.  b.  wiegmann and v.  a.  kazakov , `` large n chiral field in two - dimensions , '' phys .
* 73 * , 1750 ( 1994 ) .
j.  ecalle ,  les fonctions resurgentes " ,
i - iii ( publ . math .
orsay , 1981 ) .
m.  marino , r.  schiappa and m.  weiss , `` nonperturbative effects and the large - order behavior of matrix models and topological strings , '' commun .
* 2 * , 349 ( 2008 ) [ arxiv:0711.1954 [ hep - th ] ] .
r.  couso - santamaria , j.  d.  edelstein , r.  schiappa and m.  vonk , `` resurgent transseries and the holomorphic anomaly : nonperturbative closed strings in local @xmath425 , '' commun .
phys .   * 338 * , no . 1 , 285 ( 2015 ) [ arxiv:1407.4821 [ hep - th ] ] .                  e.  b.  bogomolny , `` calculation of instanton - anti - instanton contributions in quantum mechanics , '' phys .
b * 91 * , 431 ( 1980 ) .
j.  zinn - justin , `` multi - instanton contributions in quantum mechanics , '' nucl .
b * 192 * , 125 ( 1981 ) .
j.  zinn - justin , `` multi - instanton contributions in quantum mechanics . 2 . , '' nucl .
b * 218 * ( 1983 ) 333 .
j.  zinn - justin and u.  d.  jentschura , `` multi - instantons and exact results i : conjectures , wkb expansions , and instanton interactions , '' annals phys .   *
313 * , 197 ( 2004 ) [ quant - ph/0501136 ] .
j.  zinn - justin and u.  d.  jentschura , `` multi - instantons and exact results ii : specific cases , higher - order effects , and numerical calculations , '' annals phys .   * 313 * , 269 ( 2004 ) [ quant - ph/0501137 ] .
u.  d.  jentschura and j.  zinn - justin , `` multi - instantons and exact results .
iv : path integral formalism , '' annals phys .   * 326 * , 2186 ( 2011 ) .
m.  a.  escobar - ruiz , e.  shuryak and a.  v.  turbiner , `` three - loop correction to the instanton density .
i. the quartic double well potential , '' arxiv:1501.03993 [ hep - th ] ; `` three - loop correction to the instanton density .
ii . the sine - gordon potential , '' arxiv:1505.05115 [ hep - th ] .
m.  eto , t.  fujimori , y.  isozumi , m.  nitta , k.  ohashi , k.  ohta and n.  sakai , `` non - abelian vortices on cylinder : duality between vortices and walls , '' phys .
d * 73 * , 085008 ( 2006 ) [ hep - th/0601181 ] ; m.  eto , t.  fujimori , m.  nitta , k.  ohashi , k.  ohta and n.  sakai , `` statistical mechanics of vortices from d - branes and t - duality , '' nucl .
b * 788 * , 120 ( 2008 ) [ hep - th/0703197 ] .
m.  eto , y.  isozumi , m.  nitta , k.  ohashi and n.  sakai , `` solitons in the higgs phase : the moduli matrix approach , '' j.  phys .
a * 39 * , r315 ( 2006 ) [ hep - th/0602170 ] .
f.  bruckmann , `` instanton constituents in the o(3 ) model at finite temperature , '' phys .
lett .   * 100 * , 051602 ( 2008 ) [ arxiv:0707.0775 [ hep - th ] ] ; w.  brendel , f.  bruckmann , l.  janssen , a.  wipf and c.  wozar , `` instanton constituents and fermionic zero modes in twisted cp**n models , '' phys .
b * 676 * , 116 ( 2009 ) [ arxiv:0902.2328 [ hep - th ] ] ; d.  harland , `` kinks , chains , and loop groups in the cp**n sigma models , '' j.  math .
phys .   * 50 * , 122902 ( 2009 ) [ arxiv:0902.2303 [ hep - th ] ] .
f.  bruckmann and t.  sulejmanpasic , `` nonlinear sigma models at nonzero chemical potential : breaking up instantons and the phase diagram , '' arxiv:1408.2229 [ hep - th ] .
e.  r.  c.  abraham and p.  k.  townsend , `` q kinks , '' phys .
b * 291 * , 85 ( 1992 ) ; e.  r.  c.  abraham and p.  k.  townsend , `` more on q kinks : a ( 1 + 1)-dimensional analog of dyons , '' phys .  lett .
b * 295 * , 225 ( 1992 ) ; m.  arai , m.  naganuma , m.  nitta and n.  sakai , `` manifest supersymmetry for bps walls in n=2 nonlinear sigma models , '' nucl .
b * 652 * , 35 ( 2003 ) [ hep - th/0211103 ] ; m.  arai , m.  naganuma , m.  nitta and n.  sakai , `` bps wall in n=2 susy nonlinear sigma model with eguchi - hanson manifold , '' in * arai , a. ( ed . )
et al . : a garden of quanta * 299 - 325 [ hep - th/0302028 ]
. j.  p.  gauntlett , d.  tong and p.  k.  townsend , `` multidomain walls in massive supersymmetric sigma models , '' phys .
d * 64 * , 025010 ( 2001 ) [ hep - th/0012178 ] ; y.  isozumi , m.  nitta , k.  ohashi and n.  sakai , `` construction of non - abelian walls and their complete moduli space , '' phys .
lett .   *
93 * , 161601 ( 2004 ) [ hep - th/0404198 ] ; y.  isozumi , m.  nitta , k.  ohashi and n.  sakai , `` non - abelian walls in supersymmetric gauge theories , '' phys .  rev .
d * 70 * , 125014 ( 2004 ) [ hep - th/0405194 ] ; m.  eto , y.  isozumi , m.  nitta , k.  ohashi , k.  ohta and n.  sakai , `` d - brane construction for non - abelian walls , '' phys .
d * 71 * , 125006 ( 2005 ) [ hep - th/0412024 ] ; m.  eto , y.  isozumi , m.  nitta , k.  ohashi , k.  ohta , n.  sakai and y.  tachikawa , `` global structure of moduli space for bps walls , '' phys .
d * 71 * , 105009 ( 2005 ) [ hep - th/0503033 ] .
m.  nitta , `` josephson vortices and the atiyah - manton construction , '' phys .
d * 86 * , 125004 ( 2012 ) [ arxiv:1207.6958 [ hep - th ] ] ; m.  kobayashi and m.  nitta , `` sine - gordon kinks on a domain wall ring , '' phys .
d * 87 * , no . 8 , 085003 ( 2013 ) [ arxiv:1302.0989 [ hep - th ] ] .
y.  tanizaki , `` lefschetz - thimble techniques for path integral of zero - dimensional @xmath426 sigma models , '' phys .
d * 91 * , no . 3 , 036002 ( 2015 ) [ arxiv:1412.1891 [ hep - th ] ] ; t.  kanazawa and y.  tanizaki , `` structure of lefschetz thimbles in simple fermionic systems , '' jhep * 1503 * , 044 ( 2015 ) [ arxiv:1412.2802 [ hep - th ] ] ; y.  tanizaki , h.  nishimura and k.  kashiwa , `` evading the sign problem in the mean - field approximation through lefschetz - thimble path integral , '' phys .
d * 91 * , no .
10 , 101701 ( 2015 ) [ arxiv:1504.02979 [ hep - th ] ] .