Resurgence in η-deformed Principal Chiral Models

We study the SU(2) Principal Chiral Model (PCM) in the presence of an integrable η-deformation. We put the theory on ℝ×S1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathbb{R}\times {S}^1 $$\end{document} with twisted boundary conditions and then reduce the circle to obtain an effective quantum mechanics associated with the Whittaker-Hill equation. Using resurgent analysis we study the large order behaviour of perturbation theory and recover the fracton events responsible for IR renormalons. The fractons are modified from the standard PCM due to the presence of this η-deformation but they are still the constituents of uniton-like solutions in the deformed quantum field theory. We also find novel SL2ℂ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathrm{S}\mathrm{L}\left(2,\mathbb{C}\right) $$\end{document} saddles, thus strengthening the conjecture that the semi-classical expansion of the path integral gives rise to a resurgent transseries once written as a sum over Lefschetz thimbles living in a complexification of the field space. We conclude by connecting our quantum mechanics to a massive deformation of the N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 4-d gauge theory with gauge group SU(2) and Nf = 2.


The resurgence paradigm
The calculation of the anomalous magnetic moment of the electron to ten significant figures stands as testament to efficacy of perturbation theory in quantum field theory (QFT) [1]. The fact that perturbative expansions, even in the simpler setting of quantum mechanics (QM), are quite generally divergent may then be a cause of some cognitive dissonance to the practicing physicist. Resurgence theory is perhaps the therapy that perturbation theory demands. Consider, for example, a ground state energy obtained as a perturbative expansion in some real and positive coupling g 2 , (1.1) Due to the proliferation of Feynman diagrams it is quite generic to find factorial growth c n ∼ A n n! and the perturbative expansion, whilst initially providing increasingly accurate results, will begin to diverge. Borel summation can provide a way to attach a meaning to such asymptotic series. In this method the factorial growth is essentially traded for an integral via n! = There is some danger here! The integral (1.3) may not be well defined; the Borel transformed series can have poles along the integration path. It is oft said that this occurs when the coefficients c n are non-alternating but even an alternating series can be non-Borel summable and indeed we will encounter this. In this case it is natural to deform the integration contour by extending t to the complex plane and integrating along a ray i.e., which is called directional Borel resummation. If B[E pert ](t) is regular for arg(t) = θ (and decays sufficiently fast at infinity), then S θ [E pert ](g 2 ) defines an analytic function in the wedge of the complex coupling plane −π/2 + θ ≤ arg( g 2 ) ≤ π/2 + θ. If B[E pert ](t) has singularities along the direction arg(t) = θ, this is called a Stokes direction. For example if the original direction of integration (1.3), arg(t) = 0, is a Stokes JHEP07(2016)088 direction we can dodge the singularities by picking θ = > 0 small and positive, in this way we avoid poles in the original integration cycle and we can calculate the well defined S [E pert ](g). However we could equally choose θ = − < 0 and calculate S − [E pert ](g).
If B[E pert ](t) has singularities in the complex wedge − ≤ arg( t) ≤ , these two lateral summations, S ± [E pert ](g), are generically different but still yield the same asymptotic perturbative expansion (1.1) once expanded at weak coupling. This jumping as one deforms an integration cycle is known as the Stokes phenomenon.
Generically, the original real positive coupling arg(g 2 ) = 0 lies on a Stokes line -the location at which such a jump takes place. Thus we have traded our initial asymptotic series expression for the ground state energy for a finite but ambiguous result. Evidently this is still unsatisfactory; physical observables surely shouldn't be blighted with such ambiguity.
We have of course missed the vital part of the physics; the non-perturbative sector. In these cases one should include in the path integral contributions from non-perturbative saddle points and the fluctuations around them. These non-perturbative contributions quite certainly correct the real part of the energy as in quantum mechanics where instanton effects can lift perturbative degeneracies or give band structures. These non-perturbative contributions are, however, by themselves ambiguous. As shown by Bogomolny [2] and Zinn-Justin [3], even the leading contribution to a path integral from integration over the quasi-zero-mode associated to separated instanton I and anti-instantonĪ pairs has an ambiguous imaginary part. Far from compounding our problems, this is just the medicine needed -the [IĪ] ambiguity exactly cancels that of the perturbative sector. This is remarkable: the divergent behaviour of late terms in the perturbative sector is compensated by early terms in a non-perturbative sector. This is the first example of resurgent behaviour.
The story runs deeper; the perturbative series around the one instanton sector [I] will itself also have factorial growth and gives rise, in a similar fashion, to an ambiguous imaginary part that will cancel against early terms in the [I IĪ] sector and so on. The resurgence conjecture is then that the unambiguous physical observable can be understood as an example ofÉcalle's resurgent trans-series [4] of the form, 2 (1.5) in which in addition to the perturbative sector one should sum over all non-perturbative saddles labelled by {i} with corresponding actions S i . From a path-integral perspective the transseries representation of physical observables looks exactly like the semi-classical approximation. The perturbative expansion around the standard vacuum configuration gives rise to an asymptotic series without any exponentially suppressed factor sitting in front, while from each non-trivial saddle point we obtain an exponentially suppressed factor, coming from the classical action of the saddle, times its associated asymptotic perturbative expansion.
A not so widely appreciated fact about the semiclassical expansion is that the amplitudes associated with certain saddle points, for example the aforementioned instantonanti-instanton events, are not well-defined precisely along the Stokes line g 2 real. This

JHEP07(2016)088
has a beautiful geometric origin, we can decomposeà la Morse [6] the original contour of integration over the space of fields into the sum of certain privileged contours of integrations called Lefschetz thimbles. As we vary the argument of the coupling constant two things will happen: first the semi-classical expansion will generically receive contributions from all finite action complex saddles living in a complexification of the space of fields and secondly precisely at Stokes line this thimble decomposition will jump discontinuously, that is, in addition to the jump in the resummation of the perturbative series, the sum over non-perturbative saddle points will also jump and the two "ambiguities" will precisely cancel.
After we sum over all the finite action complexified saddles, the semiclassical expansion (1.5) should not be thought of as an approximation, a resurgent trans-series is an exact coded representation of the physical observable, i.e. it really is an analytic function in a certain wedge of the complex coupling plane.

Resurgence in QFT and the PCM
The general philosophy above is expected to hold in quantum field theories in which both perturbative and non-perturbative contributions play a role. Factorially divergent series abound in 4-d Yang-Mills or QCD arising both from growth in diagrammatica and from the so-called renormalons [7] coming from IR and UV contributions to loop momenta integration. This standard renormalons diagrams argument suggest that generically, in asymptotically free theory, we should expect singularities in the Borel transform along the real axis. However, once we consider these theories on R 3 × S 1 , due to the presence of non-trivial holonomies around the circle, we see [8] that renormalon type divergences are actually reproduced by the combinatorics of diagrams getting mixed with the discretized momenta along the S 1 . We refer to IR/UV renormalons as singularities in the Borel transform along the positive/negative real axis that remain at finite locations once we take the approriate large-N limits, without making any reference to their diagrammatic origins or otherwise. 3 As a first step towards full QCD, it is natural to study the ideas and implications of resurgence in a simpler setting of two-dimensions [9]. The SU(N ) Principal Chiral Model (PCM), whose action is given by, with g ∈ SU(N ), provides a rather effective toy model to QCD since it displays asymptotic freedom, a dynamically generated mass gap and confinement. Being an integrable theory in two-dimension [10,11] the PCM has the advantage of also being sufficiently simple to remain tractable; for instance the mass gap can be exactly determined [12,13]. Although the SU(N ) PCM does not possess instantons (in the sense that π 2 [SU(N )] = 0 and there is no topological charge), it does have a non-perturbative sector whose constituents, known as unitons [14], are finite action field configurations that solve second order Euclidean equations of motion, i.e. they are not BPS objects but genuine saddles [15].

JHEP07(2016)088
In [9,16] the resurgent structure of the two-dimensional SU(N ) PCM was studied. The technique employed was to study the spatially compactified Euclidean theory with twisted boundary conditions corresponding to turning on a background gauge field for the diagonal SU(N ) ⊂ SU(N ) L × SU(N ) R global symmetry of the theory. This regime is adiabatically connected to the PCM defined on R 2 and many features of the original theory remain visible in an effective one-dimensional quantum mechanics that comes when taking the radius of spatial compactification very small. In particular, the large order behaviour of this QM exhibits non-alternating factorial growth with poles along the positive real half-line of the Borel plane. These poles can be matched precisely to certain instantonanti-instanton configurations in the QM. What is quite astonishing is that the origin of these non-perturbative saddles can be exactly identified in the two-dimensional quantum field theory -they correspond to "fractons". When the PCM is put on R × S 1 the unitons fractionalize and break up into these fractons when their size becomes comparable to the radius of compactification. In this way, the fractons lead to a semi-classical realisation of IR renormalons and the mass gap in the circle-compactified theory. This is a generic feature of 2-d (and perhaps 4-d) asymptotically free QFTs as for example the CP N model [17,18], the Grassmanian models and the O(N ) model [19]. In particular, the mass gap of the O(6) sigma model directly relates to non-perturbative scales in AdS/CFT [20] and more specifically it translates to non-perturbative effects in the 4-d N = 4 Super-Yang-Mills cusp anomalous dimension, whose resurgence properties were the subject of [21,22]. We should emphasise resurgence is not the statement that the perturbative sector encodes all non-perturbative information. Instead resurgence can allow sectors of a theory of a given fixed topological charge to communicate when need to cancel ambiguities -this idea is called the resurgence triangle [17]. In some cases extra symmetries can mean certain sectors simply do not exist and the resurgence triangle can truncate [17].
Whilst our main focus in the present context is in 2-d QFTs, important progress has already been made in higher dimensional QFTs in the context of 3d Chern-Simons matter (ABJM) and N = 2 supersymmetric 4d Yang-Mills gauge theories [23] building on [24], where supersymmetric localisation allows for comparison to exact quantities. Recently this was further developed in [25] and extended to 5d N = 1 supersymmetric theories. The examples of 4d N = 2 [23,24] and 5d N = 1 [25] supersymmetric gauge theories show that the perturbative expansion is Borel summable within a given topological sectorpresumably as a consequence of the extended supersymmetry -and the Borel resummation does not require communication between different sectors.

