Skip to content

Stage C: Instanton Transseries, Large Order, and Ambiguity Cancellation

Page 6 turned two admissible endpoint lines into an exact boundary determinant. When that determinant contains an exponentially small Voros symbol, solving it displaces a perturbative root by an exponentially small amount. Repeating the implicit solution produces an instanton transseries for the spectral parameter.

That formal enlargement is not optional decoration. On a singular Borel ray, the perturbative series has two lateral sums. Their difference is of the same exponential size as a neighboring instanton sector. A physical answer becomes independent of the lateral choice only after the transseries parameter, or more generally the full collection of sector parameters, transforms with the Stokes automorphism.

This page derives that mechanism from the quantization function, relates neighboring sectors to large-order growth, and then explains the double-well selection rule. Complex turning points and genuinely complex spectra are deferred to Page 8.

An exponential in the determinant becomes a transseries variable

Section titled “An exponential in the determinant becomes a transseries variable”

Let the lateral exact quantization function be organized near one perturbative branch as

Q^(E,;q)=r0qrQ^r(E,),\widehat{\mathcal Q} (E,\hbar;q) = \sum_{r\geq0} q^r\widehat{\mathcal Q}_r(E,\hbar),

where

q():=αexp ⁣(A)q^().q(\hbar) := \hbar^\alpha \exp\!\left(-\frac{A}{\hbar}\right) \widehat q(\hbar).

Here AA is the leading action with Re(A/)>0\operatorname{Re}(A/\hbar)>0, while q^\widehat q is a normalized formal fluctuation series with nonzero constant term. In exact WKB, qq usually comes from a small cycle or path Voros symbol after its classical exponential has been separated. The power α\hbar^\alpha records a one-loop or endpoint normalization and must not be discarded.

For the recursion below, qq is temporarily an independent formal variable. If A=A(E)A=A(E) is substituted too early, differentiating exp[A(E)/]\exp[-A(E)/\hbar] produces additional A(E)/A'(E)/\hbar terms. One can instead promote the relevant Voros symbols and logarithmic partners to independent graded variables, perform the implicit recursion, and then restore their declared EE-dependence. If one works directly with q(E)q(E), those derivative terms must be included explicitly.

For example, set q0=q(E0,)q_0=q(E_0,\hbar) and write E=E0+q0e1+q02e2+E=E_0+q_0e_1+q_0^2e_2+\cdots. Then the second coefficient becomes

e2=1Q0[Q2+e1(Q1+Q1Elogq)+12Q0e12]E=E0.\begin{aligned} e_2 = -\frac1{\mathcal Q_0'} \biggl[ &\mathcal Q_2 +e_1 \left( \mathcal Q_1' +\mathcal Q_1\partial_E\log q \right) \\ &+ \frac12\mathcal Q_0''e_1^2 \biggr]_{E=E_0}. \end{aligned}

This moving-action term is one reason a higher instanton sector is not merely a power of the first.

Assume that the perturbative equation

Q^0(E0,)=0\widehat{\mathcal Q}_0(E_0,\hbar)=0

has a simple formal root:

EQ^0(E0,)0.\partial_E\widehat{\mathcal Q}_0(E_0,\hbar) \neq0.

Seek a root in the form

E^=E0+qE1+q2E2+.\widehat E = E_0+qE_1+q^2E_2+\cdots.

Taylor expansion gives, with every quantity on the right evaluated at E=E0E=E_0,

E1=Q^1EQ^0,E_1 = -\frac{\widehat{\mathcal Q}_1} {\partial_E\widehat{\mathcal Q}_0},

and

E2=1EQ^0[12E2Q^0E12+EQ^1E1+Q^2].\begin{aligned} E_2 = -\frac1{\partial_E\widehat{\mathcal Q}_0} \biggl[ &\frac12 \partial_E^2\widehat{\mathcal Q}_0\,E_1^2 \\ &+ \partial_E\widehat{\mathcal Q}_1\,E_1 + \widehat{\mathcal Q}_2 \biggr]. \end{aligned}

Higher sectors follow recursively. Thus the first instanton coefficient is not appended by hand: it is the response of the perturbative root to the first exponentially small term in the boundary determinant.

The simple-root hypothesis is load bearing. Near a level collision, exceptional point, or threshold, EQ^0\partial_E\widehat{\mathcal Q}_0 can vanish. The correct local scale may then be q1/2q^{1/2} or another Puiseux power, and a uniform multiple-root analysis replaces the formulas above.

The elementary model

Q(E,q)=EE+qeλE\mathcal Q(E,q) = E-E_*+q\ee^{\lambda E}

has the exact root

E(q)=E1λW ⁣(λqeλE),E(q) = E_* -\frac1\lambda W\!\left( \lambda q\ee^{\lambda E_*} \right),

where WW is Lambert’s function. Its first terms are

E(q)=EqeλE+λq2e2λE32λ2q3e3λE+O(q4).\begin{aligned} E(q) =E_* &-q\ee^{\lambda E_*} \\ &+\lambda q^2\ee^{2\lambda E_*} \\ &-\frac32\lambda^2q^3 \ee^{3\lambda E_*} +O(q^4). \end{aligned}

This toy determinant will be audited below. It illustrates the implicit algebra only; its qq has not been derived from an ODE and therefore is not, by itself, an exact-WKB theorem.

“Instanton” names three objects that need a dictionary

Section titled ““Instanton” names three objects that need a dictionary”

The same word is often used at three analytic levels.

ObjectDefinitionInformation required
Euclidean instantonA finite-action saddle of a Euclidean path integralAction normalization, endpoints or topological sector, fluctuation determinant, zero and quasi-zero modes
Exact-WKB actionA period or relative period on the spectral double coverSheet, cycle or path, orientation, regularization, and chamber
TransmonomialA formal scale such as αeA/\hbar^\alpha\ee^{-A/\hbar}Action branch, power, logarithms, sector normalization, and lateral prescription

