MSc Physics Track: Theoretical Physics. Getting a grip on the Stokes Phenomenon

Transcription

1 MSc Physics Track: Theoretical Physics Master Thesis Getting a grip on the Stokes Phenomenon Stokes Constants and their relations by Bruno Eijsvoogel 6487 May 7 6 EC //6 - //7 Supervisor dr. M.L. Vonk Second Examiner dr. D.M. Hofman Instituut Theoretische Fysica Amsterdam ITFA

2 Abstract The asymptotic behaviour of a function as it approaches a limit point can depend on the direction it is taken in. The presence of a discontinuity in this dependence is called the Stokes phenomenon and its precise form can be characterized by a set of constants known as Stokes constants. This master s thesis describes the calculation of Stokes constants for the Painlevé I equation in a pedagogical manner. Starting from a treatment of some summation techniques for asymptotic series, we introduce the necessary concepts within resurgence theory in order to calculate the Stokes constants for Painlevé I. Previously uncalculated S turns out to equal 4 4 n S. This result supports the conjecture n S = n S n. The origins of this and other conjectured relations between Stokes constants remains unknown. Finding more relations could lead to a better understanding of the Stokes phenomenon. Finally we briefly consider resurgence in the context of topological string theory, to motivate Stokes constant calculation using the Holomorphic Anomaly Equations.

3 Contents Introduction. Structure of the Thesis Research Aims Physical Context Series Expansions 4. Power Series Asymptotic Expansions Divergence and non-perturbative contributions Transseries Expanding about Stokes Phenomenon & Stokes Constants 9. Borel Summation Method Singularities in the Borel Plane Alien Calculus Pointed Alien Derivative The Bridge Equation Painlevé I 4. Asymptotic Expansion Solution Transseries Solution Large Order Calculations Stokes constant calculation strategy Topological Strings 9 5. Complex Structure Linear Complex Structure Complex Structure on the Torus Cycles and Periods of the Torus Calabi-Yau manifolds Cycles and Periods on the Calabi-Yau The Holomorphic Anomaly Equations Summary and Outlook 7 7 Acknowledgements 8 Appendices 9 A Approximation Methods 9 A. Borel-Padé Summation A. Richardson Transform B More Large Order relations 4 B. Large Order B. Large Order B. Large Order

4 Introduction. Structure of the Thesis The Stokes phenomenon can be described as the sudden change in asymptotic behaviour of functions in different regions in the complex plane. This thesis has firstly the aim to give a pedagogical treatment of the Stokes phenomenon. A treatment that is both intuitive but also presents some detailed calculations. We will mainly look at the Stokes phenomenon in the context of Borel summation in which the Stokes phenomenon has a fairly clear explicit formulation in simple examples. We will be tackling the two-parameter transseries solution of the Painlevé I equation which exihibits the Stokes phenomenon in a way that is still unclear. The second objective is to present this difficulty and the research that was done in the name of this master project, aimed at further understanding the Stokes phenomenon in this context. Since the majority of the thesis will be this first pedagogical review, the present introduction will be a sketch of the second part: the research question. This way we can keep in mind what the review is all for. The research done in this thesis can be viewed as a small addition to the body of work which has [] as its foundation. The work done is however more closely related to the progress made in []. We have tested the methods used in [] pertaining to sections regarding Painlevé I and used them to make further calculations. Initially the goal of this project included conducting calculations of Stokes constants in the context of the Holomorphic Anomaly Equations. Time proved to be too short, but the literature review persists near the end of the thesis in section 5.. Research Aims We will be looking at the Stokes phenomenon in the context of the so-called transseries solution of the Painlevé I equation. They will be introduced in sections. and 4 respectively. Such a solution is divergent, and to make sense of it we need to perform a Borel summation, see section.. To do this we will need to evaluate an integral through the so-called Borel-plane. The point is that there are singularities in the Borel-plane resulting in more than one way to do this integral. This ambiguity is just the Stokes phenomenon in a different guise. In the simplest cases one-parameter transseries we can explicitly describe how the solution changes if we integrate around the other side of a singularity. But for an even slightly more complicated case two-parameter transseries the change is not clear at all. What we do know is that these changes are encoded by the so-called Stokes constants. Let s sketch what this looks like. Denote a one-parameter transseries as Φz, σ where z is the expansion parameter and σ is the integration constant that parametrizes the family of solutions and denote the Borel summation operator along a path θ as S θ. In case we have a singularity on the positive real axis and θ + and θ are paths passing above and below the singularity respectively, the two different Borel summations are related by S θ +[Φ]z, σ = S θ [Φ]z, σ + S. where S is a Stokes constant. For a two-parameter transseries Φz, σ, σ so far we only know that it should look like S θ +[Φ]z, σ, σ = S θ [Φ] z, ˆσ σ, σ, ˆσ σ, σ,. where the ˆσ i are power series in σ and σ and whose coefficients contain the Stokes constants. A list of these constants S i, S i have been calculated before in [], and in this thesis we calculate one more. These Stokes constants have various relations amongst themselves, some of which can be derived relatively straightforwardly. The origins of some other relations however are not yet understood. The upshot is that these non-trivial relations should lead us to understanding the form of these functions ˆσ i and through this provide a better understanding of the Stokes phenomenon. In this light it should be clear that it would have been exceedingly interesting to produce a similar analysis in another context than that of Painlevé I. This was the original plan for this project. Namely to see if similar nontrivial relations existed between the Stokes constants of for example the Holomorphic Anomaly Equations. For the tilde. The fact that we are dealing with a two-parameter transseries gives us two classes of Stokes constants that we distinguish with the

5 Figure : The two integration contours around a singularity corresponding to different Borel summations. research done in this area see [8] [9] []. Within the context of Painlevé I however, we did work out some long calculations and produced a new Stokes constant S. The result supported a conjectured relation between Stokes constants: n n S = n S n... Physical Context Throughout most of this thesis the objectives might often seem very far from any physical application. The Stokes phenomenon is a mathematical phenomenon but it occurs in approximations to solutions of many physical problems. Airy s equation that describes the diffraction of light see [] has an asymptotic series solution that exihibts the Stokes phenomenon, as does a geometric model that describes crystal growth [4] and the first of the Painlevé transcendents, the last of which will be our main focus in section 4. Painlevé I shows up among other things in D quantum gravity [6] while Painlevé II can be used to mode the electric field inside a semiconductor [5]. Aside from physical problem that are directly effected by the Stokes phenomenon, analyzing this problem is physically interesting because of its link to resurgence theory: a mathematical framework that offers a bridge between perturbative and non-perturbative information. We will however not allude to any of these physical examples. This thesis will be about mathematical methods. That this is related to physical problems is only in the back of our mind. Series Expansions In this chapter we will start by discussing power series as representations of functions. More specifically, the method of assuming solutions to differential equations are of the form of a power series. We will encounter power series that diverge and this will lead us to a more general assumption for solutions to differential equations called a transseries. After that we will specifically look at solutions of this form that are resurgent. Unfortunately a full introduction to the intricacies of resurgence goes beyond the scope of this thesis. So instead of starting out with proper definitions, we will introduce the concepts we will need as we go along by their most important properties and refer to more advanced texts for interested readers.