Resurgence puzzles
The study of resurgence has thus far proven to be very profitable but let us mention two puzzling features that the present study will help to shed light on: 1. Complex Saddles and QFT. The situation described above in which the cuts in the Borel plane lie on either the positive or the negative real axis is not at all generic.
In fact we should expect that in general the singularities might be scattered across JHEP07(2016)088 the entire complex plane. To understand how resurgence works in these cases will require a better understanding of the non-perturbative saddles that give rise to such singularities. These will typically be associated to complex field configurations that will need to be understood in a field theory context.

2.
Multiple couplings and Borel flow. The Lagrangian defining the QFT may have several bare parameters and complicated RG flows. In this case it is not obvious how to understand physical observables as resurgent trans-series. Furthermore the semiclassical regime in which we will apply resurgence is only adiabatically connected to the original quantum field theory. The structure of the Borel plane does not remain invariant under this adiabatic flow and having multiple couplings, together with RG-invariant quantities, might tell us how this Borel flow takes place not only qualitatively but also quantitatively.

The η-deformed Principal Chiral Model
To advance the study of resurgence in QFT and target the puzzling features, we would like to find a candidate theory that is rich enough to include a variety of non-perturbative sectors and multiple couplings but yet remain tractable. A class of deformations of the PCM chiral model that has received substantial recent interest are the Yang-Baxter deformed σ-models introduced by Klimcik [26,27]. These theories, which exactly meet our requirements, are specified by a certain matrix R, that we shall detail later, and a parameter η (which a priori should be taken to be real) and are defined by the action, Whilst the deformation breaks the global symmetry down to SU(N ) L × U(1) R , the theory remains integrable (at least classically) and thus amenable to study. The broken SU(N ) R symmetry is not completely lost; its avatar is present in the structure of non-local charges whose Poisson bracket relations give rise to a classical version of a quantum group U q (sl(N )) with parameter q = exp (πηt) [28]. This quantum group symmetry is expected to be manifest in the exact S-matrix of the theory. This Yang-Baxter deformation, after suitable modifications, can be equally well applied to the PCM on symmetric or semi-symmetric spaces (supercosets) [29] including the wellknown example of PSU(2, 2|4)/SO(4, 1) × SO(5) relevant to superstrings in AdS 5 × S 5 [30]. This offers intriguing prospects as a Lagrangian description for quantum group deformations of holography. In the later context the Yang-Baxter σ-model is often called an η-deformation and we adopt this nomenclature.
There is much left to be understood about η-deformed AdS 5 × S 5 [31,32] and in particular the status of the deformed σ-model as a CFT is questionable but that concern is irrelevant to our current investigation. It would be remiss not to mention that a related class of integrable theories [33][34][35] known as λ-deformations has recently been developed that extend the idea of a current-current perturbation of a (gauged-)Wess-Zumino-Witten model. These λ-deformations are classically related to (1.7) via a Poisson-Lie T-duality JHEP07(2016)088 together with some analytic continuations [36][37][38]. This class of theories is expected to display a similar quantum group symmetry but with the parameter q = exp iπ k a root of unity. There is by now good evidence that λ-deformed σ-models based on semi-symmetric spaces have target space times that solve the equations of type II supergravity and may define CFTs [39][40][41]. For a recent review see [42].

Summary of results
In a fortuitous case of symbiosis, in this work we apply the ideas of resurgence to the η-deformed σ-model and draw concurrent insight about the technique of resurgence itself and the structure of the η-deformed theory. In this work we focus on the simplest case of the SU(2) η-deformed principal chiral model -the extension to higher rank groups should be quite similar and will be pursued elsewhere. Let us summarise some of the key findings: • Deformed unitons. We demonstrate explicitly that the known uniton solutions of the PCM persist in a suitably deformed sense as classical saddles of the action (1.7). These give rise to finite, and quantised, actions which are not classified by the usual topological instanton charge.
• Complex saddles. Intriguingly we find new complex saddles of (1.7) whose actions are finite but as η → 0 diverge and so do not contribute to the undeformed PCM.
• Dimensional reduction. After applying a twisted spatial compactification the η-deformed theory gives rise to a quantum mechanics with a two-term periodic potential of the Whittaker-Hill type (known to be a quasi-integrable QM [43,44]). This theory has both instanton configurations and complex field configurations whose actions match respectively those of the fractons constituents of the unitons and the complex saddles in the two-dimensional QFT.
• Large order behaviour and resurgence. The perturbative expansion of the ground state energy gives an asymptotic series whose Borel transform has cuts in both the positive and the negative half line. The leading and first sub-leading terms of the large order perturbation theory are precisely tied to the above non-perturbative saddles in line with expectations from resurgence.
• Critical coupling η c . For η > η c ≈ 0.27 the dominant contribution to the large order behaviour of the perturbative energy coefficients, and the singularity nearest the origin (but on the negative real axis) in the Borel plane, comes not from instantonanti-instanton events but from the complex saddle. The large order perturbation theory for η > η c is alternating in sign but, due to the sub-dominant contribution coming from the singularity in the positive half line associated with instanton-antiinstanton events, this series remains non-Borel summable along the positive real line.
• Analytic continuation in η. We consider also the case of pure imaginary η = iη R , with η R ∈ R. The pole structure on the Borel plane now includes similar cuts in the positive real axis and also a complex conjugate pair of cuts (with positive real part)

JHEP07(2016)088
that correspond to complex saddles. There are three regimes: 0 ≤ η R 0.4 in which the singularity along the real axis dominate; 0.4 η R ≤ 1 √ 2 in which the dominating large order behaviour is due to the complex conjugate pairs of singularities; 1 √ 2 < η R ≤ 1 in which the imaginary part of the complex conjugate pair vanishes and the relevant quantum mechanics has two different real types of instanton corresponding to tunnelling between vacua separated by large and small barriers.
• Stokes jumps and Seiberg-Witten theory. A Stokes graph consists of trajectories with constant phase of a WKB functional. 4 As the coupling is analytically continued over a Stokes line the topology of this graph will change. The same construction arises in Seiberg-Witten theory where the topology jumps are a manifestation of Wall Crossing. We find that the QM of the η-deformed theory gives a Stokes graph for an N f = 2 mass deformation of an SU(2) gauge theory with equal masses for the flavours. See also the very recent papers [45,46].

Outline
We proceed in section 2 where we shall first recapitulated the requisite details of the η-deformed Principal Chiral Model before discussing its non-perturbative semi-classical configurations of the deformed uniton and the complex saddle. In section 3 we then consider the theory on R × S 1 and illustrate the fractionalisation of the saddle point configurations and give details of the reduction to quantum mechanics. In section 4 we examine the large order behaviour of perturbation theory of this reduced theory paying attention to the location of the poles in the Borel plane. In section 5 we demonstrate the resurgent aspects of this system by comparison of the large order perturbative behaviour to low orders of perturbation theory around non-perturbative saddles. In section 6 we discuss the Stokes graphs and match to Seiberg-Witten theory. We then return to a discussion and interpretation of our results.

PCM background
The Principal Chiral Model (PCM) for the group G is given by, in which Σ is the world sheet taken to be R (1,1) and the right-invariant forms are R a ± = R a µ ∂ ± X µ = Tr(T a ∂ ± gg −1 ) for generators T a of the algebra 5 g whose dimension is |G|. The PCM is asymptotically free and strongly coupled in the IR with a mass scale generated through dimensional transmutation. As is well known, this theory has a G L ⊗G R symmetry under which g transform as g → U L g U R with U L/R ∈ G L/R two constant matrices, and is JHEP07(2016)088 both classically and quantum mechanically integrable [10,11]. The S-matrix is two-particle factorisable and hence multi-particle scattering processes can be obtained as a succession of two-particle scattering events. The S-matrix can be determined from the factorised bootstrap and the theory solved using Bethe-ansatz. [2,3] Although these theories, after analytic continuation to an Euclidean signature, have no topological sector, π 2 [SU(N )] = 0, they do posses non-perturbative saddle configurations known as unitons introduced in the seminal work of Uhlenbeck [14]. A uniton is a finite action field configuration solving the Euclidean equations of motion which are a harmonic condition, and obey, modulo a constant rotation, Let us specialise to the case of the G = SU(2). We adopt a Hopf parametrisation for the group element, with ranges θ ∈ [0, π], φ 1 ∈ [0, π] and φ 2 ∈ [0, 2π], such that the PCM (2.1) defines a σ-model whose target space is a round S 3 with metric, ds 2 = 1 t dθ 2 + cos 2 θdφ 2 1 + sin 2 θdφ 2 2 . (2.5) The uniton U : R 2 → SU(2) is specified by a holomorphic function f (z) (z = t + ix) , When the function f (z) is a polynomial of degree k we speak of the k-uniton and we find that the uniton action is, It may seem counter intuitive to have such a quantised action but if we write the uniton in Hopf coordinates, it becomes evident that we are dealing not with general maps into the whole of S 3 but rather with maps that are varying over an S 2 since the uniton can be seen as the embedding of a CP 1 lump into SU (2). More precisely the quantisation of the action can be related to a non-contractible loop in field space, i.e. F = {g : S 2 → S 3 } has π 1 (F ) = π 3 (SU(2)) = Z, see the beautiful orange book by Manton and Sutcliffe [47].