In a standard one-dimensional polynomial oscillator these descriptions can often be matched. The match still has to be demonstrated: a WKB cycle may give 2SI2S_I when the Euclidean convention calls a single crossing SIS_I, and endpoint factors can change the algebraic power of \hbar. For complex saddles there may be no real Euclidean trajectory at all.

The invariant content is the exponential scale and the normalized connection data in the chosen problem. A name such as “one instanton” does not fix either.

The general spectral ansatz needs actions, powers, and logarithms

Section titled “The general spectral ansatz needs actions, powers, and logarithms”

For rr independent small actions, a useful formal template is

E^(σ;)=nNrσnexp ⁣(nA)×nβp=0P(n)(log)pE^n,p().\begin{aligned} \widehat E(\boldsymbol\sigma;\hbar) = \sum_{\boldsymbol n\in\mathbb N^r} &\boldsymbol\sigma^{\boldsymbol n} \exp\!\left( -\frac{\boldsymbol n\mathbin{\cdot}\boldsymbol A}{\hbar} \right) \\ &\times \hbar^{\boldsymbol n\mathbin{\cdot}\boldsymbol\beta} \sum_{p=0}^{P(\boldsymbol n)} (\log\hbar)^p \widehat E_{\boldsymbol n,p}(\hbar). \end{aligned}

Each E^n,p\widehat E_{\boldsymbol n,p} is ordinarily an asymptotic power series. The data in this display have different origins:

  • the actions A\boldsymbol A come from critical paths, cycles, or saddle differences;
  • the powers β\boldsymbol\beta come from Gaussian fluctuations, endpoint normalizations, and zero modes;
  • logarithms arise when action sectors resonate or quasi-zero-mode integrals collide;
  • the allowed multi-indices and parameter powers are restricted by symmetry, topology, and the boundary condition.

Even when only one action AA occurs, the transseries need not reduce to nσnenA/E^n\sum_n\sigma^n\ee^{-nA/\hbar}\widehat E_n. A symmetric double well has parity information in odd sectors and genuine Stokes freedom in neutral even sectors. A periodic potential adds a Bloch angle. Several equal or commensurate actions can force logarithmic sectors and a nontrivial parameter lattice.

The Stokes automorphism acts on the family, not just one sector

Section titled “The Stokes automorphism acts on the family, not just one sector”

For a minimal one-parameter family, write

E^(σ;)=n0σnenA/E^n().\widehat E(\sigma;\hbar) = \sum_{n\geq0} \sigma^n \ee^{-nA/\hbar} \widehat E_n(\hbar).

Assume the action ray is θ=argA\theta=\arg A and the family is closed under the corresponding Stokes automorphism. In the simplest normalization,

SθE^(σ)=E^(σ+S),\mathfrak S_\theta \widehat E(\sigma) = \widehat E(\sigma+\mathsf S),

where S\mathsf S is the normalized Stokes constant. Combining this with the book’s operational convention gives

Sθ+E^(σ)=SθE^(σ+S).\mathcal S_{\theta+}\widehat E(\sigma) = \mathcal S_{\theta-} \widehat E(\sigma+\mathsf S).

To first exponential order, suppose

DiscθanE^0=SeA/SθE^1+O(e2A/).\operatorname{Disc}^{\mathrm{an}}_\theta \widehat E_0 = \mathsf S\, \ee^{-A/\hbar} \mathcal S_{\theta-}\widehat E_1 +O(\ee^{-2A/\hbar}).

Then

Sθ+E^(σ+)=SθE^0+(σ++S)eA/×SθE^1+.\begin{aligned} \mathcal S_{\theta+} \widehat E(\sigma_+) = \mathcal S_{\theta-}\widehat E_0 &+ (\sigma_++\mathsf S) \ee^{-A/\hbar} \\ &\times \mathcal S_{\theta-}\widehat E_1 +\cdots. \end{aligned}

Equality with the lower lateral representative requires

σ=σ++S.\sigma_-=\sigma_++\mathsf S.

The jump belongs to the coordinate σ\sigma, while the resummed spectral value is unchanged.

Median summation uses half the full automorphism

Section titled “Median summation uses half the full automorphism”

The parameter translation makes the balanced choice transparent:

σ+=σS2,σ=σ+S2.\sigma_+ = \sigma-\frac{\mathsf S}{2}, \qquad \sigma_- = \sigma+\frac{\mathsf S}{2}.

For a general closed transseries, the invariant definition is

Sθmed:=Sθ+Sθ1/2=SθSθ1/2.\mathcal S_\theta^{\mathrm{med}} := \mathcal S_{\theta+} \circ\mathfrak S_\theta^{-1/2} = \mathcal S_{\theta-} \circ\mathfrak S_\theta^{1/2}.

The half-powers are formal exponentials of 12logSθ\tfrac12\log\mathfrak S_\theta in the completed action filtration. This construction keeps all coupled sectors and composite alien operations. It is stronger than taking the arithmetic mean of two numbers or the real part of each sector separately.

If the ODE, contour, and boundary conditions are real, complex conjugation interchanges the two laterals, and the balanced parameter is compatible with that conjugation, the median sum is real. None of these hypotheses follows from the word “median.”

The Euler pole shows cancellation with every sign visible

Section titled “The Euler pole shows cancellation with every sign visible”

Return to Page 1’s nonalternating Euler series,

Φ^0():=n=0n!n+1,BΦ^0(ξ)=11ξ.\widehat\Phi_0(\hbar) := \sum_{n=0}^{\infty}n!\hbar^{n+1}, \qquad \mathcal B\widehat\Phi_0(\xi) = \frac1{1-\xi}.

For >0\hbar>0, the book’s upper-minus-lower convention gives

S0+Φ^0=ΦPV+πie1/,S0Φ^0=ΦPVπie1/,\begin{aligned} \mathcal S_{0+}\widehat\Phi_0 &= \Phi_{\mathrm{PV}} +\pi\ii\ee^{-1/\hbar}, \\ \mathcal S_{0-}\widehat\Phi_0 &= \Phi_{\mathrm{PV}} -\pi\ii\ee^{-1/\hbar}, \end{aligned}