6 . Power Series Let s take an ordinary differential equation ODE as an example to see how we can use power series to describe its solutions. df = fx. dx This ODE is very simple and the solution is obvious: fx = Ce x. But what if we hadn t thought of this? We could have made the assumption that the solution was of the form: fx = n a n x n.. Inserting this form into the ODE gives us an expression for the thus far undetermined coefficients: na n x n = a n x n n n na n x n = a n x n n n + a n+ x n n n = a n x n. n. Equating the equal powers of x we get resulting in the following series n + a n+ = a n for n,.4 fx = a + x + x + x 6 + x This gives us exactly the expected solution with a = C. This is the power series method at its best and easiest. The expansion describes the solution completely and the power series converges for all x both real and complex. This is due to a n shrinking extremely fast. We will make the concept of fast growing and shrinking of coefficients more precise in section..... Asymptotic Expansions Now let s look at an example that doesn t go as smoothly. x df = fx.6 dx Following the same procedure of inserting a power series solution, we find the recursion for the coefficients a n+ = a n..7 So our solution is the geometric series times an arbitrary factor: fx = a x n..8 This sum converges to a function of closed form fx = a x, but only for values x <. The sum does fully describe the solution within this region, but outside the sum diverges and we cannot assign values to it in the usual fashion. However the function fx is valid in the whole complex plane except at x =. It also remains a solution to our ODE for these values. Expanding the domain of fx outside the unit disc is an example of analytic continuation. It is an important result in elementary complex analysis that the divergence of the sum f at x = is entirely due n

7 to the singularity at the point x =. So far we have considered series that, when they converge, are expressible in a closed form. This isn t always necessarily the case. For a converging series without such an expression, one can calculate more and more coefficients in order to achieve better precision. Alternatively one can also take values of x closer and closer to to make the approximation better. We obviously prefer to add terms, since this leaves us free to keep x where it is. There are however expansions for which adding terms doesn t improve the approximation because they don t converge! But if for a fixed number of terms the approximation does improve as we take x closer to the origin it is called an asymptotic expansion. Putting this more formally: fx and gx are asymptotic as x x if fx lim =,.9 x x gx which we denote as fx = gx x x.. An asymptotic expansion is then just a series that is asymptotic to the function under consideration with a specific limit. We will often omit this limit because it will be obvious from the context what it is most often or. The direction in the complex plane from which we take this limit is of crucial importance! In fact for a given function fx one can generally find different asymptotic expansions when taking the limit to x from different directions. This is called the Stokes phenomenon. We will come back to what this effect looks like in the first section of chapter. For now let s do a few numerical calculations with the geometric series to see in practice that divergent series can be asymptotic. Dropping a and starting at x = our true value will be f = =,. and truncating our expansion at n = 5 we approximate it by 5 n f 5 = = with an error ɛ 5 =.5.. n= If we add another term n = 6 the error will be smaller ɛ 6 =.565, but we could also reduce the error by taking a smaller x like f = 4 4 = 4,. and truncating our expansion at n = 5 we approximate it by 5 n f 5 = = with an error ɛ 5 = n= But why should we consider this second method of minimizing the error? It alters which values we can look at of our original function, while adding extra terms doesn t. The important part is that adding extra terms doesn t make your approximation better when the expansion is divergent! Divergent expansions however can still be asymptotic. Let s take our truncated geometric series outside its radius of convergence f = =,.5 and truncating our expansion at n = 5 we approximate it by f 5 = 5 n = 64 with an error ɛ 5 = n= This means that every convergent expansion is also asymptotic. But since this property isn t very relevant in convergent cases one usualy only includes divergent expansions when referring to asymptotic expansions.

8 Looks pretty bad, but it gets better when we go slightly closer to the origin and truncating our expansion at n = 5 we approximate it by f 5 = f = =,.7 5 n = 6 with an error ɛ 5 = n= At the risk of overemphasizing, we repeat that adding extra terms here is going to make our approximation worse. We have the number of terms fixed and bring x closer to. This example may seem a little silly since the values the truncated expansion gives seem completly unrelated to the true values. There are however examples of expansions that for small enough x seem to converge, but after a number of terms diverge anyway. Truncating these expansions at the seemingly converging part ie the term with the smallest absolute value. gives good approximations and is called optimal truncation. In the next section we will look more closely at the role of the expansion coefficients, how their increase and decrease affects the behaviour of the truncated asymptotic expansions and in which cases the use of optimal truncation is a justified method. There are however more sophisticated ways to extract values from asymptotic expansions. One of them is called Borel summation which we will discuss in section.... Divergence and non-perturbative contributions Let s take a closer look at divergence. Whenever a series diverges we can try to put a bound on the growth the sequence of its coefficients. In the last example take f = a We can bound the sequence of coefficients by a large enough constant C n= since the coefficients do not grow. In another case we have fx = = a + a + a a n < C n N. n x n a n = n. n= These coefficients do grow, so there is no constant large enough to bound all coefficients but we can bound the growth as a n < CA n n N.. The minus sign in the exponent is just a convention because of the interpretation A will have later. We will be most interested in series whose coefficients are bounded by a n < CA n n!. These are called Gevrey- formal series. Formal because the series don t necessarily equal a number in the usual sense right now for any x. This is because we are allowing the coefficients grow so fast that the radius of convergence could be. In [] it is shown that in these cases the value at which you should optimally truncate is n = N opt A x. This seems reasonable since the larger A, the slower the coefficients grow and the smaller x, the slower the terms grow as well, so the larger N opt should be. At the end of section will show that the error in this procedure is proportional to a non-analytic function They are definitely not unrelated as we will see later in.4. error e A/x.4

9 This is very important because an expression like this exponential is not expressible as a power series around. If you try Taylor expand this expression about the origin each coefficient will turn out to be. It is a contribution lurking in the function that an asymptotic series is trying to represent, but can t. Because asymptotic series are ubiquitous in perturbation theory, we call this error a non-perturbative contribution. This intuitively leads us to a more general form for a solution of an ODE than our original power series. We could add this non-perturbative contribution with a power series to get a n x n + e A/x n= m= b m x m,.5 in an attempt to better represent the solution. In the very simplest case we might assume that these coefficients b n grow like the a n so the that the total error our new expression would be proportional to e A/x. Now we could add another series with this exponential as a factor in front and continue to do this. We are now building the simplest case of a transseries.. Transseries Instead of assuming the solution to an ODE is a power series, we could use a more general expression called a transseries. Transseries are very general but we will start of with an extremely simple form just to get an idea. Φx, σ = σ n e An/x x nβ n m a n m x m.6 where σ is called the instanton counting parameter and should be thought of as the integration constant that parametrizes a family of solutions to an ODE. A and β are constants that depend on the ODE. The significance of these constants will be made more clear in section 4. It is important to note that the innermost sum is just a power series. The n = term of the outer sum is really just a power series with nothing added. This sum is called the perturbative expansion or perturbative contribution. Each of the other sums is called the n-instanton contribution. This is also due to that instanton contributions to a quantum mechanical amplitude are of this form. The reason for naming things after instantons is that again in the context of perturbation theory the exponentially suppressed contributions are attributed to instantons. The exponential is the main ingredient here, since it is not expressible as a power series. The importance of this object in the theory of resurgence is that a so-called resurgent transseries has the crucial but not defining property that its coefficients are asymptotically related for large values of m! Such a relation is called a large order relation and we will derive the simplest kind for the Painlevé I equation in section 4 and two more complicated examples in appendix B. This is a very interesting property since it implies we can calculate parts of the non-perturbative contributions from just the perturbative expansion. This is the most important property of resurgence and also the origin of its name. Contributions to one power series resurge in another.. Expanding about So far we have been looking at expansions about the origin. Taking powers of x a with a =. We could have just as well taken a. Perturbative expansions about are studied when the true parameter value eg coupling will be small. If however our relevant parameter is actually very big we might want to expand about. Of course we can t plug in infinity for a, but we can take powers of x instead of x see.7. That way if we take x all the terms of the expansion except the constant term vanish, just like setting x = would in an expansion about. For a mathematically rigorous treatment see [9] a n x n.7 n=

