$\lambda \phi^4$ Theory II: The Broken Phase Beyond NNNN(NNNN)LO

We extend the study of the two-dimensional euclidean $\phi^4$ theory initiated in ref. [1] to the $\mathbb Z_2$ broken phase. In particular, we compute in perturbation theory up to N$^4$LO in the quartic coupling the vacuum energy, the vacuum expectation value of $\phi$ and the mass gap of the theory. We determine the large order behavior of the perturbative series by finding the leading order finite action complex instanton configuration in the $\mathbb Z_2$ broken phase. Using an appropriate conformal mapping, we then Borel resum the perturbative series. Interestingly enough, the truncated perturbative series for the vacuum energy and the vacuum expectation value of the field is reliable up to the critical coupling where a second order phase transition occurs, and breaks down around the transition for the mass gap. We compute the vacuum energy using also an alternative perturbative series, dubbed exact perturbation theory, that allows us to effectively reach N$^8$LO in the quartic coupling. In this way we can access the strong coupling region of the $\mathbb Z_2$ broken phase and test Chang duality by comparing the vacuum energies computed in three different descriptions of the same physical system. This result can also be considered as a confirmation of the Borel summability of the theory. Our results are in very good agreement (and with comparable or better precision) with those obtained by Hamiltonian truncation methods. We also discuss some subtleties related to the physical interpretation of the mass gap and provide evidence that the kink mass can be obtained by analytic continuation from the unbroken to the broken phase.