where

ΦPV=e1/Ei(1/).\Phi_{\mathrm{PV}} = \ee^{-1/\hbar} \operatorname{Ei}(1/\hbar).

Adjoin the flat solution

Φ^(σ):=Φ^0+σe1/.\widehat\Phi(\sigma) := \widehat\Phi_0+\sigma\ee^{-1/\hbar}.

The Stokes translation is S=2πi\mathsf S=2\pi\ii. Therefore

S0+Φ^(σπi)=ΦPV+σe1/,S0Φ^(σ+πi)=ΦPV+σe1/.\begin{aligned} \mathcal S_{0+} \widehat\Phi(\sigma-\pi\ii) &= \Phi_{\mathrm{PV}} +\sigma\ee^{-1/\hbar}, \\ \mathcal S_{0-} \widehat\Phi(\sigma+\pi\ii) &= \Phi_{\mathrm{PV}} +\sigma\ee^{-1/\hbar}. \end{aligned}

For real σ\sigma, the common value is real. The cancellation is not a mysterious disappearance of imaginary parts: the perturbative lateral ambiguity and the jump of the flat-sector coefficient are the same Stokes datum written in two coordinates.

This is an exact linear calibration, not a spectral model. In an oscillator, the neighboring sector has its own divergent fluctuation series, and cancellation continues recursively through higher action grades.

Large order reads the first loops around a neighboring sector

Section titled “Large order reads the first loops around a neighboring sector”

Write the perturbative tail in the shifted convention of Page 1 as

E^0()=e+g0eg(0)g+1.\widehat E_0(\hbar) = e_* + \sum_{g\geq0} e_g^{(0)}\hbar^{g+1}.

Suppose the nearest relevant resurgent relation on the chosen determination is

DiscθanE^0=S01eA/βSθE^1()+larger actions,\operatorname{Disc}^{\mathrm{an}}_\theta \widehat E_0 = \mathsf S_{01} \ee^{-A/\hbar} \hbar^\beta \mathcal S_{\theta-} \widehat E_1(\hbar) +\text{larger actions},

with

E^1()=k0ek(1)k.\widehat E_1(\hbar) = \sum_{k\geq0}e_k^{(1)}\hbar^k.

Under the required continuation, growth, and dominance assumptions, Cauchy’s coefficient formula deformed onto the cut at AA gives

eg(0)S012πik=0K1ek(1)Γ(g+1βk)Ag+1βk,\begin{aligned} e_g^{(0)} \sim \frac{\mathsf S_{01}}{2\pi\ii} \sum_{k=0}^{K-1} e_k^{(1)} \frac{ \Gamma(g+1-\beta-k) }{ A^{g+1-\beta-k} }, \end{aligned}

up to the retained loop order and contributions from other singularities. The leading term is

eg(0)S01e0(1)2πiΓ(g+1β)Ag+1β.e_g^{(0)} \sim \frac{\mathsf S_{01}e_0^{(1)}}{2\pi\ii} \frac{\Gamma(g+1-\beta)} {A^{g+1-\beta}}.

For the Euler calibration, A=1A=1, β=0\beta=0, S01=2πi\mathsf S_{01}=2\pi\ii, and e0(1)=1e_0^{(1)}=1. The formula gives eg(0)=Γ(g+1)=g!e_g^{(0)}=\Gamma(g+1)=g! exactly.

Thus the perturbative coefficients know the action AA, the exponent β\beta, and the products S01ek(1)\mathsf S_{01}e_k^{(1)}. They do not separate a Stokes constant from the arbitrary normalization of E^1\widehat E_1.

If one positive action dominates and the leading amplitude is nonzero,

Ag:=(g+1β)eg(0)eg+1(0)A.A_g := (g+1-\beta) \frac{e_g^{(0)}}{e_{g+1}^{(0)}} \longrightarrow A.

After estimating AA and β\beta,

Cg:=Ag+1βeg(0)Γ(g+1β)S01e0(1)2πi.C_g := \frac{ A^{g+1-\beta}e_g^{(0)} }{ \Gamma(g+1-\beta) } \longrightarrow \frac{\mathsf S_{01}e_0^{(1)}}{2\pi\ii}.

These limits can be accelerated, but acceleration cannot repair a wrong singularity model. Equal-modulus actions must be summed before taking a ratio. A conjugate pair produces oscillations. A negative action produces alternating signs. A vanishing leading Stokes-weighted coefficient exposes a more distant singularity or a subleading loop.

For several isolated actions the leading expansion is a sum,

eg(0)jS0j2πik0ek(j)Γ(g+1βjk)Ajg+1βjk.e_g^{(0)} \sim \sum_j \frac{\mathsf S_{0j}}{2\pi\ii} \sum_{k\geq0} e_k^{(j)} \frac{ \Gamma(g+1-\beta_j-k) }{ A_j^{g+1-\beta_j-k} }.

The closest singularity controls the exponential envelope in gg, not necessarily every coefficient. More distant actions appear as exponentially small corrections such as (A1/A2)g(A_1/A_2)^g after the nearest contribution is subtracted.

Logarithmic sectors leave logarithms of the order

Section titled “Logarithmic sectors leave logarithms of the order”

Because

βpβ=β(log)p,\partial_\beta^p\hbar^\beta = \hbar^\beta(\log\hbar)^p,

a companion sector containing (log)p(\log\hbar)^p differentiates the factorial-over-power expression with respect to β\beta. The result contains logg\log g, polygamma functions, and inverse powers of gg. Fitting a pure gamma law to a resonant logarithmic sector therefore misidentifies both the amplitude and the exponent.

Large order is consequently a stringent audit of a proposed transseries, but not a substitute for its global boundary condition. It cannot fix a transseries parameter that is exponentially invisible to the formal perturbative branch.

A quartic thimble supplies a nontrivial companion series