10 Stokes Phenomenon & Stokes Constants Whenever we approximate a function we would like to be able to categorize the subsequent contributions to this approximation, most preferably from most important to continually less important. For example take a power series in a small parameter x: or in a large parameter z: x + x + x z + 6z We might call the first term the leading contribution and the second term the subleading contribution. Or in the case of a transseries like in.6 we might, for small x, call the perturbative series leading and the other series proportional to a power of e A/x. subleading. It does however matter how we take this limit of x to. Just looking at the exponentials below we already notice lim e A/x = x lim e A/x x.4 that we can only call the first exponential a subleading contribution. And if we approach the origin from the positive and negative imaginary axis we also get very different results. In fact those limits don t exist since they oscillate at an ever increasing rate. lim e A/ix x.5 lim e A/ix x In the complex plane we have a whole π of different directions in which to approach. If we have a large parameter z we look at e Az which has different behaviour depending in which direction we go to infinity. The formal way as presented in [4] to compare contributions in asymptotic analysis is as follows: fx is negligible compared to gx as x x if fx lim =,.6 x x gx and is denoted fx << gx x x. To make this slightly more tangible we briefly look at the asymptotic expansions for large values of the Airy function Aix. But for a more thorough treatment the author strongly recommends examples concerning this function in []. The function Aix has a asymptotic expansion in x Aix = x / e x /4 π x/ 96x x /..7 We have here an asymptotic relation between a function and a series, but we haven t been specific about the limit. Large z implies the limit should be x, but does the direction matter? In fact the expansion does agree fairly well both along the positive and negative real axis for an arbitrary truncation. However this is not what it means for to functions to be asymptotic. Taking the limits we have that f N x lim =.8 x + Aix

11 while the limit f N x lim x Aix does not exist, where we have denoted f N x as the truncated sum. The true Airy function oscillates towards while the asymptotic expansion does not. We can find another expansion that is asymptotic to Aix at along the negative real axis. Ai x π = sin πx /4 4 + x/ 85 π 468x cos 4 + x/ 5 48x x where we We then say that the Airy function has different asymptotic behaviour in different sectors of the complex plane. By which we mean that the limiting behaviour towards complex infinity of the function is captured by different power series for different directions. And the most striking part of it is that the change is discontinuous in that direction! In a more in depth analysis we can find specific angles at which the asymptotic behaviour suddenly changes. The rays at these angles are called Stokes line and anti Stokes lines. A Stokes line is a boundary between two directions where a a contribution that is negligible on one side that is no longer negligible on the other, while an anti-stokes line is a direction where the asymptotic behaviour of the function is oscillatory.. Borel Summation Method Whenever the coefficients of a power series a n grow as n! the radius of convergence is. However there is a method to still extract values from this series called Borel summation. The essence of this method is outlined by the following steps.. Transform the divergent sum to a different sum that has a finite radius of convergence.. Analytically continue the converged sum outside this radius.. Transform the analytic continuation back using the inverse transformation. The transformation in step is called the Borel Transform B. For a formal series about infinity without constant term φz = a n z n B[ φ]s n= n= a n n! sn.9. A shorthand for B[ φ] is ˆφ. Then the inverse transformation is the generalized Laplace transform L θ, which we first only apply on a single term of the Borel transform series L θ [ a n n! sn] z = e iθ an = a n n! z n n! sn e sz ds. e iθ s n e s ds. = a n z n,.4 because the inverse transformation might not be well defined for the whole Borel transform. We have to analytically continue the Borel transform outside its radius of convergence in a certain direction in the s-plane, also referred to Notice that we don t mean that positive and negative infinity are different limit points. In the complex plane there is one unique complex infinity. But it is one point that can be approached from different directions.

12 as the Borel plane. We have a certain freedom to choose θ for our inverse transformation since we get the same value for the following integral n! = e iθ s n e s ds for θ < π,.5 but to use L θ on the Borel transform we also need to avoid its singularities! A formal series φz is called finesummable in the direction θ if it satisfies the following requirements: Its Borel transform can be analytically continued to complex infinity in the direction θ The series of the inverse transforms of the Borel transform terms [ L θ a n n! sn] z,.6 converges. n= The final result is then valid for { z Reze iθ > c }, for some constant c. Note that this is a half plane centered around the ray in the θ direction. For a formal approach to the question when this is allowed, for which θ s it is allowed and in which part of the z-plane the result holds see [8]. It is a mathematically rigorous treatment that is still very legible for a mathematically inclined physics graduate. Now that everything is defined we can walk through the steps of the Borel summation method.. We first have a Gevrey- formal series φz = n= a nz n. We take the Borel Transform ˆφs = B[ φ]s = So that we now have a finite radius of convergence. n=. We analytically continue ˆφs outside its radius of convergence out to e iθ along a path that avoids any singularities, see figure. We discuss possible singularities in the next section.. We preform the Laplace transformation along this path. The whole procedure can be contained inside the Borel summation operator S θ. It takes an asymptotic series and maps it to a function. This isn t some sort of magic trick that randomly asigns values to asymptotic series. If for example the series under consideration was the formal solution to a differential equation, then the Borel summation will also be a solution. We will sometimes use similar notation like S C to signify analytically continuing along a specified path C that connects the origin to infinity. And in case there are singularities in the direction θ we have the so-called lateral Borel Summations S θ + and S θ to indicate a path in the θ direction that avoids them on the left or on the right respectively, see figure page 4. We have so far excluded formal sums with a constant term. This causes some difficulty at the Borel transform step, but not in the entire procedure of Borel summation. So we can just extend the definition. a n n! s n S C [a + φ] a + S C [ φ].7. Singularities in the Borel Plane If we now have a look at an asymptotic series φx of the Gevrey- kind that has coefficients that satisfy., we will now take its Borel transform. B[ φ] = n= a n n! sn with a n < CA n n!.8 We can then also extend the definition of the Borel transform by B[a + φ]s = a δs + B[ φ]s

13 Figure : We can analytically continue the Borel transform outside the radius of convergence as long as we aviod singularities. The new coefficient an n! grows less fast and gives us a radius of convergence of at least A. If the radius is not infinite this must mean that we have at least one singularity in the Borel plane a distance A away and that it is the closest singularity to the origin. If for example a n = A n n! φz = B[ φ] = A n n!z n.9 n= s n. A Our Borel transform is just the geometric series, which we saw above has a finite radius of convergence n= B[ φ]s = s A for s < A.. We can analytically continue B[ φ]s outside the unit disc. The function has a singularity at s = A though. For now let s say A R >. Then if we want to transform back by integrating over the positive real axis we will have to choose a direction in which to avoid this point. If we take the difference of these lateral Borel summations we get a clockwise closed contour integral that encloses s = A. S +[ φ] S [ φ] e sz = s ds A e Asz = A ds = πiae Az s. So Borel summation ambiguity tells us that the same asymptotic expansion can be related to different functions and the Stokes phenomenon tells us that the the same functions can be related to different asymptotic expansions. It is also notable that this ambiguity in Borel summations is of the same form of the minimum error we found handwavingly in.4. Another reminder that divergence of asymptotic series is directly linked to non-perturbative contributions. Ultimately we would like to apply this summation method to the asymptotic series in a transseries solution and

14 see the Stokes phenomenon there. In section 4 we will consider the specific case Painlevé I but we will encouter an obstacle. Not only is the Borel transform not expressible in closed form, the expansion coefficients aren t even expressible in closed form. We only have a recursion relation. This is why we wanted the Stokes constants in the first place: To get a grip on the Stokes phenomenon. The analysis we have described so far is considered to be classical asymptotics []. The theory of resurgence goes beyond this. To do the central calculation of this thesis in section 4. we need some of its tools. The definition of a resurgent function is quite involved so we refer the interested readers to [], [], [] for introductions within a physical perspective and to [8] for a more mathematically formal introduction which is still readable for the mathematically inclined physicist. Jean Écalle, the inventor of resurgence theory, is the primary source for the complete theory of resurgence, so the truly invested mathematician is referred to []. Instead this chapter will be devoted to displaying a couple of necessary tools for the calculations of section 4. and to give some indication of what they represent.. Alien Calculus The objective of alien calculus is to quantify and compare singular behaviour of functions. There are a few different types of Alien Operators but we will only be discussing the alien derivative ω. The index ω refers to a point in the Borel plane which will singular or else the operator gives zero. In [] this operator is implicitly defined as S θ + = S θ exp e ωz ω. ω Γ θ where Γ θ is the set of all singularities on the ray in the θ direction. Notice the left hand side is just like the calculation of the non-perturbative contribution in.. Hence apparently the Stokes phenomenon can be quantified by these alien derivatives. So we would very much like to calculate them. Remark Before looking at an example let s display two key objects in. that we will encounter later. First of all the outermost exponential S θ = exp e ωz ω.4 ω Γ θ is called the Stokes Automorphism. As pointed out in [7], this can be intuitively compared to generating finite translations of a function from inifitesimal ones. a d dx fx e ωz ω φz.5 fx + a = expa d fx exp e ωz ω φz.6 dx So where the Stokes Automorphism encodes the change in asymptotic behaviour across a direction θ in the complex plane, the alien derivative is like an infinitesimal change in asymptotic behaviour due to a single singular point in that direction. As such we will also be interested in the analogy to fx fx + a = expa d dx fx: The Discontinuity Disc θ = S θ.7 For complete and explicit definitions see [8] and for a shorter summary see []. ω