!
Unbroken Broken g c 0 0 1 1 g < l a t e x i t s h a 1 _ b a s e 6 4 = " E i J K / E 3 t r R L z 9 T A m n U y a j Z 6 C u P w = " > A A A B 6 H i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 l E q M e i F 4 8 t 2 A 9 o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g o b m 1 v b O 8 X d 0 t 7 + w e F R + f i k r e N U M W y x W M S q G 1 C N g k t s G W 4 E d h O F N A o E d o L J 3 d z v P K H S P J Y P Z p q g H 9 G R 5 C F n 1 F i p O R q U K 2 7 V X Y C s E y 8 n F c j R G J S / + s O Y p R F K w w T V u u e 5 i f E z q g x n A m e l f q o x o W x C R 9 i z V N I I t Z 8 t D p 2 R C 6 s M S R g r W 9 K Q h f p 7 I q O R 1 t M o s J 0 R N W O 9 6 s 3 F / 7 x e a s I b P + M y S Q 1 K t l w U p o K Y m M y / J k O u k B k x t Y Q y x e 2 t h I 2 p o s z Y b E o 2 B G / 1 5 X X S v q p 6 b t V r X l f q t 3 k c R T i D c 7 g E D 2 p Q h 3 t o Q A s Y I D z D K 7 w 5 j 8 6 L 8 + 5 8 L F s L T j 5 z C n / g f P 4 A z H W M 6 w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " E i J K / E 3 t r R L z 9 T A m n U y a j Z 6 C u P w = " > A A A B 6 H i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 l E q M e i F 4 8 t 2 A 9 o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g o b m 1 v b O 8 X d 0 t 7 + w e F R + f i k r e N U M W y x W M S q G 1 C N g k t s G W 4 E d h O F N A o E d o L J 3 d z v P K H S P J Y P Z p q g H 9 G R 5 C F n 1 F i p O R q U K 2 7 V X Y C s E y 8 n F c j R G J S / + s O Y p R F K w w T V u u e 5 i f E z q g x n A m e l f q o x o W x C R 9 i z V N I I t Z 8 t D p 2 R C 6 s M S R g r W 9 K Q h f p 7 I q O R 1 t M o s J 0 R N W O 9 6 s 3 F / 7 x e a s I b P + M y S Q 1 K t l w U p o K Y m M y / J k O u k B k x t Y Q y x e 2 t h I 2 p o s z Y b E o 2 B G / 1 5 X X S v q p 6 b t V r X l f q t 3 k c R T i D c 7 g E D 2 p Q h 3 t o Q A s Y I D z D K 7 w 5 j 8 6 L 8 + 5 8 L F s L T j 5 z C n / g f P 4 A z H W M 6 w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " E i J K / E 3 t r R L z 9 T A m n U y a j Z 6 C u P w = " > A A A B 6 H i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 l E q M e i F 4 8 t 2 A 9 o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g o b m 1 v b O 8 X d 0 t 7 + w e F R + f i k r e N U M W y x W M S q G 1 C N g k t s G W 4 E d h O F N A o E d o L J 3 d z v P K H S P J Y P Z p q g H 9 G R 5 C F n 1 F i p O R q U K 2 7 V X Y C s E y 8 n F c j R G J S / + s O Y p R F K w w T V u u e 5 i f E z q g x n A m e l f q o x o W x C R 9 i z V N I I t Z 8 t D p 2 R C 6 s M S R g r W 9 K Q h f p 7 I q O R 1 t M o s J 0 R N W O 9 6 s 3 F / 7 x e a s I b P + M y S Q 1 K t l w U p o K Y m M y / J k O u k B k x t Y Q y x e 2 t h I 2 p o s z Y b E o 2 B G / 1 5 X X S v q p 6 b t V r X l f q t 3 k c R T i D c 7 g E D 2 p Q h 3 t o Q A s Y I D z D K 7 w 5 j 8 6 L 8 + 5 8 L F s L T j 5 z C n / g f P 4 A z H W M 6 w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " E i J K / E 3 t r R L z 9 T A m n U y a j Z 6 C u P w = " > A A A B 6 H i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 l E q M e i F 4 8 t 2 A 9 o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g o b m 1 v b O 8 X d 0 t 7 + w e F R + f i k r e N U M W y x W M S q G 1 C N g k t s G W 4 E d h O F N A o E d o L J 3 d z v P K H S P J Y P Z p q g H 9 G R 5 C F n 1 F i p O R q U K 2 7 V X Y C s E y 8 n F c j R G J S / + s O Y p R F K w w T V u u e 5 i f E z q g x n A m e l f q o x o W x C R 9 i z V N I I t Z 8 t D p 2 R C 6 s M S R g r W 9 K Q h f p 7 I q O R 1 t M o s J 0 R N W O 9 6 s 3 F / 7 x e a s I b P + M y S Q 1 K t l w U p o K Y m M y / J k O u k B k x t Y Q y x e 2 t h I 2 p o s z Y b E o 2 B G / 1 5 X X S v q p 6 b t V r X l f q t 3 k c R T i D c + j E + M E u J u L E 0 J j a f q 7 4 m M R E q N o s B 0 R k T 3 1 b w 3 E f / z 2 q n u X n o Z E 0 m q Q d D Z o m 7 K s Y 7 x J A o c M g l U 8 5 E h h E p m b s W 0 T y S h 2 g R W N C G 4 8 y 8 v k s Z Z x X U q 7 u 1 5 q X q V x 1 F A R + g Y l Z G L L l A V 3 a A a q i O K H t E z e k V v 1 p P 1 Y r 1 b H 7 P W J S u f O U B / Y H 3 + A A y p l P E = < / l a t e x i t > < l a t e x i t s h a + j E + M E u J u L E 0 J j a f q 7 4 m M R E q N o s B 0 R k T 3 1 b w 3 E f / z 2 q n u X n o Z E 0 m q Q d D Z o m 7 K s Y 7 x J A o c M g l U 8 5 E h h E p m b s W 0 T y S h 2 g R W N C G 4 8 y 8 v k s Z Z x X U q 7 u 1 5 q X q V x 1 F A R + g Y l Z G L L l A V 3 a A a q i O K H t E z e k V v 1 p P 1 Y r 1 b H 7 P W J S u f O U B / Y H 3 + A A y p l P E = < / l a t e x i t > < l a t e x i t s h a  a connection has been found recently between the Lefschetz thimble decomposition of path integrals and Borel summability of perturbative series [24,25], that allowed us to show the Borel summability of a broad class of Euclidean 2d and 3d scalar field theories [1].These include the 2d φ 4 theory with both m 2 > 0, known already to be Borel resummable for parametrically small couplings [26], and m 2 < 0.
The aim of this paper is to extend to the Z 2 broken phase the study of the 2d Euclidean φ 4 theory of ref. [1], where the unbroken phase was analyzed.In this way we will also provide numerical evidence of the Borel summability of the theory and the first applications of the methods introduced in refs.[24,25] in QFT.The Euclidean Lagrangian density (modulo counterterms) reads with m 2 > 0. We denote the parameters in the broken phase with a tilde to distinguish them from the ones appearing in the Lagrangian density used in ref. [1] to analyse the unbroken phase (m 2 > 0).The effective expansion parameter in the Lagrangian (1.1) is the dimensionless coupling We start in section 2 by reviewing Chang duality and the phase structure of the theory, as expected by the duality.This is summarized in fig. 1.We see that starting from a given value of the coupling, g B in the unbroken phase, the φ 4 theory admits three equivalent descriptions: one in terms of a theory with tree-level mass term m 2 > 0 and coupling g = λ/m 2 (black line) and two in terms of a theory with tree-level mass term m 2 = − m 2 /2 < 0 and coupling g (red lines).We call weakly and strongly coupled branches the regions 0 ≤ g ≤ g B and g ≥ g B , respectively.The points g c , g (w) c and g (s) c in fig. 1 denote the critical couplings in the three descriptions where the theory has a second order phase transition and flows to the 2d Ising model.In the Z 2 broken phase the only description that can be explored by resumming the perturbative series and without encountering phase transitions is the part of the weakly coupled branch with 0 ≤ g ≤ g (w) c .This is the region we will mostly focus on, although we will also explore other sectors of the phase diagram.In ref. [1] we instead focused on the region 0 ≤ g ≤ g c in the black line of fig. 1.
In section 3 we set the stage for the methods we will use in the paper to Borel resum the perturbative series.In subsection 3.1 we look for finite action solutions to the complexified Euclidean equations of motion (complex instanton configurations).As well-known [27], the configurations with the smallest action determine the leading large-order behavior of the perturbative expansion of a given observable and the leading singularity of its Borel transform function.This is useful information, that we exploit to set up a suitable conformal mapping for the numerical Borel resummation of the perturbative series.In subsection 3.2 we consider a modified perturbative series, where the trilinear coupling defined below in eq.(1.4) gets replaced by λ 3 → √ 2λ m 2 / √ λ 0 .It has been shown in refs.[24,25], in the context of quantum mechanical models, that modified perturbative expansion of this kind, where one expands in g with g 0 ≡ λ 0 / m 2 held fixed, and sets back g 0 = g after the resummation has been performed, leads to a significant improvement of the Borel resummation of perturbative series at strong coupling.Expansions of this kind were dubbed in refs.[24,25] Exact Perturbation Theory (EPT). 1 EPT will allow us to compute the vacuum energy in the Z 2 broken phase at strong coupling, namely in the region g ≥ g 0 in fig. 1 (strongly coupled branch).Note that for g = g 0 the Z 2 symmetry is explicitly broken, so EPT is equivalent to adding an explicit symmetry breaking term which is eventually set to zero.
In section 4 we compute the perturbative series expansion up to order g 4 of the 0-point, 1-point and 2-point Schwinger functions.This is the maximal order that can be reached by computing Feynman diagrams with up to eight interaction vertices (the maximum we could reach), because of the presence of the trilinear interaction proportional to We use a renormalization scheme which is equivalent to normal ordering, but perform our computations in an intermediate auxiliary scheme which allows us to efficiently treat the corrections to the 1-point tadpole that are generated order by order in perturbation theory.The main results of this section are the perturbative expressions for the vacuum energy Λ, the 1-point function φ and the physical mass M 2 in eqs.(4.7), (4.9) and (4.12), respectively.In addition to that, we is the one found in this paper from eq. (1.5), g Chang c is g expressed in terms of the unbroken variables using Chang relation (2.9) and g c is the value found in the unbroken phase in ref. [1].provide in eq.(4.8) the perturbative expression of Λ up to g 8 in EPT.
We report in section 5 the final results of our investigation.We start by discussing the main differences in the numerical resummation procedures with respect to ref. [1].In particular, we point out that the shortness of the ordinary perturbative series, together with the slow convergence of the series after the conformal mapping,2 leads to results which are much less accurate than those found by Borel resumming in the unbroken case.On the other hand, we find that in the entire regime 0 ≤ g ≤ g (w) c , the O( g4 ) perturbative series is reliable for Λ and φ while the series for M 2 breaks down slightly before g (w) c .For all the three observables we find that the central values of the Borel resummed results are in very good agreement with the perturbative result, see figs.5, 8 and 9.We determine the critical coupling g ( We report in table 1 the value of g (w) c .For convenience, we also report its value g Chang c in terms of the unbroken variables, using Chang duality, and compare it with the value g c found in ref. [1].Chang duality predicts g Chang c = g c , compatible with our results.As discussed in ref. [1], the Borel resummed Schwinger functions are expected to coincide with the exact ones in a given phase of the theory.At g = g (w) c the two vacua of the broken phase collide in the unique Z 2 invariant vacuum and for g the Z 2 symmetry is restored.Consequently our perturbative results for φ and M 2 for g > g (w) c do not have an obvious direct physical interpretation. 3n the other hand, Λ should be continuous along the transition and its value past g (w) c should still be identified with the vacuum energy in the symmetric phase and then after g c with the vacuum energy in the strongly coupled branch of the broken phase.This comment applies also for the vacuum energy Λ for g > g c computed starting from the unbroken phase with m 2 > 0. We then compute Λ using the O( g 8 ) series associated to EPT, which allows us to explore the strongly coupled branch of the theory with a better accuracy than the ordinary O( g 4 ) series.We determine Λ in the broken phase for a certain range of g in all three descriptions, see fig. 7. We consider the agreement of the results as a convincing numerical check of Chang duality, 4 of the Borel summability of the theory and of the usefulness of EPT in QFT.
The vacuum energy and the mass gap of the 2d φ 4 theory in the broken phase have also been computed using Hamiltonian truncation methods in refs.[10,11].We compare our findings with those of the above works in section 6, finding very good agreement.In particular, we confirm that both the results for Λ of ref. [10] and the ones for M 2 in ref. [11] are within the regime of validity of perturbation theory!The φ 4 theory considered in refs.[10,11] is not however in a genuinely broken phase.This makes a comparison with refs.[10,11] non-trivial in a range of the coupling where the elementary particle is supposed to decay in pair of kink and anti-kink, since single kink states decouple in a theory where the Z 2 symmetry is spontaneously broken.In subsection 6.1 we provide numerical evidence that the kink sector of the theory can be accessed starting from the unbroken phase considered in ref. [1].More precisely, we show that the value of the mass gap |M | computed in the unbroken phase for g ≥ g c is in agreement with the mass of the kink state computed in ref. [11], see fig.11.
We conclude in section 7.In appendix A the coefficients for the series expansion of the 0, 1 and 2-point function are reported.The coefficients have been determined for independent cubic and quartic couplings λ 3 and λ and could be used also for theories with explicit breaking of the Z 2 symmetry φ → −φ.We would finally like to emphasize that while several works along more than forty years have analyzed the 2d φ 4 theory in the unbroken case by means of various resummation procedures [19-23, 28, 29, 1], as far as we know this is the first paper addressing the Z 2 broken phase.

Chang Duality
In the φ 4 theory with m 2 > 0, aside from the normalization of the free theory path integral, there are only two superficially divergent one-particle irreducible (1PI) diagrams: the two loop number eight graph occurring in the 0-point function and the one-loop tadpole of the 2-point function.Correspondingly, the counterterms δΛ and δm 2 contain up to O(λ) terms to all orders in perturbation theory, with no need of a wave function and coupling constant counterterms.Such simple renormalization property are at the base of one of the simplest strong-weak dualities in QFT: Chang duality [4].The original derivation worked with normal ordering prescriptions (recently nicely reviewed in ref. [10]), but the same analysis can be repeated in other regularizations.We choose here dimensional regularization (DR).This formulation allows for a straightforward generalization of the duality in the 3d λφ 4 theory [30], where normal ordering is no longer enough to cancel all divergencies.The theory in d dimensions defined by the (Euclidean) Lagrangian where = 2 − d and µ is the usual RG sliding scale, has a dual description in terms of the theory L with negative squared-mass term which is derived below using dimensional regularization and a modified Minimal Subtraction (MS) renormalization scheme.The counterterms for L are fixed by requiring that for µ = m the 2-point tadpole diagram at one loop is exactly canceled and that the vacuum energy vanishes up to O(λ 2 ). 5 We get where we defined γ ≡ γ E − log(4π), with γ E ≈ 0.577 . . . the Euler-Mascheroni constant.Note that in the above expression the 1/ term of δΛ depends on the squared mass defined at the scale µ at which we are evaluating the theory (whose explicit expression is obtained below).This is needed for the counterterm to cancel the divergent part at any scale.The finite contribution to δΛ depends instead on the arbitrary renormalization point chosen, i.e. m 2 , and hence the squared mass terms that appear in the O( 0 ) term are defined as m 2 (m) ≡ m 2 .To all orders in perturbation theory the β functions for the mass, the coupling and the vacuum energy read and hence Renormalizing with a normal ordering mass µ is equivalent to use m 2 (µ) in the Lagrangian (2.1).
The counterterms for L are chosen in the same modified MS renormalization scheme as before, i.e. by requiring that for µ = m the 2-point tadpole diagram at one loop and the divergent 0-point terms are exactly canceled.We obtain6 Using eq.(2.5) and setting µ = m, we find that the theory (2.1) is equivalent to (2.2) provided the following identities hold (2.7) We have then established a relation between two theories in different phases: the theory in the broken phase with squared mass term − m 2 /4 (giving rise to a state of mass squared m 2 when expanded around the tree-level VEV) in the modified MS scheme with µ = m is identical to the theory in the unbroken phase in the same analogous modified MS scheme with µ = m.This is Chang duality [4].The solutions to the Chang equation (2.7) have been already reviewed in ref. [10], but for completeness we will again briefly summarize them here.Defining the dimensionless coupling constants we can rewrite eq.(2.7) as where At fixed g, we look for solutions in g of eq.(2.9).Since f ( g) > 0 ∀ g, no solutions can evidently exist for sufficiently small g, where f (g) < 0. The minimum of f ( g) occurs at g B = π/6 and a solution exists for where W is the Lambert function (also known as product logarithm or omega function).For g ≥ g B there are two solution branches g w,s (g).We label the two branches as strong (s) and weak (w) branches according to their behavior: (2.12) The existence of the weak branch allows to prove the existence of a phase transition in the φ 4 theory in the classically unbroken phase with m 2 > 0 at sufficiently strong coupling g.Indeed, for parametrically small g the theory is well described by its classical potential and is in the unbroken phase, while at parametrically large couplings the duality implies it is in the broken phase, since it can be described by a weakly coupled theory with m 2 = − m 2 /2 < 0. By continuity there should exist a critical coupling g c where the phase transition occurs.The value of the critical coupling in the normal ordering scheme has been computed by different methods [1, 5-8, 11, 14, 15, 23] to be g c ≈ 2.76.Correspondingly, we can predict the value of the two critical couplings: The phase structure of the theory, as predicted from Chang duality, is then the following.Starting from a perturbative description in the unbroken phase, m 2 > 0, the theory develops a (secondorder) phase transition at g = g c and above g c remains in the Z 2 broken phase.Starting from a perturbative description in the broken phase, m 2 = − m 2 /2 < 0, we first encounter a phase transition at g = g c , the Z 2 symmetry is restored for g c , and at g = g (s) c we have another phase transition.For g > g (s) c the theory remains in the broken phase.Chang duality predicts that the three regimes c , and g ≤ g ≤ ∞ are different descriptions of the same physical theory in the broken phase.Similarly, the three regimes g B ≤ g ≤ g c , g c are different descriptions of the same physical theory in the unbroken phase.In particular the three critical points represent the very same transition in different descriptions.The region 0 ≤ g ≤ g B admits a single description in terms of the unbroken theory with m 2 > 0. See fig. 1 for a summary.In this paper, we will mostly focus on the weak branch regime 0 ≤ g ≤ g (w) c .However, by using different techniques, we will be able to compute the vacuum energy in the broken phase in all the three descriptions.

Borel Summability in the Broken Phase
As discussed in ref. [1] the φ 4 theory is Borel resummable to the exact result also in the broken phase, when the perturbative expansion is performed around either one of the two degenerate vacua at infinite volume.Indeed, thanks to Derrick's theorem [31], the non-trivial real saddles connecting the two vacua have infinite action at infinite volume and do not obstruct the Borel transformation procedure.Besides, the vacuum selection, performed by introducing a fictitious Z 2 breaking parameter removed only after taking the infinite volume limit, decouples any contribution from the Borel resummable perturbative expansion around the 'other' vacuum.The perturbative expansion around one of the two Z 2 -breaking vacua is therefore Borel resummable to the exact result of the symmetry broken theory, similarly to the unbroken case (and with the same caveats associated to phase transitions and manipulations of the N -point functions).
Borel summability implies that a generic observable7 F ( g) function of the coupling g (with a divergent series expansion ∞ n=0 F n g n ) can be recovered by performing the integral where the Borel function B(t) is the analytic continuation of the function defined by the convergent series n F n t n /n!.While Borel summability corresponds to B(t) being regular over the positive real axis, singularities (in general an infinite number of them) are present in the complex t plane.Their position corresponds to the value of the action on all possible (complex) solutions of the classical equations of motion.The saddles with the minimum absolute value for the action determine the radius of convergence of the perturbative series of B(t), thus the leading growth of the perturbative coefficients: if we call t ± 1 the position of such saddles (real action complex saddles come in pair, t − i = t + i * ) the leading asymptotic growth of the coefficients is of the form F n ∝ a n Γ(n + b + 1) + h.c.(with a ≡ 1/t + 1 and b a constant).Because of its finite radius of convergence, the Borel function can be well approximated by a truncated perturbative expansion only in the interval 0 < t < |t ± 1 |, this determines an intrinsic error in the reconstruction of F ( g) from eq. (3.1), which at small coupling is O(e −|t ± 1 |/ g ), reproducing the accuracy limit of the original divergent asymptotic series for F ( g).In order to really improve over the truncated series, B(t) must be analytically continued beyond its radius of convergence |t ± 1 |.Two commonly used techniques to achieve this are the method of Padè approximants and the conformal mapping.To be effective, the first method in general requires the knowledge of a large number of coefficients, unfortunately this is not our case so the improvement obtained with this method is limited.The second technique exploits the knowledge of the position of the singularities in the complex t plane to perform a clever change of variable t(u) after which the integration in eq.(3.1) is mapped over the interval u ∈ [0, 1) and all the singularities of B on the unit circle |u| = 1 of the complex u-plane: The new series expansion of B(u) ≡ B[t(u)] in powers of u is therefore convergent over the entire range of integration-the original divergent truncated series has been transformed into a convergent one!Of course the method requires the knowledge of all singularities of B, i.e. all the finite action complex saddle points of the classical action.Except in very special cases such information is not available.However, even the mapping of the sole leading singularity t ± 1 to the unit circle in the u-plane represents a big improvement: it effectively enlarges the radius of convergence of the original Borel function to the next-to-leading singularities t ± 2 improving the accuracy of the series to O(e −|t ± 2 |/ g ).We discuss this technique more in detail in the next subsection, while in section 3.2 we will describe a different approach based on EPT introduced in refs.[24,25].

Weak Coupling: Conformal Mapping of Complex Saddles
In order to implement the conformal mapping method we must find the saddle points of the classical action (at least the leading ones) and the required mapping t = t(u).Saddles with finite action configurations must have the field flowing to the minimum of the potential at infinity in all Euclidean directions.We concentrate on SO(2) symmetric configurations which are expected to be the dominant ones.In polar coordinates the problem reduces to a system of 2nd order For the unbroken phase the problem could be simplified by focusing on purely imaginary solutions, since in such case the system collapses to a single differential equation corresponding to a simple shooting problem (see ref. [32]).The value of the position of the singularity in this case is real and negative (t u. 1 = −1.4626...).In the broken phase, non-trivial saddle points necessarily require both the real and the imaginary parts of the field to vary.With some work the complex trajectories for the field and the corresponding values of the action can be obtained by numerical integration, following the complex solutions from the unbroken case for increasing values of the φ 3 deformation that convert the unbroken potential to the double-well one.After rescaling the field as φ = (8 g) −1/2 ϕ and the coordinate as r = √ 2ρ/ m , the spherically symmetric solutions to the equation of motion satisfy the following differential equation The trajectories for the first three leading saddle points (among the SO(2) invariant ones) are shown in fig. 2. They always come in complex conjugate pairs ϕ ± i and are quite non-trivial.Interestingly enough, the complex trajectory of the leading saddle exactly matches a circular arc in the complex ϕ-plane.The corresponding values for the action computed numerically lead to the following values for the position of the leading singularities in the complex t plane: where the numerical error should be smaller than the last digit reported, where As mentioned before the value of the leading singularities determines the leading growth of the perturbative coefficients.In particular |t ± 1 | = 1.6181... is the radius of convergence of the Borel transform and roughly measures for what value of the coupling the theory turns completely non-perturbative.The fact that the second saddle is far at |t ± 2 | = 8.89... means that, for our purposes, it can be neglected in the conformal mapping.On the other hand π/arg(t + 1 ) 3.8 measures the half-period of oscillation of the sign of the coefficients, in particular the vicinity of the saddle to the real axis decreases the oscillation of the perturbative coefficients affecting its convergence. 8We do not know the position (if present) of SO(2) non-invariant complex saddles, however they are expected to be subleading with respect to the leading SO(2) invariant one.
Having identified the position of the leading singularities we now turn to the identification of the suitable change of variable t = t(u) for the resummation.This should correspond to the (conformal) map which sends the positive semi axis t ∈ [0, ∞) to the segment u ∈ [0, 1), all singularities t ± i to the unit circle |u| = 1 and be regular around the origin t = u = 0. Schwarz-Christoffel transformations, which can map (degenerate) polygons to the unit disc, exactly achieve this.For the simple case of the mapping of a single couple of complex singularities located at the complex conjugate points t = |t ± 1 |e ±iπα 1 (with α 1 ∈ (0, 1]) the mapping takes the simple form Notice that for a real negative singularity (α 1 = 1) the usual conformal mapping, used also in the unbroken phase (see e.g.ref. [1]), is recovered.For generic complex instantons eq.(3.6) maps the singularities and the associated rays on the unit circle |u| = 1 as in fig. 3.
The factor e −t(u)/ g in eq.(3.2) exponentially suppresses the integrand as u → 1.In the absence of any other singularity besides those mapped on the unit circle, B(u) has unit radius of convergence, and its truncated power series in u can be used to approximate B(u) with arbitrary precision on any point of the integration interval.Hence we can legitimately exchange the integration sign in eq.(3.2) with the summation sign of the power series defining B(u).In this case the conformal mapping has turned the original asymptotic series in g into a new convergent series for every finite value of g.This is not true in the presence of other singularities (as in our case) since the latter reduce the radius of convergence of B(u and the above exchange of integration and summation signs is no longer legitimate.The conformally mapped perturbative expansion can therefore approximate the integrand of eq.(3.2) only up to |u ± 2 |, which at small coupling corresponds to an irreducible error of O(e −|t 2 |/ g ).In this case the original series is mapped into another asymptotic one but with a smaller degree of divergence.Similarly to the hyperasymptotics techniques of refs.[34,35] the knowledge of nearby complex instantons can be used to improve the degree of convergence of the original series.In our specific case the known subleading singularities t ± i>1 are far enough (see fig. 4) and given the number of coefficients available and the range of couplings probed, the map (3.6) is not yet limited by the presence of these subleading instantons.
While in principle any singularity can be moved away by an appropriate choice of t(u), the efficiency of the mapping critically depends on the original position of the singularities, in particular on the angular distance from the positive real axis.For singularities on the negative real axis (α = 1), the region of points at t > |t ±  1 | (which carry the non-trivial information required to improve beyond the original truncated expansion) is mapped into the region u > 3 − 2 √ 2 0.17, well inside the radius of convergence of B(u).As the singularities move closer to the positive real axis in the t-plane, such region is pushed more and more towards the unit circle in the u-plane, in particular for α 1 → 0, the region with When the pair of complex singularities pinch the positive real axis the region t > |t 1 | becomes inaccessible and the series non-Borel summable.As a result the closer the singularities are to the real axis the less performant the conformal mapping is, and the larger is the number of coefficients requested to reach a certain accuracy.As we will discuss in section 5 this is the main limiting factor of this method in our computation.
Note that the conformal mapping (3.6) would have mapped the original asymptotic series into a convergent one if the other saddles were aligned among them (including the origin) along the cut t-plane depicted in fig. 4.This is the situation expected in the unbroken Z 2 -symmetric phase, where the known singularities are aligned over the real negative t axis and are mapped at the boundary of the unit u-disc.As we have seen, this is not the case in the spontaneously broken phase, where the saddles are not aligned.

Strong Coupling: Exact Perturbation Theory
In the context of quantum mechanics it has recently been shown that it is possible to define modified perturbative expansions that are Borel resummable in theories where ordinary perturbation theory is not [24,25].For this reason such modified perturbative series were denoted Exact Perturbation Theory (EPT) in refs.[24,25].EPT can also be useful when the ordinary expansion is Borel resummable to start with, by improving the Borel resummation of perturbative series.Namely, it can improve the accuracy of results obtained by the approximate knowledge of the Borel function with a finite number of perturbative terms, especially at strong coupling.We now show that EPT can similarly be applied in QFT.We will focus in what follows to the 2d φ 4 theory, though most considerations apply also in more general settings.

Consider a n-point Schwinger function
where N is an irrelevant constant factor, φ is the quantum fluctuation around the classical minimum φ cl = + m/ √ 8λ and L is the Lagrangian (2.2) expanded around φ cl : where For simplicity we have omitted to write in eq.(3.8) the counterterms δ m 2 and δ Λ in eqs.(2.6).
The first is higher order in λ, while the last is field-independent, so they can be neglected when establishing the classical finite action field configurations.It is straightforward to see that the expansion in λ is equivalent to a loopwise expansion in and that all the terms appearing in the Lagrangian (3.8) are of the same order in .Consider now the following n-point Schwinger function: where ) For simplicity, as before we omit to write the counterterms necessary to make the theory UV finite.The Lagrangian (3.11) is identical to that in eq.(3.8) except for the cubic coupling.At fixed λ 0 , we have effectively turned the classical cubic term into a quantum one.Hence the classical finite action configurations of eq.(3.11) coincide with those of a theory in the Z 2 unbroken phase with mass term m 2 and quartic coupling λ, although of course the perturbative expansions in the two theories are different because of the cubic term.At strong coupling the Borel resummability of the modified perturbative series of Ĝ(n) is expected to be similar to the one of the unbroken φ 4 theory and hence better than that of the original expansion in λ of G (n) .We can then Borel resum the modified perturbative series of Ĝ(n) and, after that, recover the original Schwinger function by setting λ 0 = λ: Note that whenever λ 0 = λ the Lagrangian L breaks explicitly the Z 2 symmetry.We have then explicitly broken the symmetry, resummed, and switched off the breaking term after the resummation.This is precisely what we are supposed to do any time a vacuum should be non-perturbatively selected in presence of spontaneous symmetry breaking.While the Green functions G (n) should have singularities at the phase transition points, the Green functions Ĝ(n) are smooth for λ 0 = λ and should become singular only in the limit λ 0 → λ.
However, at weak coupling EPT is not as good as ordinary perturbation theory and requires more terms to correctly reproduce the weak coupling expansion of G (n) .In particular, we have verified that with the number of perturbative terms at our disposal, EPT can be reliably used only for the vacuum energy.
4 Perturbative Coefficients up to g 8  3 Order In the 2d φ 4 theory the vacuum energy and the mass are the only terms that require the introduction of counterterms δΛ and δm 2 .In the broken phase, when we perform the vacuum selection by shifting the field φ → φ cl + φ, a divergent one-loop 1-point tadpole is also generated (but the corresponding counterterm is fixed in terms of δm 2 ) as well as a cubic interaction term.Because of the cubic term in the Lagrangian, at a given order in the perturbative expansion the number of topologically distinct diagrams for the broken symmetry phase is much higher than for the symmetric phase.Moreover, diagrams with different number of cubic and quartic vertices contribute at each order in the effective coupling g = λ/ m 2 , making cancellations possible and lowering the numerical accuracy.The task of computing the perturbation series in this theory is then much more challenging.It should also be noted that, in contrast to the unbroken phase, this renormalization scheme is not optimal for computations because of the presence of a radiatively generated 1-point tadpole that should be taken into account order by order in perturbation theory.The computation of the n-point functions will in fact involve diagrams decorated with 1point tadpole terms-i.e.sub-diagrams with zero net momentum flow.A better renormalization scheme can be found by considering the auxiliary Lagrangian L obtained by expanding eq. ( 2.2) around an arbitrary constant configuration φ → (φ cl + φ 0 ) + ϕ and choosing the scale µ such that the divergent 2-pt tadpole is exactly canceled.Choosing φ cl = + m/ √ 8λ and normalizing the vacuum energy so that ∆ Λ = 0 in eq.(2.2), we get where we defined with Z = (4π) −1 log m 2 /m 2 and λ 3 = √ 2λ m.The term 3λ 3 Z for λ 1 and the term 3λ 3 Zφ 0 for ∆Λ come from the running of the linear term due to the one-loop tadpole divergence.The counterterms for the mass and the vacuum energy have the same form as in eq.(2.3) with m replaced by m and the counterterm for the linear term is given by δm 2 λ 3 /(4λ).In this theory all the divergent diagrams are exactly canceled and leave no finite part.If we additionally choose λ 1 such that the 1PI one-point function Γ 1 vanishes, 10 the Lagrangian in eq.(4.1) is the most convenient for computations since we can now forget about both 1-point and 2-point tadpole terms.
Practically the auxiliary theory L can be used to determine the perturbative series in the 2d φ 4 theory L (in the renormalization scheme defined after eq.(2.5)) for arbitrary cubic coupling λ 3 .We first fix the perturbative series for λ 1 in terms of λ 3 and λ by requiring Γ 1 = 0. Up to three loops the non-vanishing diagrams are where we have omitted the multiplicities.Therefore we get where c ij are the dimensionless coefficients obtained by computing all the Feynman diagrams with i cubic and j quartic vertices.We then use the set of eqs.The expansion for λ 3 is trivially obtained by the third eq. in (4.2).At this point we can obtain the perturbative series of our original theory L by re-expanding the parameters m and λ 3 in the series for the auxiliary theory L.Moreover, since ϕ = 0, we have that the VEV is simply given by φ = φ cl + φ 0 .From eq. ( 4.5) we see that φ is picking additional contributions with respect to the 1PI diagrams.At order λ 3 3 λ and λ 5 3 they are Similarly, using the second eq. in (4.5), for the two point function in momentum space Γ 2 = p 2 + m 2 + . . .we see that the expansion of m 2 provides to lowest order the contribution of the following two graphs for the theory L which have to be added to the 1PI graphs at the same order (see eq. ( 4.11) below).Proceeding at higher order one must be careful in expanding also all the higher order terms of the expansion.
In order to verify that the above procedure was correctly implemented, we checked that the series obtained for the VEV and for the 2-pt function matched a direct computation in the theory L up to five loops.
In the following we focus on the 0-, 1-and 2-point functions.Multi-loops computations have been addressed as in ref. [1] using the Montecarlo VEGAS algorithm [36].We refer the reader to ref. [1] for further details.We have computed Feynman diagrams up to the insertion of eight vertices, independently of the nature of the vertex, cubic or quartic.Since λ 3 ∼ m √ λ, this implies that our ordinary perturbative series can reach O( g 4 ).Using instead EPT as described in section 3.2, where parametrically λ EPT 3 ∼ λ, we effectively can reach O( g 8 ).

Vacuum Energy
We have computed all the vacuum energy 1PI graphs with up to eight vertices in the auxiliary theory.The number of topologically distinct graphs (in the chosen scheme) as a function of the number of cubic and quartic diagrams is reported in tab. 2. For illustration, the non-vanishing diagrams with up to three loops are Expanding the parameters as explained above we get the following expression for the series of the vacuum energy is the polygamma function and the numbers in parenthesis indicate the error in the last two digits due to the numerical integration.The coefficients of the vacuum energy Λ for generic values of the couplings λ 3 and λ are reported in the appendix (tab.6).In order to access the strong coupling regime of Λ we will use EPT as described in sec.3.2, therefore we report here the series obtained by setting

1-Point Tadpole
The series coefficients for the VEV have been obtained from Γ 1 as explained above.We have computed all 1PI 1-pt graphs with up to eight vertices in the auxiliary theory.The number of topologically distinct graphs (in the chosen scheme) as a function of the number of cubic and quartic diagrams is reported in tab.Following the method described at the beginning of sec. 4 we get the series for the vacuum expectation value of φ as where φ cl is the tree-level value.The coefficients of φ for generic values of the couplings λ 3 and λ are reported in the appendix (tab.7).

Physical Mass
We define the physical mass as the smallest zero of the 1PI two-point function in momentum space for complex values of the Euclidean momentum: By a perturbative expansion of M 2 in powers of g, the zero of Γ 2 (p 2 = − M 2 ) is determined in terms of Γ 2 (p 2 = − m 2 ) and its derivatives with respect to p 2 (see ref. [1] for further details), which are in turn determined from Γ 2 (p 2 = −m 2 ) and derivatives in the auxiliary theory.The number of topologically distinct 1PI graphs (in the chosen scheme) as a function of the number of cubic and quartic diagrams is reported in tab. 4. For illustration, the non-vanishing diagrams with up to two loops are In tab.8 in the appendix we report the coefficients for Γ 2 (p 2 = − m 2 ) and its derivatives for generic values of the couplings λ 3 and λ.For λ 3 = √ 2λ m we get the following expression for the physical mass 1529(18) g 2 − 14.886(30) g 3 − 50.62(99) g 4 + O( g 5 ) . (4.12)

Large Order Behavior
The large order behavior of the perturbative expansion of n-point Schwinger functions G n in the 2d φ 4 theory in the broken phase has not been determined before.On general grounds, we expect that the coupling expansion n g k behaves, for k 1, as where c n and b n are n-dependent constants, while a is an n-independent constant given by the inverse of t − 1 in eq. ( 3.4): The half-period of oscillation of the large order coefficients is given by π/arg(a) 3.8.The evaluation of the coefficients b n and c n require a detailed analysis of small fluctuations around the instanton configuration which we have not attempted to perform.However, the knowledge of a is enough to allow us to use a conformal mapping, as discussed in sec.3.1.
The shortness of our perturbative series does not allow for a reliable comparison with the asymptotic large order behavior.We just observe that the period of oscillation given by eq.(4.13) is consistent with the same-sign behavior of the first terms in eqs.(4.7), (4.9) and (4.12).

Results
We report in this section the numerical results obtained from the truncated perturbative series and from its Borel resummation.We have considered conformal mapping and Padé-Borel approximants methods to obtain a numerical estimate of the Borel function.Most of the details of our numerical implementation, as well as a short introduction to these resummation methods, has been given in ref. [1], so here we will only focus on some characteristic features of the broken phase.
The generalized conformal mapping method explained in section 3.1 allows us to write any observable F ( g) using eqs.(3.2) and (3.6), with |t ± 1 | and α 1 the modulus and phase of the inverse of the leading (complex) instanton action in eq.(3.4).In terms of the Le Roy -Borel function we have In order to not clutter the notation, we omit to write the (b, s) dependence of F B ( g).As in ref. [1], we introduced two summation variables denoted by b and s to further improve the behavior of the u-series and to have more control on the accuracy of the results.Due to the presence of additional singularities within the unit u-disc, after the conformal mapping (3.6) the series is still asymptotic.The conformal mapping in this case is supposed to extend the region in coupling space where the asymptotic series behaves effectively as a convergent one, giving a more reliable estimate of the observable with respect to optimal truncation of the perturbative series.It is useful to estimate the rate of "convergence" of the series after conformal mapping to have a rough expectation on the accuracy of the results obtained with finite truncations.
In order to simplify our discussion, we might consider the ideal situation where no additional singularities are present.In this case the series in the last line of eq. ( 5.2) is truly convergent for any value of the coupling constant and its rate of convergence is governed by the large n behavior of the u integral.For n 1 a saddle point approximation gives where α 1 is the angle between the real positive axis and the leading singularity in the Borel t-plane appearing in eq.(3.6), C is a b-, s-and α 1 -dependent coefficient which is polynomial in n, and c is a smooth function of α 1 of order one for all values of |α 1 | ≤ 1.We see that the exponential convergence in n of the u integral sensitively depends on α 1 .It is maximal for α 1 = 1, it monotonically decreases for smaller values of α 1 and eventually vanishes for α 1 = 0, when the series is no longer Borel resummable.In the Z 2 unbroken and broken theories we have α 1 = 1, α 1 0.260 respectively.At fixed number of orders we then expect a slower convergence of the conformal mapping method in the broken case with respect to the unbroken one.
Of course the overall convergence of the series in n in eq. ( 5.2) is also determined by the large order behavior of the coefficients B (b,s) n , which is governed by the other singularities in the Borel plane.In the unbroken theory the known singularities are all aligned with the leading one along the negative real axis, the best case scenario, while we have explicitly seen in section 3.1 that this is not the case for the broken case.We then expect that next to leading singularities would most likely further increase the gap in accuracy between resummations in the unbroken and broken phases, though we believe that for g ≤ g (w) c this is a sub-dominant effect.
Due to the shortness of the series we obtained in sec.4, the error minimization procedure of ref. [1] that is used to select the central values for the parameters b and s is sometimes unstable (especially for the 2-point function).We also note that in the broken phase the error minimization procedure tends to select lower values of b than in the unbroken phase.In order to avoid the dangerous b = −1 point at which the Γ-function in the Borel-Le Roy transform diverges, in this paper we fix ∆b = 1, where ∆b is the semirange used to scan the parameter b for the error estimation as in ref. [1].The resulting error estimate typically yields large errorbars but this is expected since we are only resumming few perturbative terms.Tests on simple toy models showed that the error estimate usually correctly represents the difference between the resummed and true values.
In contrast to the unbroken phase, Padé-Borel approximants are not very useful because the same-sign form of the first available series coefficients is typically responsible for spurious unphysical poles that hinder a proper use of this method.Moreover, the shortness of the series does not allow us to systematically select the "best" Padé-Borel approximant as explained in ref. [1].Nevertheless, we report for completeness the results obtained using Padé-Borel approximants, when available.The shortness of the series does not allow us to estimate the contribution to the error coming from the convergence.
We have also computed the vacuum energy using EPT as explained in section 3.2.In this case the leading singularity of the Borel function B(t, g 0 ), at fixed g 0 , is the same as in the unbroken phase, with α 1 = 1, |a| 0.683708.We now have The value of an observable F B ( g) is recovered by eventually setting g 0 = g in eq.(5.4): The different analytic structure of the Borel function and a longer perturbative series allows us to use Padé-Borel approximants in EPT.Unless stated differently, we set in what follows m 2 = 1.

Vacuum Energy: Weak Coupling
The perturbative expression for Λ up to order g 4 is reported in eq.(4.7).We show in the left panel of fig. 5 Λ( g) as a function of g.Surprisingly enough, the vacuum energy series is within the perturbative regime up to the critical coupling g c and well beyond, as evident from the fact that optimal truncation and untruncated perturbation theory coincides for all the couplings shown in the figure and in this whole regime they well approximate the central value of the Borel resummed result (blue dashed line).For g 0.2, the error associated to the resummation rapidly increases, despite the central values remain quite close.
In the right panel of fig. 5 we compare Λ( g) computed using conformal mapping and Pad-Borel resummation techniques.The results with the conformal mapping (5.2) are obtained using resummation parameters s = 1/5 and b = 0.The approximant used in the Padè-Borel method is [1/2] with parameter b = 1.The results are all well compatible with each other, confirming that the vacuum energy is perturbatively accessible in the whole range of the weakly coupled branch.

Vacuum Energy: Strong Coupling and Chang Duality Checks
The perturbative expression for Λ up to order g 8 in EPT is reported in eq.(4.8).As discussed in detail in ref. [24], EPT works at its best at strong coupling.More quantitatively, at fixed number of loops N , we expect that EPT improves over ordinary perturbation theory when g 2 > 1/N . 11t is then the ideal tool to compute the vacuum energy at strong coupling (strong branch in the broken phase).We report in fig.6 Λ as a function of g and compare the results obtained using conformal mapping and Padé-Borel approximants.The light blue line corresponds to the conformal mapping at N = 8.In order to avoid dangerous poles, in the Padé-Borel method (light red) we have removed the vanishing O( g 0 ) and O( g) coefficients from the series and effectively resummed Λ( g)/ g 2 .The approximant shown is  in the unbroken, and weak/strong branches of the broken phases for values of the couplings associated to the same physical theory (points in the same vertical line in fig.1).Indeed, it has been argued in ref. [1] that the Borel resummation of perturbation theory around the unbroken vacuum, when applied beyond the phase transition point and without a proper selection of the vacuum, reconstructs the correlation functions F in a vacuum where cluster decomposition is violated. 12However, since the vacua |± are degenerate, the vacuum energy Λ computed for g ≥ g c starting from the unbroken phase, coincides with the vacuum energies as computed from the broken phase in the weak and the strong branches.We summarize our findings in fig. 7. We report the vacuum energy as a function of the coupling constants in the various phases, g in the unbroken phase (blue), g (w) in the weakly branch of the broken phase (green) and g (s) in the strong branch of the broken phase (red).In comparing the vacuum energy in different phases one has to pay attention to the different units of mass and vacuum energy normalizations in the three descriptions, which are related as in eq.(2.7).In fig.7 we have set to unity the squared mass term m 2 (so that m 2 = 1 in both weak and strong branches) and normalized the vacuum energy Λ to be zero for g = 0 in the unbroken phase.The three vacuum energies are consistent with each other, as expected from Chang duality.We consider this result a numerical check of the Borel summability of the φ 4 theory in the broken phase and an example of the use of EPT in QFT.
Figure 7: The vacuum energy Λ for the theory described by the Lagrangian (2.1) as a function of the coupling constants g (blue), g (w) (green) and g (s) (red).The blue band is computed by resumming the series in the unbroken phase (conformal mapping as in ref. [1]).The green band is obtained via Chang duality from the weak branch of the broken phase (conformal mapping in eq.(3.6))).The red band is obtained via Chang duality from the strong branch using the EPT series in eq.(4.8) (conformal mapping as in ref. [1]).

