Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions

We analyze the asymptotically free massless scalar ³ quantum field theory in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer-Connes Hopf-algebraic Dyson-Schwinger e...

Повний опис

Збережено в:

Бібліографічні деталі
Опубліковано в: :	Symmetry, Integrability and Geometry: Methods and Applications
Дата:	2021
Автори:	Borinsky, Michael, Dunne, Gerald V., Meynig, Max
Формат:	Стаття
Мова:	Англійська
Опубліковано:	Інститут математики НАН України 2021
Онлайн доступ:	https://nasplib.isofts.kiev.ua/handle/123456789/211440
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions. Michael Borinsky, Gerald V. Dunne and Max Meynig. SIGMA 17 (2021), 087, 26 pages

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

_version_	1860196575087689728
author	Borinsky, Michael Dunne, Gerald V. Meynig, Max
author_facet	Borinsky, Michael Dunne, Gerald V. Meynig, Max
citation_txt	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions. Michael Borinsky, Gerald V. Dunne and Max Meynig. SIGMA 17 (2021), 087, 26 pages
collection	DSpace DC
container_title	Symmetry, Integrability and Geometry: Methods and Applications
description	We analyze the asymptotically free massless scalar ³ quantum field theory in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer-Connes Hopf-algebraic Dyson-Schwinger equations for the anomalous dimension. This scalar conformal field theory is asymptotically free and has a real Lipatov instanton. In the Hopf-algebraic approach, we find a trans-series having an intricate Borel singularity structure, with three distinct but resonant non-perturbative terms, each repeated in an infinite series. These expansions are in terms of the renormalized coupling. The resonant structure leads to powers of logarithmic terms at higher levels of the trans-series, analogous to logarithmic terms arising from interactions between instantons and anti-instantons, but arising from a purely perturbative formalism rather than from a semi-classical analysis.
first_indexed	2026-03-20T15:48:17Z
format	Article
fulltext	Symmetry, Integrability and Geometry: Methods and Applications SIGMA 17 (2021), 087, 26 pages Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson–Schwinger Equations: ϕ3 QFT in 6 Dimensions Michael BORINSKY a, Gerald V. DUNNE b and Max MEYNIG b a) Nikhef Theory Group, Amsterdam 1098 XG, The Netherlands E-mail: michael.borinsky@nikhef.nl b) Department of Physics, University of Connecticut, Storrs CT 06269-3046, USA E-mail: gerald.dunne@uconn.edu, max.meynig@uconn.edu Received April 07, 2021, in final form September 16, 2021; Published online September 23, 2021 https://doi.org/10.3842/SIGMA.2021.087 Abstract. We analyze the asymptotically free massless scalar ϕ3 quantum field theo- ry in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer–Connes Hopf-algebraic Dyson–Schwinger equations for the anomalous dimension. This scalar conformal field theory is asymptotically free and has a real Lipatov instanton. In the Hopf-algebraic approach we find a trans-series having an intricate Borel singular- ity structure, with three distinct but resonant non-perturbative terms, each repeated in an infinite series. These expansions are in terms of the renormalized coupling. The resonant structure leads to powers of logarithmic terms at higher levels of the trans-series, analo- gous to logarithmic terms arising from interactions between instantons and anti-instantons, but arising from a purely perturbative formalism rather than from a semi-classical ana- lysis. Key words: renormalons; resurgence; non-perturbative corrections; quantum field theory; renormalization; Hopf algebra; trans-series 2020 Mathematics Subject Classification: 81T15; 81Q15; 34E10 1 Introduction The seminal work of Kreimer and Connes showed that there is an underlying Hopf-algebraic structure to the renormalization of quantum field theory (QFT) [28, 29, 61]. This new per- spective has led to deep insights into QFT, and also to novel computational methods that have enabled significant progress in higher order perturbative computations [14, 16, 17, 20, 22, 23, 27, 60, 62, 63, 64, 81, 83, 84, 85, 88, 89, 90, 91]. The Hopf-algebraic formulation of QFT is inherently perturbative in nature, so an important open question is to understand how the non-perturbative features of QFT arise naturally within the perturbative Hopf algebra struc- ture. In a recent paper [18] we showed how this works for 4 dimensional massless Yukawa theory, using Écalle’s theory of resurgent trans-series and alien calculus [3, 33, 44, 48, 80, 82]. Here we extend this analysis to a conformal field theory: massless scalar ϕ3 theory in six dimensional space-time. This QFT has been studied extensively from numerous directions, and has many interesting features, both perturbative and non-perturbative. The theory is This paper is a contribution to the Special Issue on Algebraic Structures in Perturbative Quan- tum Field Theory in honor of Dirk Kreimer for his 60th birthday. The full collection is available at https://www.emis.de/journals/SIGMA/Kreimer.html mailto:michael.borinsky@nikhef.nl mailto:gerald.dunne@uconn.edu mailto:max.meynig@uconn.edu https://doi.org/10.3842/SIGMA.2021.087 https://www.emis.de/journals/SIGMA/Kreimer.html 2 M. Borinsky, G.V. Dunne and M. Meynig asymptotically free for real coupling g [26, 30, 67, 68], and has a Yang–Lee edge singularity when g is imaginary [50]. The perturbative beta function and anomalous dimensions have been computed to 4 loop order [53] (and very recently to 5 loop order [19, 20]). The per- turbative Hopf algebra structure of Dyson–Schwinger equations of this model was formulated in the pioneering papers [22, 23]. On the non-perturbative side, this QFT has a real Lipa- tov instanton when g is real, for which the conventional one-instanton semi-classical analy- sis [21, 66, 78, 93] of the fluctuation determinant has been studied [76, 77]. Further extensions to multi-dimensional cubic interactions have many interesting applications and implications for conformal quantum field theories in general [11, 19, 38, 39, 49, 52, 53, 54, 55, 56, 57]. For other analyses of resurgence properties of renormalization group and Dyson–Schwinger equa- tions see [5, 6, 7, 8, 9, 13]. Our technical analysis is based on the fundamental result [22, 23, 62, 63] that the Dyson– Schwinger equations have a recursive Hopf-algebraic structure which, when combined with the renormalization group equations describing the anomalous scaling under re-scaling of parameters and in the absence of vertex renormalization, reduces the problem to a non-linear ordinary dif- ferential equation (ODE), where the variable is the renormalized coupling. This Hopf-algebraic approximation goes well beyond the familiar rainbow [42] and chain [22, 23] approximations to the Dyson–Schwinger equations. These results cast the Hopf algebra renormalization prob- lem in a form in which very high orders of perturbation theory become accessible, and as we show here it also enables direct access to the associated non-perturbative structure. We employ the trans-series approach to the resurgence properties of non-linear differential equations, along the lines of [31, 32, 33]. Our main new result is that the perturbative Hopf algebra formula- tion encodes a non-perturbative trans-series that involves powers of all three trans-monomial elements: x, e−1/x, and log(x), all expressed in terms of the renormalized coupling. More- over, this trans-series has the form of an all-orders multi-instanton expansion, and the loga- rithms appear with the characteristic structure of logarithmic terms arising from the interaction of instantons and anti-instantons.1 Logarithmic terms are familiar in semi-classical computa- tions [1, 40, 41, 45, 46, 47, 65, 71, 79, 92, 94, 95], and have been studied in differential equations where resonant Borel singularities ±A interact [4, 37, 51], but here we find a quite different res- onant Borel structure, with three resonant singularities of the same sign yet in integer multiples. All this non-perturbative information is encoded in the original perturbative Hopf-algebraic for- mulation, which at first sight makes no explicit mention of instantons, let alone interactions between instantons and anti-instantons. 2 Perturbative Hopf-algebraic analysis of massless ϕ3 theory in 6 dimensions In this paper we analyze the massless scalar ϕ3 theory in 6 dimensional spacetime. This is the critical dimension in which the theory is asymptotically free [68] and in which it has a Lipatov instanton [66, 76, 77]. We analyze the non-perturbative features arising in the Hopf-algebraic approach of [22, 23, 62, 63]. The Lagrangian density is L = 1 2 (∂µϕ) 2 − g 3! ϕ3. 1This logarithmic structure does not occur for the 4 dimensional Yukawa model studied in [18]. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 3 = 1 2 + 1 2 + 1 2 + 1 2 + · · · (a) Rainbow approximation = 1 2 + 1 2 + 1 4 + 1 8 + · · · (b) Chain approximation = 1 2 + 1 2 + 1 2 + 1 4 + 1 2 + + 1 8 + 1 4 + 1 4 + 1 4 + · · · (c) Hopf approximation Figure 1. A comparison of the diagrams included in various approximations for the self-energy of the massless ϕ3 theory: (a) the rainbow approximation; (b) the chain approximation; and (c) the Hopf approximation studied in this paper. The relevant symmetry factors are also indicated. Note that the three approximations agree up to two-loop order, but differ at higher orders. At the third loop order, the Hopf approximation is the sum of the other two. Beyond third order, the Hopf approximation includes new classes of diagrams which are not present in either the rainbow or the chain approximation. As in [22, 23] we consider the renormalized scalar self-energy Π ( q2 ) := and take all propagator self-insertions into account. This Hopf-algebraic approach is depicted by the Dyson–Schwinger equation = 1 2 + + + · · · − subtractions (2.1) with the appropriate BPHZ subtractions indicated. Another way to describe the relevant set of graphs is to start with the one-loop graph and add all possible iterated and multiple insertions of this graph into one of the propagators. Figure 1 shows the resulting low order dia- grams and compares this Hopf expansion with two other common approximations to the Dyson– Schwinger equations: the rainbow approximation [42] and the chain approximation [22, 23]. The Hopf expansion includes a much larger class of diagrams than either the rainbow or the chain approximation, and leads to a much richer non-perturbative structure. The differences between these approximations is discussed further below, in Sections 3 and 4.2 2The effect of including also the vertex corrections will be addressed in future work. 4 M. Borinsky, G.V. Dunne and M. Meynig The pictorial Dyson–Schwinger equation (2.1) corresponds to the integral equation Π ( q2 ) = a π3 ∫ d6k 1 (q + k)2 ( 1 k2 + 1 k2 Π ( k2 ) 1 k2 + 1 k2 Π ( k2 ) 1 k2 Π ( k2 ) 1 k2 + · · · ) − subtractions, where the natural coupling expansion parameter is a := g2 (4π)3 . It is convenient to extract a factor of q2 by defining the function Π̃ ( q2 ) , such that Π ( q2 ) = q2Π̃ ( q2 ) , and the integral equation reduces to q2Π̃ ( q2 ) = a π3 ∫ d6k 1 k2(q + k)2 ( 1− Π̃ ( k2 )) − subtractions. The BPHZ subtractions are chosen such that the momentum subtraction renormalization con- dition Π̃ ( µ2 ) = 0 is satisfied. The anomalous dimension is defined in the momentum subtraction scheme as γ(a) := d d ln q2 ln ( 1− Π̃ ( q2 ))∣∣∣∣ q2=µ2 . Broadhurst and Kreimer [22, 23] showed that the Hopf-algebraic anomalous dimension satisfies the following nonlinear ordinary differential equation: 0 = 8a3γ(a) ( γ(a)2γ′′′(a) + γ′(a)3 + 4γ(a)γ′(a)γ′′(a) ) + 4a2γ(a) ( 2(γ(a)− 3)γ(a)γ′′(a) + (γ(a)− 6)γ′(a)2 ) + 2aγ(a) ( 2γ(a)2 + 6γ(a) + 11 ) γ′(a)− γ(a)(γ(a) + 1)(γ(a) + 2)(γ(a) + 3)− a. (2.2) This equation is quartic in the anomalous dimension γ(a), and third order in derivatives with respect to the renormalized coupling a. Contrast this with the massless Yukawa theory in 4 dimensions where the corresponding nonlinear equation for the anomalous dimension is quadratic in γ(a), and first order in derivatives with respect to a [18, 22, 23]. We therefore expect the ϕ3 theory to have a richer perturbative and non-perturbative structure, as we demonstrate explicitly in this paper. In this massless theory, the full self-energy can be expanded formally in powers of L ≡ ln q2 µ2 , or in powers of the renormalized coupling a [62, 91]: Π̃ ( q2 ) = − ∞∑ j=1 γj(a)L j = − ∞∑ j=1 cj(L)a j , L ≡ ln q2 µ2 . The first term, γ1(a), is just the anomalous dimension γ(a), and all subsequent higher coefficients are expressed recursively in term of γ1(a) [62, 91]: γk(a) = 1 k γ1(a)(1− 2a∂a)γk−1(a), k ≥ 2. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 5 Therefore the basic trans-series structure of γ(a) is inherited by all the subsequent γk(a), and hence also by the self-energy Π̃ ( q2 ) . The unique perturbative solution in R[[a]] to the nonlinear ODE (2.2) can be generated straightforwardly with a formal perturbative ansatz [22, 23] γ(a) := ∞∑ n=1 (−1)n An 62n−1 an = −a 6 + 11 a2 63 − 376 a3 65 + 20241 a4 67 − 1427156 a5 69 + 121639250 a6 611 − 12007003824 a7 613 + · · · . (2.3) The coefficients are all rational, and the normalization of the coefficients was chosen in [22, 23] to make the An integer-valued. The first few integers An are{ 1, 11, 376, 20241, 1427156, 121639250, 12007003824, 1337583507153, 165328009728652, 22404009743110566, 3299256277254713760, . . . } . This sequence is listed in the OEIS [86] as entry A051862. There is currently no known com- binatorial interpretation of the integer An coefficients, in contrast to the 4d massless Yukawa model analyzed in [18, 22, 23], where the corresponding perturbative expansion is the generating function for connected chord diagrams [36, 69, 75, 91]. 3 Asymptotics of the perturbative solution of the Hopf-algebraic Dyson–Schwinger equation Given the ODE (2.2), it is straightforward to generate recursively very high orders of the per- turbative expansion (2.3).3 Via a simple ratio test combined with high-order Richardson extra- polation [10] we can experimentally deduce that the leading large order growth is given by the classical gamma function, multiplied by a power of 12: An = S1 6 12nΓ ( n+ 23 12 )( 1 +O ( n−1 ) ) , for n → ∞. (3.1) Therefore the perturbative expansion in (2.3) is a factorially divergent series. The overall coef- ficient S1, the Stokes constant, can be determined to very high precision: S1 ≈ 0.08759555290917912448379544742126299062738801740682153692058109 . . . . (3.2) The normalization factor 1 6 was chosen for later convenience. S1 does not appear to be a simple recognizable number. The leading large n factorial dependence in (3.1) is close to, but importantly different from, the apparent An ∼ 12nΓ(n + 2) large order growth estimated in [22, 23] based on the first 30 terms. The rational offset 23/12 of the argument of the Gamma function in (3.1) follows analytically from a trans-series analysis (see Section 5) of the Dyson–Schwinger equation (2.2).4 Figure 2 shows a numerical illustration of this offset parameter. It is interesting to compare the Hopf-algebraic perturbative expansion (2.3) with two simpler perturbative approximations to the Dyson–Schwinger equations: the rainbow and the chain approximations, as depicted in Figure 1. The rainbow approximation [22, 23, 42] yields a simple 3We thank David Broadhurst for also providing a list of the first 2000 integer coefficients An. 4See also [5]. https://oeis.org/A051862 6 M. Borinsky, G.V. Dunne and M. Meynig 0 100 200 300 400 500 n 0.005 0.010 0.015 0.020 0.025 growth ratio Figure 2. The ratio of the perturbative An coefficients in (2.3) to the leading order growth 12nΓ ( n+ 23 12 ) (blue), compared with other nearby offsets: 12nΓ ( n+ 24 12 ) (gold), and 12nΓ ( n+ 22 12 ) (green). The origin of the exact factorial offset 23 12 is explained analytically in Section 5. closed-form algebraic expression for the anomalous dimension, which has a convergent expansion (for ease of comparison, we use the same normalization convention as in (2.3)): γrainbow(a) := ∞∑ n=1 (−1)n Arainbow n 62n−1 an = −a 6 + 11 a2 63 − 206 a3 65 + 4711 a4 67 − 119762 a5 69 + 3251262 a6 611 + · · · = 1 2 ( 3− √ 5 + 4 √ 1 + a ) . On the other hand, the chain approximation [22, 23] yields a divergent perturbative expansion of the anomalous dimension (once again with the same normalization convention as in (2.3)): γchain(a) := ∞∑ n=1 (−1)n Achain n 62n−1 an (3.3) = −a 6 + 11 a2 63 − 170 a3 65 + 3450 a4 67 − 87864 a5 69 + 2715720 a6 611 + · · · (3.4) = 3 ∞∑ n=1 (−1)nΓ(n) ( 1 6n − 2 12n + 1 18n ) an. (3.5) Note that the first two terms of the rainbow, chain and Hopf approximations all agree, since they involve the same diagrams, but that the a3 coefficient differs. This is because at this order the relevant diagrams differ, as illustrated in Figure 1. Furthermore, note that the Hopf-algebraic a3 coefficient in (2.3) is −376/65 = −206/65 − 170/65, equal to the sum of the rainbow and chain contributions at this order, as can be understood diagrammatically from Figure 1. This also illustrates the fact that the Hopf-algebraic analysis incorporates a larger class of diagrams than either the rainbow or chain approximations at this order and beyond. It is also interesting to note that for the first 6 terms it appears that the rainbow approximation coefficients are growing more rapidly than those of the chain approximation, but eventually the factorially divergent growth of the chain approximation coefficients overtakes the growth of the coefficients of the convergent rainbow approximation. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 7 We use 500 Hopf expansion coefficients in (2.3) to experimentally extract subleading power- law corrections to the leading factorial growth in (3.1). The methodology is explained in Appen- dix A: An ∼ S1 6 12nΓ ( n+ 23 12 )( 1− 97 48( n+ 11 12 )− 53917 13824( n− 1 12 ) ( n+ 11 12 )− 3026443 221184( n− 13 12 ) ( n− 1 12 ) ( n+ 11 12 ) − 32035763261 382205952( n− 25 12 ) ( n− 13 12 ) ( n− 1 12 ) ( n+ 11 12 ) − · · · ) + · · · , n → ∞. (3.6) See Figure 3 for numerical illustrations of the precision of these subleading power-law corrections to the leading factorial large-order growth. The first few subleading coefficients in (3.6) are{ 1,−97 48 ,−53917 13824 ,−3026443 221184 ,−32035763261 382205952 , . . . } . (3.7) The physical significance of these subleading coefficients in (3.6)–(3.7) is discussed below – see Section 5.1. 100 200 300 400 500 n 0.0144 0.0145 0.0146 0.0147 growth ratio Figure 3. The ratio of the perturbative Hopf An coefficients in (2.3) to the leading order growth, 12nΓ ( n+ 23 12 ) (blue). This blue curve is a zoomed-in view of the blue curve in Figure 2. The other curves in this figure include successively the first 3 subleading power-law corrections in (3.6) (gold, green, red). The asymptotic value is equal to the constant S1/6 in (3.2). The final ellipsis in (3.6) refers to further exponentially suppressed subleading corrections, be- yond the power-law subleading corrections indicated inside the parentheses. These exponentially small corrections cannot be resolved by simple ratio tests and Richardson extrapolations. How- ever, these exponentially small corrections are resolved by the Borel analysis in Section 4, and by the trans-series analysis in Section 5. Physically, they correspond to “higher instanton” non- perturbative terms, while the expansion in (3.6) corresponds to just the leading “one-instanton” term: the leading factor characterizes the one-instanton, while the subleading power-law correc- tions in (3.6) encode the perturbative fluctuations about the single instanton. Note that while the chain approximation expansion (3.3)–(3.5) is also factorially divergent, there is a simple closed-form expression for the expansion coefficients. Using the same normalization of the coef- ficients as in the Hopf expansion (2.3), we see that the chain coefficients defined in (3.3)–(3.5) also grow factorially fast: Achain n = 1 2 · 6nΓ(n) ( 1− 2 2n + 1 3n ) . (3.8) 8 M. Borinsky, G.V. Dunne and M. Meynig This makes it clear that the growth of the chain coefficients is slower than that of the Hopf coefficients in (3.6), because of the factor 6n instead of 12n. The argument of the leading factorial factor, Γ(n), is also different. The exact expression (3.8) also shows that there are exponentially smaller corrections to the leading factorial growth, Achain n ∼ −1 26 nΓ(n), encoded in the 1/2n and 1/3n correction terms in (3.8). The significance of these exponentially smaller corrections will become clear in the Borel analysis of the next section. 4 Borel analysis of the perturbative expansion of the anomalous dimension In this section we use Borel methods, combined with conformal maps, to analyze in more detail the structure of the formal perturbative solution (2.3) to the Hopf-algebraic Dyson–Schwinger equation (2.2). Remarkably, the complicated-looking nonlinear equation (2.2) can be factored:[ G(x) ( 2x d dx − 1 ) − 1 ][ G(x) ( 2x d dx − 1 ) − 2 ][ G(x) ( 2x d dx − 1 ) − 3 ] G(x) = −3x. (4.1) Here we have defined G(x) := γ(−3x) (4.2) rescaling the variable as x := −a 3 , to account for the alternating sign and the power-law growth factor, 12n/62n = 1/3n, coming from (2.3) and the leading growth in (3.1).5 Then the formal perturbative series in (2.3) becomes Gpert(x) := 6 ∞∑ n=1 An 12n xn (4.3) ∼ x 2 + 11x2 24 + 47x3 36 + 2249x4 384 + 356789x5 10368 + 60819625x6 248832 + · · · . (4.4) With this scaling, the coefficients of xn in Gpert(x) have leading growth that is purely factorial, ∼ Γ ( n+ 23 12 ) , with the exponential factor 12n in (3.1) scaled out. We define the corresponding Borel transform: Bpert(t) := 6 ∞∑ n=1 An 12n tn n! . (4.5) The formal perturbative series for G(x) is recovered by the Laplace transform:6 Gpert(x) = 1 x ∫ ∞ 0 dt e−t/xBpert(t). With this choice of scaling, the Borel transform (4.5) has radius of convergence equal to 1. To probe the leading Borel singularity more closely, a natural first step is to use Padé approxi- mants [10]. Figure 4 shows the poles of the diagonal Padé approximant to Bpert(t), truncated 5In terms of the physical expansion variable a = −3x = g2 (4π)3 , the change of sign from a to x is of interest in the context of the Lee–Yang edge singularity for the theory where g is pure imaginary [39, 50, 76, 77]. It is also of interest for the inclusion of vertex diagrams, as will be discussed in future work. For these reasons we choose to work here in terms of the variable x, where exponential terms are most easily identified. 6This formal integral expression is of course to be understood in the Borel–Écalle sense [33, 44, 48, 80, 82]. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 9 1 2 3 4 5 6 7 Re[t] -2 -1 0 1 2 Im[t] Figure 4. The blue points show the poles of the diagonal Padé approximant of the Borel transform (4.5) of the Hopf-algebraic perturbative series (4.3). There is a clear accumulation of Padé poles to the point t = 1+, indicating the presence of a branch point there. There are also hints of further singularities at integer multiples of the leading singularity, with diminishing resolution due to the finite order of Padé. This structure is resolved more clearly in Figures 5 and 6 after using a conformal map. after 500 terms. Recall that a Padé approximant, being a rational approximation of the function as a ratio of polynomials, has only pole singularities. But a Padé approximant of a function with branch points represents a branch point as an accumulation point of poles [34, 35, 87]. Figure 4 suggests that the leading Borel singularity is at t = 1, consistent with the radius of convergence being 1. Furthermore, the accumulation of Padé poles to t = 1+ suggests that this leading singularity is a branch point rather than a pole. This reveals a drawback of the Padé approximant: in attempting to represent a branch cut by a line of poles along the interval t ∈ [1,∞), accumulating to t = 1, it obscures the possible existence of other branch points along this same line. Since the Borel singularities correspond physically to non-perturbative “instanton” terms, with the leading singularity at t = 1 being the “one-instanton” singularity, in a nonlinear problem such as this we expect this leading singularity to be repeated at integer multiples, corresponding to the “multi-instanton” terms. A closer look at Figure 4 hints at the possible existence of other singularities at integer values along the positive t axis, but is unable to resolve them clearly. Fortunately there is a simple way to cure this problem. We can resolve these higher Borel singularities using a conformal mapping method that has been widely used (albeit for different reasons) in the physics literature [24, 25, 59, 93]. The idea is to make a conformal map7 from the cut Borel t plane, based on the leading singularity at t = +1, to the unit disk in the z plane: z = 1− √ 1− t 1 + √ 1− t , t = 4z (1 + z)2 . (4.6) Note that this conformal map does not require knowledge of the nature of the leading Borel singularity, just its location. In fact, the high precision of the conformal mapping step means that it can be used to iteratively refine the location of the leading singularity if it is only known approximately [34, 35]. By construction, any further singularities on the positive Borel t axis, beyond the leading one at t = +1, are mapped to the unit circle in the z plane. The conformal map takes t = +1 to z = +1, and t = +2 to z = ±i, representing both sides of the first cut, and so on. 7Significantly higher precision can be obtained with a uniformization map [35], which is particularly useful if fewer perturbative terms are available. 10 M. Borinsky, G.V. Dunne and M. Meynig This Padé–conformal–Borel procedure is as follows: 1. Re-expand the Borel transform B ( 4z/(1 + z)2 ) about z = 0 inside the unit disk of the conformally mapped z plane, to the same order as the original expansion in the Borel t plane. 2. Make a diagonal Padé approximant to the resulting truncated series in z. 3. Find the singularities (poles) in z of this Padé approximant. Higher branch points in the t plane will now be separated as points of accumulation of poles to the unit circle in the z plane. 4. To obtain an accurate analytic continuation in the original Borel t plane, especially near the branch points, map this Padé approximant in z back to the Borel t plane with the inverse conformal map in (4.6): the resulting analytic continuation of the Borel transform is denoted PCB(t). Figure 5 plots the z poles of the Padé approximant of the conformally mapped series (in step 3, above), demonstrating that multiple branch points, and their associated branch cuts, have been separated and resolved by the conformal map to the conformal z plane [34, 35]. All these branch points lie on the real positive Borel axis, in the interval t ∈ [1,∞). The singularity at z = +1 is the conformal map image of the leading singularity at t = +1. The singularities at z = ±i correspond to the conformal map image of the two-instanton Borel singularity at t = +2, on either side of the branch cut. We can further resolve a third singularity at z = −1 3± 2 √ 2 3 i, which is the conformal map image of the three-instanton Borel singularity at t = +3, once again on either side of the branch cut. There is also weaker evidence of a further singularity at the conformal image of t = +4, corresponding to z = e±2πi/3. This can be resolved by taking more terms in the original perturbative expansion, therefore allowing a higher order Padé approximation. -1.0 -0.5 0.5 1.0 Re[z] -1.0 -0.5 0.5 1.0 Im[z] Figure 5. The blue points indicate poles in the complex z plane of the diagonal Padé approximant of the Borel transform (4.5) after the conformal mapping (4.6), followed by re-expansion about z = 0. The poles accumulate to the conformal map images (black points) of the points t = +1,+2,+3,+4 in the Borel t plane. This illustrates how the conformal map separates and resolves the genuine physical Borel singularities. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 11 This structure of higher Borel singularities can also be seen after mapping the Padé approx- imant back to the Borel t plane (step 4 of the PCB algorithm above). Figure 6 shows a log plot of the imaginary part of the Borel transform just above the real t axis, revealing singu- larities at t = +1,+2,+3,+4 in the Borel t plane, beyond the leading singularity at t = +1. This Padé–conformal–Borel analysis confirms our physical expectation that there should be -2 2 4 6 t 10 20 30 Figure 6. Plot of log (∣∣Im (PCBHopf ( t+ i 50 ))∣∣) for the Padé–conformal–Borel extrapolation. The loga- rithm is plotted in order to show all the singularities on the same figure. The peaks at t = +1,+2,+3,+4 indicate the location of a branch point in the Borel t plane. higher Borel singularities at integer multiples of the leading one at t = +1, corresponding to a “multi-instanton expansion”. Furthermore, these higher Borel singularities all appear to be branch point singularities. These numerical results, and further extensions thereof, are derived analytically using the method of trans-series in Section 5. To conclude this Borel analysis section, we comment further on the comparison with the chain approximation in (3.5), which also produces a divergent formal perturbative expansion. Adopting the same scaling as in the Hopf expansion in (4.2), in terms of the variable x = −a/3, the corresponding formal chain approximation expansion is Gchain(x) := γchain(−3x) ∼ 3 ∞∑ n=1 Γ(n) 2n ( 1− 2 2n + 1 3n ) xn = −24 ∫ ∞ 0 dt (t− 2)(t− 4)(t− 6) e−t/x. (4.7) The final expression here has the form of a Borel integral,8 with Borel singularities consisting of just three simple poles: at t = +2, t = +4 and t = +6. The singularities are poles rather than branch points, and there are only three of them, not an infinite “multi-instanton” tower of integer-spaced singularities. Also note that the leading Borel singularity for the chain expansion is at twice the distance from the origin compared to the leading singularity of the Borel transform for the Hopf expansion. This corresponds directly to the fact that the chain expansion is less divergent than the Hopf expansion: compare (2.3) and (3.4), or (3.1) and (3.8). However, despite these significant differences between the chain and Hopf approximations, we note one interesting similarity: in both cases the first three Borel singularity locations are in the relative proportions 1 : 2 : 3. In the next section a trans-series analysis reveals analytically a similar feature of the Borel transform of the Hopf expansion, and explains the physical significance of this fact. 8This formal integral expression is of course to be understood in the Borel–Écalle sense [33, 44, 48, 80, 82]. 12 M. Borinsky, G.V. Dunne and M. Meynig 5 Trans-series analysis of the Hopf-algebraic Dyson–Schwinger equation The perturbative solution for G(x) in (4.3)–(4.4) cannot be the general solution of the Dyson– Schwinger equation (4.1), because the general solution of this third order equation involves 3 boundary condition parameters. The perturbative expansion in (4.3) has no boundary condi- tion parameters. The resolution of this mis-match is to generalize the divergent formal pertur- bative series in (4.3)–(4.4) to a trans-series expansion, which may involve exponentially small non-perturbative corrections and also possibly logarithmic terms. The trans-series involves three trans-series parameters, which we denote as σ1, σ2, σ3, and which characterize the boundary conditions associated with a three-parameter family of solutions, each of which shares the same divergent formal power series perturbative expansion in (4.3). The full structure of the trans- series can be extracted from the differential equation (4.1), and exhibits a rich network of resurgence relations connecting different non-perturbative sectors [31, 32, 33]. 5.1 Identifying the “seed” exponential terms: the linearized equation The first step in constructing the trans-series from the differential equation (4.1) is to identify the basic exponential “building blocks”, or “seeds”, of the trans-series. To do this, we linearize (4.1) with the ansatz form [31, 32, 33] G(x) ∼ Gpert(x) +Gnon-pert(x), x → 0+. (5.1) Here Gnon-pert(x) is exponentially small, beyond all orders in perturbation theory as x → 0+, also with an accompanying prefactor power of x: Gnon-pert(x) ∼ xβe−λ/x(1 +O(x)), x → 0+. (5.2) We substitute the ansatz (5.1)–(5.2) into the Dyson–Schwinger equation (4.1) and linearize the equation for Gnon-pert(x). This results in the following third order linear and homogeneous equation for Gnon-pert(x): 8x3Gpert(x) 3G′′′ non-pert(x) + ( 32x3Gpert(x) 2G′ pert(x) + 8x2Gpert(x) 3 − 24x2Gpert(x) 2 ) G′′ non-pert(x) + ( 24x3Gpert(x)G ′ pert(x) 2 + 8x2Gpert(x) 2G′ pert(x)− 48x2Gpert(x)G ′ pert(x) + 32x3Gpert(x) 2G′′ pert(x)+4xGpert(x) 3+12xGpert(x) 2+22xGpert(x) ) G′ non-pert(x) + ( 64x3Gpert(x)G ′ pert(x)G ′′ pert(x) + 8x3G′ pert(x) 3 + 8x2Gpert(x)G ′ pert(x) 2 − 24x2G′ pert(x) 2 + 24x2Gpert(x) 2G′′ pert(x)− 48x2Gpert(x)G ′′ pert(x) + 24x3Gpert(x) 2G′′′ pert(x) + 12xGpert(x) 2G′ pert(x) + 24xGpert(x)G ′ pert(x) + 22xG′ pert(x)−4Gpert(x) 3−18Gpert(x) 2− 22Gpert(x)−6 ) Gnon-pert(x) = 0. (5.3) This linear equation for Gnon-pert(x) determines three possible pairs of solutions (since it is a third order equation) for the parameters β and λ in the non-perturbative ansatz (5.2). It is convenient to express these in a vector notation: λ⃗ = (1, 2, 3), β⃗ = ( −23 12 ,+ 1 6 ,−11 4 ) . (5.4) Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 13 The existence of three different values of the parameter λ corresponds to the fact that for the perturbative Hopf expansion (4.2) there are three Borel singularities at different locations: at t = 1, 2, 3. This explains analytically the existence of Borel singularities at t = 1, 2, 3, which were found numerically by the conformal mapping analysis of the Borel transform in the previous section. Recall Figures 5 and 6. Note also the similarity to, and difference from, the Borel singularities of the chain approximation expansion, which appear at t = 2, 4, 6 (in the same normalization: see equation (4.7)). The linearized equation (5.3) for Gnon-pert(x) also generates the corresponding subleading power-law factor, xβ, characterized by the β parameter in (5.2). The fact that the three values of β in (5.4) are all different and are all non-integer implies that the corresponding Borel singu- larities at t = 1, 2, 3 are branch point singularities, each with a different exponent [31, 32, 33]. We denote the solutions of the linearized equation as: Gnon-pert, linearized k⃗ (x) := ( k⃗ · σ⃗ ) xk⃗·β⃗e−(k⃗·λ⃗)/xF k⃗ (x). (5.5) Here we label the three different solutions (5.5) to the third order linearized equation by the three basis unit vectors k⃗ = (1, 0, 0), k⃗ = (0, 1, 0), k⃗ = (0, 0, 1). (5.6) Thus, β1 = (1, 0, 0) · β⃗ = −23 12 , and λ1 = (1, 0, 0) · λ⃗ = 1, and so on. Notice that each of the three “seed” solutions in (5.5) has a free multiplicative parameter, σ1, σ2, σ3, since each solves the third order homogeneous linear equation (5.3). For these seed solutions of the linearized equation, the final factor in (5.5), denoted as F k⃗ (x), is a formal fluctuation series (also labeled by k⃗):9 F k⃗ (x) := ∞∑ n=1 ak⃗nx n−1. (5.7) The coefficients ak⃗n in (5.7) are generated recursively by simple substitution of the ansatz (5.5) into the differential equation (4.1){ a(1,0,0)n } = { −1, 97 48 , 53917 13824 , 3026443 221184 , 32035763261 382205952 , . . . } , (5.8) { a(0,1,0)n } = { −1, 151 24 ,−63727 3456 , 7112963 82944 ,−7975908763 23887872 , . . . } , (5.9) { a(0,0,1)n } = { −1, 227 48 , 1399 4608 , 814211 73728 , 3444654437 42467328 , . . . } . (5.10) � We observe that the a (1,0,0) n coefficients in (5.8) coincide with (up to an overall minus sign) the coefficients of the subleading corrections to the large-order growth of the perturbative coefficients An in (3.6)–(3.7). This can be confirmed to very high subleading order. This is a clear example of the generic low-order/large-order resurgent behavior connecting the large order growth of the perturbative series to the low orders of the fluctuations around the first instanton term [12]. � We also observe that the pre-factor exponent, (1, 0, 0) · β⃗ = −23 12 , associated with the leading exponential term, e−((1,0,0)·λ⃗)/x = e−1/x, coincides with the offset of the argument of the leading factorial growth of the perturbative An coefficients in (3.1). This is also an example of generic behavior relating a formal perturbative series with the leading exponential correction. 9We index the coefficients ak⃗ n beginning with n = 1 to match the indexing convention for the perturbative coefficients An in (2.3). 14 M. Borinsky, G.V. Dunne and M. Meynig We record the leading large order behavior of the expansion coefficients in (5.8)–(5.10). We generated 500 coefficients (5.8)–(5.10) (all rational numbers) for each of the fluctuation expansions in (5.7). Each series is factorially divergent. The leading large-order growth of the expansion coefficients is10 a(1,0,0)n = 4S1Γ ( n+ 23 12 )( 1 +O ( n−1 )) , n → ∞, (5.11) a(0,1,0)n = 4S1Γ ( n+ 23 12 )( 1 +O ( n−5/6 )) , n → ∞, (5.12) a(0,0,1)n = 20 3 S1Γ ( n+ 23 12 )( 1 +O ( n−1 )) , n → ∞. (5.13) Note that the factorial growth factor is the same for each series, and also agrees with that of the perturbative series; recall (3.1). Also note that the overall Stokes constants are expressed as sim- ple rational multiples of the Stokes constant S1 for the perturbative series: recall (3.1) and (3.2). Thus, at this leading order level, no new independent Stokes constant is generated. Corrections to these leading large order growth expressions in (5.11)–(5.13) are discussed in Section 5.3. 5.2 Beyond the linearized equation: resonant trans-series and logarithms The three exponential terms, e−λ1/x, e−λ2/x and e−λ3/x, are just the “seed” non-perturbative terms coming from the linearized equation (5.3) for Gnon-pert(x). But the full Dyson–Schwinger equation is nonlinear, so each of these three exponential terms will re-appear in all integer powers of the seed term, generating a trans-series expansion that includes both perturbative and non-perturbative terms to all orders. For a generic non-resonant trans-series [31, 32, 33], in which there are no integer relations between the different λj values in (5.4), the trans-series would be generated as a three dimensional infinite sum generating all powers of the basic seed factors, xk⃗·β⃗e−(k⃗·λ⃗)/x, each being further multiplied by a formal fluctuation series. The terms in such a non-resonant trans-series would all have the form indicated in (5.5), but now labeled by vectors k⃗ being integer linear combinations of the basis vectors in (5.6) [31, 32, 33]. However, our case here is resonant. This is because the Dyson–Schwinger equation (4.1) has the special resonant property that the three different possible values for the exponent λ in (5.2) satisfy integer relations: λ1 : λ2 : λ3 = 1 : 2 : 3. This resonant property has profound implications for the structure of the trans-series beyond the leading exponential order, leading to an even richer structure.11 In a resonant case, the exponent coefficient (k⃗ · λ⃗) in (5.5) may take the same value for different integer-valued vectors k⃗. For example, e−2/x appears through one power of the seed term with λ2 = 2, but it also appears via the square of the seed term with λ1 = 1. Therefore, when we grade the solution by its exponential order, e−(integer)/x, a given order can have contributions from different k⃗ vectors, and when they mix there can also appear logarithmic terms in the solution to the Dyson–Schwinger equation. An analogous feature is familiar for linear ODEs, but for a nonlinear ODE there will appear higher powers of logarithms as we go higher in the exponential order of the trans-series (i.e., to higher “instanton” order). 10Notice that the coefficients a (0,1,0) n alternate in sign at low order, but eventually settle down to have the same sign. 11A familiar and illustrative example is the Painlevé I equation, a second-order nonlinear equation which has two resonant values λ = ±1, suitably normalized [4, 51]. Here the resonant structure is quite different. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 15 To characterize this resonant structure we write the trans-series expansion in an exponentially graded form G(x) ∼ G(0)(x)+ e−1/xG(1)(x)+ e−2/xG(2)(x)+ e−3/xG(3)(x)+ e−4/xG(4)(x) + · · · . (5.14) Here G(0)(x) ≡ Gpert(x), is the formal perturbative series in (4.2)–(4.3). By construction, each G(m)(x), for m ≥ 1, satisfies a third order linear differential equation. For the first expo- nential term, G(1)(x), this equation is both linear and homogeneous, and the solution is: G(1)(x) = σ1x −23/12F(1,0,0)(x). (5.15) This corresponds to the first of the “seed” solutions (5.5) from the previous section, with the exponential factor separated out as in (5.14). Up to this first exponential order there are no logarithmic terms. At second exponential order, and beyond, the equation for G(m)(x), with m ≥ 2, is linear but inhomogeneous. We thus need a particular solution in addition to the homogeneous solution. For m = 2 the inhomogeneity comes from a term involving ( G(1)(x) )2 . This means that the inhomogeneity vanishes if G(1)(x) vanishes. This is the case if the first trans-series parameter σ1 is chosen to vanish. In this case G(2)(x) satisfies a linear homogeneous equation, and there exists a solution, which is multiplied by a new constant parameter, the second trans-series parameter σ2: Ghomogeneous (2) (x) = σ2x 1/6F(0,1,0)(x), with σ1 = 0. (5.16) Noting that (0, 1, 0) · β⃗ = 1 6 , we recognize this as corresponding to the second “seed” solution in (5.5), once again with the exponential factor separated out. However, if σ1 ̸= 0, the full solu- tion to the inhomogeneous equation for G(2)(x) requires also a particular solution. To generate this solution we substitute a Frobenius ansatz: G(2)(x) ∼ xδ [∑ n ∆nx n + (∑ n µnx n ) log(x) ] . The power parameter δ, and the expansion coefficients ∆n and µn are determined recursively by the linear inhomogeneous equation for G(2)(x). This leads to the solution G(2)(x) ∼ ( 1 x23/12 )2[ −2σ2 1 1 x F(2,0,0)(x) + x4 ( σ2 1 21265 2304 log(x) + σ2 ) F(0,1,0)(x) ] , (5.17) where we find a new fluctuation series F(2,0,0)(x) ( the labeling notation refers to the two powers of σ1 multiplying this term: σ2 1 = σ2 1σ 0 2σ 0 3 := σ⃗(2,0,0) ) : F(2,0,0)(x) ∼ 1− 49 12 x− 13235 3456 x2 − 43049 3456 x3 − 2496477497 23887872 x4 − 0 · x5 − 3315185066507813 247669456896 x6 − · · · . (5.18) Several comments are in order concerning the structure of the expression for G(2)(x) in (5.17)– (5.18): 1. The overall rational-exponent prefactor in G(2)(x) is naturally expressed in terms of the square of the corresponding rational-exponent prefactor for G(1)(x) in (5.15). 2. Noting that ( 1 x23/12 )2 x4 = x1/6, we see that when σ1 = 0 we recover the second “seed” solution Ghomogeneous (2) (x) in (5.16), with fluctuation series F(0,1,0)(x), and multiplied by its new trans-series parameter σ2 ≡ (0, 1, 0) · σ⃗. 16 M. Borinsky, G.V. Dunne and M. Meynig 3. When σ1 ̸= 0 there appears a log(x) term, proportional to σ2 1, and multiplied by the same “seed” fluctuation series F(0,1,0)(x) mentioned in the previous item, and with a specific fixed coefficient: 21265 2304 x 4. 4. When σ1 ̸= 0 there is also a new series, once again proportional to σ2 1, denoted as F(2,0,0)(x), and whose first terms are listed in (5.18). 5. The coefficient of x5 in F(2,0,0)(x) vanishes. Note that because the value of σ2 is not fixed and behaves like an integration constant, the fact that this coefficient vanishes is not a coincidence, but a convenient choice. In fact, we could have fixed this coefficient to any value while appropriately changing higher order terms of F(2,0,0)(x). At third exponential order in the graded trans-series (5.14), the function G(3)(x) satisfies a linear and inhomogeneous equation, and the inhomogeneity involves both (G(1)(x)) 3 and G(1)(x)G(2)(x). Thus, setting σ1 = 0, we obtain a homogeneous linear equation for G(3)(x), with solution: Ghomogeneous (3) (x) = σ3x −11/4F(0,0,1)(x), with σ1 = 0. (5.19) We recognize this as corresponding to the third “seed” solution in (5.5), with k⃗ = (0, 0, 1) and with the third independent factor σ3 (and once again with the exponential factor separated out). But when σ1 ̸= 0, resonant logarithmic terms appear again. For G(3)(x) we find an exp- ression involving two new series, denoted F(3,0,0)(x) and F(1,1,0)(x), multiplied by σ3 1 and σ1σ2, respectively: G(3)(x) ∼ ( 1 x23/12 )3[ σ3 1 ( − 6 x2 ) F(3,0,0)(x) + x3 ( −5σ3 1 21265 2304 log(x) + σ3 ) F(0,0,1)(x) + 12x4 ( σ3 1 21265 2304 log(x) + σ1σ2 ) F(1,1,0)(x) ] . (5.20) The two new series, F(3,0,0)(x) and F(1,1,0)(x), appearing in (5.20) have initial terms: F(3,0,0)(x) ∼ 1− 103x 16 + 8821x2 4608 − 454379x3 221184 − 1344528799x4 14155776 + 0 · x5 − 9242013290874467x6 587068342272 + · · · , (5.21) F(1,1,0)(x) ∼ 1− 131x 16 + 153997x2 4608 − 12555605x3 73728 + 484910403x4 524288 − 4580515441493x5 679477248 + 28651912194246539x6 587068342272 − 30305985730060738841x7 65751654334464 + · · · . (5.22) We comment on the structure of the expression for G(3)(x) in (5.20)–(5.22): 1. The overall rational-exponent prefactor in G(3)(x) is naturally expressed in terms of the cube of the corresponding rational-exponent prefactor for G(1)(x) in (5.15). 2. Noting that ( 1 x23/12 )3 x3 = x−11/4, we see that when σ1 = 0 we recover the third “seed” solution Ghomogeneous (3) (x) in (5.19), with fluctuation series F(0,0,1)(x), and multiplied by its new trans-series parameter σ3 ≡ (0, 0, 1) · σ⃗. 3. When σ1 ̸= 0 there appears a log(x) term, proportional to σ3 1. Remarkably, this logarithmic term is multiplied by a rational-exponent factor x−11/4 = x(0,0,1)·λ⃗, and by a fluctuation series that naturally splits into two pieces. One piece coincides with the fluctuation series F(0,0,1)(x) appearing in the homogeneous solution Ghomogeneous (3) (x) in (5.19), while the other piece is x times the fluctuation factor F(1,1,0)(x) which multiplies the σ1σ2 ≡ σ⃗(1,1,0) factor in the full solution (5.20). Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 17 4. When σ1 ̸= 0 there is also a new series, proportional to σ3 1, denoted as F(3,0,0)(x), and whose first terms are listed in (5.21). Note that the coefficient of x5 in F(3,0,0)(x) vanishes, just like for F(2,0,0)(x) in (5.18). Proceeding to higher exponential orders of the graded trans-series in (5.14), one finds that there are no new homogeneous solutions, even after setting σ1 = σ2 = σ3 = 0. This is because the three independent homogeneous solutions were generated by the linearized equation discussed in Section 5.1. Hence no new independent trans-series parameters are generated, which is consistent with the interpretation of (σ1, σ2, σ3) as the three boundary condition parameters of the third order Dyson–Schwinger equation (4.1). The inhomogeneity in the linear equation for the 4th order graded term G(4)(x) involves various combinations of lower-order functions: (i) 4 factors involving G(1)(x); (ii) 2 factors involving G(1)(x) and 1 factor involving G(2)(x); (iii) 1 factor involving G(1)(x) and 1 factor involving G(3)(x); (iv) 2 factors involving G(2)(x). Since one power of log(x) appears in G(2)(x), we find that the solution for G(4)(x) involves the first appearance of a log2(x) term. This can be confirmed from the ODE by direct substitution. In general, continuing to all orders in the exponential grading we see that eventually all powers of log(x) are generated. However, at any given exponential order only a finite number of powers of log(x) appear. This trans-series structure is a direct consequence of the resonant character of the Dyson–Schwinger equation (4.1). Thus, the full trans-series solution involves all powers of the three basic trans-monomial elements: x, e−1/x, and log(x),12 with log(x) terms appearing only at the second exponential order and beyond. The physical significance of this fact is discussed in the conclusions. 5.3 Logarithmic behavior in large-order growth The appearance of log(x) terms in the trans-series starting at the second exponential order of the graded trans-series (5.14) is reflected in the appearance of log(n) corrections appearing in the subleading large-order growth of the coefficients of the fluctuation series beyond the original perturbative series G(0)(x) ≡ Gpert(x). 13 To illustrate this phenomenon, consider the coefficients a (1,0,0) n of the fluctuation series F(1,0,0)(x), which appears in the first exponential term G(1)(x) in (5.15). The first few coefficients were already listed in (5.8) and their leading factorial growth was shown in (5.11). To probe the subleading corrections to this leading factorial growth, it is instructive to begin with a conventional resurgent ansatz involving power-law subleading corrections: Γ ( n+ 23 12 )[ 1 + ∞∑ k=1 bk∏k l=1 ( n+ 23 12 − l )]. (5.23) Here the bk are expected to be rational numbers related to coefficients of fluctuation expansions higher in the trans-series. We will soon see that this ansatz (5.23) is in fact insufficient, but remarkably it is extremely precise for the first 4 subleading power-law corrections. Using 500 terms of the a (1,0,0) n coefficients, we use high-order Richardson extrapolation meth- ods to identify the first few such rational coefficients: a(1,0,0)n ∼ 4S1Γ ( n+ 23 12 )[ 1− 49 12( n+ 11 12 ) − 13235 3456( n− 1 12 )( n+ 11 12 ) − 43049 3456( n− 13 12 )( n− 1 12 )( n+ 11 12 ) − 2496477497 23887872( n− 25 12 )( n− 13 12 )( n− 1 12 )( n+ 11 12 ) − · · · ] . (5.24) 12The fact that only the three basic trans-monomial elements, x, e−1/x, and log(x), suffice to construct the trans-series solution of the ODE follows from [31, 32, 33]. 13Analogous logarithmic large-order growth is seen in the Painlevé I equation [4, 51]. 18 M. Borinsky, G.V. Dunne and M. Meynig Figure 7 illustrates the extra precision gained by including these subleading corrections; compare with the analogous plot in Figure 3 for the perturbative series coefficients. The precision is extremely good for the first four subleading corrections displayed in (5.24). 100 200 300 400 n 0.340 0.345 0.350 0.355 ratio Figure 7. Plot of the ratio of the coefficients a (1,0,0) n to the leading growth Γ ( n+ 23 12 ) in (5.11) (blue curve), and then subsequent curves include the further subleading large n corrections shown in (5.24). The large n values tend to 4S1 = 0.3503822 . . . , in terms of the Stokes constant of the large order growth of the original perturbative coefficients in (3.1)–(3.2). Notice that the first coefficients of the subleading correction terms for large-order growth of the coefficients a (1,0,0) n , as shown in (5.24), coincide with the low order coefficients of the fluctua- tion series F(2,0,0)(x) in (5.18). Thus, superficially this looks like a generic large-order/low-order resurgence relation connecting the large orders of the series coefficients entering the e−1/xG(1)(x) term in the graded trans-series (5.14) and the e−2/xG(2)(x) term in (5.14). However, this cor- respondence breaks down at the next subleading order. This can be traced to the fact that at this order in G(2)(x) we find a logarithmic term x5 log(x). Moreover, the coefficient of this logarithmic term is fixed in (5.17). This suggest that we extend the large-order growth ansatz in (5.23) to incorporate a log(n) growth factor: a(1,0,0)n ∼ 4S1Γ ( n+ 23 12 )[ 1− 49 12( n+ 11 12 ) − 13235 3456( n− 1 12 ) ( n+ 11 12 ) − 43049 3456( n− 13 12 ) ( n− 1 12 ) ( n+ 11 12 ) − 2496477497 23887872( n− 25 12 ) ( n− 13 12 ) ( n− 1 12 ) ( n+ 11 12 ) + d− 21265 4608 log(n)( n− 37 12 ) ( n− 25 12 ) ( n− 13 12 ) ( n− 1 12 ) ( n+ 11 12 ) + · · · ] . (5.25) The logarithmic term in front of log(n) can be found using a modified Richardson acceleration method, as described in Appendix A. The coefficient is found to be rational and in correspon- dence with the rational factor 21265 2304 multiplying the log(x) term in G(2)(x) in (5.17). The con- stant coefficient d in (5.25) can similarly be determined to extremely high precision but does not appear to be rational: d = −846.3135357762392334586752585470377294931556414335975009793245 . . . . (5.26) Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 19 6 Conclusions In this paper we have shown that the non-perturbative completion of the perturbative Hopf- algebraic Dyson–Schwinger equation (4.1) for the anomalous dimension of the massless ϕ3 scalar QFT in 6 dimensions leads to a trans-series of the following form: G(x) ∼ ∞∑ n=0 ∞∑ k=0 Lk∑ l=0 cn,k,lx n+1−k ( e−1/x x23/12 )k( log(x) )l . (6.1) We stress that this trans-series was derived directly from the Dyson–Schwinger equation (4.1), which was itself derived from the manifestly perturbative Hopf-algebraic renormalization for- malism. No semi-classical non-perturbative argument was invoked; just a resurgent asymp- totic analysis of the Dyson–Schwinger differential equation. Nevertheless, we observe that the resulting trans-series has the structure of a multi-instanton expansion, typical of semi- classical computations in quantum mechanics and quantum field theory. Each multi-instanton term, e−k/x, is multiplied by a linear combination of formal asymptotic perturbative series in powers of x, and polynomial factors in powers of logarithms of x. The first logarithmic term arises at two-instanton order, and at a given instanton order there is a maximal power Lk = [k/2] of log(x) factors. Remarkably, this structure matches the trans-series structure found using instanton-calculus methods in quantum mechanical models and also in certain quantum field theories [45, 47, 51, 93, 94, 95]. We note that other infinitely-iterated approximations to the Dyson–Schwinger equations, such as the rainbow or chain approximation, do not lead to a multi- instanton trans-series structure. In semi-classical instanton-calculus computations the logarithmic terms arise from quasi-zero modes associated with bion solutions, which can be interpreted as molecules of instantons and anti-instantons [47, 93]. By contrast, our starting point was purely perturbative, with the non- perturbative features being recovered from a trans-series analysis of the Hopf-algebraic Dyson– Schwinger equation. The sum over instanton-like exponential terms arises from the generic integer-spaced repetition of Borel singularities in the solution of the nonlinear Dyson–Schwinger differential equation for the anomalous dimension, while the logarithmic terms in (6.1) come from the special resonant structure of the nonlinear differential equation. This nonlinear equation is third order, so generically there would be three independent Borel singularities [31, 32, 33], therefore associated with three independent “instantons”. However, the Dyson–Schwinger equa- tion has extra symmetry that results in the three Borel singularities being collinear, and also in the integer proportions 1 : 2 : 3. This resonant form can be traced to the particular iterative structure of the Hopf-algebraic Dyson–Schwinger equations, and has interesting implications for the resurgent large-order/low-order relations. Resonance produces log(n) factors in the large- order growth of the coefficients of the fluctuations about the instanton factors, but absent in the large order growth of the original formal perturbative series. Similar log(n) behavior is seen in the large order growth of fluctuation coefficients at higher orders of a semiclassical instanton calculus computation. While there is a long history of studying logarithmic terms in semiclassical approximations in quantum mechanical and quantum field theoretics models, and in resonant nonlinear differential equations, it is much more rare to be able to identify such behavior in terms of the renormalized coupling, and especially in an asymptotically free QFT. Renormalons, which are also associated with iterative perturbative structures, have been studied recently using ideas from resurgence, in a wide variety of theories: see for example [2, 43, 58, 70, 72, 73, 74], and references therein. There are interesting analogies between the iterative Hopf Dyson–Schwinger formalism and renormalon and large N computations in QFT, and these issues are currently under investigation. It would also be valuable to understand in detail the relation to a quite different Hopf approximation for the ϕ3 6 theory, which includes some aspects of vertex corrections [8]. Ultimately the goal is 20 M. Borinsky, G.V. Dunne and M. Meynig to understand more deeply how the non-perturbative trans-series structure fits naturally within the underlying perturbative Hopf algebra structure of QFT renormalization. The Dyson–Schwinger ODE (4.1) could also be analysed using Écalle’s alien differential cal- culus [48, 80]. For the 4 dimensional Yukawa theory analysed in [18], an alien calculus approach led to a closed form expression for the all orders trans-series for the anomalous dimension. In this Yukawa model the trans-series contained no logarithmic terms, so the alien derivative formalism of [15] was well adapted to derive the closed form solution. For the 6 dimensional ϕ3 model, the appearance of resonant logarithmic terms produces a structure that lies beyond the simplified alien calculus framework of [15]. This motivates a generalization of [15] that uses more aspects of Écalle’s full theory of resurgence, a project beyond the scope of this article, and left for future work. Such a formalism would have many potential applications. It could enable a systematic evaluation of the large order behaviour at all orders, as was done in the Yukawa model [18]. An alien calculus treatment should also prove to all orders the correspondence between the trans-series fluctation coefficients and the large order behaviour of the expansion coefficients. A rigorous treatment of the log term asymptotic behaviour would be especially interesting as there are intriguing parallels to resonance phenomena in quantum mechanics, where the appear- ance of the log terms can be associated with the inherent stability properties of the underlying quantum mechanical system. A Modified Richardson method for probing logarithmic large-order growth In this appendix we explain an effective general method to obtain high-precision numerical results for the coefficients of high-order growth terms involving log(n) and powers of log(n). For applications to log(n) growth see, for example, [51, 92]. We first summarize the conven- tional Richardson extrapolation method [10] in a form that makes the generalization simple to formulate, and also simple to implement. As input data, we have a sequence of numbers fn, whose asymptotic behaviour is to be determined experimentally. Even though it is not technically necessary, it increases the efficacy of the method dramatically if the numbers fn are given as explicit rational numbers. If the numbers are not represented as rational numbers, they have to be available up to very high precision for the method to give reliable results. To apply the original version of Richardson extrapolation, we need to assume that the numbers fn have the asymptotic behaviour, fn ∼ K∑ k=0 akn −k +O ( n−K−1 ) as n → ∞, (A.1) for some arbitrary K ≥ 0 with certain (not necessarily known) coefficients ak. This condition is satisfied by the subleading correction factor in the conventional resurgent large-order growth ansatz (5.23) truncated to order K, and matches, for example, the subleading large-order growth in (3.6). In determining the coefficients of such subleading corrections, the technical task is to extract effectively the coefficients an in (A.1), starting with a0 and then proceeding to higher coefficients. The precision should be high enough to be able to recognize rational values of the coefficients, whose rationality can then be tested by probing further corrections. It follows immediately from (A.1) that fn ∼ a0 + O ( n−1 ) . Therefore, if we naively use the value of fn for some large n as an estimate for an, the error will be of the order 1 n . Due to this slow rate of convergence, we would need a huge number of coefficients fn to get a sufficiently accurate estimate. Richardson extrapolation provides a more efficient way to determine the value of a0 at higher precision with fewer input terms. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 21 Let ∆n be the (forward) difference operator, i.e., for any sequence fn ∆nfn = fn+1 − fn. (A.2) Clearly, the operator ∆n is linear and can be iterated. We treat it as a derivative operator that acts on everything on its right. For instance, ∆2 nfn = ∆n(∆nfn) = ∆nfn+1 −∆nfn = fn+2 − 2fn+1 + fn. With this operator the K-th order Richardson extrapolation of a sequence fn can be defined as, RK [fn] := 1 K! ∆K n nKfn. (A.3) This form of Richardson extrapolation is well adapted to straightforward implementation beca- use the difference operator is a natural computer operation. The key observation concerning extrapolation is the fact that if the sequence fn has the asymptotic form in (A.1), then RK [fn] = 1 K! ∆K n nKfn ∼ a0 +O ( n−K−1 ) as n → ∞. (A.4) Thus, the Richardson extrapolated sequence RK [fn] converges much more rapidly, with order O ( n−K−1 ) . If K is chosen sufficiently large, the Richardson extrapolation (A.3) can therefore be used to estimate a0 to very high accuracy. To verify (A.4), observe that by (A.1) the sequence nNfn has the asymptotic form, nKfn ∼ p(n) + ∞∑ k=1 aK+kn −k as n → ∞, (A.5) where p(n) = a0n K+O ( nK−1 ) is a polynomial of orderK. In order to prove (A.4) from (A.5), we need to use the definition (A.3) and also the following elementary properties of the ∆n-operator: ∆K n nK = K!, ∆K n nk = 0 for all 0 ≤ k < K, ∆K n n−k ∼ O ( n−K−k ) for all k > 0 as n → ∞. (A.6) Having formulated conventional Richardson extrapolation in terms of the forward difference operator (A.2), it is now straightforward to construct generalizations that can be applied to sequences with a more general asymptotic behaviour than that in (A.1). For example, suppose a sequence gn has the following asymptotic behaviour with a subleading log n contribution: gn ∼ K∑ k=0 akn −k + log n K∑ k=1 bkn −k +O ( n−K−1 log n ) as n → ∞. Naive application of the Richardson extrapolation operator does not lead to an accelerated convergence, due to the extra log-terms. To derive the variant of Richardson extrapolation that can also deal with this more general asymptotic behaviour, we first evaluate some identities for the difference operator applied to log-terms. We have, ∆K n nk−1 log n = ∆K−1 n (( (n+ 1)k−1 − nk−1 ) log n+ (n+ 1)k−1 log ( 1 + 1 n )) ∼ K∑ ℓ=K−k+1 c̃ℓn −ℓ +O ( n−K−1 ) for all 1 ≤ k ≤ K as n → ∞, 22 M. Borinsky, G.V. Dunne and M. Meynig for some coefficients c̃ℓ whose specific values are not important for this discussion, and where we used the expansion log ( 1 + 1 n ) = − ∑∞ k=1 (−1)k k n−k. Similarly, we obtain ∆K n n−k log n ∼ O ( n−k−K log n ) for all k > 0 as n → ∞. Using these observations together with the previous ones in (A.6), we find that the application of the Richardson extrapolation operator ∆K n nK on the sequence gn leads to a sequence with the following asymptotic behaviour: RK [gn] = 1 K! ∆K n nKgn ∼ a0 + K∑ k=1 ckn −k +O ( n−K−1 log n ) as n → ∞, with some coefficients ck. The log terms are suppressed now, so we can apply the normal Richardson RK operator in (A.3) again, in order to get rid of the subleading non-log-terms with equation (A.4), and obtain a rapidly converging sequence, RK [RK [gn]] = 1 K!2 ∆K n nK∆K n nKgn ∼ a0 +O ( n−K−1 log n ) as n → ∞. This can be used for efficient extraction of the a0 term. Note that this double Richardson extrapolation operator is equivalent to a difference operator of order 2K. Analogously, we can derive related difference operators that produce rapidly converging sequ- ences for other kinds of asymptotic behaviours. If, for instance, the sequence hn is expected to behave asymptotically as follows with a leading log-term, hn ∼ log n K∑ k=0 akn −k + K∑ k=0 bkn −k +O ( n−K−1 log n ) as n → ∞, and we are again interested in the leading a0 coefficient, then the following 2K + 1-th order difference operator applied to hn leads to the desired result: 1 K!2 ∆K n nK+1∆K+1 n nKhn ∼ a0 +O ( n−K−1 log n ) as n → ∞. The proof works analogously to the previous derivations. This is the operator that we used to verify the expected asymptotic behaviour in (5.25), obtaining sufficient precision to clearly identify the rational coefficient 21265 4608 , which is associated with the coefficient of the log(x) term in G(2)(x) in (5.17). Difference operators that can deal with sequences involving higher power logℓ-terms can be easily developed in a similar fashion. Acknowledgements This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of High Energy Physics under Award Number DE-SC0010339 (GD, MM) and by the NWO Vidi grant 680-47-551 “Decoding Singularities of Feynman graphs” (MB). This work was begun during visits by the first two authors to Humboldt University in 2018 and 2019, and at the Les Houches Summer School in 2018, and we thank these institutions for hospitality. We are grateful to Marc Bellon, David Broadhurst, Ovidiu Costin, John Gracey, Dirk Kreimer, Enrico Russo and Karen Yeats for discussions. We also want to thank David Broadhurst for helping with the estimation of the constant d in equation (5.26), and for pointing out some typos in a previous version. Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 23 References [1] Álvarez G., Howls C.J., Silverstone H.J., Anharmonic oscillator discontinuity formulae up to second- exponentially-small order, J. Phys. A: Math. Gen. 35 (2002), 4003–4016. [2] Anber M.M., Sulejmanpasic T., The renormalon diagram in gauge theories on R3×S1, J. High Energy Phys. 2015 (2015), no. 1, 139, 34 pages, arXiv:1410.0121. [3] Aniceto I., Başar G., Schiappa R., A primer on resurgent transseries and their asymptotics, Phys. Rep. 809 (2019), 1–135, arXiv:1802.10441. [4] Aniceto I., Schiappa R., Vonk M., The resurgence of instantons in string theory, Commun. Number Theory Phys. 6 (2012), 339–496, arXiv:1106.5922. [5] Bellon M.P., An efficient method for the solution of Schwinger–Dyson equations for propagators, Lett. Math. Phys. 94 (2010), 77–86, arXiv:1005.0196. [6] Bellon M.P., Clavier P.J., Alien calculus and a Schwinger–Dyson equation: two-point function with a non- perturbative mass scale, Lett. Math. Phys. 108 (2018), 391–412, arXiv:1612.07813. [7] Bellon M.P., Russo E.I., Resurgent analysis of Ward–Schwinger–Dyson equations, SIGMA 17 (2021), 075, 18 pages, arXiv:2011.13822. [8] Bellon M.P., Russo E.I., Ward–Schwinger–Dyson equations in ϕ3 6 quantum field theory, Lett. Math. Phys. 111 (2021), 42, 31 pages, arXiv:2007.15675. [9] Bellon M.P., Schaposnik F.A., Renormalization group functions for the Wess–Zumino model: up to 200 loops through Hopf algebras, Nuclear Phys. B 800 (2008), 517–526, arXiv:0801.0727. [10] Bender C.M., Orszag S.A., Advanced mathematical methods for scientists and engineers. I. Asymptotic methods and perturbation theory, Springer-Verlag, New York, 1999. [11] Benedetti D., Delporte N., Remarks on a melonic field theory with cubic interactions, J. High Energy Phys. 2021 (2021), no. 4, 197, 30 pages, arXiv:2012.12238. [12] Berry M.V., Howls C.J., Hyperasymptotics for integrals with saddles, Proc. Roy. Soc. London Ser. A 434 (1991), 657–675. [13] Bersini J., Maiezza A., Vasquez J.C., Resurgence of the renormalization group equation, Ann. Physics 415 (2020), 168126, 15 pages, arXiv:1910.14507. [14] Bogner C., Borowka S., Hahn T., Heinrich G., Jones S.P., Kerner M., von Manteuffel A., Michel M., Panzer E., Papara V., Loopedia, a database for loop integrals, Comput. Phys. Commun. 225 (2018), 1–9, arXiv:1709.01266. [15] Borinsky M., Generating asymptotics for factorially divergent sequences, Electron. J. Combin. 25 (2018), 4.1, 32 pages, arXiv:1603.01236. [16] Borinsky M., Graphs in perturbation theory: algebraic structure and asymptotics, Springer Theses, Springer, Cham, 2018. [17] Borinsky M., Tropical Monte Carlo quadrature for Feynman integrals, arXiv:2008.12310. [18] Borinsky M., Dunne G.V., Non-perturbative completion of Hopf-algebraic Dyson–Schwinger equations, Nuc- lear Phys. B 957 (2020), 115096, 17 pages, arXiv:2005.04265. [19] Borinsky M., Gracey J.A., Kompaniets M.V., Schnetz O., Five-loop renormalization of ϕ3 theory with applications to the Lee–Yang edge singularity and percolation theory, Phys. Rev. D 103 (2021), 116024, 35 pages, arXiv:2103.16224. [20] Borinsky M., Schnetz O., Graphical functions in even dimensions, arXiv:2105.05015. [21] Brezin E., Le Guillou J.C., Zinn-Justin J., Perturbation theory at large order. I. The ϕ2N interaction, Phys. Rev. D 15 (1977), 1544–1557. [22] Broadhurst D.J., Kreimer D., Combinatoric explosion of renormalization tamed by Hopf algebra: 30-loop Padé–Borel resummation, Phys. Lett. B 475 (2000), 63–70, arXiv:hep-th/9912093. [23] Broadhurst D.J., Kreimer D., Exact solutions of Dyson–Schwinger equations for iterated one loop integrals and propagator coupling duality, Nuclear Phys. B 600 (2001), 403–422, arXiv:hep-th/0012146. [24] Caliceti E., Meyer-Hermann M., Ribeca P., Surzhykov A., Jentschura U.D., From useful algorithms for slowly convergent series to physical predictions based on divergent perturbative expansions, Phys. Rep. 446 (2007), 1–96, arXiv:0707.1596. https://doi.org/10.1088/0305-4470/35/18/302 https://doi.org/10.1007/JHEP01(2015)139 https://arxiv.org/abs/1410.0121 https://doi.org/10.1016/j.physrep.2019.02.003 https://arxiv.org/abs/1802.10441 https://doi.org/10.4310/CNTP.2012.v6.n2.a3 https://doi.org/10.4310/CNTP.2012.v6.n2.a3 https://arxiv.org/abs/1106.5922 https://doi.org/10.1007/s11005-010-0415-3 https://doi.org/10.1007/s11005-010-0415-3 https://arxiv.org/abs/1005.0196 https://doi.org/10.1007/s11005-017-1016-1 https://arxiv.org/abs/1612.07813 https://doi.org/10.3842/SIGMA.2021.075 https://arxiv.org/abs/2011.13822 https://doi.org/10.1007/s11005-021-01377-2 https://arxiv.org/abs/2007.15675 https://doi.org/10.1016/j.nuclphysb.2008.02.005 https://arxiv.org/abs/0801.0727 https://doi.org/10.1007/978-1-4757-3069-2 https://doi.org/10.1007/JHEP04(2021)197 https://arxiv.org/abs/2012.12238 https://doi.org/10.1098/rspa.1991.0119 https://doi.org/10.1016/j.aop.2020.168126 https://arxiv.org/abs/1910.14507 https://doi.org/10.1016/j.cpc.2017.12.017 https://arxiv.org/abs/1709.01266 https://doi.org/10.37236/5999 https://arxiv.org/abs/1603.01236 https://doi.org/10.1007/978-3-030-03541-9 https://arxiv.org/abs/2008.12310 https://doi.org/10.1016/j.nuclphysb.2020.115096 https://doi.org/10.1016/j.nuclphysb.2020.115096 https://arxiv.org/abs/2005.04265 https://doi.org/10.1103/physrevd.103.116024 https://arxiv.org/abs/2103.16224 https://arxiv.org/abs/2105.05015 https://doi.org/10.1103/PhysRevD.15.1544 https://doi.org/10.1103/PhysRevD.15.1544 https://doi.org/10.1016/S0370-2693(00)00051-4 https://arxiv.org/abs/hep-th/9912093 https://doi.org/10.1016/S0550-3213(01)00071-2 https://arxiv.org/abs/hep-th/0012146 https://doi.org/10.1016/j.physrep.2007.03.003 https://arxiv.org/abs/0707.1596 24 M. Borinsky, G.V. Dunne and M. Meynig [25] Caprini I., Fischer J., Abbas G., Ananthanarayan B., Perturbative expansions in QCD improved by confor- mal mappings of the Borel plane, in Perturbation Theory: Advances in Research and Applications, Nova Science Publishers, Inc., Hauppauge, 2018, 211–254, arXiv:1711.04445. [26] Cardy J.L., High-energy behavior in ϕ3 theory in six-dimensions, Nuclear Phys. B 93 (1975), 525–546. [27] Clavier P.J., Borel–Écalle resummation of a two-point function, Ann. Henri Poincaré 22 (2021), 2103–2136, arXiv:1912.03237. [28] Connes A., Kreimer D., Renormalization in quantum field theory and the Riemann–Hilbert problem, J. High Energy Phys. 1999 (1999), no. 9, 024, 8 pages, arXiv:hep-th/9909126. [29] Connes A., Kreimer D., Renormalization in quantum field theory and the Riemann–Hilbert problem. II. The β-function, diffeomorphisms and the renormalization group, Comm. Math. Phys. 216 (2001), 215–241, arXiv:hep-th/0003188. [30] Cornwall J.M., Morris D.A., Toy models of nonperturbative asymptotic freedom in ϕ3 in six-dimensions, Phys. Rev. D 52 (1995), 6074–6086, arXiv:hep-ph/9506293. [31] Costin O., Exponential asymptotics, transseries, and generalized Borel summation for analytic, nonlin- ear, rank-one systems of ordinary differential equations, Int. Math. Res. Not. 1995 (1995), 377–417, arXiv:math.CA/0608414. [32] Costin O., On Borel summation and Stokes phenomena for rank-1 nonlinear systems of ordinary differential equations, Duke Math. J. 93 (1998), 289–344, arXiv:math.CA/0608408. [33] Costin O., Asymptotics and Borel summability, Chapman & Hall/CRC Monographs and Surveys in Pure and Applied Mathematics, Vol. 141, CRC Press, Boca Raton, FL, 2009. [34] Costin O., Dunne G.V., Physical resurgent extrapolation, Phys. Lett. B 808 (2020), 135627, 7 pages, arXiv:2003.07451. [35] Costin O., Dunne G.V., Uniformization and constructive analytic continuation of Taylor series, arXiv:2009.01962. [36] Courtiel J., Yeats K., Next-tok leading log expansions by chord diagrams, Comm. Math. Phys. 377 (2020), 469–501, arXiv:1906.05139. [37] Damburg R.J., Propin R.K., Graffi S., Grecchi V., Harrell E.M., Č́ıžek J., Paldus J., Silverstone H.J., 1/R expansion for H+ 2 : analyticity, summability, asymptotics, and calculation of exponentially small terms, Phys. Rev. Lett. 52 (1984), 1112–1115. [38] de Alcantara Bonfim O.F., Kirkham J.E., McKane A.J., Critical exponents to order ϵ3 for ϕ3 models of critical phenomena in six ϵ-dimension, J. Phys. A: Math. Gen. 13 (1980), L247–L251, Erratum, J. Phys. A: Math. Gen. 13 (1980), 3785–3785. [39] de Alcantara Bonfim O.F., Kirkham J.E., McKane A.J., Critical exponents for the percolation problem and the Yang–Lee edge singularity, J. Phys. A: Math. Gen. 13 (1981), 2391–2413. [40] Delabaere E., Dillinger H., Pham F., Exact semiclassical expansions for one-dimensional quantum oscillators, J. Math. Phys. 38 (1997), 6126–6184. [41] Delabaere E., Pham F., Resurgent methods in semi-classical asymptotics, Ann. Inst. H. Poincaré Phys. Théor. 71 (1999), 1–94. [42] Delbourgo R., Elliott D., McAnally D.S., Dimensional renormalization in ϕ3 theory: ladders and rainbows, Phys. Rev. D 55 (1997), 5230–5233, arXiv:hep-th/9611150. [43] Dondi N.A., Dunne G.V., Reichert M., Sannino F., Towards the QED beta function and renormalons at 1/N2 f and 1/N3 f , Phys. Rev. D 102 (2020), 035005, 15 pages, arXiv:2003.08397. [44] Dorigoni D., An introduction to resurgence, trans-series and alien calculus, Ann. Physics 409 (2019), 167914, 38 pages, arXiv:1411.3585. [45] Dunne G.V., Ünsal M., Resurgence and trans-series in quantum field theory: the CPN−1 model, J. High Energy Phys. 2012 (2012), no. 11, 170, 85 pages, arXiv:1210.2423. [46] Dunne G.V., Ünsal M., Generating nonperturbative physics from perturbation theory, Phys. Rev. D 89 (2014), 041701, 5 pages, arXiv:1306.4405. [47] Dunne G.V., Ünsal M., New nonperturbative methods in quantum field theory: from large-N orbifold equivalence to bions and resurgence, Ann. Rev. Nuclear Part. Sci. 66 (2016), 245–272, arXiv:1601.03414. [48] Écalle J., Les fonctions résurgentes, Tome I, Les algèbres de fonctions résurgentes, Publications Mathématiques d’Orsay 81, Vol. 5, Université de Paris-Sud, Orsay, 1981. https://arxiv.org/abs/1711.04445 https://doi.org/10.1016/0550-3213(75)90518-0 https://doi.org/10.1007/s00023-021-01057-w https://arxiv.org/abs/1912.03237 https://doi.org/10.1088/1126-6708/1999/09/024 https://doi.org/10.1088/1126-6708/1999/09/024 https://arxiv.org/abs/hep-th/9909126 https://doi.org/10.1007/PL00005547 https://arxiv.org/abs/hep-th/0003188 https://doi.org/10.1103/PhysRevD.52.6074 https://arxiv.org/abs/hep-ph/9506293 https://doi.org/10.1155/S1073792895000286 https://arxiv.org/abs/math.CA/0608414 https://doi.org/10.1215/S0012-7094-98-09311-5 https://arxiv.org/abs/math.CA/0608408 https://doi.org/10.1016/j.physletb.2020.135627 https://arxiv.org/abs/2003.07451 https://arxiv.org/abs/2009.01962 https://doi.org/10.1007/s00220-020-03691-7 https://arxiv.org/abs/1906.05139 https://doi.org/10.1103/PhysRevLett.52.1112 https://doi.org/10.1088/0305-4470/13/7/006 https://doi.org/10.1088/0305-4470/13/12/529 https://doi.org/10.1088/0305-4470/13/12/529 https://doi.org/10.1088/0305-4470/14/9/034 https://doi.org/10.1063/1.532206 https://doi.org/10.1103/PhysRevD.55.5230 https://arxiv.org/abs/hep-th/9611150 https://doi.org/10.1103/physrevd.102.035005 https://arxiv.org/abs/2003.08397 https://doi.org/10.1016/j.aop.2019.167914 https://arxiv.org/abs/1411.3585 https://doi.org/10.1007/JHEP11(2012)170 https://doi.org/10.1007/JHEP11(2012)170 https://arxiv.org/abs/1210.2423 https://doi.org/10.1103/PhysRevD.89.041701 https://arxiv.org/abs/1306.4405 https://doi.org/10.1146/annurev-nucl-102115-044755 https://arxiv.org/abs/1601.03414 Semiclassical Trans-Series from the Perturbative Hopf-Algebraic DSEs 25 [49] Fei L., Giombi S., Klebanov I.R., Critical O(N) models in 6−ϵ dimensions, Phys. Rev. D 90 (2014), 025018, 19 pages, arXiv:1404.1094. [50] Fisher M.E., Yang–Lee edge singularity and ϕ3 field theory, Phys. Rev. Lett. 40 (1978), 1610–1613. [51] Garoufalidis S., Its A., Kapaev A., Mariño M., Asymptotics of the instantons of Painlevé I, Int. Math. Res. Not. 2012 (2012), 561–606, arXiv:1002.3634. [52] Giombi S., Huang R., Klebanov I.R., Pufu S.S., Tarnopolsky G., O(N) model in 4 < d < 6: instantons and complex CFTs, Phys. Rev. D 101 (2020), 045013, 25 pages, arXiv:1910.02462. [53] Gracey J.A., Four loop renormalization of ϕ3 theory in six dimensions, Phys. Rev. D 92 (2015), 025012, 31 pages, arXiv:1506.03357. [54] Gracey J.A., Asymptotic freedom from the two-loop term of the β function in a cubic theory, Phys. Rev. D 101 (2020), 125022, 13 pages, arXiv:2004.14208. [55] Gracey J.A., Ryttov T.A., Shrock R., Renormalization-group behavior of ϕ3 theories in d = 6 dimensions, Phys. Rev. D 102 (2020), 045016, 9 pages, arXiv:2007.12234. [56] Grinstein B., Stone D., Stergiou A., Zhong M., Challenge to the a theorem in six dimensions, Phys. Rev. Lett. 113 (2014), 231602, 5 pages, arXiv:1406.3626. [57] Houghton A., Reeve J.S., Wallace D.J., High order behavior in ϕ3 field theories and the percolation problem, Phys. Rev. B 15 (1977), 2956–2964. [58] Ishikawa K., Morikawa O., Shibata K., Suzuki H., Takaura H., Renormalon structure in compactified space- time, Prog. Theor. Exp. Phys. (2020), 013B01, 14 pages, arXiv:1909.09579. [59] Jentschura U.D., Soff G., Improved conformal mapping of the Borel plane, J. Phys. A: Math. Gen. 34 (2001), 1451–1457, arXiv:hep-ph/0006089. [60] Kompaniets M.V., Panzer E., Minimally subtracted six-loop renormalization of O(n)-symmetric ϕ4 theory and critical exponents, Phys. Rev. D 96 (2017), 036016, 26 pages, arXiv:1705.06483. [61] Kreimer D., On the Hopf algebra structure of perturbative quantum field theories, Adv. Theor. Math. Phys. 2 (1998), 303–334, arXiv:q-alg/9707029. [62] Kreimer D., Yeats K., An étude in non-linear Dyson–Schwinger equations, Nuclear Phys. B Proc. Suppl. 160 (2006), 116–121, arXiv:hep-th/0605096. [63] Kreimer D., Yeats K., Recursion and growth estimates in renormalizable quantum field theory, Comm. Math. Phys. 279 (2008), 401–427, arXiv:hep-th/0612179. [64] Krüger O., Log expansions from combinatorial Dyson–Schwinger equations, Lett. Math. Phys. 110 (2020), 2175–2202, arXiv:1906.06131. [65] Lapedes A., Mottola E., Complex path integrals and finite temperature, Nuclear Phys. B 203 (1982), 58–92. [66] Lipatov L.N., Divergence of the perturbation theory series and the quasiclassical theory, Sov. Phys. JETP 45 (1977), 216–223. [67] Ma E., Asymptotic freedom and a ‘quark’ model in six-dimension, Progr. Theoret. Phys. 54 (1975), 1828– 1832. [68] Macfarlane A.J., Woo G., ϕ3 theory in six dimensions and the renormalization group, Nuclear Phys. B 77 (1974), 91–108, Erratum, Nuclear Phys. B 86 (1975), 548–548. [69] Mahmoud A.A., Yeats K., Connected chord diagrams and the combinatorics of asymptotic expansions, arXiv:2010.06550. [70] Maiezza A., Vasquez J.C., Renormalons in a general quantum field theory, Ann. Physics 394 (2018), 84–97, arXiv:1802.06022. [71] Mariño M., Lectures on non-perturbative effects in large N gauge theories, matrix models and strings, Fortschr. Phys. 62 (2014), 455–540, arXiv:1206.6272. [72] Mariño M., Reis T., Renormalons in integrable field theories, J. High Energy Phys. 2020 (2020), no. 4, 160, 36 pages, arXiv:1909.12134. [73] Mariño M., Reis T., A new renormalon in two dimensions, J. High Energy Phys. 2020 (2020), no. 7, 216, 34 pages, arXiv:1912.06228. [74] Mariño M., Reis T., Resurgence and renormalons in the one-dimensional Hubbard model, arXiv:2006.05131. [75] Marie N., Yeats K., A chord diagram expansion coming from some Dyson–Schwinger equations, Commun. Number Theory Phys. 7 (2013), 251–291, arXiv:1210.5457. https://doi.org/10.1103/PhysRevD.90.025018 https://arxiv.org/abs/1404.1094 https://doi.org/10.1103/PhysRevLett.40.1610 https://doi.org/10.1093/imrn/rnr029 https://doi.org/10.1093/imrn/rnr029 https://arxiv.org/abs/1002.3634 https://doi.org/10.1103/physrevd.101.045013 https://arxiv.org/abs/1910.02462 https://doi.org/10.1103/PhysRevD.92.025012 https://arxiv.org/abs/1506.03357 https://doi.org/10.1103/physrevd.101.125022 https://arxiv.org/abs/2004.14208 https://doi.org/10.1103/physrevd.102.045016 https://arxiv.org/abs/2007.12234 https://doi.org/10.1103/PhysRevLett.113.231602 https://doi.org/10.1103/PhysRevLett.113.231602 https://arxiv.org/abs/1406.3626 https://doi.org/10.1103/PhysRevB.17.2956 https://doi.org/10.1093/ptep/ptz147 https://arxiv.org/abs/1909.09579 https://doi.org/10.1088/0305-4470/34/7/316 https://arxiv.org/abs/hep-ph/0006089 https://doi.org/10.1103/physrevd.96.036016 https://arxiv.org/abs/1705.06483 https://doi.org/10.4310/ATMP.1998.v2.n2.a4 https://arxiv.org/abs/q-alg/9707029 https://doi.org/10.1016/j.nuclphysbps.2006.09.036 https://arxiv.org/abs/hep-th/0605096 https://doi.org/10.1007/s00220-008-0431-7 https://doi.org/10.1007/s00220-008-0431-7 https://arxiv.org/abs/hep-th/0612179 https://doi.org/10.1007/s11005-020-01288-8 https://arxiv.org/abs/1906.06131 https://doi.org/10.1016/0550-3213(82)90477-1 https://doi.org/10.1143/PTP.54.1828 https://doi.org/10.1016/0550-3213(74)90306-X https://doi.org/10.1016/0550-3213(75)90361-2 https://arxiv.org/abs/2010.06550 https://doi.org/10.1016/j.aop.2018.04.027 https://arxiv.org/abs/1802.06022 https://doi.org/10.1002/prop.201400005 https://arxiv.org/abs/1206.6272 https://doi.org/10.1007/jhep04(2020)160 https://arxiv.org/abs/1909.12134 https://doi.org/10.1007/jhep07(2020)216 https://arxiv.org/abs/1912.06228 https://arxiv.org/abs/2006.05131 https://doi.org/10.4310/CNTP.2013.v7.n2.a2 https://doi.org/10.4310/CNTP.2013.v7.n2.a2 https://arxiv.org/abs/1210.5457 26 M. Borinsky, G.V. Dunne and M. Meynig [76] McKane A.J., Vacuum instability in scalar field theories, Ph.D. Thesis, University of Southampton, 1978. [77] McKane A.J., Vacuum instability in scalar field theories, Nuclear Phys. B 152 (1979), 166–188. [78] McKane A.J., Perturbation expansions at large order: results for scalar field theories revisite, J. Phys. A: Math. Theor. 52 (2019), 055401, 23 pages, arXiv:1807.00656. [79] Misumi T., Nitta M., Sakai N., Resurgence in sine-Gordon quantum mechanics: exact agreement be- tween multi-instantons and uniform WKB, J. High Energy Phys. 2015 (2015), no. 9, 157, 41 pages, arXiv:1507.00408. [80] Mitschi C., Sauzin D., Divergent series, summability and resurgence. I. Monodromy and resurgence, Lecture Notes in Math., Vol. 2153, Springer, Cham, 2016. [81] Panzer E., Algorithms for the symbolic integration of hyperlogarithms with applications to Feynman inte- grals, Comput. Phys. Commun. 188 (2015), 148–166, arXiv:1403.3385. [82] Sauzin D., Introduction to 1-summability and resurgence, arXiv:1405.0356. [83] Schnetz O., Graphical functions and single-valued multiple polylogarithms, Commun. Number Theory Phys. 8 (2014), 589–675, arXiv:1302.6445. [84] Schnetz O., The Galois coaction on the electron anomalous magnetic moment, Commun. Number Theory Phys. 12 (2018), 335–354, arXiv:1711.05118. [85] Schnetz O., Numbers and functions in quantum field theory, Phys. Rev. D 97 (2018), 085018, 20 pages, arXiv:1606.08598. [86] Sloane N.J.A., The on-line encyclopedia of integer sequences, http://oeis.org/. [87] Stahl H., The convergence of Padé approximants to functions with branch points, J. Approx. Theory 91 (1997), 139–204. [88] van Baalen G., Kreimer D., Uminsky D., Yeats K., The QED β-function from global solutions to Dyson– Schwinger equations, Ann. Physics 324 (2009), 205–219, arXiv:0805.0826. [89] van Baalen G., Kreimer D., Uminsky D., Yeats K., The QCD β-function from global solutions to Dyson– Schwinger equations, Ann. Physics 325 (2010), 300–324, arXiv:0906.1754. [90] Yeats K., Growth estimates for Dyson–Schwinger equations, Ph.D. Thesis, University of Waterloo, 2008, arXiv:0810.2249. [91] Yeats K., A combinatorial perspective on quantum field theory, SpringerBriefs in Mathematical Physics, Vol. 15, Springer, Cham, 2017. [92] Zinn-Justin J., Expansion around instantons in quantum mechanics, J. Math. Phys. 22 (1981), 511–520. [93] Zinn-Justin J., Quantum field theory and critical phenomena, International Series of Monographs on Physics, Vol. 77, The Clarendon Press, Oxford University Press, New York, 1989. [94] Zinn-Justin J., Jentschura U.D., Multi-instantons and exact results. I. Conjectures, WKB expansions, and instanton interactions, Ann. Physics 313 (2004), 197–267, arXiv:quant-ph/0501136. [95] Zinn-Justin J., Jentschura U.D., Multi-instantons and exact results. II. Specific cases, higher-order effects, and numerical calculations, Ann. Physics 313 (2004), 269–325. https://doi.org/10.1016/0550-3213(79)90086-5 https://doi.org/10.1088/1751-8121/aaf768 https://doi.org/10.1088/1751-8121/aaf768 https://arxiv.org/abs/1807.00656 https://doi.org/10.1007/JHEP09(2015)157 https://arxiv.org/abs/1507.00408 https://doi.org/10.1007/978-3-319-28736-2 https://doi.org/10.1016/j.cpc.2014.10.019 https://arxiv.org/abs/1403.3385 https://arxiv.org/abs/1405.0356 https://doi.org/10.4310/CNTP.2014.v8.n4.a1 https://arxiv.org/abs/1302.6445 https://doi.org/10.4310/cntp.2018.v12.n2.a4 https://doi.org/10.4310/cntp.2018.v12.n2.a4 https://arxiv.org/abs/1711.05118 https://doi.org/10.1103/physrevd.97.085018 https://arxiv.org/abs/1606.08598 http://oeis.org/ https://doi.org/10.1006/jath.1997.3141 https://doi.org/10.1016/j.aop.2008.05.007 https://arxiv.org/abs/0805.0826 https://doi.org/10.1016/j.aop.2009.10.011 https://arxiv.org/abs/0906.1754 https://arxiv.org/abs/0810.2249 https://doi.org/10.1007/978-3-319-47551-6 https://doi.org/10.1063/1.524919 https://doi.org/10.1016/j.aop.2004.04.004 https://arxiv.org/abs/quant-ph/0501136 https://doi.org/10.1016/j.aop.2004.04.003 1 Introduction 2 Perturbative Hopf-algebraic analysis of massless phi3 theory in 6 dimensions 3 Asymptotics of the perturbative solution of the Hopf-algebraic Dyson–Schwinger equation 4 Borel analysis of the perturbative expansion of the anomalous dimension 5 Trans-series analysis of the Hopf-algebraic Dyson–Schwinger equation 5.1 Identifying the "seed" exponential terms: the linearized equation 5.2 Beyond the linearized equation: resonant trans-series and logarithms 5.3 Logarithmic behavior in large-order growth 6 Conclusions A Modified Richardson method for probing logarithmic large-order growth References
id	nasplib_isofts_kiev_ua-123456789-211440
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn	1815-0659
language	English
last_indexed	2026-03-20T15:48:17Z
publishDate	2021
publisher	Інститут математики НАН України
record_format	dspace
spelling	Borinsky, Michael Dunne, Gerald V. Meynig, Max 2026-01-02T08:34:15Z 2021 Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions. Michael Borinsky, Gerald V. Dunne and Max Meynig. SIGMA 17 (2021), 087, 26 pages 1815-0659 2020 Mathematics Subject Classification: 81T15; 81Q15; 34E10 arXiv:2104.00593 https://nasplib.isofts.kiev.ua/handle/123456789/211440 https://doi.org/10.3842/SIGMA.2021.087 We analyze the asymptotically free massless scalar ³ quantum field theory in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer-Connes Hopf-algebraic Dyson-Schwinger equations for the anomalous dimension. This scalar conformal field theory is asymptotically free and has a real Lipatov instanton. In the Hopf-algebraic approach, we find a trans-series having an intricate Borel singularity structure, with three distinct but resonant non-perturbative terms, each repeated in an infinite series. These expansions are in terms of the renormalized coupling. The resonant structure leads to powers of logarithmic terms at higher levels of the trans-series, analogous to logarithmic terms arising from interactions between instantons and anti-instantons, but arising from a purely perturbative formalism rather than from a semi-classical analysis. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Officeof High Energy Physics under Award Number DE-SC0010339 (GD, MM) and by the NWO Vidi grant 680-47-551 “Decoding Singularities of Feynman graphs” (MB). This work was begun during visits by the first two authors to Humboldt University in 2018 and 2019, and at the Les Houches Summer School in 2018, and we thank these institutions for their hospitality. We are grateful to Marc Bellon, David Broadhurst, Ovidiu Costin, John Gracey, Dirk Kreimer, Enrico Russo, and Karen Yeats for discussions. We also want to thank David Broadhurst for helping with the estimation of the constant in equation (5.26) and for pointing out some typos in a previous version. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions Article published earlier
spellingShingle	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions Borinsky, Michael Dunne, Gerald V. Meynig, Max
title	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions
title_full	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions
title_fullStr	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions
title_full_unstemmed	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions
title_short	Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions
title_sort	semiclassical trans-series from the perturbative hopf-algebraic dyson-schwinger equations: ³ qft in 6 dimensions
url	https://nasplib.isofts.kiev.ua/handle/123456789/211440
work_keys_str_mv	AT borinskymichael semiclassicaltransseriesfromtheperturbativehopfalgebraicdysonschwingerequations3qftin6dimensions AT dunnegeraldv semiclassicaltransseriesfromtheperturbativehopfalgebraicdysonschwingerequations3qftin6dimensions AT meynigmax semiclassicaltransseriesfromtheperturbativehopfalgebraicdysonschwingerequations3qftin6dimensions

Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: ³ QFT in 6 Dimensions

Репозитарії

Схожі ресурси