15 The use of the discontinuity will be our means of introducing resurgence into the calculations in section 4. We will then need the bridge equation from section.5 to compute the alien derivatives. Example Let s start out by looking at an easy example. Starting by immediately restricting ω to be a so-called simple singularity. A function ˆφs has a simple singularity at s = ω if close to ω it looks like a πis ω + ˆψs logs ω ω + ˆχs ω,.8 πi for some functions ˆψs and ˆχs that are holomorphic near the origin and a constant a. In other words our singularities are only allowed to be poles and logarithmic singularities. No log s ω or crazier things. We will make our example even easier by requiring ω be the only singularity in the argω direction. This last assumption makes the exponential on the r.h.s. of. much more manageable exp e ωz ω = + n= e nωz n! ω n..9 Remember that ω commutes with any function that doesn t have a singularity at ω. So plugging this into. we can now rewrite this as our familiar difference in lateral Borel summations. e nωz S θ + S θ = S θ ω n. n! Now lets take an asymptotic series φz whose Borel transform looks like.8 and let the l.h.s. of. work on it, S θ +[ φ]z S θ [ φ]z = ˆφse sz ds C ω n= a = C ω πis ω + ˆψs logs ω ω πi = ae ωz logs ω + ˆψs ω e sz ds C ω πi + ˆχs ω e sz ds. where C ω is the path that goes from to infinity along θ + and back along θ so that we effectively have a closed path around the singularity ω. To simplify the last integral it is useful to separate it into two integrals from to e iθ. And while we re at it we can change shift the integration variable s s + ω e iθ+ɛ = ae ωz + e ωz ˆψs logs e iθ ɛ πi e sz ds e ωz ˆψs logs πi e sz ds. As ɛ goes to the only difference between these two integrals is their logarithms have are on different branches. Their difference being πi. So we can combine them as e iθ ae ωz + e ωz Now compare with r.h.s. of.9 which seems to imply that ˆψse sz ds = ae ωz + e ωz S θ [ ψ]z. S θ [ ω φ]z = a + Sθ [ ψ]z.4 The mathematically precise way to formulate this is that we can find a small enough open neighbourhood around ω where ˆφ coincides with the above expression.8. Note we don t have to change the integration interval since these contours are deformable anyway as long as we don t cross the only singularity ω.

16 and that all higher powers of ω vanish. This last fact should come as no surprise since we have taken ˆψ to be holomorphic at ω. If we now asymptotically expand the Borel summations again we have a more explicit expression This example closely followes [] where more examples are treated. ω φz = a + ψz.5.4 Pointed Alien Derivative One very important property of the alien derivative is how it commutes with the normal z-derivative. We will derive its commutation relation for this specific example. It holds true in general and one can find a proper derivation in [8]. After that we will introduce the pointed alien derivative that does commute with the z-derivative. We want to calculate the commutator [ ] d dz, ω = d dz d ω ω.6 dz The first term is easy since we can just take the z-derivative of.5. But for the second term we need show how B Let s take [ d φ dz d dz φ ω = d ψ dz ] [ ] relates to B φ..7 φz = n= [ ] a n z n with B φ = n= a n n! sn.8 d φ dz = n + a n z n = na n z n.9 n= n= [ d B φ ] na n = s n a [ n = dz n! n! sn+ = s B φ].4 n= So taking the derivative in the z-plane corresponds to multiplication by s in the s-plane. So since we are considering its Borel Transform near ω will behave as.8 times s d φ dz n= as πs ω + s ˆψs logs ω ω πi So we again operate with the l.h.s. of. to get sˆχs ω.4

17 S θ + [ d φ ] [ d S φ ] dz θ = dz C ω ase sz πis ω ds + = aωe ωz + e ωz C ω = aωe ωz + e ωz e C ω s ˆψs ωe sz logs s + ω ˆψse πi ds iθ s + ω ˆψse sz ds sz logs ω πi [ ] e iθ = aωe ωz ωs θ ψ e ωz + e ωz s ˆψse sz ds [ ] [ e iθ = ω a + S θ ψ e ωz + e ωz d B ψ ] e sz ds dz [ [ = ω S θ a + ψ ] e ωz d + S ψ ] θ e ωz dz [ ] [ ] d = ω S θ ω φ e ωz + S θ dz φ ω e ωz ds.4 Where in the last line we used.5 and.7 to rewrite everything in terms of alien derivatives. Now comparing to the r.h.s. of. we must have that all derivatives higher than first order must vanish. The only term that can contribute is [ d S θ φ ] ω e ωz..4 dz So now equating the above and dropping the Borel summation and exponential we finally get our commutator. d φ ω dz = ω φ ω + d dz φ ω [ ] d dz, ω φ = ω ω φ.44 With this commutation in mind we can define the pointed alien derivative which exactly commutes with d dz since ω e ωz ω.45 [ ] [ ] d dz, ω φ = ω ω φ + e ωz d dz, ω φ =.46 The reason this commuting operator is important is that we can use it to derive a specific case of an important equation: The bridge equation. This allows us to express alien calculus in ordinary calculus. We will take some time to show its derivation in the next subsection so we can obtain some simple expressions for Alien derivatives in the case of Painlevé I in section.5..5 The Bridge Equation Now we bring the transseries back to the discussion. We will start with the one-parameter transseries in.6. But our main goal is the two-parameter transseries because the general solution to Painlevé I will be of this form.

18 We assume the one-parameter transseries uz, σ is the solution to some possibly non-linear and inhomogeneous first order ODE in z, but that has coefficients that have vanishing l derivative. F u z, σ, uz, σ, z =.47 If we now operate on this equation with the l, thanks to its commutation with d dz we will get an linear, homogeneous, first order ODE in l uz, σ. The coefficients might depend on uz, σ and that s fine. For example if we have then uz, σ + du cosz =.48 dz l uz, σ + du dz cosz = l uz, σ uz, σ + d uz, σ = dz.49 where we have used that the alien derivative satisfies the Leibniz rule and that the cosine term vanishes because its Borel transform doesn t have singularities for finite l. Now we re not going to solve this equation for l uz, σ. Instead we notice that this reasoning holds for any other derivatives as well, as long as it commutes with d dz. So if we do this with u σ we would get exactly the same ODE for σ. u uz, σ + d σ dz u σ =.5 Since a linear homogeneous first order ODE has a one dimensional solution space, this must mean that they are proportional l uz, σ = S l σ u.5 σ but of course this proportionality is only independent of z. The transseries we have used thusfar is called one-parameter transseries because of the single parameter σ. It also only has one Instanton action A. We now introduce a specific kind of two-parameter transseries uz, σ, σ = n= m= σ n σ m e n maz Φ n m z.5 There is a lot to unpack here. We discuss everything one at a time. First of all and least surprising, we now have two parameters σ and σ. Second our exponential looks different. It is basically the multiplication of e naz, which we had in our old transseries, with a new e maz which is again a non-perturbative contribution with a new instanton action A. Remark One may be concerned about this new instanton action having the opposite sign. The one-parameter transseries converged for z along the positive real axis. Now we don t have this anymore, but this is not a problem. This is a formal expansion. We are not interested in convergence right now and if we are, we will use methods like Borel summation to produce finite numbers. But for now we only interested in the expansion coefficients themselves. Then finally we have Φ n m z which in some cases eg n = m is just an asymptotic series as in our old transseries. Strictly speaking this is a PDE since u also depends on σ, eventhough there are no derivatives with respect to σ. Making this distinction is a bit of a nuisance because it implies that our equation is either ODE or PDE depending on our ansatz. An ODE if we solve for an asymptotic series and a PDE if we solve for a transseries. We will thus call it an ODE and even leave the derivatives with respect to z as non-partial derivatives.