Section titled “A quartic thimble supplies a nontrivial companion series”

The Euler pole has a constant companion sector. A zero-dimensional quartic integral gives an exact next step in which both sectors carry nontrivial divergent fluctuations. It is a resurgence calibration, not a Schrödinger spectral problem.

For g>0g>0, define the convergent stable integral

Zst(g):=12πgRexp ⁣[1g(x22+x424)] ⁣dx.Z_{\mathrm{st}}(g) := \frac1{\sqrt{2\pi g}} \int_{\mathbb R} \exp\!\left[ -\frac1g \left( \frac{x^2}{2}+\frac{x^4}{24} \right) \right] \dd x.

The analytically continued unstable action profile

U(x)=x22x424U_-(x) = \frac{x^2}{2}-\frac{x^4}{24}

has critical points

x0=0,x±=±6,x_0=0, \qquad x_\pm=\pm\sqrt6,

and the common saddle-action difference

A=U(x±)U(x0)=32.A = U_-(x_\pm)-U_-(x_0) = \frac32.

The real unstable integral diverges. Its perturbative expansion is defined through analytic continuation or lateral thimbles, not by an ordinary real-axis integral.

Introduce the two formal sectors

Φ^0(g)=n0angn,Φ^1(g)=n0bngn,\begin{aligned} \widehat\Phi_0(g) &= \sum_{n\geq0}a_ng^n, \\ \widehat\Phi_1(g) &= \sum_{n\geq0}b_ng^n, \end{aligned}

with

an=(4n)!96n(2n)!n!,bn=(1)nan.\begin{aligned} a_n &= \frac{(4n)!} {96^n(2n)!\,n!}, \\ b_n&=(-1)^na_n. \end{aligned}

Thus

Φ^0(g)=1+g8+35g2384+,Φ^1(g)=1g8+35g2384.\begin{aligned} \widehat\Phi_0(g) &= 1+\frac g8+\frac{35g^2}{384}+\cdots, \\ \widehat\Phi_1(g) &= 1-\frac g8+\frac{35g^2}{384}-\cdots. \end{aligned}

To use the book’s shifted transform, apply it to gΦ^jg\widehat\Phi_j. The two Borel germs are

B(gΦ^0)(ξ)=2F1 ⁣(14,34;1;ξA),B(gΦ^1)(ξ)=2F1 ⁣(14,34;1;ξA).\begin{aligned} \mathcal B(g\widehat\Phi_0)(\xi) &= {}_2F_1\!\left( \frac14,\frac34;1;\frac{\xi}{A} \right), \\ \mathcal B(g\widehat\Phi_1)(\xi) &= {}_2F_1\!\left( \frac14,\frac34;1;-\frac{\xi}{A} \right). \end{aligned}

On the first cut, for ξ>A\xi>A, their boundary values obey

B(gΦ^0)(ξ+i0)B(gΦ^0)(ξi0)=i2B(gΦ^1)(ξA).\begin{aligned} & \mathcal B(g\widehat\Phi_0)(\xi+\ii0) - \mathcal B(g\widehat\Phi_0)(\xi-\ii0) \\ &\qquad= \ii\sqrt2\, \mathcal B(g\widehat\Phi_1)(\xi-A). \end{aligned}

Laplace transformation therefore gives

Disc0anΦ^0=i2eA/gS0Φ^1.\operatorname{Disc}_0^{\mathrm{an}} \widehat\Phi_0 = \ii\sqrt2\, \ee^{-A/g} \mathcal S_0\widehat\Phi_1.

The companion-sector formula predicts

an1π2k=0K1bkΓ(nk)Ank.a_n \sim \frac1{\pi\sqrt2} \sum_{k=0}^{K-1} b_k \frac{\Gamma(n-k)} {A^{n-k}}.

The first loop corrections are

anΓ(n)π2An[1316(n1)+105512(n1)(n2)+].\begin{aligned} a_n \sim \frac{\Gamma(n)} {\pi\sqrt2\,A^n} \biggl[ 1 &-\frac{3}{16(n-1)} \\ &+ \frac{105} {512(n-1)(n-2)} +\cdots \biggr]. \end{aligned}

In particular,

nanan+132,AnanΓ(n)1π2.n\frac{a_n}{a_{n+1}} \longrightarrow \frac32, \qquad \frac{A^na_n}{\Gamma(n)} \longrightarrow \frac1{\pi\sqrt2}.

The stable companion has an independent exact evaluation,

Zst(g)=32πge3/(4g)K1/4 ⁣(34g).\begin{aligned} Z_{\mathrm{st}}(g) = \sqrt{\frac{3}{2\pi g}}\, \ee^{3/(4g)} K_{1/4}\!\left(\frac{3}{4g}\right). \end{aligned}

It equals the regular positive-ray sum of Φ^1\widehat\Phi_1. For the completed unstable family

Z^(σ;g)=Φ^0(g)+σeA/gΦ^1(g),\widehat Z(\sigma;g) = \widehat\Phi_0(g) +\sigma\ee^{-A/g}\widehat\Phi_1(g),

the Stokes constant is S=i2\mathsf S=\ii\sqrt2. Hence the upper parameter i/2-\ii/\sqrt2 and lower parameter +i/2+\ii/\sqrt2 give the same balanced real value. In this linear example that value is also the arithmetic mean of the two perturbative laterals; that simplification is not the general definition of median summation.

This laboratory separates three independent checks: saddle geometry fixes AA, the hypergeometric cut fixes S\mathsf S, and the exact coefficients test the resulting large-order prediction.

The double well separates splitting from ambiguity

Section titled “The double well separates splitting from ambiguity”

Consider a real symmetric double well with weak positive coupling gg, two degenerate minima, and one-instanton action SI>0S_I>0. For a level label NN and parity ε=±1\varepsilon=\pm1, the model-dependent expansion has the schematic but structurally important form