Tadpole
The perturbative expression for φ , normalized to its classical value, up to order g 4 is reported in eq.(4.9).We have resummed the expression φ /φ cl to the eighth power, because it shows better convergence properties than φ /φ cl .This is not surprising.We know that in the 2d Ising model M ∝ | g c − g| (critical exponent ν = 1) and φ ∝ M 1/8 (critical exponent β = 1/8), so that T ≡ ( φ /φ cl ) 8 ∝ | g c − g| approaches the critical coupling in an analytic way. 13e show in fig.8 T as a function of g in the weak coupling regime.The series for T is within the perturbative regime up to the critical coupling g c and soon after breaks down.The value of g (w) c perturbatively found using the 4-loop series is surprisingly closer to the value expected from Chang duality g The physical mass M as a function of the coupling constant g in the weakly coupled branch using the coefficients up to the g 4 order.We compare the Borel resummation result using conformal mapping (blue dashed line) with optimal truncation (red dotted line).The vertical dotted line at g = g k k ≈ 0.125 is the coupling (taken from ref. [11]) where the mass of the elementary particle is twice the mass of the kink and might become unstable.
lowest energy state in the Hilbert space corresponds to the ground state of the periodic sector, while the ground state in the antiperiodic sector is identified with the kink state in the infinite length limit L → ∞.No spontaneous symmetry breaking can occur at finite L, since the vacuum is a unique state linear combination of |+ and |− , where |± denote the two vacua where at tree-level φ = ±v.The periodic Hilbert space sector is characterized by a quasi degenerate spectrum of states, whose energy splitting decreases exponentially with L and is governed by the energy of the antiperiodic vacuum, i.e. the kink mass.In particular, for L large enough, the lowest energy state beyond the vacuum is the combination of |+ and |− orthogonal to the vacuum, which becomes degenerate with it for L → ∞.In this limit all states in the periodic sector become exactly degenerate and the Hilbert space is expected to also contain states that can be seen as composed of an even number of kink and anti-kink states. 14Analogously, states in the antiperiodic sector can be interpreted as composed of an odd number of kink and anti-kink states.A superselection rule forbids transitions that do not preserve a Z 2 topological charge, the kink number operator.Semi-classical arguments [37] suggest that the elementary φ-particle excitation decays into a pair of kink anti-kink states at some value of the coupling g k k < g c .The presence of such decay has been checked numerically in ref. [11] (see also ref. [38]), where the absence of single particle states in the periodic sector for g k k ≈ 0.125 is interpreted as its decay in a pair of kink and anti-kink states.
The φ 4 theory discussed in refs.[10,11], in the L → ∞ limit, becomes a non-compact theory with two degenerate vacua connected by topological kink excitations.No vacuum selection has  [10,11] of the physical mass as a function of g.The blue dashed line are our results obtained by Borel resummation with conformal mapping, the black points are the results of ref. [11] extrapolated at infinite volume, the red points are the results of ref. [10] at finite volume L = 20.The vertical dotted line at g = g k k ≈ 0.125 is the coupling (taken from ref. [11]) where the mass of the elementary particle is twice the mass of the kink and might become unstable.
been performed, and therefore φ = 0 and no spontaneous symmetry breaking can occur.In contrast, our results are based on ordinary perturbation theory in non-compact space, expanding around |+ or |− , where spontaneous symmetry breaking occurs.This point should be taken into account when comparing our results with those of refs.[10,11].We do not expect subtleties related to the choice of vacuum for the vacuum energy Λ, since this is a continuous and smooth function at least up to g (w) c .Similarly, M is a smooth function as long as the particle is stable and should coincide with the mass of the lightest single particle state for g < g k k.
We compare in fig. 10 the values of Λ and M computed respectively in refs.[10] and [10,11] with our results. 15Note that the results of ref. [10] have not been extrapolated at infinite volume: for Λ we plot their points for two different volumes L = 12 and L = 20, while the points for M are at L = 20.In tab. 5 we make the comparison explicit for some values of the coupling g.The values we report for Λ as computed by ref. [10] are obtained as the means of the values at two renormalization scales µ = 0.9 m and µ = 1.1 m and the reported error is the semidifference.The values of M taken from ref. [11] have been normalized accordingly to our definitions.As it can be seen, we get more accurate results than those of refs.[10,11] and they are all in good agreement among themselves. 16or g k k > g, there might be subtleties in the interpretation of M related to spontaneous symmetry breaking.Indeed, the proper way to select a vacuum when the Z 2 symmetry is spontaneously broken is achieved by adding a small explicit breaking term such as φ and take the   [10] and [11] respectively.The results of ref. [10] are at finite volume L and the reported error is obtained as the semidifference between the results at two renormalization points, see the text for additional details.
limit → 0 only after L → ∞.The situation considered in refs.[10,11] corresponds to the opposite order of limits, → 0 first and L → ∞ after, since no breaking term was present to begin with.As usual in perturbative QFT, we performed a selection of the vacuum "by hand", by choosing to expand around any of the two vacua, neglecting the effect of the other.There is no need to add an explicit breaking term, so our configuration is equivalent to having taken L → ∞ first, since we have infinite volume to start with, and → 0 after.Since the two limits do not commute, the results for the particle decay found in ref. [11] do not obviously apply in our context.When L → ∞, at finite , the non-trivial topological Hilbert space sector containing an odd number of kink and anti-kink states decouple.While topologically trivial states of kink anti-kink do not decouple, the fundamental particle can no longer kinematically decay into freely moving kink and anti-kink states, since single kink and anti-kink states are no longer in the spectrum.In other words, the φ-particle for any finite would behave as a spatially extended but stable bound state.As becomes smaller and smaller, this state becomes less and less bound and for = 0 and g ≥ g k k it unbounds to a pair of essentially free kink-anti-kink states.In this case our results for M do not have a clear interpretation, since they have been obtained assuming the existence of a pole of the two-point function.But in the topological trivial sector no single particle state would remain, and the operator φ would only create multi-particle states.In other words, the pole of the two-point function would dissolve in a branch-cut singularity for g ≥ g k k.
In this case M would be related to the threshold energy for the multi particle production, or perhaps it would simply be an unphysical analytic continuation with no obvious significance.As we will see in the next subsection, the second hypothesis is favored by our analysis.