19 For n m it can contain powers of logarithms which we will discuss in the next section, but that we don t have to worry about just yet. Just like we called the subsequent asymptotic series of the one-parameter transseries n-instanton series, we will call the Φ n m the generalized n m-instanton series. For a two-parameter transseries satifying a second order ODE in z, it can be shown that its pointed derivative will satisfy linear, homogeneous, second order ODE. Again its ordinary derivatives with respect to the two parameters σ, will satisfy that same ODE. We will show this explicitly for Painlevé I in section 4. Since we have three distinct functions satisfying the same second order, linear ODE, we can therefore write the pointed alien derivative as a linear combination. la uz, σ, σ = S l σ, σ uz, σ, σ + σ S l σ, σ uz, σ, σ.5 σ This is then the bridge equation for the two-parameter transseries. We now define the degree of an expression as degσ n σ m e laz = n m + l.54 so that deg uz, σ, σ = la deg uz, σ, σ = l deg uz, σ, σ = σ deg uz, σ, σ = + σ.55 If we want the degree of both sides of.5 to match we need deg S l σ, σ = l deg Sl σ, σ = l.56 So if we assume that S l and S l are expressible as power series S l σ, σ = S l σ, σ = k = k = k = k = S k,k l σ k σk S k,k l σ k σk.57 then with the requirement on their degree, they must be S l σ, σ = S l σ, σ = k=max, +l k=max, l S k l+,k l σ k l+ σ k S k,k+l+ l σ k σ k+l+.58 where the constants no longer depend on two indices so we can redefine S k l+,k l S k l+ l S k,k+l+ k+l+ l S l.59 We use la in the alien derivative this time because the singularities will be at multiples of A, but this doesn t make our equation less general. We also take l because that alien derivative should be anyway. The Borel transform of u should have a finite radius of convergence Gevrey- see..

20 to get what are called Stokes constants. Now we can start collecting the expressions that we will want to fill in the bridge equation.5. la u = n,m σ n σ m e n m+laz la Φ n m z u = nσ n σ m e n maz Φ n m z σ n,m u = mσ n σ m e n maz Φ n m z σ n,m.6 S l σ, σ u σ = n,m,k S l σ, σ u σ = n,m,k nσ n+k l σ m+k e n maz S k l+ l Φ n m z.6 mσ n+k σ m+k+l n maz e Sk+l+ l Φ n m z.6 So far our indices n and m run from to. But we want to compare the bottom two lines to the one of the alien derivative. So we will shift the indices of the former. For.6: ˆn = n + k l ˆm = m + k.6 For.6: Making S l σ, σ u σ = = ˆn, ˆm,k ˆn, ˆm= ˆn = n + k ˆm = m + k + l ˆn k + l σˆn σ ˆm e ˆn ˆm+lAz S k l+ l Φ ˆn k+l ˆm k z minˆn+l, σˆn σ ˆm ˆn ˆm+lAz e ˆm k=max, +l ˆn k + l S k l+ l Φ ˆn k+l ˆm k z We have manipulated index values a bit. First off notice that since from the old indices satisfying n and m we can infer that k ˆn + l and k ˆm. This gives us the upper bound for k. Further our new indices ˆn and ˆm have starting values that can be greater than depending on k and l, so this should mean we can t pull those sums outside the k-sum. But we can set them to start at anyway since the extra terms we get are the only ones with Φ x y with x, y. We thus require those Φ to be. Similarly for the other sum. S l σ, σ u σ = = ˆn, ˆm,k ˆn, ˆm= ˆm k l σˆn σ ˆm e minˆn, σˆn σ ˆm ˆn ˆm+lAz e ˆn ˆm+lAz Sk+l+ l ˆm l k=max, l Φ ˆn k ˆm k l z ˆm k l Other natural conventions for labeling S and S are also possible. We follow the conventions of []. k+l+ S l Φ ˆn k ˆm k l z.66

21 These have an identical form as the expression for the pointed alien derivative, so this finally gives us an easily calculable expression for alien derivatives of the generalized n m-instanton series. la Φ n m = + minn+l,m k=max, +l minn,m l k=max, l n k + l S k l+ l Φ n k+l m k z m k l k+l+ S l Φ n k m k l z.67 4 Painlevé I Now we move on to introduce the equation which will be the main focus of this thesis. Painlevé I: uz d u 6 dz = z 4. In this section we will briefly show how one goes about calculating the asymptotic series solution as well as the two-parameter transseries solution. The focus however will lie on the so-called large order relations. These are approximate expressions for the transseries coefficients that contain the Stokes constants. The large order relations are of interest to us because they are our means to calculate Stokes constants. 4. Asymptotic Expansion Solution This is an equation whose asymptotic series solution exhibits Gevrey- behaviour. We can calculate it with the following powerseries ansatz u z = z α u g z αg with α i R, and α > 4. g where we have expanded about infinity instead of this time. This is because in the context of D quantum gravity a small coupling corresponds to a large z []. First off we might want to determine α. Let s look at the forms u and u take u = z α g u = z α g u = z α g α g + α u g z αg α g + α α g + α u g z αg c g z αg 4. where c g = g m= u g m. First off notice that to equal the z on the r.h.s. of 4. we need the u to provide it for us. Since u only has more negative powers of z. Therefore we should set α = and require the succesive powers of z from u and u to cancel. If we want these contributions to cancel we need their powers in z to be comparable. This determines α by equating the powers in z m u Note that we take α > because we want to expand in /z. α α g = α α g 4.4

22 where we need g = g +. This leads to α = 5. Now that our parameters are fixed we can start calculating the series coefficients. We can use our expressions from 4. to plug in 4. and compare equal powers. Firstly we will have u = 4.5 which we normalize to u =. But subsequently the u term enters and the loose z on the r.h.s. leaves, giving us g m= g u u g + m= c g = 6 α g + α α g + α u u m u g m = 6 u m u g m = 5 g + Rewritting a few things and replacing u = we get our recursion u g for g g + u g 4.7 5g 65g 4 u g g = 5n u g 48 g u m u g m 4.9 m= When we calculate these coefficients we see that they diverge factorially. Such that for large g u g CA g g! 4. In other words the expansion is asymptotic and of Gevrey- type. We know from section.. that this constant A is related to the error of optimal truncation. We have approximated it by taking the ratio of coefficients lim g u g g + u g+ = A 4. This convergence is relatively slow. A useful trick to speed up this procedure is the Richardson transform for which we have provided a short manual in appendix A.. We can find the exact value of A in a similar way to finding the parameters α i and β of the asymptotic series. By looking at the transseries ansatz for Painlevé I, one finds that we have see [] []. A = Transseries Solution The Painlevé I equation also admits a one-parameter and two-parameter transseries solution. For the procedures needed to calculate series coefficients for these series we refer to [] []. Here we just present and discuss the form of the two-parameter solution. We first change the notation a bit according to the conventions of [] x = z 5/4 ûx = uz z z=x 4/5, 4. We can make this choose it to be instead of - because of the ambiguity in the sign of the leading z term. It could have just as well been z.