The η-deformed PCM generalities
An interesting class of theories called variously as Yang-Baxter σ-models and η-deformations (we prefer the later term) has been obtained by Klimcik [26,27] as integrable deformations of the PCM. To introduce these we require an operator R acting on the algebra g of a compact bosonic group G, that solves the modified Yang-Baxter equation, 6 Given a solution R of the modified YB equation we can consider the general action, This is integrable -at least classically. Lets look at the equations of motion which, with the introduction of, can be written as, In addition we have the Bianchi identity, The classical integrability is then established [27] since the equations of motion and Bianchi identity follow as a flatness conditions for a spectral parameter dependent Lax connection L(µ), (2.14) In the simplest rank 1 example, G = SU(2), there is a unique R that acts on generators normalised as T i = 1 √ 2 σ i by, 15) and the matrix entering the deformed PCM is, (2.16) Then the deformed action is simply,

JHEP07(2016)088
As a σ-model we have a target space whose metric is a squashed three-sphere equipped with a pure gauge NS two-form, Since the B-field is pure gauge in this case there is no danger in allowing η to become pure imaginary. In this case one see that η 2 = −1 corresponds to the O(3) sigma model (the action will develop a gauge invariance allowing one to reduce the degrees of freedom to just an S 2 's worth), whilst η = 0 is the SU(2) PCM, or equivalently, the O(4) sigma-model. This interpolating theory was used in the classic work of Polyakov and Wiegmann [10].
Using TBA, the work [48] established the exact mass gap of the theory. The classical integrability of the σ-model on the squashed three-sphere was established long ago [49] and an understanding of the symmetries as q-deformations was developed in [50]. Quantum mechanically, of the two couplings η and t, there is an RG invariant parameter ζ = t η and a single RG flow equation chosen between, where c G is the quadratic Casimir in the adjoint. A more convenient combination, that is closer to an overall coupling that sits outside the action eq. (2.10), is given by, in terms of which ζ is given as, The solution in parametric form (for ζ > 0) is, (2.23) In figure 1 we plot the phase portrait of the RG flow of the couplings g 2 and η 2 allowing η 2 to take both positive and negative values. For η 2 ≤ 0 the RG invariant quantity ζ 2 ≤ 0 and the theory is asymptotically free [48] in the coupling g 2 with the coupling η 2 flowing to −1 in the UV. There are two fixed points, one at η 2 = −1 which corresponds to the O(3) σ-model, and the second one for η 2 = 0 which is the undeformed PCM. Note that the η = 0 fixed line is "unstable" under this deformation, any non-zero ζ will automatically make us flow away from it. Although the figure illustrates also the flow for η 2 < −1 one should keep in mind that the target space metric will not have a positive definite signature. For the case of η 2 > 0 we have that the RG invariant quantity ζ 2 > 0 and the situation is different; the theory in the IR is strongly coupled in g 2 as before however as we approach the UV this coupling will decrease to a minimum value g 2 min = 2ζ which occurs when η = 1 JHEP07(2016)088 Figure 1. A phase portrait of the RG evolution of the couplings g 2 (horizontal) and η 2 (vertical). The arrows point towards the UV and the colouring is such that from blue to red the absolute value of the RG invariant parameter ζ 2 increases in magnitude (ζ 2 > 0 for η 2 > 0). at a scale log µ/µ 0 ≈ 1/ζ. This minimum value of g 2 is not a fixed point of the RG flow since the running of η 2 continues at this point and then for higher energies the theory becomes once again strongly coupled in g 2 . At large energies, for ζ > 0 target space of deformed σ-model will again become extremely curved rendering the σ-model description incomplete. Note however that the ratio between the two strong coupling scales Λ UV and Λ IR is parametrically large Λ 2 UV /Λ 2 IR = e π/(2ζ) , so by picking a very small RG invariant ζ we can obtain an RG flow for which the coupling g 2 ∼ 2ζ remains small and η 2 ∼ O(1) for a long RG flow time thus justifying later on our weak coupling expansion g 1 and η fixed.
Our philosophy here is that we will view the deformed σ-model being an effective theory up to some UV scale. For all values of ζ, the Ricci scalar of the target space of deformed σ-model is diverging in the IR indicating that non-perturbative effects will take hold. In general our goal in the resurgence analysis that follows is to understand this IR behaviour and thus the issue of the UV behaviour, and its resolution, will not have an impact on our analysis. We will use resurgence for the weak coupling expansion in g 1 while η will be a fixed, possibly complex, number. It is very likely that, to make a connection (i.e. the Borel flow mentioned in the Introduction) between the Borel plane of the QFT and that one of the reduced model that we will shortly introduce, we will have to perform some appropriate double scaling limit in g and η or possibly fix η to be close to 1 or i modulo some small, order g corrections. We will study this precise connection in future works.

Deformed-uniton solution
We would like to find uniton solutions, U η of the SU(2) η-deformed PCM 7 i.e. maps that obey the Euclidean equations of motion and obey, To ensure that the first of these conditions is solved let us make an assumption that in Hopf coordinates the solution for φ 1 and φ 2 is the same as in the undeformed case given in eq. (2.8). This is consistent; upon plugging this into the full second order equation of motions in the deformed theory one finds that the equations of motion for φ i reduce to an η-independent condition on θ, which is simply solved by θ = θ(|f | 2 ). The second order equations of motion for θ(|f | 2 ) now dramatically simplify to, This has a solution that matches to the undeformed uniton in the limit η → 0 given by, We can also rewrite equation (2.26) by changing variables to x = log |f | 2 , obtaining, which is precisely the equation of motion of a particular quantum mechanical model that we will shortly study in detail.
Although the full uniton is not a BPS configuration, note that the solution for the function θ obeys a rather suggestive 1 st order equation, Indeed this equation arises from a Bogomolny rewriting of the action of the η-deformed PCM evaluated on the ansatz θ = θ(|f | 2 ) and φ i given by (2.8). In summary we have a deformed uniton solution, Substitution of the uniton ansatz into the action (2.17) making use of the pseudo-BPS condition eq. (2.29) shows that, In [27] a bijection was set up between classical solutions of the undeformed and η-deformed PCM.
However, this argument relied on Minkowski signature and the solutions we present in Euclidean do not immediately follow. We thank C. Klimcik for correspondence on this. from which it is easy to see, after changing variables w = f (z) and making further use of eq. (2.29) that the action for a k-uniton (i.e. for f (z) a polynomial of degree k) evaluates to,

JHEP07(2016)088
where, will correspond later to the action of a quantum mechanical instanton (i.e. the fracton). Note that as η → 0 the uniton action S η [U η ] reduces to the undeformed case (2.7) and the instanton action S I → 2. The minimal uniton solution is obtained by plugging into eq. (2.30), where λ i ∈ C are the moduli of the solution, with λ 2 = 0. The modulus λ 1 can be thought as the centre of the uniton while λ 2 is a size modulus. The real uniton on R 2 behaves like a lump and can be best visualised by plotting its Lagrangian density, as shown in figure 2. The effect of the deformation is rather mild, as we increase η ∈ R the Lagrangian density becomes slightly punched out. We can easily consider the case η = iη R with η R ∈ R just by analytical continuation of the deformed uniton solution (2.30) and the uniton action (2.32) remains real. As we increase η R , the action density for the deformed uniton in this case steepens and its support gets smaller as we can see in figure 2. In both cases, η real and purely imaginary, we can vary over the moduli space of λ i and no fractionalisation occurs, the lump remains coherent.

Complex uniton solution
Conventionally one would consider only saddle points of PCM for which the field g is indeed an element of G. However, there is by now some substantial evidence [6,9,[51][52][53][54] that suggests that complex field configurations could play an important role in determining the resurgent structure of the theory and that we should also look for finite action solutions to the equations of motion when g ∈ G C . To this end, in the η-deformed theory we can search for solutions in which, as in the uniton, φ i are real and obey (2.8) but θ(|f | 2 ) is allowed to become complex but still satisfy the equation of motion eq. (2.26). Remarkably there is indeed such a solution, which corresponds to a complex-uniton U c η ∈ SL(2, C) given by, We call this solution to the complexified η-deformed PCM equation of motion a complex-uniton because the matrix U c η lives in the natural complexification SL(2, C) of the original target space SU (2).
In this case, the action when f (z) has degree k is given by, where, will correspond later to a complex-instanton solution in the quantum mechanically reduced deformed PCM. The reason quite evidently that these solutions have not been, to the best of our knowledge, considered before in the literature is that in the limit η → 0 the action diverges and can not give any contributions to the PCM path integral. This is readily seen by looking at the Lagrangian density for η → 0 and f (z) = λz (the constant term in f (z) just corresponds to a shift in the centre of the uniton and can be set to zero here) which reads, and evidently diverges along a circle in the complex plane where the denominator vanishes. The effect of the η-deformation is to smooth out this singularity giving a finite value to the action. The minimal complex-uniton is obtained from eq. (2.37) by choosing,

JHEP07(2016)088
where λ i ∈ C are the moduli of the solution.
In contrast, for the complex-uniton on R 2 the effect of the η-deformation is drastic, it has the effect of smoothing out a singular Lagrangian density to give a finite action configuration for η > 0 (see figure 3).
Note also that if we analytically continue the complex-uniton solution (2.37) to η = iη R with η R ∈ R we obtain a Lagrangian density which has a non-integrable singularity. The reason is that our solution (2.36) becomes singular when analytically continued to purely imaginary η. The complex-uniton action in this regime can be obtained via analytic continuation of (2.38) from arg(η) = π/2 ± (see later on section 3.3 as well as [54]) obtaining two types of actions, (2.43)

JHEP07(2016)088
There is another way of smoothening these field configuration for η purely imaginary. If we allow, as we should, all the field configuration to become complex, then we can replace in (2.27)-(2.36) |f (z)| 2 → Z 0 |f (z)| 2 with Z 0 ∈ C. This corresponds to the choice of a complex centre x 0 = log Z 0 for the instanton or complex-instanton solution to (2.28). The modulus of Z 0 acts simply as a scale modulus but its phase allow us to move the singularities of the Lagrangian density away from the domain of the spatial integration. In particular the action density for the complex-uniton for Z 0 = e iα in the undeformed theory becomes, to compare against (2.40). For α ∈ (0, 2π) the above action density has no singularities over z ∈ C and it integrates precisely to the uniton action but we stress again that the corresponding complex-uniton is a genuine SL(2, C) matrix even in the undeformed theory. This issue of complex moduli is at the very core of the Morse-Picard-Lefschetz approach to the path integral of quantum field theories as discussed in section 5.2.

Twisted spatial compactification
Recent works [17,18,[55][56][57] have emphasised the idea of adiabatic continuity as a useful technique in the study of the IR physics of QCD-like gauge theories on R 4 (or R 2 ). The idea is to place the theory on R 3 ×S 1 L (or R×S 1 L ) such that when the radius L is sufficiently small the theory enters a weakly coupled domain.
Of course, the physics in this regime may be quite different to that of the theory on R 4 ; for a thermal compactification we have a deconfined phase at small L and a confined phase at large L. In such a case calculations performed at small L will not tell us a great deal about the long distance physics relevant for the theory at large L. If instead one takes a spatial compactification, the small L behaviour can be adiabatically connected to the large L-regime, note that the fermions need a spatial compactification as well to prevent any Hagedorn instability [58]. For 2-d theories on R × S 1 L a thermal compactification would only give rise to a cross-over between the deconfined regime at small L and the confined one at large L, it is only in the N = ∞ limit that we obtain a true phase transition as we dial the radius of the S 1 .
This idea has been applied to the PCM eq. (1.6), which as discussed in the introduction, can serve in many ways as a good prototype for QCD but in the easier scenario of twodimensions. In [9] it was studied in detail how the properties of compactification of the SU(N ) PCM on R × S 1 L are related to a choice of boundary conditions, in which ψ denote any fermions in the game and the − sign corresponds to anti-periodic thermal boundary conditions and the + sign to periodic spatial boundary conditions. These boundary conditions can equally well be thought of as turning on background gauge fields

JHEP07(2016)088
for H L × H R inside the global SU(N ) L × SU(N ) R symmetry. Two conditions are required of the choice of H L and H R for adiabaticity. First is that the free energy F should be at a stationary point as we vary the boundary conditions and second that F /N 2 → 0 for any L. The eigenvalues of the axial H A = 1 L (H L − H R ) obey a tree level effective potential which is extremised only when H A = 0 and thus attention can be restricted to vectorial gauge field for H V = 1 L (H L + H R ). By integrating out the KK-modes one finds an effective potential V (Ω) for the Wilson loop Ω = exp i dxH V = exp iLH V which depends crucially on the sign choice of fermionic boundary conditions. For thermal boundary conditions this potential is minimised by, leading to a free energy that does not satisfy the adiabaticity requirements as expected. For QCD-like theories on R 3 × S 1 L , this holonomy would correspond to a VEV for the Polyakov line that breaks the Z N centre symmetry.
Conversely, the minimisation of the effective potential for spatial compactifications is achieved by, For QCD-like theories on R 3 × S 1 L this would correspond to a VEV for the Polyakov line that preserves the Z N centre symmetry since tr(Ω) = 0.
For SU(2) this has the simple interpretation of turning on a background field for the U(1) V , In the case at hand, the η-deformations the situation is slightly subtle. For −1 < η 2 < 0, as discussed previously, the theory is asymptotically free and the story outlined above transfers directly. Indeed one can calculate the effective potential in this case and establish that the same boundary conditions, eq. (3.4), are required. For the case of η 2 > 0 the theory on R 2 is both IR and UV strongly coupled. We can, and shall, of course still study the spatial compactification but should be appropriately more cautious in interpreting the implications of the results of the reduced theory to that on R 2 . Nonetheless what we shall find is a rather consistent story that is valid for both −1 < η 2 < 0 and η 2 > 0.
Furthermore, since we are not working at N = ∞, our choice of twisted boundary condition on R × S 1 L is by no means a guarantee of adiabatic continuity. There is no phase transition as we shrink the radius of S 1 , nonetheless, as noted in [59], this choice of maximal twist is very particular even at finite N hinting towards volume independence also when N = 2. These caveats aside, we will study the compactification of the η-deformed SU(2) PCM in the presence of the maximal twist (3.4) and leave for further studies the precise connection with the theory on R 2 .
Rather than working with twisted boundary conditions it is expedient to work instead with periodic fields,g

JHEP07(2016)088
Following the discussion above, we will turn on only a vectorial twist, where ξ = π/2 will correspond to the maximal twist (3.4). We calculate the Lagrangian for this periodic field given in terms of right-invariant forms as, using the coordinates x ± = 1 2 (t ± ix) resulting in, In the above, we have introduced the coupling g 2 = t(1 + η 2 ) and we have dropped the pure gauge B-field of the η-deformation (the term in the action is a total derivative and imaginary in Euclidean signature) and note that the final imaginary term coming from the twisted boundary conditions is also a total derivative and will be henceforth discarded. The U(1) L -and U(1) R -currents, corresponding to the symmetry g → e i L σ 3 ge −i R σ 3 of the untwisted theory, are given by, (3.9) Reflecting the isometries of the squashed sphere target space, we note that the left acting symmetry is insensitive to the deformation while the right acting symmetry is modified.
We further see that the twisting used is equivalent to turning on a spatial component of the background gauge field of the vectorial action since, We will use now a small-L effective quantum mechanics by performing a KK-reduction and dropping all the x dependence from the fields. Clearly such a reduction will be forgetful regarding the UV-renormalon singularities, which can be only extracted from the microscopic theory. An important subtlety in performing the KK-reduction is that it is not quite true that all KK-momentum carrying states can be discarded. The issue is that states that carry non-zero winding number can contribute to the low-energy dynamics on the same footing as states that carry zero winding number, as already noted in [9].

JHEP07(2016)088
In particular we have to consider when φ 2 carries a non-trivially winding around the S 1 , with n ∈ Z.
Taking into account the factor of L that comes from performing the integral around the S 1 one finds, after rotating back to Lorentzian signature, an Hamiltonian of the form, where the momenta P I = {p θ , p φ 1 , p φ 2 } are coupled via the inverse metric on the squashed sphere namely, As explained above we see that, after picking the maximal twist ξ = π/2, the states with zero winding number for φ 2 will have exactly the same energy as the states with winding number −1.
We are interested in the ground state of this system and to this end we employ the Born-Oppenheimer approximation which separates the dynamics of the light degrees of freedom p θ with the heavy ones p φ i as in [9]. Remembering the subtleties about winding modes we can also set n = 0 from now on and choose the maximal twist ξ = π/2. Thus we arrive at a one-dimensional quantum mechanics of interest given by eq. (3.12) with p φ i = 0, n = 0 and ξ = π/2, (3.14) We can rescale the time variable to remove the energy scale 1/L and some numerical factor to put the above Hamiltonian in the form, with the corresponding Lagrangian given by, The corresponding Schrödinger equation can be put in the form, with the parameters,

JHEP07(2016)088
This differential equation is called the three-term Whittaker-Hill equation and some of its properties are known. 8 Although it will not be directly relevant for what follows, it is notable that whilst this is not an integrable QM (in the sense that its spectrum can't be given in closed form in terms of known functions), it is an example of a quasiexactly-solvable QM (for a recent comprehensive treatment [61]) in which some -but not all -eigenstates for particular values of the parameters can be found exactly. This Whittaker-Hill equation also emerges in certain inflationary scenarios as describing entropy fluctuations during reheating [62] (this paper also contains a useful exploration of stability properties of the equation).
We notice that our Hamiltonian can be directly related (for η 2 < 0) to the supersymmetric quantum mechanics studied in [54,63] with a superpotential of the form W (θ) = cos(θ) after the fermions are integrated out.
That our deformed bosonic theory reduces in one-dimension to the supersymmetric completion of the undeformed theory seems, at first sight, rather peculiar. One anticipates that the twisted reduction of the two-dimensional supersymmetric Principal Chiral Model would be the natural origin for the supersymmetric quantum mechanics studied in [54,63] and that this would be a rather different theory from the bosonic η-deformed PCM. However there is one sense in which the η-deformation and the supersymmetric deformation are similar; they both modify the S-matrix of the theory by means of a quantum group. Indeed the supersymmetric part of the S-matrix of the SUSY PCM for spinors in the O(4) model is identical to the soliton S-matrix of the supersymmetric sine-Gordon theory evaluated at a value of coupling such that scattering is reflectionless [64,65] . This in turn is related to an affine SU(2) quantum group with parameter q a root of unity. This was extended in [66] to the SU(N ) PCM wherein the supersymmetric part of the S-matrix is a q root unity affine quantum group deformation of SU(N ). It would be interesting to try and clarify this possible relation at the level of Lagrangians perhaps by directly integrating out the fermions from the PCM.
As expected, for η → 0 the QM defined by eq. (3.17) reduces to the well known Mathieu equation studied in this PCM context in [9,16]. More generally the Mathieu equation has been carefully studied in a number of works as an example of a QM where degenerate minimum give rise to a banded spectrum [67,68]. This system also plays a key role in the study of 4-dimensional N = 2 supersymmetric pure SU(2) gauge theory and its wall crossing phenomenon [69,70].
Before proceeding, let us make a small digression concerning the reduced quantum mechanics eq. (3.12). If instead of adopting Euler angles we use R 4 embedding coordinates to parametrise the group element as,

JHEP07(2016)088
then the Lagrangian corresponding to eq. (3.12) can be expressed compactly as, For η 2 = 0 this quantum mechanics on an S 3 with a quadratic potential in the embedding coordinates has been well studied under the name of the C. Neumann model [71] and is known to be integrable both classically [14,72] and quantum mechanically [73]. 9 It seems quite likely then, given the integrable 2-d origin of eq. (3.20) that the full η-deformed quantum mechanics is also integrable and will be pursued elsewhere [74]. Intriguingly an integrable deformation of the Neumann model has already arisen in relation to spinning strings in η-deformed AdS 5 × S 5 in [75] -though the structure there seem more involved than the simple Lagrangian of eq. (3.20).

Unitons on R × S 1
To put the minimal unitons and the minimal complex-unitons on R × S 1 L , and to take care of the twisted boundary conditions we use the results of [76,77] and we simply make the replacement, with z = t + ix such that, where H V is the maximal twist matrix of (3.4).
Here something quite different happens as we vary over the moduli-space, both the real and complex-uniton solutions break-up or fractionalize. In the regime λ 1 1 λ 2 there is no fractionalisation however when λ 1 ∼ λ 2 1, regardless of the value of η these semi-classical solutions break up into widely separated lumps of Lagrangian density. This is illustrated in figure 4 for the real uniton and figures 5 and 6 for the real and imaginary contributions to the Lagrangian density of the complex-uniton.
As well as seeing the fractionalisation in the Lagrangian density, we can see it directly at the level of field space. Let us focus on the field θ(x, t) defined via eq. (2.27) for the real uniton and (2.36) for the complex-uniton but with the function f (z) appropriate for R × S 1 L with for concreteness real moduli. For λ 1 λ 2 1, which holds in the fractionalized regime, the derivative with respect to x -i.e. the dependence in the S 1 L -is suppressed and θ(x, t) ∼ θ(t) effectively truncates to the zero-KK mode sector. The field φ 2 given by equation (2.8) does instead depend on x and interpolates between a configuration with 0 winding and a configuration with winding −1. For the real uniton this gives a correlated fracton-anti-KK-fracton profile so schematically we write the fractionalisation as U → [F 0 F −1 ] where the subscript are to remind that the field φ 2 has different winding numbers.
As we will see later on, the action of each fracton event can be computed from our quantum mechanic model and it is precisely equal to (2.33) so that the uniton, even in the JHEP07(2016)088 deformed theory, is indeed composed by a fracton and a KK-anti-fracton and its action is given by (2.32).
For the complex-uniton, we find that the profile of θ(t) is such that the real part remains constant however the imaginary part is non-constant and it fractionalizes into two equal constituents. As for the uniton, even in this case the field φ 2 is given by equation (2.8) so it does depend on x and interpolates between a configuration with 0 winding and a configuration with winding −1. For these reasons we denote the fractionalisation of the complex-uniton as U C → [CF 0 CF −1 ], i.e. the complex-uniton is formed by a complexfracton and a KK-complex-fracton. We will shortly show that the action of a complexfracton can be computed from our quantum mechanical model and it is precisely equal to (2.39) explaining why the complex-uniton action is given by (2.32).
The profiles of θ for both these cases are plotted in figure. 7.

Non-perturbative saddles in the reduced quantum mechanics
We can understand the origin of the fracton and the complex-fracton from the reduced quantum mechanics. Let us rewrite here the Hamiltonian of interest after a rescaling 10 becomes H = p 2 θ + sin 2 (θ) 1 + η 2 sin 2 (θ) = p 2 θ + V , whose variation gives the second order equation of motion To classify all the non-perturbative solutions to the equation of motion we first have to find the complex critical points of the potential V (θ) = sin 2 (θ) (1 + η 2 sin 2 (θ)) (3.26) in the strip Re θ ∈ [0, π] of the complex plane θ ∈ C. We can easily solve for V (θ) = 0 for arbitrary η ∈ C with |η| ≤ 1, finding generically 7 critical points: where θ cr ∈ C is given by, (3.28) If we start with η = 0 we immediately see that θ cr → ∞ so we are only left with the three critical points {0, π/2, π} of the undeformed Mathieu case, as discussed in [9]. As we turn on η = 0 with arg(η) = 0 we see that θ cr becomes finite and purely imaginary, of the 7 critical points 3 are real while 4 are complex, see figure 8.
If we increase the argument of η, θ cr acquires a real part that becomes exactly equal to π/2 for η = iη R with 0 < η R < 1/ √ 2 so generically we have 7 critical points, 3 real and 4 complex. When η is purely imaginary and in modulus less than 1/ √ 2 we have that θ cr = π −θ cr so we are left with only 5 critical points, 3 real and 2 imaginary, see figure 8. As η approaches i/ √ 2 the two critical points θ cr andθ cr approach one another, and for η = i/ √ 2 we have that θ cr = π/2 so we are left with only 3 real critical points, see figure 9. Im Θ η = 2 3 e iπ/2 Figure 8. The red dots correspond to the critical points of the potential (3.26) for different values of η. As we increase the argument of η while keeping its modulus fixed, the two critical points θ cr and π −θ cr come closer together moving on the blue curve, whileθ cr and π − θ cr move closer along the green curve. For arg(η) = π/2 we are left with 5 critical points. If we keep increasing the modulus of η while keeping its argument to π/2 we see that θ cr remains real and moves on the real axis towards the origin, while θ − θ cr moves on the real axis towards π. In this regime we have once again 5 critical points all reals, see figure 9. For η purely imaginary and in modulus larger than 1/ √ 2, the critical points can be easily understood by plotting the potential V (θ) for real θ as in figure 10. In this regime there are two different types of instanton events corresponding to tunnelling between the tall barrier or the short one.
The most general solution to the equation of motion (3.25) can be obtained by making use of the Weierstrass elliptic functions [54] but we will only be interested in finite action solutions to (3.25). The first of such solutions is a typical instanton event interpolating between θ → 0 as t → −∞ and θ → π as t → +∞, with profile given by, where t 0 ∈ R denotes the instanton centre zero-mode. The action can be computed, and it is precisely equal to the fracton action (2.33). In the η → 0 limit the instanton reduces to the standard instanton of the reduced PCM case of [9] and its action stays finite S I → 2.
Clearly an anti-instanton solution can also be found and it has exactly the same action while its profile is given by, For η real the instanton and the anti-instanton are real solutions (if the modulus t 0 is chosen to be real) while for generic complex η the profile θ I and the action S I have a nonzero imaginary part. As η becomes purely imaginary, η = iη R with η R ∈ R, the function θ I becomes real once again (for t 0 ∈ R) and the instanton action can be rewritten as, which is clearly real.

JHEP07(2016)088
Note that for any η ∈ C both the instanton and the anti-instanton interpolates between different vacua, namely θ = 0 and θ = π. This tells us that only correlated event IĪ interpolating between the same vacuum θ = 0 can (and will) communicate with standard perturbation theory around the vacuum θ = 0.
Perhaps more interesting are the complex-instanton solutions, called complex-bions in [54], together with {π − θ CI , θ CI , π − θ CI }, all with action, As we see in figure 11, for η real θ CI (t) is purely imaginary (for t 0 ∈ R). The complexinstanton solution tunnels from θ → 0 as t → −∞ to θ → θ cr given in (3.28) and then tunnels back to the same vacuum θ → 0 as t → +∞, and similarly for all the other complex-instanton solutions written above. Note that despite the complex-instanton being a solution living in the complexification of the field space its action is real. In the limit η → 0 the complex-instanton θ CI (t) becomes singular for t = t 0 and its action (3.34) diverges but this only happens if we keep the centre modulus real! If we allow ourself to consider t 0 ∈ C, and we should since the field space should be complexified in its entirety, we can see the presence of this additional saddle even in the η = 0 case (see the discussion on complex-unitons of section 2).
If we increase the argument of η the solution θ CI (t) will stop being purely imaginary and the action 3.34 will become complex, but the complex-instanton will remain a configuration interpolating between the vacuum θ → 0 for t → −∞, then off to the complex critical point θ cr and back again for t → +∞ to the same vacuum θ → 0.
When η = iη R with 0 < η R < 1/ √ 2 the complex-instanton (3.33) develops two nonintegrable singularities as displayed in figure 11. To obtain a smooth profile function we can either reach η purely imaginary from arg(η) = π/2 ± or we can consider a complex centre t 0 modulus. As explained in detail in [54], when η = iη R with 0 < η R < 1/ √ 2 the branch cut of arctanh produces two types of complex-instanton solutions with actions: one complex conjugate of the other. Note that for any η ∈ C the complex-instanton interpolates between the same vacuum. This tells us that a single complex-instanton event CI can (and will) communicate with standard perturbation theory around the vacuum θ = 0.

Perturbation theory in the reduced quantum mechanics
Now that we have a better knowledge of the non-perturbative sector of our η-deformed PCM and of the reduced quantum mechanics, we want to understand if these instanton JHEP07(2016)088 (fracton) and complex-instanton (complex-fracton) events play any role in the semi-classical expansion of physical quantities. To achieve this we consider a very simple quantity: the ground state energy of the quantum mechanics (3.15) obtained from the maximally twisted reduction of our QFT.

WKB Perturbation theory and Borel singularities
We will compute the ground state energy in perturbation theory and we will see how this will be insufficient to obtain an unambiguous result because we will be missing some crucial non-perturbative effect. To compute the large orders perturbative coefficients, we found that using a standard WKB approach was numerically faster than using uniform WKB.
We consider the time independent Schrödinger equation associated with the Hamiltonian (3.23), with E the rescaled energy and V (x) = sin x + η 2 sin 2 x. We can use a WKB ansatz for the wave-function, reducing the time-independent Schrödinger equation to a Riccati equation for S(x): We are interested in studying the resurgence properties of the energy levels of the quantum mechanical system above in the limit g → 0 with η ∈ C fixed. In particular for the ground state energy we can make the ansatz:

JHEP07(2016)088
with E 0 = 1, while the coefficients E n will be some factorially growing, η dependent, complex numbers. The auxiliary functions S n (x) can be expanded as for some complex, η dependent, coefficients S n,m , since we have assumed that the functions S n (x) are regular 11 at x = 0. We plugged the ansatz (4.5) in (4.3), expanded as a power series in g 2 , and imposed that the Riccati equation is satisfied order by order in g 2 . These are the first few terms in the small g 2 expansion of the Riccati equation At this point we can use the small x expansion (4.6) and solve for the unknown coefficients {E n , S n,m } order by order in g and in x.
Note that this procedure is completely independent from the argument of η, in particular we can find the perturbative coefficients {E n } for the energy levels of the quantum mechanical system (4.1) for all values of η ∈ C. Obviously, the Hamiltonian associated with (4.1) is an Hermitian operator only for arg(η) = 0 or arg(η) = π/2. Nonetheless, the study of the large order behaviour of the coefficients {E n } will be extremely interesting if we allow arg(η) to vary in [0, π/2], despite the quantum mechanical model (4.1) not describing a physical system for generic η ∈ C.
We can compare the first few perturbative coefficients E n for different values of η obtained with standard WKB analysis as described above, against the coefficients (4.20), for B = 1/2 (i.e. the ground state), that we will shortly obtain using uniform WKB approach and the two match perfectly for all values of η ∈ C.
As anticipated, we expect these coefficients {E n } to grow factorially with n, hence all the series in (4.5) are only asymptotic. The first thing we want to understand is the structure of non-perturbative contributions that one must add to our ansatz (4.5) to construct a proper transseries of the form (1.5). 11 For excited states this ansatz has to be changed since the wave-function ψ(x) must have some nodes. In particular if we look for the N -th energy level, with N = 2ν +1 and ν a half integer, the energy will be E (ν) = 1+2ν+O(g 2 ) while the WKB wave-function must take the form S(x) = S   We first fix η ∈ C and then generate, using the procedure outlined above, the first N p = 150 perturbative coefficients {E n }, while keeping 200 significant precision digits. We Borel transform the ground state energy, truncated to order 2N p : E (0) Since we only kept a finite number N p of terms, this Borel transform is an entire function in the complex t-plane. For this reason we study the pole structure of the Borel-Padé approximants (see e.g. [84]), , (4.9) where the polynomials P n (t) and Q m (t) can be constructed from B T (t) and have degrees n and m respectively, with N p = n+m. The idea of the Borel-Padé technique is certainly not new and its applications to strongly coupled particle physics were summer school material in 1970 [78]. 12 The use of Borel-Padé approximants are supported by their strong convergence properties. For instance taking a diagonal approximant m = n, we can compare the first predicted value for the series with the known result as a function of n and find an exponential convergence, This is illustrated by the log plot of figure 12 in which we find σ ∼ 1.69 for η = 0.2. Similar exponential efficacy of the Borel-Padé method was conjectured in [79] and the method was used to predict quite successfully the then unknown 4-loop contribution to the βfunction of QCD. More recently the exponential convergence of the Borel-Padé method was verified against explicit calculations using localisation techniques of correlation functions of Coulomb branch operators in four-dimensional N = 2 superconformal field theories in [80]. The Borel-Padé approximants, being rational functions, will have poles in the complex t-plane. As we increase the number of coefficients N p we can obtain better and better Borel-Padé approximants and their pole structure is supposed to converge to the branchcut structure of the complete Borel transform, We can see from figures 13-14 that the singularity structures of the Borel-Padé approximants change dramatically as we dial η ∈ C.
For arg(η) = 0, the only two Stokes directions, i.e. singular directions, for B(t) are θ = arg(t) = 0 and θ = π. In this regime we know that there are two non-perturbative field configurations which can mix with perturbation theory, an instanton-anti-instanton configuration with action (3.30), and a complex-instanton configuration with action (3.34), Note that for η ∈ R both actions are real, S IĪ > 0 while S CI < 0. We see in figure 13 that the Borel transform has two branch-cuts starting precisely at t = S IĪ and t = S CI . If we try to go back to the undeformed case, η = 0, we can see that the branch cut along the Stokes lines θ = π tries to disappear, moving towards −∞, since S CI ∼ − π 2η + O(η 0 ) as η → 0, while the branch cut along θ = 0 starts precisely at t = S IĪ ∼ 4 + O(η), as expected from the undeformed Principal Chiral Model analysis carried out in [9] and shown in figure 13.
As we increase the argument of η, we see in figure 14 that the two Stokes lines move from θ = 0, π to the directions θ = arg(S IĪ ) and θ = arg(S CI ) and the two branch cuts start precisely at t = S IĪ and t = S CI .
When we reach arg(η) = π/2, i.e. η = iη R with η R ∈ R and |η R | < 1/ √ 2, the Borel transform of the ground state energy has three branch cut. The reason is that the complex-instanton action (3.34) has a branch cut discontinuity in the complex η plane running precisely between −i and i. This gives rise to a third type of non-perturbative saddle (see the discussion in section 3.3) which can mix with perturbation theory. The actions (3.32)-(3.35) of the three relevant non-perturbative field configurations for η = iη R with η R ∈ R and |η R | < 1/ √ 2 are: 14) Note that S IĪ is once again real, while S CI and S CI are complex and furthermore they are complex conjugates to one another: S CI = S CI . We have thus three Stokes lines at θ = 0 and θ = ±arg(S CI ), and the branch cuts of the Borel transform start precisely at t = S IĪ , t = S CI and t = S CI respectively, as shown in figure 14. Finally, if we crank up the modulus of η such that |η| ≥ 1/ √ 2, while keeping arg(η) = π/2, we know from the analysis of section 3.3 that the complex-instantons cease to exist and we only have two types of non-perturbative saddles, both with real actions and the only Stokes line is for θ = 0.

Uniform WKB
To analyse the perturbative expansion in detail we follow the uniform WKB method that is well-suited to potentials with degenerate minima [67]. This consists of two steps, first one is to solve the Schrödinger equation about a given minimum perturbatively and then one invokes a global boundary condition that will fully determine the energy levels.
Since the potential in (3.23) is locally harmonic, a sensible ansatz for the wavefunction is, in which D ν (z) is a parabolic cylinder function and ν = B − 1 2 is parameter that will be determined via the global boundary condition but is exponentially close to an integer. Substitution of this ansatz into the Schrödinger equation yields a non-linear equation, which is studied perturbatively with, The first few terms of the energy series are for the potential of equation (3.23) are given by, (4.20) We can set η = 0 and read the coefficients of the undeformed case already computed in [67] and, for B = 1/2 we can directly compare against the numerical coefficients (4.7) found using standard WKB for different values of η with perfect agreement. We compute the wave-function in perturbation theory and at zero'th order the wavefunction is given by,

JHEP07(2016)088
The first and second order wave-functions are given in terms of u 0 by, 22) in which the function χ is quite unpleasant to write in full but is obtained by solving, with the constant of integration such that u 2 (0) = 0. The values of these wave-functions at the mid-point θ = π 2 will presently become useful and they are given by, (4.24) in which we recall that S I is precisely the instanton action defined in equation (3.30).

Resurgence analysis
We are now in position to obtain a transseries ansatz for the ground state energy of the quantum mechanics (3.23) and obtain some novel predictions.
As explained in the Introduction and as we have just seen in practice, the Borel transform of the perturbative ground state energy (4.5) will generically have branch cuts along some Stokes directions arg(t) = θ. The two lateral resummations S θ+ (E pert ) and S θ− (E pert ) will differ from one another and furthermore their difference will be exponentially suppressed, with S θ ∈ C and arg(S θ ) = θ and Φ θ (g 2 ) a new function with an asymptotic expansion for g 2 ∼ 0. Thanks to a standard dispersion-like argument [81,82] we can relate the large order behaviour of the perturbative coefficients E n in (4.5) to all the discontinuities of the Borel transform B[E pert ]: where θ * runs over all the Stokes directions while the discontinuity is precisely related to the jump in the lateral resummations Disc θ (E pert (w)) = S θ+ (E pert ) − S θ− (E pert ).
Let us start with the η ∈ R case. A careful analysis of equation (5.2) together with our numerical studies of the perturbative coefficients E n for numerous values of η, reveals that the large order perturbative coefficients behave for n 1 as, where we omitted higher instanton contributions. In the above equation A(η) and B(η) correspond to the Stokes constants in the singular directions θ = 0 and π respectively, while S I and S CI are precisely the instanton and complex-instanton actions (3.31)-(3.34). The coefficients a For η = iη R with 0 < η R < 1/ √ 2 the situation is very similar and the large order perturbative coefficients behave for n 1 as, 4) where the actions S CI and S CI are precisely given by (4.15)-(4.16) and we omitted again higher instanton contributions. Note that in the above equation the Stokes constants A(η) and B(η), the instanton-anti-instanton perturbative coefficients a (i) IĪ (η), and the complexinstanton coefficients a (i) CI are exactly the same as in (5.3), simply evaluated at purely imaginary η = iη R . We have checked this also for generic complex values of η and it remains correct.
Making use of the uniform WKB approach 13 [67] we can compute the Stokes constant A(η) and the perturbative coefficients on top of the instanton-anti-instanton solution a (i) IĪ (η). We first write the energy transseries for the N th level,

5)
13 Roughly we need to use the resurgent asymptotic behaviour of the parabolic cylinder functions in the complex plane, we refer to [67] for all the technical details.