Kink States from the Unbroken Phase?
As we mentioned in the last subsection, we cannot have a direct access to single kink and antikink states starting from a vacuum where spontaneous symmetry breaking occurs.Interestingly enough, however, we might possibly access the kink sector of the theory starting from the unbroken theory with m 2 > 0! Indeed, it has been conjectured in ref. [1] that the Borel resummation of correlation functions starting from the unbroken phase reconstructs for g > g c the correla- as a function of g in the weakly coupled branch.In brown we have |M | as obtained from the Borel resummation of the perturbative series in the unbroken phase up to order g 8 using conformal mapping.The black dashed line is the kink mass as given by the semiclassical equation (6.1).The green and red points are the results obtained by ref. [11] from the splitting and antiperiodic sector respectively.See the main text for further details.
tion functions in a vacuum linear combination of |+ and |− connected by kink configurations, where φ = 0 and cluster decomposition is violated, that is the configuration in refs.[10,11].A consequence of this conjecture is that the vacuum energy Λ for g > g c should be identified with the vacuum energy Λ as computed in the weakly and strongly coupled branches.More interestingly, |M (g)| as computed in ref. [1] for g g c should be identified with the mass gap in the non-clustered vacuum, which is given by the kink mass.The kink is indeed the lightest single particle excitation close to the phase transition (in fact, the only one when g is sufficiently close to g c ).The kink mass as a function of g has been numerically studied in ref. [11].It turns out that the semi-classical kink mass formula, including one-loop corrections, is in very good agreement with the numerical results for all values of g in the weakly coupled branch. 17In fig.11 we compare the results of ref. [11] with eq. ( 6.1) and our results for |M | as computed in ref. [1] and analytically continued beyond the phase transition.The red and green points are the results of ref. [11].The former are directly obtained by computing the vacuum energy in the anti-periodic sector, while the latter are obtained by computing the energy splitting between the first two states in the topological trivial sector. 18Note that the region of |M (g)| beyond the transition for g > g c correspond to g < g  11, while the latter using conformal mapping in the broken phase, as in fig.9.The vertical dotted line at g = g k k ≈ 0.125 is the coupling (taken from ref. [11]) where the mass of the elementary particle is twice the mass of the kink and becomes unstable.The color coding is the same as in figs.9 and 11.
m 2 = 1, while in ref. [1] we had m 2 = 1.The results are in very good agreement with each other, providing evidence to our conjecture.The relation between perturbative and non-perturbative states in different phases of a theory is typical in theories enjoying duality symmetries, but we are unaware of relations of this sort that involve correlation functions in a vacuum that does not satisfy cluster decomposition.
We conclude by comparing M and 2M kink as a function of g, see fig.12.This is useful because if M still describes the mass gap in the theory it should follow the curve 2M kink .Though not conclusive, fig.12 does not support this hypothesis and seems to suggest instead that for g ≥ g k k M is an analytic continuation with no obvious interpretation.