23 and to keep the notation less messy, we won t use the hat on the û. transseries it will be the rescaled one. Now our ansatz for the solution is So whenever we use the two-parameter where Φ n m x = ux, σ, σ = k nm k= n= m= σ n σ m e n ma/x Φ n m x 4.4 log k x k Φ [k] n mx with Φ[k] n m x = g= u n m[k] xg+β[k] g+β nm [k] nm 4.5 There is quite a lot of notation to be explained here. Let s start with Φ [k] n mx. This is a power series in x with coefficients u. The coefficients obviously have labels n, m and k for the series they belong to and g as the basic summation index. β nm [k] plays the same role as α of determining the starting power in asymptotic series in 4.. We also shift the index by the amount β nm [k] so that index and the power of x are the same besides the factor of. The expression for β nm [k] is given below in 4.9. Then Φ n m x which is a sum of power series with powers of logx that goes up until k nm minn, m mδ nm. 4.6 These log-terms are necessary due to a well known effect in power series solutions called Resonance. Below in 4.7 we will show a recursion for all the transseries coefficients. If you try to derive it for a transseries ansatz without log terms this leads to inconsistencies. In 4.7 for certain values of n and m the leading coefficient u n m g will cancel for all g. Sometimes this is fine and we can isolate u n m g, but this can also lead to equations like =. If we had not taken the log terms this would have also led to equations of the form const=, which is why we need them. This is will not cause any other trouble though, because the log term coefficients can be easily calculated from the others see 4.. Now we can calculate as many coefficients as we want of each n m-instanton series, by inserting our transseries ansatz back into the differential equation. This gives us an implicit recursion relation for the coefficients. This is a long and tedious calculation so we just give the result presented in []. δ g δn δ m δ k = g n ĝ= ˆn= m k u ˆm= ˆk= ˆn ˆm[ˆk] ĝ n ˆn m ˆm[k ˆk] u g ĝ 5 96 n m A u n m[k] g + 5 m nk + Aun m[k+] g m ng Aun m[k] g 5 k + k + un m[k+] g k + g 4un m[k+] g 4 5g 65g 4un m[k] g When we want to compute the coefficients for the n m-instanton series for specific n and m, we need to rewrite this equation to get something of the form u n m[] g = where the r.h.s. should only contain terms of lower generalized instanton coefficients. The recursion is correct for any combination of four integers n, m, k, g. We just have to set the coefficients with negative g and k > k nm to zero. When using this equation to calculate coefficients one will notice that for higher n and m that first non-zero coefficient has a ever increasing index. This starting point is encoded by β nm [k] = n + m [ ] knm + k 4.9 The factor of is again just a convention that comes from another choice of variables x = w that is also used in the literature. The power of in x = w shows up here again. This time as a factor of in front of the logarithm. I

24 where the brackets [α] I stand for the integer part of α. From 4.7 we can derive some important relations between coefficients such as u n m[k] g+β u n m[k] g+β = k! = g+β[k] nm n+m/ u m n[k] g+β 4m n u n k m k[] g+β. 4. The second equation holds similarly for the whole corresponding series Φ [k] n m = k! since the factor in front does not depend on the summation index g. 4. Large Order Calculations 4m n Φ [] n k m k, 4. Now that we have the full two-parameter transseries solution at our disposal or at least arbitrarily many coefficients, we will derive a few relations between these coefficients for large values of g, called large order relations. Our tools will be alien calculus and the bridge equation that we discussed in section.5. The important part is that these relations will include the Stokes constants we are after and are therfore a means to calculate them. In the present section we will show how this is done for u [], but u [] and u [] are presented in appendix B. From these examples it should become clear how to produce similar calculations for any n m-instanton coefficients. In Table below we show the large order relations that we have calculated for this thesis project. u [] u [] u [] u [] u [] u [] u 5 [] Table : Calculated large order Relations The initial trick that will pave our path back to resurgence theory is one that works for any complex function with a branch-cut along the θ direction, but that is analytic everwhere else. We rewrite the function with a slightly altered version of Cauchy s Integral Theorem. See figure for a graph of the integration contour. πifz = = e iθ+ɛ e iθ+ɛ R ɛ = e iθ dw fw w z + dw fw w z dw Disc θfw w z θ ɛ+π θ+ɛ e iθ ɛ dφ ire iφ fre iφ Re iφ z + dw fw w z + θ ɛ+π θ+ɛ e iθ ɛ dw fw w z dφ ire iφ fre iφ Re iφ z 4. In the last line we took the arc part of the contour integral to vanish for R. We can generalize this for a function with several branch-cuts such that πifz = θ e iθ dw Disc θfw w z 4. The branch-cut in our case will be a Stokes line, the discontinuity is the one in.7 and the functions we will study are none other than the Φ n m s of the Painlevé I two-parameter transseries solution. Fortunately we know how to calculate their alien derivatives from.67.

25 Figure : The blue line represents the integration contour in the w-plane with an aribitrary point z within the contour of integration. The solid black line is f s branch-cut. Remark We cannot use the expression for the Discontinuity in.7 as is. Its form depends on the expansion variable. We are now expanding in x about, but in section. we were effectively expanding in z = x about. Let us briefly sketch some of the consequences of using different expansion variables. All we really need to look at is how the complex plane maps onto itself by these variable transformations. z = x 4.4 If we look at this map in polar coordinates it should be clear that this inverts every point in the plane about the unit circle and then mirrors it in the real axis. The mirroring is the important part for us, because it flips orientation. Concretely, going counterclockwise in the z-plane coincides with going clockwise in the x-plane. So if we do go counterclockwise in the x-plane and want to cross the branch-cut, that is the same as crossing the z-plane branch-cut clockwise ie with an overall minus sign. S θ + S θ z-plane S θ S θ + x-plane 4.5 Disc θ z-plane Disc θ x-plane Finally let s also look at a different variable since it is also used sometimes used []. x = ω 4.6 Here the new variable has the same orientation but a complete cycle around the origin of the ω-plane is equal to two such cycles in the x-plane. If we where to use the expression.7 only replacing z = ω we would get a result with an incorrect sign and with an additional factor of, because we have crossed the branch-cut twice. We are now ready to derive the large order relation of the perturbative coefficients u [] of the two-parameter transseries solution. Using 4.5 and 4. Φ x = πi + πi Disc Φ w dw w x Disc π Φ w dw, w x Not the same z as in 4.. It can be checked that the inversion about the unit circle leaves the discontinuity invariant. 4.7

26 we can asymptotically expand this as g= u [] g x g + = πi g= x g w g Disc Φ wdw + πi x g w g Disc π Φ wdw 4.8 where we have simlutaneously expanded w x = g= xg w g about x = aswell. Now we can simply equate equal powers of x to get an expression for the coefficient u [] g = πi + πi g= w g Disc Φ wdw w g Disc π Φ wdw 4.9 Now we want to calculate the discontinuities using.7 but with the opposite sign as explained in the remark above. We will take the discontinuities up to first order in exponentials to keep things as simple as possible for now and illustrate the main idea. with Φ x = Disc Φ x = S e A/x Φ x + O[e A/x ] Disc π Φ x = S ea/x Φ x + O[e A/x ] h= u [] h+ xh+ Φ x = where we have used 4. to switch the instanton indices Plugging all of this back into 4.9 we get 4. h u [] h+ xh+, 4. h= u [] h+ = h u [] h+. 4. u [] g = S πi + S πi = S πi h= u [] h+ h u [] h= h= g S i πi h+ w g+h e A/w dw u [] Γg h h+ A g h h= = S i g S πi u [] Γg h h+ A g h h= w g+h e A/w dw u [] Γg h h+ A g h 4. Last thing to note about this is that from 4.7 we can infer that u [] g must be for odd g. This implies that we must have S = is. 4.4 For every Stokes constant we will calculate there will be another Stokes constant that is constrained by the others. The equation above is the simplest case but we will see this happen as well in appendix B. In fact all Stokes constants with negative lower index are constrained in the same way by using the symmetry of 4. on large order That the higher order exponentials are indeed negligible in the corresponing directions to.