JHEP07(2016)088
where the asymptotic series E(B = N + 1 2 , g 2 ) is precisely the object we computed in (4.19) in terms of the coefficients (4.20). At leading non-perturbative order, the parabolic cylinder function index ν in (4.17) is exponentially close to the level number N and δν = ν − N as explained in [67].
As mentioned in section 4.2 to determine fully the energy levels we have to impose a global boundary condition. This global boundary condition is an extra constraint that the uniform WKB function (4.17) must satisfy and it arises from the Bloch phenomenon. The eigenfunctions of a periodic potential must satisfy: ψ(x + π) = e iθ B ψ(x), where θ B is the Bloch angle θ B ∈ [0, π] that labels states in a given band of the spectrum. We can compute δν by making use of this global boundary condition expressed by equation (63) of [67] and if we consider the lowest part of the band and set the Bloch angle to θ B = 0 we obtain, where u(x) is the uniform WKB wave function, solution to (4.18), evaluated at mid turning point. We can expand the wave-function at mid-turning point as an asymptotic power series in g and the first few coefficients are given in (4.24). At leading order we have, where S I is once again given by (3.31). By substituting (5.7) in our transseries (5.5) we see that the splitting of the lowest band is, as expected, a one instanton effect. By expanding the global boundary condition to order e −2S I /g 2 , we obtain an imaginary part of δν, where the sign ± is correlated with arg(g 2 ) ≷ 0. The leading imaginary part of the energy transseries (5.5) coming from the two-instanton sector, including the perturbative fluctuations around it, can be found from, (5.9) By imposing that this purely imaginary term in the transseries (5.5) is cancelled exactly by the discontinuity of the resummation of the perturbative expansion Disc 0 (E pert (g)) translates immediately into the large-order behaviour of the perturbative coefficients and we can read the Stokes constant A(η) and the first perturbative correction  [67]. A path integral derivation, as opposed to our approach based on the Schrödinger equation, of (5.11) would be an amazing check of the resurgent program. For the sine-Gordon model this agreement has been shown up to three loops [83].
Note that the uniform WKB is completely oblivious of the argument of η ∈ C so we use (5.10) and (5.11) for generic η. It is also interesting to notice that from the uniform WKB is not at all obvious how to extract the Stokes constant B(η) and the perturbative coefficients a (i) CI (η) associated with the complex-instanton. A superficial expansion of δν seems to contain only powers of exp(−S I /g 2 ) while we know from the singularities of B(E (pert) ) that we should also find terms of the form exp(−S CI /g 2 ). Clearly these terms must be hiding in the asymptotic nature of the uniform WKB wave-function expansion (4.19)-(4.24) but for the moment we do not know how to extract them. 14 Thanks to our study of the large order behaviour of the perturbative coefficients we predict CI (η) = 1 48 −29 + 95η 2 + 8 1 + η 2 .