E^ε,N(g)=0eN,(0)g+n1[CN(g)eSI/g]n×p=0n1[log ⁣(κg)]p0eε,N;n,p,g.\begin{aligned} \widehat E_{\varepsilon,N}(g) ={}& \sum_{\ell\geq0} e_{N,\ell}^{(0)}g^\ell \\ &+ \sum_{n\geq1} \left[ \mathcal C_N(g)\ee^{-S_I/g} \right]^n \\ &\quad\times \sum_{p=0}^{n-1} \left[ \log\!\left(-\frac{\kappa}{g}\right) \right]^p \sum_{\ell\geq0} e_{\varepsilon,N;n,p,\ell}g^\ell. \end{aligned}

The prefactor CN(g)\mathcal C_N(g), the positive constant κ\kappa, and the coefficient normalizations depend on the potential convention. The maximum logarithmic power n1n-1 is characteristic of the standard multi-instanton interaction analysis; it is not a universal law for every transseries.

Two selection rules prevent a common misreading.

  1. A single instanton changes wells. It produces the leading parity splitting

    E+,NE,N=O(eSI/g).E_{+,N}-E_{-,N} = O(\ee^{-S_I/g}).
  2. Perturbation theory about one specified well begins and ends in the same vacuum. Its first neutral nonperturbative neighbor is an instanton–anti-instanton sector with action 2SI2S_I:

    DiscE^N(0)=O(e2SI/g).\operatorname{Disc} \widehat E_N^{(0)} = O(\ee^{-2S_I/g}).

Consequently, the one-instanton splitting need not be the singularity seen by the leading large order of the common perturbative series. The neutral two-event sector is.

One precise realization uses the normalization

[22 ⁣d2 ⁣dx2+x2(1+x)22]ψ=Eψ.\left[ -\frac{\hbar^2}{2}\frac{\dd^2}{\dd x^2} + \frac{x^2(1+x)^2}{2} \right]\psi = E\psi.

For energies below the barrier, let

t:=ΠA2π,tD:=iΠB,t:=\frac{\Pi_A}{2\pi}, \qquad t_D:=-\ii\Pi_B,

in the cycle convention of van Spaendonck and Vonk. Their two lateral quantization conditions can be written

Dε±=1+e±2πit/iεetD(t;)/(2)=0,\begin{aligned} D_\varepsilon^\pm = 1 +\ee^{\pm2\pi\ii t/\hbar} \mp\ii\varepsilon \ee^{-t_D(t;\hbar)/(2\hbar)} =0, \end{aligned}

where ε=±1\varepsilon=\pm1 is parity. This kinetic normalization differs from the unit-kinetic convention used elsewhere in the book, so its periods must not be imported without translation.

The perturbative root is

t0=(N+12).t_0=\hbar\left(N+\frac12\right).

Set

q(t,):=etD(t;)/(2)q(t,\hbar) := \ee^{-t_D(t;\hbar)/(2\hbar)}

and t=t0+Δt±t=t_0+\Delta t^\pm. The exact implicit displacement is

Δt±=i2πlog ⁣[1iεq(t0+Δt±,)].\begin{aligned} \Delta t^\pm = \mp\frac{\ii\hbar}{2\pi} \log\!\left[ 1\mp\ii\varepsilon q(t_0+\Delta t^\pm,\hbar) \right]. \end{aligned}

Retaining the tt-dependence of qq gives

Δt±=ε2πq+[i4π8π2ttD]q2+O(q3),\begin{aligned} \Delta t^\pm ={}& -\varepsilon\frac{\hbar}{2\pi}q \\ &+ \left[ \mp\frac{\ii\hbar}{4\pi} -\frac{\hbar}{8\pi^2} \partial_t t_D \right]q^2 \\ &+ O(q^3), \end{aligned}

with qq and ttD\partial_t t_D evaluated at t0t_0. The first correction is real and parity dependent. Lateral imaginary ambiguity first appears at q2q^2, while q2=exp(tD/)q^2=\exp(-t_D/\hbar) is the neutral action scale seen by perturbative large order. This explicit calculation is the promised counterexample to the slogan that the first spectral splitting must control the perturbative coefficients.

For g>0g>0, the logarithm has two lateral values,

log ⁣(κg)±=log ⁣(κg)±πi.\log\!\left(-\frac{\kappa}{g}\right)_{\pm} = \log\!\left(\frac{\kappa}{g}\right) \pm\pi\ii.

The imaginary ambiguity carried by the neutral two-instanton term cancels the ambiguity in the lateral Borel sum of the perturbative sector. At the next grades, the ambiguity of a one-instanton fluctuation series communicates with three-instanton sectors, the two-instanton series with four-instanton sectors, and so on. The exact quantization condition fixes the coefficients and signs of this cancellation ladder.

Balanced upper and lower transseries representatives meet at the median sum beside a double-well event lattice separating parity splitting from neutral ambiguity cancellation.

Balanced upper and lower representatives meet after opposite half-Stokes shifts. In the symmetric-double-well event lattice, the vertical coordinate counts instanton events and the horizontal coordinate is the net well-changing charge. Odd sectors control parity splitting, whereas the neutral column communicates with the perturbative vacuum and carries its ambiguity-canceling even sectors. Lattice arrows add an instanton I\mathcal I or anti-instanton I\overline{\mathcal I}; their coefficients are model-dependent Stokes data.

The triangle records action and topological selection rules. It does not compute fluctuation determinants, quasi-zero-mode integrals, Stokes constants, or boundary-condition weights. Nor does it imply that every allowed node is nonzero. Symmetry or a vanishing connection coefficient can remove a sector.

The neutral pair is sometimes called a bion. Depending on the model, its semiclassical representative can involve a correlated saddle rather than a real exact two-event solution. The resurgent statement concerns the normalized sector and its Stokes relation, not the existence of a particular real trajectory.

Exact WKB makes the cancellation a covariance statement

Section titled “Exact WKB makes the cancellation a covariance statement”