27 Figure 4: A plot of the number of decimals that S for different g agrees with S for g = 7. The last few dots grow faster than linear because S for g close to 7, the approximation that S for g = 7 is the true value does not hold anymore. relations. We will elaborate on this topic at the end of the chapter. We can now replace every g with g and get our final expression for the large order relation. u [] 4g = S πi h= u [] Γg h h+ A g h 4.5 We can calculate all the coefficients for arbitrarily high g and h using 4.7 and we can approximate the asymptotic series using the Borel-Padé procedure see appendix A.. So now we can calculate although it is the only Stokes constant for Painlevé I that has actually been calculated exactly []. S = i /4 π 4.6 So far this section was just to illustrate the large order calculation in its simplest form for Painlevé I. In appendix B we extend this procedure to find further Stokes constants. In [] this method was used to calculate a list of Stokes constants not only for Painlevé I, and for this project we have calculated an additional one, namely S = i 4.7 We have calculated this number to what seems to be an accuracy of 9 decimal places at g = 7. The way we estimate this error is firstly by comparing much lower values of g and seeing that the number of decimal places of that agree with S for g = 7 rises linearly with g. Extrapolating this line gives us that estimate see figure 4. The linear behaviour only holds aslong as g is much lower than 7. That way the value at g = 7 is effectively the true value. One of the tantalizing results of [] was that many Stokes constants seemed to be related and a couple of conjectured relations were suggested such as which the result of 4.7 affirms. In [] the related Stokes constant n n S = n S n 4.8 S 4 = i, 4.9 was calculated to 8 decimals! Our calculation converges to the conjectured value S 4 4 S. 4.4

28 Therefore the relation is confirmed to 9 decimals for n = 4. It is important to clarify there are relations between Stokes constants which immediately follow from the large order derivation, such as 4.4 that are entirely different from relations like 4.8. The conjectured relations are done only on the basis of the numerical results. Their origin is not explicitly understood, but it is exactly those non-trivial relations that we hope will lead to a better understanding of the Stokes phenomenon. The form of the Stokes phenomenon we initially set out to understand looks like this S θ +[u]x, σ, σ = S θ [u] x, ˆσ σ, σ, ˆσ σ, σ, 4.4 where u is the full two-parameter transseries solution to Painlevé I. Recall the goal was to find hints on how to describe the ˆσ i. This relation is non trivial for directions θ =, π. For simplicity let s consider the θ = case. We can rewrite this relation in terms of alien calculus using.. This then gives us exp e ω/x ω ux, σ, σ = u x, ˆσ σ, σ, ˆσ σ, σ 4.4 ω Γ which to first order in the outer exponential becomes I + e la/x la ux, σ, σ = u x, ˆσ σ, σ, ˆσ σ, σ. 4.4 l= Now we make the same approximation as before exactly because we made it before too. We only calculated first order Stokes constants so it only makes sense for now to calculate the Stokes phenomenon to first order in alien derivatives. I + e A/x A ux, σ, σ = u x, ˆσ σ, σ, ˆσ σ, σ We plug in the transseries for both u s now σ n σ m e n ma/x Φ n m + e A/x A Φ n m = n,m= Now we use our alien derivative rules from.67 to get A Φ n m = minn+,m k= minn,m + k= We can now use the conjectured relation 4.8 that reads in this case, to simplify this to A Φ n m = minn+,m k= n,m= n k + S k Φ n k+ m k m k S k+ Φ n k m k. ˆσ n ˆσ m e n ma/x Φ n m S k+ = k + k + Sk n k + where we have also shifted the summation index of the second sum. This gives us minn+,m σ n σ m e n ma/x Φ n m + e A/x fn, m, ks k n,m= k= k m k S k k + Φ n k+ m k, 4.48 Φ n k+ m k = n,m= ˆσ n ˆσ m e n ma/x Φ n m, 4.49

29 with fn, m, k n k + k m k. 4.5 k + Equation 4.49 implicitly describes the behaviour of the ˆσ i. We have it to first order in Stokes constants so far. Using more eventual relations between Stokes constants would make it possible to encorporate more contributions to the l.h.s. of We conclude this section about the large order calculations with some practical notes concerning carrying out further Stokes constant calculations. 4.. Stokes constant calculation strategy Much of the very long calculations have been postponed to the appendices. These are surely insightful for anyone attempting to do similar calculations. But a sense of direction is also valuable. This is what the present section will be about. This master s project has mainly focussed on calculations of Stokes constants of the form S p ± and S p ±. There is a recurring pattern in the specifics of the large order calculations that are necessary for these specific constants. Not every large order relation produces new Stokes constants. In general the large order relation of u n n[] will enable one to calculate S n, but will also give a constraint for S n by observing that un n[] g+β = for odd g, just like in 4.4. Large order relations of u n+ n[] in combination with those of u n n+[] will allow one to calculate S n+ and constrain S n+. The constraint is derived from comparing the large order relations of u n+ n[] and u n n+[], by means of the switching the instanton indices as in 4. u n+ n[] = n+[n/] I + u n n+[]. 4.5 The other constrained Stokes constants that we will come across during the appendices are S = is + 4π S S = i π S + S = is π S + S = i S π S S π is S + π i S + S 6π 9 S 4.5 S = i S + π S + 4π i S. The careful reader notices that there seem to be no constraints for S n for n =,. That is because we don t have n these at all. It is purely due to the labeling convention taken in.59. The same is true for the S with positive lower index. These also start at n =.

30 5 Topological Strings A different field that has been the focus of resurgence research is that of topological strings. Since time did not permit it, no analysis of related Stokes constants was conducted. We do however present a pedestrian approach to complex structure, meant as an intuitive handle for students who are interested in complex structure moduli and topological strings. The link being that the free energy of the topological string depends on the complex structure moduli of the manifold in which the topological string lives. The free energy can be calculated by the Holomorphic Anomaly Equations. The resurgent analysis of this problem is indeed very different from the Painlevé I case. Firstly because the coefficients of the asymptotic series and transseries for the free energy are not constants! They each are functions of the complex structure moduli. The instanton actions turn out to also depend on them and it could quite possibly be that the corresponding Stokes constants do too. Another big difference is that the Holomorphic Anomaly Equations are a set of differential equations recursively determining the the series coefficients. The expansion parameter is thus the index of this recursion and not the variable of the differential equations. This means we cannot, as in Painlevé I, derive an analogous bridge equation. The goal of the present chapter is to give an explicit description of what complex structure moduli are in the cases of vector spaces and the torus, and a qualitative description for more general spaces called Calabi-Yau manifolds. Finally we will briefly mention the setting in topological string theory in which complex structure and resurgence play a role. 5. Complex Structure Before we go on to describe what complex structure means on a Calabi-Yau manifold, let s start out simpler and discuss the complex structure of a vector space. In fact we will start off by examining the correspondence between R and the complex plane C. The main difference between these two spaces is that C has a multiplicative operation defined on it. When a complex number is multiplied by another complex number, it is effectively rotated and scaled in the complex plane by the rules of this multiplicative operation. What we will see in this examination of complex structure on the vector space R, is that we can define rules that achieve exactly the same operation effectively reconstructing C but that we can also define other rules than these standard ones. This is what it means to have different complex structure. In section 5.. we will apply this concept to a less trivial space: the torus. 5.. Linear Complex Structure Take an even-dimensional real vector space V. We can equip it with a so-called complex structure by defining a real linear transformation J : V V with J = I 5. Where I is the identity map. To make things as simple as possible, let s consider R with the basis {ˆx, ŷ} = {, }. The standard complex structure is defined as follows: J ˆx = ŷ J ŷ = ˆx J = 5. Notice that this is equivalent to the way the complex plane C works. Where J represents the multiplication by i the imaginary unit and the basis vectors represent addition by and i respectively. The main idea of complex structure is that it gives the real vector space a way to multiply its elements together with ˆx as its identity ˆx ˆx def = ˆx ˆx ŷ = J ˆx = ŷ, ŷ ŷ = J ˆx = ˆx. 5. Let s write this out explicitly for C and R and compare. Actually this is only true for so called B-model topological string theories. a + ibc + id = ac db + iad + bc 5.4