(5.12)
A proper complexified path integral derivation of (5.12) would be an amazing check of the resurgent program.

Large order behaviour in quantum mechanics
Let us analyse more in detail the large order behaviour of the perturbative coefficients (5.3)-(5.4). Firstly, for η ∈ R, we have that |2S I | ≤ |S CI | for 0 ≤ η ≤ η c ≈ 0.274 so that the leading asymptotic behaviour of the perturbative energy coefficients (5.3) is given by, We plugged in the above equation the Stokes constant (5.10) and the first instanton coefficient (5.11) obtained via resurgence theory and we can see from figure 15 the perfect agreement with the numerical coefficient E n evaluated in section 4.1.
We can push it even further and consider the difference,  Figure 15. Ratio between the exact perturbative coefficients E n and their asymptotic large order behaviours E As n for different η. On the left plot η < η c and the instanton-anti-instanton action is the dominant effect, while on the right plot η > η c and the complex-instanton is dominant.
We applied Richardson extrapolation (see e.g. [84]) and obtained predictions for a IĪ (η) can be in principle computed by expanding (5.9) to higher order in g 2 making use of the global boundary condition, although it is necessary to know the uniform WKB wave-function u 3 (π/2) at mid turning point as in (4.24).
For η > η c we have that |2S I | > |S CI | so that the leading asymptotic behaviour of the perturbative energy coefficients (5.3) is, We substituted in the above equation the predicted (5.12) Stokes constant 15 and the first perturbative correction to the complex-instanton and we can see from figure 15 the perfect agreement with the numerical coefficient E n evaluated in section 4.1.
Turning to the complex sector, we can also push it further and consider the difference, .
CI (η) from the large order behaviour of the perturbative expansion. We use 150 perturbative coefficients and the 12th Richardson extrapolation. perturbative energy coefficients (5.4) is given by, We plugged in the above equation the Stokes constant (5.10) and the first instanton coefficient (5.11) obtained via resurgence theory and in figure 16 we see the perfect agreement with the numerical coefficient E n evaluated in section 4.1.
Once again we can consider the difference and after Richardson extrapolation we can obtain predictions for a Finally for η R >η c ≈ 0.402 we have that |2S I | > |S CI | = |S CI |. In this regime the leading asymptotic behaviour of the perturbative energy coefficients (5.4) is given by,  Figure 16. Ratio between the exact perturbative coefficients E n and their asymptotic large order behaviours E As n for different η = iη R . On the left plot the instanton-anti-instanton action is the dominant effect, while on the right plot the two complex-instantons are the competing saddles producing an oscillating behaviour.
From (4.15)-(4.16) we know that S CI = S CI so that for n 1 the two saddles will compete with one another producing an oscillating asymptotic behaviour on top of the factorial growth, matching beautifully the numerical coefficients E n as depicted in figure 16.