Page 5 supplied the Stokes automorphism of Voros symbols. Page 6 inserted those symbols into a boundary determinant. Combining the two gives the cleanest route to spectral ambiguity cancellation.

For one active saddle class δ\delta, orient q=VδeA/q=\mathcal V_\delta\sim\ee^{-A/\hbar} to be small. Page 5’s DDP formula gives, for a relative class β\beta,

Xβ,+=Xβ,(1+q)δ,β.\mathcal X_{\beta,+} = \mathcal X_{\beta,-} (1+q_-)^{-\langle\delta,\beta\rangle}.

On a continuous logarithmic branch,

DiscanlogXβ=δ,βlog(1+q)=δ,β1(1)+1q.\begin{aligned} \operatorname{Disc}^{\mathrm{an}} \log\mathcal X_\beta ={}& -\langle\delta,\beta\rangle \log(1+q_-) \\ ={}& -\langle\delta,\beta\rangle \sum_{\ell\geq1} \frac{(-1)^{\ell+1}}{\ell} q_-^\ell. \end{aligned}

One primitive active class therefore generates every composite action grade A\ell A. These coefficients are pieces of the full Stokes automorphism. They are not automatically distinct primitive Euclidean instantons, and raw binomial coefficients should not be relabeled as alien derivatives without taking the logarithm of the automorphism.

Let Δ^\widehat\Delta denote the formal determinant in one graph-adapted coordinate system. On a singular direction,

Sθ+Δ^=Sθ(SθΔ^).\mathcal S_{\theta+}\widehat\Delta = \mathcal S_{\theta-} \left( \mathfrak S_\theta\widehat\Delta \right).

The two sides are the same analytic boundary function expressed through the upper and lower lateral connection data. If a simple root is continued consistently, the corresponding energy transseries obeys the induced Stokes map. The root does not acquire a physical ambiguity; its formal coordinates do.

At the formal level, the same statement has a compact implicit chain-rule form. If a homogeneous pointed alien derivation is taken with the moving action lattice locally trivialized, then

0=Δ˙[Q^(E^,)],Δ˙E^=Δ˙EQ^E=E^EQ^E=E^.\begin{aligned} 0 &= \dot\Delta \left[ \widehat{\mathcal Q}(\widehat E,\hbar) \right], \\ \dot\Delta\widehat E &= -\frac{ \left.\dot\Delta_E\widehat{\mathcal Q}\right|_{E=\widehat E} }{ \left.\partial_E\widehat{\mathcal Q}\right|_{E=\widehat E} }. \end{aligned}

This identity is conditional on parametric resurgence and the same uniformity in EE needed for the analytic implicit problem.

This argument needs four separate ingredients:

  1. the DDP or relevant wall-crossing automorphism for every active Voros symbol;
  2. a boundary determinant assembled with the same path, sheet, and frame conventions on both sides;
  3. a simple root, or an appropriate multiple-root replacement;
  4. summability and analytic continuation uniform enough in EE to pass the Stokes relation through the implicit solution.

Without the fourth item, formal covariance is not yet an equality of analytic eigenvalues.

Sector normalization changes the reported Stokes constant

Section titled “Sector normalization changes the reported Stokes constant”

Rescale a companion sector by a nonzero constant,

E^1λE^1.\widehat E_1 \longmapsto \lambda\widehat E_1.

To keep both the transseries and bridge relation fixed, transform

σσλ,S01S01λ.\sigma \longmapsto \frac{\sigma}{\lambda}, \qquad \mathsf S_{01} \longmapsto \frac{\mathsf S_{01}}{\lambda}.

Therefore large order determines the invariant product S01e0(1)\mathsf S_{01}e_0^{(1)}, not either factor separately. A numerical Stokes constant is meaningful only together with the Borel convention, lateral orientation, action branch, and sector normalization.

A practical extraction and cancellation workflow

Section titled “A practical extraction and cancellation workflow”
  1. Start from a boundary determinant whose endpoint lines and normalization are already fixed.
  2. Choose a chamber and factor every exponentially small Voros symbol into its classical action, algebraic power, and normalized quantum series.
  3. Solve the zero-action equation for a simple perturbative root.
  4. Expand the determinant recursively in the small transmonomials; allow Puiseux powers or logarithms when the recursion demands them.
  5. Derive the Stokes map from the Voros-symbol automorphism and the determinant. Do not infer it only from a desired reality property.
  6. Compute many perturbative coefficients and test their large order against the closest permitted companion sectors.
  7. Subtract the leading factorial-over-power contribution before fitting a farther action.
  8. Fix the physical transseries parameters from parity, Bloch, decay, outgoing, or other global boundary data.
  9. Use matched lateral sums or the full median automorphism on a singular ray.
  10. Compare against an independent spectral calculation and vary the working precision, truncation, and matching point.

The order matters. A fitted exponential scale is evidence for a Borel singularity, not by itself a proof that a guessed saddle and boundary condition give the correct sector.

The companion script instanton-transseries-check.py audits independent representations rather than fitting a formula back to the data that defined it:

  • symbolic implicit-root coefficients against the exact Lambert-WW laboratory;
  • the moving-action term in the exact lateral double-well displacement;
  • the signed Euler lateral sums and their median cancellation at high precision;
  • quartic saddle geometry, Gaussian moments, an ODE recurrence, and the hypergeometric Borel cut;
  • companion-sector large order and the independent Bessel/quadrature evaluation of the stable quartic integral;
  • recovery of a farther signed action after exact subtraction of the dominant Borel pole;
  • the double-well event lattice rule Qn|Q|\leq n and Qn(mod2)Q\equiv n\pmod 2.

It uses explicit runtime checks rather than Python assert, so optimized mode performs the same audit.

Run it from the project root:

Terminal window
python3 public/code/advanced-ode/instanton-transseries-check.py

These checks verify algebra, signs, and numerical diagnostics. They do not prove resurgence or Borel summability for a new potential.