31 The above equation shouldn t raise any eyebrows. Below we see that we can use J to define multiplication in a way that is entirely analogous. αˆx + βŷ γˆx + δŷ = αˆx + βj ˆx γˆx + δj ˆx = α + βj γ + δj ˆx = αγ βδˆx + αδ + βγj ˆx = αγ βδˆx + αδ + βγŷ We can however define a different complex structure than the standard one. Let s look at a more general complex structure on R. f J ˆx = f J f = ˆx J = +f f 5.6 f f Where f = f ˆx + f ŷ R with f so that f is linearly independent of ˆx. So what would the multiplication of vectors look like with this complex structure? We have ŷ = f J ˆx f ˆx so working out our new multiplication we get 5.5 αˆx + βŷ γˆx + δŷ = αˆx + βj ˆx f βˆx γˆx + δj ˆx f δˆx f f f f = αγ f βδ f f αδ + βγ + βδ ˆx f f + αδ + βγ f βδ J ˆx f f = { } f αγ + βδ ˆx + αδ + βγ f f { f f } βδ ŷ 5.7 Notice first of all that when we let f = and f = so that f = ŷ we get our standard complex structure back. Further that we can modify the complex structure continuously by varying the components of this vector f. These parameters are what we will call complex structure moduli and complex structure moduli space is the corresponding parameter space. Which would be the R excluding the line f = if it werent for the fact that not every combination of f, f corresponds to a different complex structure. There is some redundancy in our description above. For example if we replace f, f by f, f we get exactly the same structure. In fact that is the only redundancy in this case. To see this notice that if we vary f and f but the expressions in the curly brackets of the last line in the above equations remain invariant, then the complex structure is the same. The first expression is constant along a hyperbola in the plane en the second is a straight line. f + = c c > 5.8 These equations solve for exactly two solutions f f f f = c c R 5.9 c f = ± c c f = ± c c 5. that differ by an overall minus sign just like we had already guessed. So our moduli space is only half of a R plane. We can choose which half but if we want it to be connected we The expressions seem to have a singularities and values for which the f i become complex but remember that we are choosing the f i first and then checking if some other f i give the same structure. We are not letting the c i just take any values. So we will never run into those problems.

32 Figure 5: Torus better choose either the lower or the upper half-plane. This is because we need f. We will find a very different, more interesting space for the complex structure of the torus in the next section. 5.. Complex Structure on the Torus Let s start by construction the torus from the regular complex plane C. The familiar representation it has as a D surface in three dimensions is shown in figure 6. We can make it from a square. Just glue together two opposing sides to get a cylinder and then glue the two circular edges together. However the surface we start out with doesn t have to be a square. It could just as well be a rectangle or even a parallelogram. You might think this will make a difference because of how you would need to stretch the shape to make it into a torus. It won t make a difference in a topological sense, because topologically a torus is a dimensional surface that allows for exactly two independent closed paths, that can t be continuously deformed into each other or to a point. So streching this surface a bit will not change that. But as we will see it will make a difference if the torus has complex structure. We now want to give this surface complex structure. We ll take the parallelogram to stay general and effectively cut out this shape from C. Imagine putting a parallelogram on C with one corner in the origin and one side on the real axis like in figure 5. Then fill up all of C with nicely stacked copies of this parallelogram. Now define an equivalence relation such that a point will be equivalent to each corresponding point on each parallelogram. That way every point in C now has a unique representative in our original parallelogram. Putting this equivalence on C effectively glues the sides together. So now we have a torus with complex structure. We could have taken C with a general structure as in the previous section to explore the moduli space. But instead we can equivalently tweak the complex structure of the torus by changing the shape of our parallelogram. As a quick example imagine the top-left corner of our parallelogram τ = + i. The standard complex structure J will map x = to J x = i, except this isn t inside the original parallelogram so we need to look for the equivalent point that is inside. It s easy to see that it turns out to be x = + i. What this x turns out to be obviously depends on the shape of our parallelogram. It determines where J maps points to. Therefore it is effectively the By independent we mean that looping around twice doens t count as a new loop. Actually if a point is on the boundary of two or more parallelograms it would have more than one representative. To solve this problem the shape to consider would only include two sides and one corner. In figure 5 we show the included sides as being darker. We don t lose any generality by assuming τ is in the upper-half plane.

33 Figure 6: Torus embedded in D. complex structure. It is useful to write x in terms of the corners τ and. x = ατ + β α, β R Re x = αreτ + β = 5. Im x = αimτ = x = Imτ τ + Reτ Imτ 5. But not every parallelogram makes for a different complex structure. You can imagine tweaking the shape of the main parallelogram such that i lies in another parallelogram that would again be equivalent to the same x. More explicitly there are specific changes to τ you can make that will keep x invariant. This would thus give us the same complex structure. So parallelograms that keep the equivalences between points the same produce the same complex structure. The two transformations below are called the generators of the modular group P SL Z and they will leave the complex structure the same: T : H H : τ τ + S : H H : τ τ, 5. where H is the upper half of the complex plane. It is quite easy to see that T will leave our complex structure the same by plugging it into 5.. But to see that S does so as well we will take a few steps that can be found in any text on modular transformations. So if we want a truly different complex structure we need a new τ that isn t related to the old one by a modular transformation. So the final question is what subset M H has the following properties? No two points related by a modular transformation. Every point in H is related to exactly one point in M by a modular transformation. A domain that satisfies this is called a fundamental domain and one is shown in figure 7. You can imagine that the shaded region has these properties if you think about how T and S map points. The strip has a width of. So only boundary points can be mapped to each other by T. S takes a point and mirrors it in the unit circle and then in the imaginary axis. So any point inside that strip is going to be mapped

34 Figure 7: A fundamental domain of the P SL Z. The three thick dots represent two of the three singular points. The two outer ones are mapped onto each other. outside by S. Again here the boundary points can be mapped to each other. So the final thing to do is to glue those boundaries together. The dotted line onto the dashed one see figure 7. These edges meet at complex infinity which is a singular point as are the others black dots in the figure. Singular in this case means that if you were to draw a circle around it, it would take you less that π to complete it. The glueing together and the addition of complex infinity make the topology of this space that of a sphere, but with singular points: the point that was on top of the circle, the two points on the sides got that got glued together to one and the point at complex infinity. This then finally is the complex structure moduli space of the torus! Each possible complex structure encoded in it and each only once! But we aren t done with the torus yet. The above approach to calculate the complex structure moduli is not feasible in any topological string theory example. In fact describing the moduli space of any such example globally is very hard []. We can however find out relatively easily how many dimensions such a space has and how to describe it locally. In our case this local description would amount to finding τ. In the next section we display how to compute complex structure moduli for the torus by integrating over its cylces. In contrast to the previous calculation for the torus, the following method can be generalized to more complicated higher dimensional examples. 5.. Cycles and Periods of the Torus As mentioned before, the torus has two independent ways in which one can draw a closed path that cannot be continuously deformed to a point. We call these two paths cycles and denote them C,C. They are drawn on both images of the torus in figure 8 and 9. Now the complex structure moduli are locally described by integrals called periods of the one-form dz over these cycles. So that we would have dz = τ = τ C 5.4 dz = =, C where we have deformed the contours to the left and bottom side of the parallelogram respectively. Not surprisingly the results are the top-left and bottom-right corner of our parallelogram. The point of this alternative view is to prepare for a dive in the deep end, when we generalize this story to higher dimensional complex spaces of the Calabi-Yau type. But first we must briefly mention a few aspects to get even the slightest feeling for what a Calabi-Yau manifold is.

35 Figure 8: Integration over the torus seen in the D embedding. Figure 9: Integration over the torus seen in the flat parallelogram.