The role of complex saddles
As already realised long time ago by Bogomolny and Zinn-Justin [2,3] (see also [85]) and recently revisited in the context of Morse-Picard-Lefschetz theory [54,86], when we consider instanton-anti-instanton events we have to be careful on the integration over quasizero modes. In a path-integral approach we would be tempted to integrate over R for the separation τ between the instanton and the anti-instanton event but this would lead to a diverging contribution to the partition function. A more careful analysis suggests that the correct contour of integration, also called Lefschetz thimble, associated to the instantonanti-instanton event lives in a complexification of the field space. Depending on arg(g) ≷ 0 the correct contour of integration in the relevant direction of the complexified field space is given by τ ∈ R ± i π.
The integral over the correct thimble for the instanton-anti-instanton contribution produces precisely the correct imaginary and exponentially suppressed term to compensate for the jump in lateral resummation of the perturbative vacuum contribution [54,86]. Since we complexify the coupling constant to avoid Stokes lines it is only natural that we complexify also the field space to allow for more general field configurations. This is also suggesting that in principle all finite action solutions to the complexified equation of motions could contribute to the transseries expansion of any physical observables. This has been seen in the context of the analytic continuation of the path integral of JHEP07(2016)088 quantum mechanics in phase space [51] and the analytic continuation of Chern-Simons [6].
These two examples are rather peculiar since they have first order in time equation of motion. For field theories with second order in time eoms the only examples we are aware of is given by the analytic continuation of Liouville theory [52].
In this work we have presented a new example of quantum field theory with second order in time equation of motions, for which there exist finite action solution living in a complexification of the field space: the complex-unitons (2.37). We believe that their fractionalize constituents, the complex-instanton, will manifest themselves in the transseries expansion of some physical observables of the η-deformed SU(2) PCM. We are currently computing such a candidate observable.