Conclusions
In this paper we have studied the 2d φ 4 theory in the broken phase, recently shown to be Borel resummable [1].We have computed the leading finite action complex instanton solutions, important to determine the large order behavior of the perturbative series and to Borel resum it using a generalized conformal mapping method.We have computed the perturbative series expansion for the first Schwinger functions up to order g 4 and Borel resummed the truncated series using our generalized conformal mapping technique.The results of the resummation are not as accurate as in the unbroken phase, but allow us to establish that the weakly coupled branch of the broken phase is almost entirely within the perturbative regime.This somewhat unexpected result is fully confirmed by comparing our perturbative (resummed or not) results with hamiltonian truncation methods.
We have also used EPT to compute the vacuum energy at strong coupling, and compared the results to the ones obtained in the weakly coupled branch and in the strongly coupled unbroken phase, proving in this way Chang duality and the Borel summability of the theory.We have finally provided a numerical evidence that the mass gap analytically continued from the unbroken to the broken Z 2 phase can be identified with a kink state.This result is in agreement with the expectation that the analytically continued Schwinger two-point function in the broken phase corresponds to the mass gap in a inhomogeneous vacuum where cluster decomposition does not hold.It would be nice to have a deeper understanding of this phenomenon and to establish if and to what extent it applies for other Schwinger functions and for other Borel resummable theories undergoing phase transitions.
7 g E D 2 p Q h 3 t o Q A s Y I D z D K 7 w 5 j 8 6 L 8 + 5 8 L F s L T j 5 z C n / g f P 4 A z H W M 6 w = = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " o f+ A 8 t H p D 1 Q H r P 7 8 9 n Y B Z l 3 6 N Z o = " > A A A B / X i c b V D L S s N A F J 3 4 r P U V H z s 3 g 0 W o m 5 K I o M u i G 5 c V 7 A P a G C a T 2 3 b o Z B J m J p Y a i r / i x o U i b v 0 P d / 6 N 0 z Y L b T 1 w 4 X D O v d x 7 T 5 B w p r T j f F t L y y u r a + u F j e L m 1 v b O r r 2 3 3 1 B x K i n U a c x j 2 Q q I A s 4 E 1 D X T H F q J B B I F H J r B 4 H r i N x 9 A K h a L O z 1 K w I t I T 7 A u o 0 Q b y b c P O 0 M W g m Y 8 B N z z 6 X 1 W H p 6 O f b v k V J w p 8 C J x c 1 J C O W q + / d U J Y 5 p G I D T l R K m 26 y T a y 4 j U j H I Y F z u p g o T Q A e l B 2 1 B B I l B e N r 1 1 _ b a s e 6 4 = " o f + A 8 t H p D 1 Q H r P 7 8 9 n Y B Z l 3 6 NZ o = " > A A A B / X i c b V D L S s N A F J 3 4 r P U V H z s 3 g 0 W o m 5 K I o M u i G 5 c V 7 A P a G C a T 2 3 b o Z B J m J p Y a i r / i x o U i b v 0 P d / 6 N 0 z Y L b T 1 w 4 X D O v d x 7 T 5 B w p r T j f F t L y y u r a + u F j e L m 1 v b O r r 2 3 3 1 B x K i n U a c x j 2 Q q I A s 4 E 1 D X T H F q J B B I F H J r B 4 H r i N x 9 A K h a L O z 1 K w I t I T 7 A u o 0 Q b y b c P O 0 M W g m Y 8 B N z z 6 X 1 W H p 6 O f b v k V J w p 8 C J x c 1 J C O W q + / d U J Y 5 p G I D T l R K m 26 y T a y 4 j U j H I Y F z u p g o T Q A e l B 2 1 B B I l B e N r 1 1 _ b a s e 6 4 = " o f + A 8 t H p D 1 Q H r P 7 8 9 n Y B Z l 3 6 N Z o = " > A A A B / X i c b V D L S s N A F J 3 4 r P U V H z s 3 g 0 W o m 5 K I o M u i G 5 c V 7 A P a G C a T 2 3 b o Z B J m J p Y a i r / i x o U i b v 0 P d / 6 N 0 z Y L b T 1 w 4 X D O v d x 7 T 5 B w p r T j f F t L y y u r a + u F j e L m 1 v b O r r 2 3 3 1 B x K i n U a c x j 2 Q q

