Life at the Landau pole

If a quantum field theory has a Landau pole, the theory is usually called 'sick' and dismissed as a candidate for an interacting UV-complete theory. In a recent study on the interacting 4d O(N) model at large N, it was shown that at the Landau pole, observables remain well-defined and finite. In this work, I study both relevant and irrelevant deformations of the said model at the Landau pole, finding that physical observables remain unaffected. Apparently, the Landau pole in this theory is benign. As a phenomenological application, I compare the O(N) model to QCD, by identifying $\Lambda_{\overline{\rm MS}}$ with the Landau pole in the O(N) model.


I. MOTIVATION
In the early days of quantum field theory, Landau and collaborators studied quantum electrodynamics in perturbation theory [1].They found that QED has positive β-function in perturbation theory, recognizing that this would lead to an uncontrolled growth of the theory's coupling constant as a function of energy.In modern notation, the QED running coupling to leading order in perturbation theory becomes where it is customary to fix Λ LP = m e e 3π 2α 0 with m e the electron mass and α 0 ≃ 1 137 .Landau noted that besides the dependence of the fine-structure constant α on the momentum scale μ, the form (1) implied that the running coupling is diverging (has a pole) at a finite momentum scale μ = Λ LP .Since the coupling diverges at this scale, it seems that one cannot meaningfully probe momentum scales μ > Λ LP in QED, so the theory does not have a well-defined continuum limit.It has even been suggested that Landau was so disturbed by this feature that he quit working on quantum field theory 12 .
Modern physics deals with the issue of the Landau pole through a mix of denial and shoulder shrugging.Denial adherents will rightly point out that (1) was derived in perturbation theory, requiring α ≪ 1, so that as a consequence (1) cannot be expected to correctly capture features such as the Landau pole, where by definition α → ∞.Shrugging adherents will (also rightly) point out that Λ LP = m e e 3π 2α 0 ≃ 10 280 MeV (2) puts the scale of the Landau pole beyond the Planck scale, so that in practice it is entirely pointless to understand QED in that regime anyway.(However, it should be noted that in the full Standard Model, Λ LP ≃ 10 34 GeV is much lower than (2), but still extremely high [6].)The prevailing dogma in both cases is that theories with a Landau pole should be 1 I thank Bill Zajc for pointing out that this statement is not true, assuming the the publication dates of papers are chronologically matched to his actual work.However, Landau's view of QFT seems to have been fairly negative after 1954, e.g. in Ref. [2] the authors write "This brings us to the conclusion that point interaction is impossible in pure electrodynamics."and in Ref. [3] Kirzhnits and Linde write about Ref. [4] as "..., which even led some authors to the conclusion that the hamiltonian quantum field theory "is dead"...". 2 See Ref. [5, chapter 8] for a first-hand historical account of Landau's approach to QED. viewed as 'UV-incomplete', or cut-off theories, that cannot be used to describe continuum physics.
In this work, I will entertain an entirely different perspective, namely that physical observables of a quantum field theory could be well-behaved even when the coupling diverges at the Landau pole, and beyond.
Unfortunately, I am unable to test my perspective in QED (yet), even though other groups have proposed similar ideas [7,8].Instead, I will focus on a quantum field theory which can be solved non-perturbatively in the limit of a large number of components, namely the O(N) model [9], using critically important input from PT -symmetric field theory [10].
I will work in 3+1 dimensions where the O(N) model is known to possess a Landau pole in the large N limit.
In a quantum field theory, the renormalized coupling is not directly observable -as is apparent through the fact that it will depend on the fictitious renormalization scale μ.For this reason, finite and well-defined physical observables at infinite renormalized coupling are certainly possible.Indeed, the idea that physical observables turn out to be finite even when the coupling diverges is well supported by several quantum field theory examples, such as N = 4 SYM in 3+1 dimensions [11,12], bosonic and fermionic large N field theories in 2+1 dimensions [13][14][15][16], as well as non-relativistic fermions at a Feshbach resonance where experimental confirmation is available [17].
In this work, I build upon and extend this idea: if physical observables remain finite when the renormalized coupling parameter diverges, maybe observables remain finite and well-defined when the renormalized coupling parameter becomes negative or even complex.After all, the renormalized coupling parameter is not directly observable, so nothing should protect it from becoming complex as long as physical observables remain well-defined.Of course this type of idea cannot be tested in perturbation theory, which is inherently a weakcoupling expansion around non-interacting field theory.For this reason, I heavily employ large N expansion techniques (which do not rely on a small perturbative coupling in order to be applicable) as well as results from PT -symmetric field theory (which allow calculations for negative or complex couplings via analytic continuation).
While this study is exploratory, I nevertheless hope that certain aspects merit further consideration when trying to interpret quantum field theory in four dimensions.