Stokes phenomena and Seiberg-Witten theory
On a slightly different note from the previous sections, we will analyse now the connection between our quantum mechanics (3.23) and N = 2 gauge theories in 4-d.
As we have seen, the presence of poles in the Borel plane means that lateral Borel resummations undergo jumping phenomena as one crosses Stokes lines. A helpful way to visualise these jumps is as topology changes of Stokes graphs (an elegant presentation with many more details than provided here can be found in [87]). As we shall detail shortly these Stokes graphs can be understood as spectral networks [88] in N = 2 gauge theories. To explain this we first return to the general Schrödinger type equation (4.1) with a coupling = g 2 > 0, and a WKB ansatz, Viewing q as a local coordinate on a Riemann surface one has that λ = χ(q)dq is a globally defined holomorphic differential one-form whose leading term is just the classical differential λ 0 = χ −1 dq = V (q) − Edq. The exact WKB quantisation condition reads, where the regular perturbative series has an integration cycle running between two real turning points (zeros of Q 0 ). More generally we can analytically continue → e −iθ such that the Schrödinger equation becomes, where now z is recognised as a complex variable. In effect, rotating θ corresponds picking out a ray in which we would like to perform a Borel resummation.

JHEP07(2016)088
In what follows a crucial role is played by Stokes curves which are trajectories of the quadratic differential, φ = Q 0 (z)dz ⊗ dz , (6.5) and its square root λ 0 = χ −1 (z)dz such that φ = λ 0 ⊗ λ 0 . A directional Stokes curves emanating from a point a is a trajectory z(t) of constant phase i.e.
A generic trajectory is asymptotic in both directions to a singular point (perhaps the same one and often at infinity). A separating trajectory connects a turning point to a singular point and a finite trajectory is asymptotic in both directions to a turning point (perhaps the same one) or is closed.
For the case at hand with V (z) = sin 2 (z)(1 + η 2 sin 2 (z)), we plot in figures 17 and 18 Stokes graphs for, respectively, η real and η pure imaginary. Both cases exhibit finite Stokes trajectories connecting real turning points indicating a Stokes direction at θ = 0. In the case of η pure imaginary there are additional Stokes directions in which finite trajectories connect real and complex turning points. As we take the energy parameter of the quadratic differential to zero, this Stokes directions occurs precisely at the angle set by the complex uniton configuration θ = arg(S CI ) and θ = arg(S CI ), as can be seen in figure 19. This is in line with the locations of poles in the Borel plane, see figure 14 and surrounding discussion.
There is a beautiful connection between Stokes graphs and N = 2 gauge theories [88] where the differential λ 0 is identified with the Seiberg-Witten differential λ SW . As the angle θ progresses these Stokes graphs can undergo topology changing morphs and for particular values of θ, where the topology jumps, a finite Stokes curve appears. These finite Stokes curves correspond exactly to the appearance of BPS states, and can be directly identified with BPS strings configuration [89] in the M-theory brane construction of the N = 2 gauge theory. What then of the full quantum differential λ? The conjecture in [90,91] is that periods give the exact prepotential of the Ω-deformed theory in the NS limit i.e. with λ = χ(q)dq. The idea is that the = expansion of the quantum mechanics exactly matches the 1 expansion of the prepotential. This connection has been most clearly exemplified by considering the quantum mechanics associated to the Mathieu equation [69,70,[91][92][93] for which the corresponding gauge theory is the SU(2) N f = 0 theory. Since the quantum mechanics associated to the undeformed Principal Chiral Model is that of the Mathieu equation and that of the η-deformed PCM is the Whittaker-Hill equation it is natural to wonder if there is some JHEP07(2016)088 similar connection at play here. Indeed this is the case; the η-deformed quantum mechanics can be related to a mass deformation of the SU(2) gauge theory. 16 For the Whittaker-Hill equation eq. (3.17) we have a quadratic differential, φ WH = (−a + 2q cos(2θ) + 2p cos(4θ))dθ ⊗ dθ . (6.9) The quadratic differential of the first realisation of the SU(2) N f = 2 gauge theory whose differential is [95], where Λ is the scale of the theory, u the parameter along the Coulomb branch and m i are masses of two flavour hypermultiplets. These two differentials coincide after the transformation z = e 2iθ and the identification, In terms of the parameters of the η PCM, eq. (3.18), we see the nice relation, It is striking that the mass parameter is the simple RG invariant combination of parameters in the η-deformed model. One interesting feature of the SU(2) N f = 2 gauge theory is that for a certain value of the bare mass, Seiberg-Witten singularities coincide and the theory enhances to an Argyres-Douglas SCFT [96]. In the conventions above this occurs at M = 2Λ [88] and in our variables this happens for η 2 = − 1 2 , precisely the value for which the complex critical points coalesce into a single real value as discussed in section 3.3.
From the quantum mechanics point of view we expect a drastic change in the transseries expansion once η 2 passes this critical value, while in the gauge theory side we have that the Argyres-Douglas point corresponds to the critical point of a phase transition [97]. In [97] the free energy of the N f = 2 SU(2) theory has been computed exactly for M < 2Λ in a suitable decompactification limit and it corresponds to the imaginary part of the Seiberg-Witten prepotential evaluated at a D = 0, i.e. at a point where the B cycle shrinks to zero (6.8). If we increase M while keeping a D = 0, we can reach the point M = 2Λ where the cycle A also shrinks to zero size. Precisely at the Argyres-Douglas point the Riemann surface defined by the Seiberg-Witten differential develops a cusp and the free energy manifests a non-analytic behaviour signalling the presence of a phase transition. It would be very interesting to understand if we can learn anything new about the two phases of the transition from the quantum mechanical model. 17 16 Whilst this manuscript was in the final stages of preparation two preprints [45] and [46] studied the Whittaker-Hill equation, in the context of the 4-d/2-d correspondence, where it emerges as a limit of the null vector decoupling equation associated to irregular conformal blocks for the N f = 2 SU(2) gauge theory (as previously noted in [94]). The preprint [46] also provides insights into the topology changes exhibited in the Stokes graphs relevant to these gauge theories. 17 We thank Jorge Russo for bringing these important aspects to our attention.  Figure 17. Stokes graphs for η = 1/2 (with E = 0.1) as we rotate the angle θ = {0, π 8 , 3π 8 , π 2 } (top left to bottom right). Turning points are indicated in purple, generic trajectories in grey and separating/finite trajectories in red. At θ = 0 there are finite Stokes trajectories connecting the real turning points (along the horizontal axis) and simultaneous Stokes trajectories connecting the real turning points to the complex turning points arising from the η-deformation.
As a final comment we note that the expansion of the quantum mechanics corresponds to the small Ω-deformation limit 1 → 0 while 2 = 0. A natural question would be to understand the connections between this expansion and the weak gauge coupling constant expansion performed in [23,25] (see also the earlier work [24]) for a similar class of N = 2 theories.

Conclusions
In this paper we analyse the η-deformed SU(2) principal chiral model. We find unitonlike solutions to the second order equation of motions that reduce to the standard uniton solutions in the PCM model when we send the deformation parameter η → 0. Furthermore we find new solitons that are solutions to the complexified equation of motions. These complex-unitons have quantized actions and live in a natural complexification of the field space Sl(2, C) ⊃ SU (2).
We put the theory on a cylinder R × S 1 L with twisted boundary conditions related to adiabaticity and subsequently performed a KK-reduction to obtain an effective quantum mechanics for the low energy degrees of freedom. The quantum Hamiltonian thus obtained Figure 19. Here we plot the critical Stokes graph for η = i/2 at the angle θ = arg(S CI ) as we take E = 0. This reflects the location of singularities in the Borel plane of the perturbation theory illustrated in figure 14.

JHEP07(2016)088
is related to the Whittaker-Hill equation and displays a large variety of non-perturbative finite action soliton solutions. These non-perturbative objects are instantons modified by the η-deformation, forming the fractionalized constituent of the previously discussed uniton solutions. Furthermore we find additional non-perturbative objects living in a complexification of the quantum mechanical model. These complex-instanton solutions are precisely the fractionalized constituent of the new complex-uniton solutions discussed.
Resurgence theory combined with a Morse-Picard-Lefschetz approach to the path integral of quantum field theory suggests that all finite action solutions living in a complexification of the field space should contribute to the semi-classical expansion of physical observables. We strengthen this claim by studying the asymptotic perturbative expansion of the ground state energy of our reduced quantum mechanics. We identify all the Stokes lines of the Borel transform of the perturbative ground state energy and relate them precisely to the instanton-anti-instanton and complex-instantons events.
We expect these fractionalized constituents of the unitons and the complex-unitons to contribute to the transseries representation of the path integral for the η-deformed principal chiral model. Note however that generically, due to the running of the coupling constant, the location of the singularities in the Borel plane might change as we perform this adiabatic continuation from R 2 to R×S 1 L as explained in [9]. In the model under consideration we have an additional RG-invariant parameter that might allow us to understand quantitatively how this Borel flow behaves. We are currently investigating in perturbation theory some bonafide quantum field theory observable to match against the quantum mechanics expectations of the Borel plane structure.
We make use of uniform WKB and resurgence theory to extract non-perturbative information out of the large order behaviour of the perturbative energy coefficients and we predict the first few non-perturbative corrections on top of these solitonic objects. A proper complexified path integral derivation of these perturbative corrections is not present at the moment and it would a beautiful confirmation of the resurgence program.
Finally we relate the Stokes phenomena associated with the jumps in lateral resummation for the perturbative expansion of the ground state energy of our quantum mechanics with the Stokes graph associated with the spectral network of a very particular 4-dimensional N = 2 supersymmetric gauge theory. We relate the Schrödinger equation associated with the Whittaker-Hill equation to the Seiberg-Witten quadratic differential of a 4-d N = 2 gauge theory with gauge group SU (2) and N f = 2. As we dial the η-deformation parameter to the particular value η = i/ √ 2 we have that multiple complexified classical vacua of the quantum mechanics coalesce and the corresponding 4-d N = 2 theory is driven to an Argyres-Douglas SCFT.

A.1 Conventions
For a compact semi-simple Lie group G corresponding to an algebra g, we parametrise a group element g ∈ G by local coordinates X µ , µ = 1, 2, . . . , dim(G). The right and left invariant Maurer-Cartan forms, as well as the orthogonal matrix (or adjoint action) relating them, are defined as, L a ± = L a µ ∂ ± X µ = Tr(T a g −1 ∂ ± g) , R a ± = R a µ ∂ ± X µ = Tr(T a ∂ ± gg −1 ) , R a µ = D ab L b µ , D ab = Tr(T a gT b g −1 ) , The generators T a obey [T a , T b ] = f ab c T c , are normalised as Tr(T a T b ) = δ ab , and with respect to the Killing metric, defined by f ac d f bd c = −c G δ ab , the structure constants with lowered indices f abc are totally antisymmetric. Group theoretic indices are frequently raised by using δ ab . World-sheet light cone coordinates are defined as σ ± = τ ± σ.

A.2 Drinfeld double and the R-matrix
It may be helpful to the reader to summarise some salient facts about the R-matrix that is used to define the η-deformations. Consider a semisimple Lie group G, a Lie algebra g, and a matrix R (an endo-morphism of g), assumed to be anti-symmetric with respect to the Killing form on g, which defines a bracket,

JHEP07(2016)088
The two Lie-brackets (the usual one [A, B] and [A, B] R ) over the same vector space define a bi-algebra. Thinking of the [A, B] R as defining an algebra g R we have a Drinfeld double defined by d = g ⊕ g R = g C viewed as a Lie algebra. Indeed, g C may be equipped with an inner product, A + iB, A + iB = Im(A + iB, A + iB ) , with respect to which g is a maximal isotropic and when R is anti-symmetric w.r.t. (·, ·) so is g R .
We can specify a Cartan basis of g C : H µ are a Hermitian basis for the Cartan sub algebra and E α and E −α are the raising and lowering operators for a root α (see e.g [27] section 2.2 for details). Then a basis for the real algebra g is given by {T µ , X α + , X α − } with α > 0) where, In this basis a canonical choice for the R matrix is given by a simple action, Although we shall not consider this in the present context, it is possible to augment this R with "homogenous" solutions to the classical c 2 = 0 Yang-Baxter Equation, i.e. rotations amongst Cartan directions associated to so-called Drinfeld-Reshetikhin twist.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.