Figure 1 :
Figure 1: Phase structure of the φ 4 theory according to Chang duality, as a function of the couplings g and g.Points in the same vertical line connecting the black and red lines correspond to different descriptions of the same theory.

Figure 2 :
Figure 2: Trajectories in the complex ϕ plane of the solutions of eq.(3.3) associated to the three leading saddles (blue) and their complex conjugates (orange).In increasing order of |t i | they are represented by continuous, dotted and dashed curves.All trajectories start for ρ = 0 in three different points of the complex plane with zero velocities ∇ϕ(0) = 0 and flow as ρ → ∞ to the same value ϕ(∞) = 1.

Figure 3 :
Figure 3: The change of variable (3.6) maps the cut t-plane (left panel) into the disk of unit radius |u| = 1 (right panel).The branch points t ± 1 are mapped into the points u ± 1 ; the real positive axis t ∈ [0, +∞) is mapped in the segment u ∈ [0, 1); the rays connecting the branch points t ± 1 to the point at complex infinity are mapped to the arcs at the boundary of the disk as shown by the colors.In drawing this figure we picked α 1 = 1/4.

Figure 4 :
Figure 4: Position of the three leading singularities of B(t) in the t plane (left) and in the u-plane (right) after the conformal mapping of eq.(3.6).The singularities u ± 2 and u ± 3 do not lie at the unit circle |u| = 1, though they are very close to it.