II. CALCULATION -THE O(N) MODEL IN 3+1 DIMENSIONS AT LARGE N
Let me first consider the case of the massless theory with quartic interaction, which is essentially a repeat of the calculation in Ref. [9], but included here for completeness.The partition function for this theory is defined through the path integral with the Euclidean action where ⃗ ϕ = (ϕ 1 , ϕ 2 , . . ., ϕ N ) is an N-component scalar field and the theory is defined on the thermal cylinder with β = 1 T the inverse temperature.The partition function may be rewritten in a more convenient form by introducing two auxiliary fields σ, ζ with (cf.Ref. [18]) The resulting path integral for σ has quadratic action, such that σ can be integrated out.
One finds Separating the auxiliary field into zero modes and fluctuations , one can verify that the path integral over fluctuations does not contribute to leading order in large N to the partition function.Since ζ 0 is a constant, the path integral over fields ⃗ ϕ is quadratic and can be done in closed form.One finds Z(λ, β) = dζ 0 e N βV p( √ 2iζ 0 ) , with the pressure per component in dimensional regularization where βV is the space-time volume, μ is the MS renormalization scale, K i (x) denotes modified Bessel functions of the second kind and I have rewritten iζ 0 = m 2 2 to simplify the appearance.
The expression (7) is divergent in the continuum lim ε → 0. However, it may be nonperturbatively renormalized as  (9) as well as for the coupling g R in the analytically continued theory (PT -symmetric theory).Adapted from Ref. [9].
which is standard procedure for large N field theories [19].The resulting running coupling is given by which has a Landau pole at μ = Λ LP , cf.Fig. 1.
In order to make sense of the theory at the Landau pole, a procedure for analytically continuing the theory beyond the Landau pole is necessary.This procedure has been provided in the form of a conjecture in Ref. [10] for so-called PT -symmetric field theory3 .Naively continuing (9) for μ > Λ LP , the sign of λ R becomes negative.So in order to make sense of the O(N) model beyond the Landau pole, one is led to consider a theory where the sign of the coupling is flipped: where the small imaginary part has been included in order to be able to 'go around' the Landau pole.Following standard nomenclature [20], the theory with flipped-sign coupling is referred to as PT -symmetric field theory, and its partition function is denoted by Z PT (g, β).
Following Ref. [10], the analytic continuation of Z(λ, β) is given by To evaluate Z PT (g, β), one may directly employ the pressure function (7), where now the sign of the coupling has been flipped.Regardless of the sign of the coupling, the expression for the pressure (7) is divergent in the 4d continuum lim ε → 0. The PT -coupling renormalization, given by 1 differs from (8) by a sign, making the PT -symmetric theory asymptotically free.Using the renormalized coupling g R (μ), one can express the renormalized pressure function as The pressure, being a physical observable, cannot depend on the choice of renormalization scale μ, so dp d ln μ = 0.This fixes the running for the renormalized coupling g R (μ) and as a consequence the form of the running coupling itself as which is the same as (9) up to a sign.Both running couplings are shown in Fig. 1, where it can be seen that they diverge at μ = Λ LP .Whereas μ = Λ LP is the Landau pole in the O(N) model, when interpreted through PT -symmetric field theory, Λ LP has all the trappings of an infrared scale parameter reminiscent of QCD (cf.section III).
For completeness, the β-function for the PT -symmetric theory in which is seen to be negative in d From the form (15) it is also easy to discuss the fix-point structure of the theory.In d = 4 − 2ε dimensions, the PT -symmetric O(N) model possesses two fixed points, one at g R = 0 and the other at g R = −4π 2 ε.From the running coupling g R (μ) shown in Fig. 1, it is clear that g R = −4π 2 ε corresponds to the IR fixed-point of the theory.Flipping the sign of both the coupling and ε, one finds that this fixed point is the same as the UV fixed point of the original theory specified by ( 4) in d = 4 + 2ε dimensions.This is consistent with the known critical dimensions 4 < d < 6 for which the O(N) model can be non-perturbatively renormalized [22], see also [23].
Inserting the form ( 14) of the running coupling into (13), one obtains At this point it is worth to pause and consider the following observations: • Any dependence of p(m) on the unphysical renormalization scale μ has dropped out • The expression ( 16) is identical to (13) evaluated at a the Landau pole, μ = Λ LP • The pressure in the (flipped sign coupling and asymptotically free) PT -symmetric theory ( 16) is identical to the the pressure in the original O(N) model ( 7), as can be seen by inserting Eqns.( 8), ( 9) into (7) • For generic values of m, the pressure at the Landau pole will be finite rather than infinite.
As a consequence of these observations, the O(N) model and the PT -symmetric O(N) model are identical at large N, even though from Fig. 1 one of these has a Landau pole and the other one is asymptotically free.For this reason I will refer to the scale μ = Λ LP as 'the Landau pole' also in the PT -symmetric theory.
One can calculate the pressure for any temperature T by noting that the remaining single integral over ζ 0 in the partition function is again dominated by the saddle points at large N, so that the physical pressure of the theory per component is given by ( 16) with m = m the solution to If (17) has more than one solution (as it generically does), then in the large N limit, the solution with the biggest real part of p( m) will dominate over all others.At zero temperature (a.k.a. the vacuum), a simple solution to ( 17) is m = 0, which corresponds to the usual starting point for perturbative calculations.However, there is a second solution to (17) located at m = Λ LP √ e, which is usually dismissed as 'being too close to the Landau pole' [19].However, the physical pressure per component for this second solution is given by which is perfectly finite.In addition, since p(m = Λ LP √ e) > p(m = 0), the solution ( 18) is thermodynamically preferred over the perturbative vacuum 4 .
Despite the presence of the Landau pole, observables in the bosonic theory seem to make physical sense.For instance, one can calculate thermodynamic properties at finite temperature by tracking solutions to (17) numerically, and evaluating p( m).As discussed in Ref. [9], at small temperature the numerical solution m of the thermodynamically preferred phase is continuously connected to m(T = 0) = √ eΛ LP .The numerical solution m becomes complex above a critical temperature T = T c ≃ Λ LP / √ e, but using results from PTsymmetric field theory [10], the analytically continued pressure is continuous across T = T c .
A plot of the pressure as a function of temperature from Ref. [9] is reproduced in Fig. 2, and results for the entropy and specific heat can be found in Ref. [9].
Besides the pressure, one may look at other observables.For instance, one can calculate the spectral function for the auxiliary field ζ ′ , finding a stable bound state in the low temperature phase [26,27].This is different from the situation in QED, where one encounters a tachyon (also sometimes called Landau's ghost) instead of a stable bound state [3].This is the main reason why this work is concerned with the Landau pole in the O(N) model rather than QED.

A. Adding Relevant Deformations
One might worry that the results from the previous sections are an artifact of tuning away all relevant and irrelevant operators.For this reason, it is useful to consider repeating the analysis for the Euclidean action Introducing the auxiliary fields as before in (5), one may again integrate out σ, and one finds in complete analogy with the previous section the PT -symmetric pressure function The explicit ε → 0 divergence can again be taken care of by using the same nonperturbative renormalization as before (12).This leads to Since the bare coupling 1 g diverges as ε → 0, there are residual divergences remaining in (21).The first one of these can be taken care of by renormalizing the bare mass parameter as The running of the renormalized mass m R (μ) is again fixed by requiring that dp(m) d ln μ = 0, which leads to with constant and finite mass scale m 0 .This leads to For the remaining term, note that in the limit ε → 0 so that the pressure function becomes For small values of m 0 , the properties of this theory are close to the unmodified version considered in the previous section.The second-order phase transition at finite temperature persists, but is pushed to higher values of Tc Λ LP .

B. Adding Irrelevant Deformations
Now let us consider what happens when adding irrelevant operators to the theory.In this case, I study where α is the bare sextic coupling parameter.Introducing the auxiliary fields as before in (5), it is possible, but not very enlightening, to integrate out σ exactly.Instead, in the large N limit it is again permissible to replace σ(x) by just its global zero mode σ 0 , so that where At large N, the integral over σ 0 is done with the saddle point method, with two saddles located at For small 3αm 2 2g 2 (justified below), one can expand the square root in this expression to obtain for the two solutions.Inserting σ 0 into (29), one finds Renormalizing the coupling g as in (12), leads to but this implies that α g 3 is divergent.Thus the bare sextic coupling parameter also needs to be renormalized as This leaves the whole tower of additional terms α 2n g 2n+3 , n ≥ 1 that are potentially divergent.However, one finds that in the ε → 0 limit because n ≥ 1.Therefore, none of these terms contribute, and one is left with where I have used the renormalization group invariance of the pressure to express with constant mass scale M .One observes that for the same reason as (35), expanding the square root in (30) is justified for σ The explicit 1 ε divergence can by renormalizing the coupling g, but the sign of the counterterm must be flipped (and as a consequence, so must the sign of the running coupling).One obtains Renormalizing the sextic coupling α as in (37), the square-root expansion in (30) is justified also for σ (2) 0 .Similar to (35), terms with positive powers of the bare coupling g in the numerator vanish, so that one finds This expression still has a divergent term ∝ 1 g 3 , but this term is independent from m.For this reason, this last divergence can be canceled by a vacuum pressure-counterterm in the Lagrangian, leading to Inspecting ( 36) and ( 41), one finds that the two saddle point solutions for σ 0 give rise to the same form for the pressure function, except for the sign of the m 6 term, which can be attributed to the fact that the sign of g R (μ) is flipped for the solution σ (1) 0 .The final integral over ζ 0 is done by finding the saddle point solution where ± corresponds to the solutions σ  (18).As M 2 decreased, the situation remains qualitatively the same until M < ∼ 84Λ LP , at which point solutions (except for m = 0) become complex-valued.
However, close to the Landau pole where g R (μ = Λ LP ) → ∞, (37) suggests that M 2 → ∞, so I will not consider small values of M in the following.
For the negative sign in (42) and large M 2 , there are two solutions m = 0 and m ≃ Λ LP √ e, where again the second solution is thermodynamically preferred.Therefore, one also recovers the unmodified theory solution (18) for the second solution σ (2) 0 .

C. Observables and Running Coupling
As mentioned above, the running coupling of a quantum field theory is not directly observable, because it is a renormalization-scheme dependent quantity.However, observables in a quantum field theory can depend in a non-trivial manner on the coupling parameter of the theory.
For instance, it is possible to calculate scattering amplitudes in the O(N) model, which are closely related to observables such as cross-sections.It turns out that the Euclidean momentum s-wave scattering amplitude to leading order in large N is given by the propagator of the auxiliary field ζ ′ introduces in (5), which for the PT -symmetric theory becomes [9] M I will use M(k) as an example of an observable quantity.
In weak coupling perturbation theory, one can expand M(k) for g R ≪ 1, so From a perturbative QFT point of view, one can therefore define an "observable" or "effective" coupling as where the identification μ ≃ k is expected to hold only approximately.This is a reasonable thing to do as long as the resulting g R,ef f is small enough in order for the observable M to be well approximated by perturbation theory.
However, since M(k) is an observable, it cannot depend on μ, and indeed it does not.
Using the explicit form of g R (μ) from ( 14), one finds This result is fully non-perturbative, and can be further simplified by noting that to this order in the 1/N expansion, the value of m is given by its saddle point value m = m = √ eΛ LP , such that For sufficiently large Euclidean momenta k ≫ m where M(k) becomes small, this leads to Comparing this to (45) and ( 14) informs the refined definition for the effective coupling I am now in a position to compare the "physical" definition of the coupling g R,ef f (μ = k) to the running coupling of the theory (14), which I do in Fig. 3. From this figure, it can be seen that the effective coupling g R,ef f (μ) is positive and finite for all energy scales μ.
Secondly, one finds that the effective coupling g R,ef f (μ) reaches a finite value in the infrared: It is trivial to see that these properties of the effective coupling g R,ef f are a direct consequence of the definition (49) in terms of a "physical" quantity (the scattering amplitude).
In particular, g R,ef f does not exhibit any remarkable feature near μ = Λ LP despite the fact that the O(N) model possesses a Landau pole at this scale.
This makes sense: since the running coupling itself is not an observable, one can always find some "physical" definition of a coupling that agrees with perturbation theory whenever the coupling is small, and has any set of properties one would want such as positiveness and finiteness in the IR.

D. Conclusions
The O(N) model in 3+1 dimensions has a Landau pole at large N. Physical observables at the Landau pole remain finite and well-behaved.This feature does not change when adding either relevant operators (e.g.mass terms) or irrelevant operators (e.g.sextic interactions) to the theory.I therefore conclude that for the O(N) model in 3+1 dimensions at large N, the Landau pole is a harmless feature of the theory, and not a sign that the theory itself is 'sick'.Most importantly, using the analytic continuation provided by PT -symmetric field theory, the O(N) model does not need to be treated as a cut-off theory.It is UV complete, and asymptotically free, despite (or perhaps because of ) the negative coupling constant.
An important omission in the present study is the question about 1 N corrections -will they destroy the features of the large N limit?While this question is without doubt important, I cannot resist pointing out that (almost) the entirety of holography is built upon the strict large N limit of field theory, without systematic discussion of 1  N corrections [30].Yet even without systematic understanding of 1 N terms, holography has indisputable been useful in building our understanding of quantum field theory, so maybe a similar attitude could be extended to the large N limit of the O(N) model.

III. A LANDAU POLE IN QCD?
Let me conclude the discussion by advancing a provocative idea: a Landau pole in Quantum Chromodynamics.According to lore, QCD does not have a Landau pole, so it seems that it should not be discussed here.But then again, there is no hard evidence against a Landau pole in QCD, so I would argue that the idea should be considered5 .QCD does Lattice QCD point at μ = 1.5 GeV is from Ref. [28], QED point at the Z-pole mass μ = M Z is from Ref. [29].See text for details.
have asymptotic freedom, and since the PT -symmetric O(N) model at large N shares this property, one can nevertheless ask how different or similar these theories are.
It would seem more appropriate to attempt a comparison to QED rather than QCD, but the experimental determination of the QED fine-structure constant does not extend much beyond the Z-pole mass μ = M Z ≃ 91 GeV.Since the QED Landau pole (2) is at vastly higher energy, the resulting QED information (shown in Fig. 4) is not particularly illuminating.
So instead of QED, in Fig. 4 I compare the running coupling from QCD to that in the PTsymmetric O(N) model.For the QCD running coupling, I am using the 3-loop perturbative QCD expression resulting from numerically integrating where a s = αs(μ) 4π and f and I am taking N f = 5 [29,33].As can be seen from Fig. 4, the QCD running coupling obtained from perturbation theory becomes very large at small scales μ.In fact, one finds that the perturbative solution for α s (μ) thus obtain diverges for a particular value of μ = Λ MS ≃ 0.3 quarks with physical masses [31], with a cross-over temperature reported as T c = 170(4)(3) MeV [32] (full vertical line).Massless O(N) model results with a second-order phase transition located at T c ≃ 0.616Λ LP (dashed vertical line) are adapted from [9].See text for details.
GeV.This value is consistent with similar values reported by other methods [34], and indeed the running coupling α s (μ) is consistent with calculations from lattice QCD [28] at values as low as μ ≃ 5Λ MS (also shown in Fig. 4).
In my opinion, Fig. 4 indicates a certain qualitative similarity between QCD and the O(N) model.Pushing the similarity further, this would lead to the interpretation that α s (μ) actually does diverge at a finite momentum scale Λ MS , and that deep in the infrared α s should be analytically continued to negative (or complex) values.
It is well-known that QCD becomes confining in the infrared, and that physical observables are well-behaved and finite for all momentum scales.In particular, thermodynamic quantities such as the pressure are continuous as a function of temperature for QCD, with a broad analytic crossover from confined to quark-gluon plasma phase around T c ≃ 170 MeV [32].Normalizing the pressure by the Stefan-Boltzmann pressure p SB (T ) = This was to be expected: the O(N) model simply is not the same as QCD, despite sharing the property of asymptotic freedom.At best, the O(N) model may serve as a crude approximation of QCD, similar to how N = 4 SYM has been used as a model for QCD [11,12].
Nevertheless, given the successes of applying results from N = 4 SYM to QCD phenomenology [35], it is possible that the O(N) model can serve as a useful phenomenological model, one which possesses a property that N = 4 SYM simply does not have: asymptotic freedom.
A. On the evidence against a Landau pole in QCD There are multiple claims in the literature that in QCD, the Landau pole in the running coupling is cured by non-perturbative effects, see e.g.Refs.[36][37][38].While some of these claims do seem not withstand rigorous scrutiny [39], others are based on careful studies of non-perturbative results in QCD, e.g. in Refs.[40,41].
As a representative example, let me consider the observation that certain quantities such as the gluon polarization tensor appear to "freeze" in the infrared for QCD [36].This is analogous to the behavior of the scattering amplitude discussed for the O(N) model in section II C. As in (49), the authors of Ref. [36] define an "effective" running coupling α s,ef f (μ) for QCD and show that it is well behaved and finite even in the zero momentum limit μ → 0.
From the O(N) model example discussed in section II C, it is clear that the approach taken in Ref. [36] is perfectly reasonable and useful as a practical definition for α s .However, from the O(N) model discussion it should also be clear that the observed infrared "freezing" of the coupling is entirely unrelated to the absence of a Landau pole.In fact, the O(N) model example shows that infrared freezing for an effective running coupling can occur in a theory with a Landau pole.To reiterate the point made in section II C, the QCD running coupling is not an observable, so it is always possible to find some "physical" definition for α s that agrees with perturbative QCD in the UV and has any set of properties in the infrared.While certain authors attribute great importance to the apparent "infrared fixed point" property based on the IR freezing of some suitably defined α s,ef f , I personally do not believe that it helps in our understanding of QCD.
Of course, it is important to state that the similarities between the O(N) model, which possesses a Landau pole, and QCD, do not constitute evidence for the presence of a Landau pole in QCD.It could simply be that observables in a theory with a "benign" Landau pole such as the O(N) model are similar to a theory without a Landau pole.
To summarize, while there is no guarantee that a Landau pole exists for QCD, I am not aware of any hard evidence against one, either.I take this to constitute evidence that at least at large N, the Landau pole in the O(N) model is harmless, and the theory does constitute a UV-complete interacting and asymptotically free theory.
Moreover, I compared results from the large N O(N) model for the running coupling and the finite temperature pressure to QCD, finding qualitative similarities between these two theories.Both the O(N) model and QCD are asymptotically free theories, but there is a key difference between these two theories.For QCD, the high energy regime possesses an intuitive interpretation as a weakly coupled equilibrium quantum field theory with welldefined particle states (quarks and gluons), whereas the infrared limit of QCD is notoriously difficult to formulate as an equilibrium continuum quantum field theory.By contrast, for the O(N) model it is the low energy regime which lends itself to an intuitive equilibrium quantum field theory formulation (with the vector and scalar bound state the perturbative degrees of freedom), whereas the high energy region is difficult to interpret.Coming full circle, while the respective "perturbative" and "non-perturbative" regimes are different for Perhaps it would be time to critically reassess the current dogma that Landau poles constitute fatal flaws of interacting quantum field theories.
FIG. 1.Running coupling in the O(N) model in 3+1 dimensions.Shown are results for the coupling λ R in the original theory (9) as well as for the coupling g R in the analytically continued theory

= 4 .
Comparison with the perturbative (weak-coupling) result for β [21, Eq.(10.54)]shows that the large N β-function is only 1/9th of the N=1 weak-coupling result.The remainder of the complete perturbative contribution arises at next-to-leading order in the large N expansion.

FIG. 2 .
FIG. 2. Pressure per component as a function of temperature for the O(N) model at large N, adapted from Ref. [9].T c ≃ 0.616Λ LP denotes the location in temperature where the solution to (17) for m becomes complex.

,
respectively.At zero temperature, where the contribution from the modified Bessel function vanishes, there is a different number of solutions depending on the sign in (42) and magnitude of M 2 .For positive sign (corresponding to solution σ (1) 0 above), and large M 2 , there are three solutions: m = 0, m ≃ Λ LP √ e and m ∝ M up to logarithmic corrections.Of these, the saddle with the largest pressure (lowest free energy) is m ≃ Λ LP √ e, hence this is the dominant saddle point at large N.One thus again recovers the solution FIG.3."Effective" running coupling (49) vs. MS coupling(14) for the PT -symmetric O(N) model at large N.The effective coupling g R,ef f is defined from the "physical" scattering amplitude M and is positive and finite for all values of μ.See text for details.
FIG. 4. Running coupling in QCD (with Λ ≡ Λ MS = 0.3 GeV) and large N PT -symmetric O(N) model.In order to have a fair comparison, I'm converting the QCD running coupling α s (μ) to the 't Hooft coupling by multiplying with a factor of N c = 3, and I am dividing the PT -symmetric O(N) model coupling by a factor of 4π.Lattice QCD point at μ = 1.5 GeV is from Ref. [28], QED

FIG. 5 .
FIG. 5. Pressure as a function of temperature for QCD (with Λ ≡ Λ MS = 0.3 GeV) and large N PT -symmetric O(N) model.The QCD pressure is from lattice QCD with N f = 2 + 1 flavors of

π 2 T 4 FIG. 6 .
FIG. 6. Speed of sound squared c 2 s (left panel) and energy density ϵ(T ) (right panel) as a function of temperature from lattice QCD [31] (with Λ ≡ Λ MS = 0.3 GeV) and large N PT -symmetric O(N) model.Unlike the pressure, both c 2 s and ϵ(T ) clearly show the second order phase transition in the O(N) model as compared to the analytic cross-over in QCD.See text for details.
In this work, I have considered the O(N) model in 3+1 dimensions at large N, which has a Landau pole.Using technology borrowed from PT -symmetric field theories, I have extended the theory beyond the Landau pole, and I have found that adding relevant and irrelevant operators do not qualitatively change the behavior of the theory close to the Landau pole.Physical observables in the O(N) model are finite and well-behaved at and close to the Landau pole.
QCD and the O(N) model, it is curious to note that both theories contain both phases, and that the Landau pole acts as a the phase boundary for the O(N) model.Based on these observations, my interpretation is as follows: Landau poles are common features for quantum field theories in 3+1 dimensions, since the O(N) model at large N, QED, and even QCD possess diverging coupling constants at a finite momentum scale μ = Λ.For two of these theories (O(N) model and QCD), we know that nothing 'bad' happens at this scale.Instead, μ = Λ merely marks the scale at which the O(N) model and QCD seem to transition from a low temperature phase to a high temperature phase.