Adding one exponential is not a transseries completion. The fluctuation series around that exponential is generally divergent, its lateral ambiguity couples to higher sectors, and resonance can introduce logarithms. Closure under the relevant Stokes automorphisms is the real test.

A small WKB cycle is not automatically one Euclidean instanton. Orientations, doubled paths, and normalization conventions can change the reported action by a sign or factor of two. Compare the actual exponent, not the label.

The closest positive action need not control every large-order coefficient. Equal-modulus or conjugate singularities can interfere, and selection rules can make the leading Stokes constant vanish. Fit the sum allowed by the geometry.

The Stokes constant depends on a sector basis. Rescaling a companion sector inversely rescales the reported constant. Quote the product tested by large order or state the normalization.

Median is not synonymous with real part. The median operation uses half of the full Stokes automorphism. A termwise principal value can miss coupled sectors and nonlinear parameter maps.

Reality is not universal. It follows for a suitable self-adjoint problem with compatible conjugation. Resonances and non-Hermitian spectra are genuinely complex and are treated on Page 8.

A transseries parameter is fixed globally. Perturbative large order can reveal Stokes-weighted neighboring sectors but not the boundary datum that selects one physical member of the family.

Let

Q(E,q)=Q0(E)+qQ1(E)+q2Q2(E)+O(q3),\mathcal Q(E,q) = \mathcal Q_0(E) +q\mathcal Q_1(E) +q^2\mathcal Q_2(E) +O(q^3),

and suppose Q0(E0)=0\mathcal Q_0(E_0)=0 with Q0(E0)0\mathcal Q_0'(E_0)\neq0. Derive E1E_1 and E2E_2 for E=E0+qE1+q2E2+E=E_0+qE_1+q^2E_2+\cdots.

Solution

Taylor expansion at E0E_0 gives