Table 3 :
Number of topologically distinct 1PI 1-pt diagrams without self-contractions with maximum eight total vertices for the auxiliary theory.

Figure 5 :
Figure 5: The vacuum energy Λ as a function of the coupling constant g using the coefficients up to the g 4 order.Left panel: the results obtained using the conformal mapping (5.2) (light blue) with parameters s = 1/5, b = 0 and the optimally truncated series (red dotted line) which coincides with the perturbative series up to g 4 for the couplings shown in the plot.Right panel: comparison between the central values of the conformal mapping (blue dashed line) and Pad-Borel (solid green) resummation techniques.The regions beyond the phase transition g > g (w) c are shaded in gray.
[3/2]  with b = −1/2.The results are in good agreement with each other.We can numerically check Chang duality in its full glory by comparing the vacuum energies 0

Figure 6 :
Figure 6: The vacuum energy Λ at strong coupling obtained by resummation of the EPT series (4.8) with the conformal mapping (light blue) and Pad-Borel (light red) methods.The region beyond the phase transition g > g (w) c is shaded in gray and the vertical dashed line signals the critical coupling in the strong branch g = g (s) c .

27 .Figure 9 :
Figure9: The physical mass M as a function of the coupling constant g in the weakly coupled branch using the coefficients up to the g 4 order.We compare the Borel resummation result using conformal mapping (blue dashed line) with optimal truncation (red dotted line).The vertical dotted line at g = g k k ≈ 0.125 is the coupling (taken from ref.[11]) where the mass of the elementary particle is twice the mass of the kink and might become unstable.

Figure 11 :
Figure 11:  Comparison between the kink mass and the mass |M | in the unbroken phase for g > g c as a function of g in the weakly coupled branch.In brown we have |M | as obtained from the Borel resummation of the perturbative series in the unbroken phase up to order g 8 using conformal mapping.The black dashed line is the kink mass as given by the semiclassical equation (6.1).The green and red points are the results obtained by ref.[11] from the splitting and antiperiodic sector respectively.See the main text for further details.

Figure 12 :
Figure 12:  Comparison between 2M kink and M as a function of g in the weakly coupled branch.The former has been obtained multiplying by a factor 2 the data in fig.11, while the latter using conformal mapping in the broken phase, as in fig.9.The vertical dotted line at g = g k k ≈ 0.125 is the coupling (taken from ref.[11]) where the mass of the elementary particle is twice the mass of the kink and becomes unstable.The color coding is the same as in figs.9 and 11.
bative series for the physical mass M 2 .

Table 8 :
The coefficients b (n) k,l relevant for the determination of the pole mass M 2 .

Table 1 :
Values of the critical coupling in the φ 4 theory.The value g

Table 2 :
Number of topologically distinct 1PI 0-pt diagrams without self-contractions with maximum eight total vertices for the auxiliary theory.

Table 4 :
Number of topologically distinct 1PI 2-pt diagrams without self-contractions with maximum eight total vertices for the auxiliary theory.
3. The non-vanishing diagrams with up to three loops have been already reported in eq.(4.3).
Figure 10: Left panel: comparison with ref. [10] of the vacuum energy as a function of g.The blue dashed line are our results obtained by Borel resummation with conformal mapping, black and green points are the results of ref. [10] at finite volume L = 12 and L = 20, respectively.Right panel: comparison with refs.

Table 5 :
The values of Λ and M for some values of g and comparison with refs.

Table 6 :
Perturbative coefficients for the vacuum energy Λ with independent cubic and quartic coupling λ 3 and λ up to eight total vertices.We set m = 1 to avoid clutter.19

Table 7 :
Perturbative coefficients for the VEV with independent cubic and quartic coupling λ 3 and λ up to eight total vertices.We set m = 1 to avoid clutter.