0=q(Q0E1+Q1)+q2(Q0E2+12Q0E12+Q1E1+Q2)+O(q3).\begin{aligned} 0 ={}& q\left( \mathcal Q_0'E_1+\mathcal Q_1 \right) \\ &+ q^2 \biggl( \mathcal Q_0'E_2 +\frac12\mathcal Q_0''E_1^2 \\ &\qquad +\mathcal Q_1'E_1 +\mathcal Q_2 \biggr) +O(q^3). \end{aligned}

Equating coefficients yields

E1=Q1Q0,E_1=-\frac{\mathcal Q_1}{\mathcal Q_0'},

and

E2=12Q0E12+Q1E1+Q2Q0.E_2 = -\frac{ \tfrac12\mathcal Q_0''E_1^2 +\mathcal Q_1'E_1 +\mathcal Q_2 }{ \mathcal Q_0' }.

Every function and derivative on the right is evaluated at E0E_0.

Starting from

EE+qeλE=0,E-E_*+q\ee^{\lambda E}=0,

derive the exact Lambert-WW solution and its first three nonzero powers of qq.

Solution

Set u=λ(EE)u=\lambda(E_*-E). The equation becomes

ueu=λqeλE.u\ee^u = \lambda q\ee^{\lambda E_*}.

Hence

E=Eλ1W ⁣(λqeλE).E = E_*-\lambda^{-1} W\!\left(\lambda q\ee^{\lambda E_*}\right).

Using

W(x)=xx2+32x3+O(x4)W(x)=x-x^2+\frac32x^3+O(x^4)

gives

E=EqeλE+λq2e2λE32λ2q3e3λE+O(q4).\begin{aligned} E =E_* &-q\ee^{\lambda E_*} \\ &+\lambda q^2\ee^{2\lambda E_*} \\ &-\frac32\lambda^2q^3 \ee^{3\lambda E_*} +O(q^4). \end{aligned}

Use the Page 1 lateral sums to show that

S0+Φ^(σπi)=S0Φ^(σ+πi)\mathcal S_{0+}\widehat\Phi(\sigma-\pi\ii) = \mathcal S_{0-}\widehat\Phi(\sigma+\pi\ii)

for Φ^(σ)=Φ^0+σe1/\widehat\Phi(\sigma)=\widehat\Phi_0+\sigma\ee^{-1/\hbar}.

Solution

The upper perturbative sum contains +πie1/+\pi\ii\ee^{-1/\hbar}, while its flat coefficient is σπi\sigma-\pi\ii. Their sum is

ΦPV+σe1/.\Phi_{\mathrm{PV}} +\sigma\ee^{-1/\hbar}.

The lower perturbative sum contains πie1/-\pi\ii\ee^{-1/\hbar}, while its flat coefficient is σ+πi\sigma+\pi\ii. It gives the same expression.

Assume

eg=CΓ(g+1β)Ag+1β.e_g = C\, \frac{\Gamma(g+1-\beta)} {A^{g+1-\beta}}.

Show that

(g+1β)egeg+1=A(g+1-\beta)\frac{e_g}{e_{g+1}}=A

exactly.

Solution

Use Γ(g+2β)=(g+1β)Γ(g+1β)\Gamma(g+2-\beta)=(g+1-\beta)\Gamma(g+1-\beta). The constant CC and the common power of AA cancel, leaving AA.

For constants rA,rBr_A,r_B and 0<A<B0<|A|<|B|, let

Bf^(ξ)=rAAξ+rBBξ.\mathcal B\widehat f(\xi) = \frac{r_A}{A-\xi} +\frac{r_B}{B-\xi}.

Find the shifted-Borel coefficients ana_n in f^=n1ann\widehat f=\sum_{n\geq1}a_n\hbar^n. Explain how subtracting the AA contribution exposes BB.

Solution

Expanding each pole at the origin gives

rAAξ=n1rAAnξn1,\frac{r_A}{A-\xi} = \sum_{n\geq1} \frac{r_A}{A^n}\xi^{n-1},

and similarly for BB. Since the shifted transform has coefficient an/Γ(n)a_n/\Gamma(n),

an=Γ(n)(rAAn+rBBn).a_n = \Gamma(n) \left( \frac{r_A}{A^n} +\frac{r_B}{B^n} \right).

Subtracting rAAnr_AA^{-n} from an/Γ(n)a_n/\Gamma(n) leaves exactly rBBnr_BB^{-n}. Ratios of that residual recover BB, including its phase or sign.

6. Locate the first neutral double-well sector

Section titled “6. Locate the first neutral double-well sector”

Let nIn_I and nIˉn_{\bar I} count instantons and anti-instantons. Define

n=nI+nIˉ,Q=nInIˉ.n=n_I+n_{\bar I}, \qquad Q=n_I-n_{\bar I}.

Show that Qn|Q|\leq n and Qn(mod2)Q\equiv n\pmod2. Which smallest nonzero nn returns to the original well?

Solution

Solving for the event numbers gives

nI=n+Q2,nIˉ=nQ2.n_I=\frac{n+Q}{2}, \qquad n_{\bar I}=\frac{n-Q}{2}.

They are nonnegative integers precisely when Qn|Q|\leq n and QQ has the same parity as nn. Returning to the original well means Q=0Q=0, so the smallest nonzero possibility is n=2n=2: one instanton and one anti-instanton.

7. Predict the effect of one logarithm on large order

Section titled “7. Predict the effect of one logarithm on large order”

Suppose the discontinuity contains

eA/βlog.\ee^{-A/\hbar} \hbar^\beta\log\hbar.

Use differentiation with respect to β\beta to determine the structure of its contribution to eg(0)e_g^{(0)}.

Solution

Differentiate the no-logarithm factor:

β[Γ(g+1β)Ag+1β].\partial_\beta \left[ \frac{\Gamma(g+1-\beta)} {A^{g+1-\beta}} \right].

Because ββ=βlog\partial_\beta\hbar^\beta =\hbar^\beta\log\hbar, the result has the required structure. Explicitly,

Γ(g+1β)Ag+1β[logAψ(g+1β)],\frac{\Gamma(g+1-\beta)} {A^{g+1-\beta}} \left[ \log A-\psi(g+1-\beta) \right],

where ψ\psi is the digamma function. Since ψ(g+1β)logg\psi(g+1-\beta)\sim\log g, the large-order sequence contains a logg\log g correction.

8. Prove the operator form of the median identity

Section titled “8. Prove the operator form of the median identity”

Starting from

Sθ+=SθSθ,\mathcal S_{\theta+} = \mathcal S_{\theta-}\circ\mathfrak S_\theta,

show that

Sθ+Sθ1/2=SθSθ1/2.\mathcal S_{\theta+}\circ\mathfrak S_\theta^{-1/2} = \mathcal S_{\theta-}\circ\mathfrak S_\theta^{1/2}.
Solution

Compose the defining relation on the right with Sθ1/2\mathfrak S_\theta^{-1/2}:

Sθ+Sθ1/2=SθSθSθ1/2=SθSθ1/2.\begin{aligned} \mathcal S_{\theta+} \circ\mathfrak S_\theta^{-1/2} &= \mathcal S_{\theta-} \circ\mathfrak S_\theta \circ\mathfrak S_\theta^{-1/2} \\ &= \mathcal S_{\theta-} \circ\mathfrak S_\theta^{1/2}. \end{aligned}

The action filtration makes the formal half-power well defined whenever logSθ\log\mathfrak S_\theta is.

9. Recover the quartic action from coefficients

Section titled “9. Recover the quartic action from coefficients”

Starting from

an=(4n)!96n(2n)!n!,a_n = \frac{(4n)!} {96^n(2n)!\,n!},

derive a first-order recurrence and use it to prove nan/an+13/2n a_n/a_{n+1}\to3/2.

Solution

Taking the ratio and canceling adjacent factors gives

an+1an=(4n+1)(4n+3)24(n+1).\frac{a_{n+1}}{a_n} = \frac{(4n+1)(4n+3)} {24(n+1)}.

Therefore

nanan+1=24n(n+1)(4n+1)(4n+3)2416=32.n\frac{a_n}{a_{n+1}} = \frac{24n(n+1)} {(4n+1)(4n+3)} \longrightarrow \frac{24}{16} = \frac32.

The recovered value agrees with the independently computed saddle-action difference.

Starting from

Δt±=i2πlog ⁣[1iεq(t0+Δt±,)],\Delta t^\pm = \mp\frac{\ii\hbar}{2\pi} \log\!\left[ 1\mp\ii\varepsilon q(t_0+\Delta t^\pm,\hbar) \right],

use

tlogq=ttD2\partial_t\log q = -\frac{\partial_t t_D}{2\hbar}

to derive the terms through q2q^2. Identify the first parity-dependent and the first lateral-ambiguous terms.

Solution

Write

q(t0+Δt)=q(tlogq)qΔt+O(q3).q(t_0+\Delta t) = q \left(\partial_t\log q\right)q\,\Delta t +O(q^3).

The coefficient of qq first gives

Δt1=ε2πq.\Delta t_1 = -\varepsilon\frac{\hbar}{2\pi}q.

Substitute this result back into the shifted qq and expand the logarithm. The result is

Δt±=ε2πq+[i4π8π2ttD]q2+O(q3).\begin{aligned} \Delta t^\pm ={}& -\varepsilon\frac{\hbar}{2\pi}q \\ &+ \left[ \mp\frac{\ii\hbar}{4\pi} -\frac{\hbar}{8\pi^2} \partial_t t_D \right]q^2 \\ &+ O(q^3). \end{aligned}

The order-qq term changes the parity levels in opposite directions. The explicit imaginary sign first appears at order q2q^2, the neutral action grade.

This page used a real double well to make the cancellation mechanism visible. The algebra of transseries and Stokes covariance survives when the relevant actions are complex, but the reality conclusion does not. Complex-conjugate Borel singularities produce oscillatory large order; outgoing boundary conditions select a nonselfadjoint spectral sheet; and turning-point motion can create Stokes transitions in the complex parameter plane.

Page 8 follows those complex actions through resonance and non-Hermitian boundary problems.