Next Article in Journal / Special Issue
Results from the Cuore Experiment
Previous Article in Journal / Special Issue
Finite-Energy Dressed String-Inspired Dirac-Like Monopoles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Four Loop Scalar ϕ4 Theory Using the Functional Renormalization Group †

by
Margaret E. Carrington
1,2,* and
Christopher D. Phillips
1
1
Department of Physics, Brandon University, Brandon, MB R7A 6A9, Canada
2
Winnipeg Institute for Theoretical Physics, Winnipeg, MB, Canada
*
Author to whom correspondence should be addressed.
This paper is based on the talk at the 7th International Conference on New Frontiers in Physics (ICNFP 2018), Crete, Greece, 4–12 July 2018.
Submission received: 28 November 2018 / Revised: 26 December 2018 / Accepted: 28 December 2018 / Published: 2 January 2019

Abstract

:
We work with a symmetric scalar theory with quartic coupling in 4-dimensions. Using a 2PI effective theory and working at 4 loop order, we renormalize with a renormalization group method. All divergences are absorbed by one bare coupling constant and one bare mass which are introduced at the level of the Lagrangian. The method is much simpler than counterterm renormalization, and can be generalized to higher order nPI effective theories.

Graphical Abstract

1. Introduction

There are many systems of physical interest that are strongly coupled and must be described with non-perturbative methods. Schwinger-Dyson (SD) equations are often used, but one problem with this approach is that the hierarchy of coupled SD equations needs to be truncated, and different truncations have been proposed. The n-particle-irreducible effective action is an alternative non-perturbative method. The action is written as a functional of dressed vertex functions, which are calculated self consistently by applying the variational principle [1,2]. A fundamental advantage of nPI is that the method provides a systematic expansion with the truncation occurring at the level of the action. Gauge invariance may be violated by the truncation [3,4], and various proposals to minimize gauge dependence have been discussed in [5,6,7,8,9]. We are primarily interested in the renormalization of nPI theories. The 2PI effective theory can be renormalized using a counterterm approach [10,11,12,13,14], but the method requires several sets of vertex counterterms and cannot be extended to the 4PI theory. It is known that higher order nPI formulations ( n > 2 ) are necessary in some situations. Transport coefficients in gauge theories (even at leading order) cannot be calculated using a 2PI formulation [15,16], and numerical calculations have shown that, for a symmetric scalar ϕ 4 theory, 4PI vertex corrections are large in three dimensions [17], and for sufficiently large coupling the 2PI approximation breaks down at the 4 loop level in four dimensions [18,19].
In this paper we work with a symmetric scalar theory, to avoid some of the complications of gauge theories, and focus on the problem of renormalizability. We use the renormalization group (RG) method that was introduced in [20]. Using this method, no counterterms are needed and the divergences are absorbed into the bare parameters of the Lagrangian, the structure of which is fixed and totally independent of the order of the approximation. In this sense, the RG method is designed to be used at any order in the nPI approximation (and at any loop order).

2. Notation

We introduce a notation that suppresses the arguments that give the space-time dependence of functions. For example, the term in the action that is quadratic in the fields is written:
i 2 d 4 x d 4 y φ ( x ) G no · int 1 ( x y ) φ ( y ) i 2 φ G no · int 1 φ ,
where G no · int is the bare propagator. The classical action is
S [ φ ] = i 2 φ G no · int 1 φ i 4 ! λ φ 4 , i G no · int 1 = ( + m 2 ) .
For notational convenience we use a scaled coupling constant ( λ phys = i λ ), and the factor of i that is introduced here will be removed when we rotate to Euclidean space for numerical calculations.
To use the functional renormalization group method, we add a non-local regulator function to the action [21]
S κ [ φ ] = S [ φ ] + Δ S κ [ φ ] , Δ S κ [ φ ] = 1 2 φ R ^ κ φ .
The scale denoted κ has dimensions of momentum. The regulator function satisfies lim Q κ R ^ κ ( Q ) κ 2 and lim Q κ R ^ κ ( Q ) 0 so that for Q κ the regulator plays the role of a large mass term which suppresses quantum fluctuations with wavelengths 1 / Q 1 / κ , while in the opposite limit fluctuations with wavelengths 1 / Q 1 / κ are unaffected. The regulated action (3) can be used to obtain the 2PI generating functionals:
Z κ [ J , J 2 ] = [ d φ ] exp i S κ [ φ ] + J φ + 1 2 φ J 2 φ , W κ [ J , J 2 ] = i ln Z κ [ J , J 2 ] .
To obtain the 2PI effective action, we take the double Legendre transform of the generating functional W κ [ J , J 2 ] with respect to the sources J and J 2 , with ϕ and G now taken as the independent variables. The resulting effective action Γ κ [ ϕ , G ] can be written
Γ κ [ ϕ , G ] = Γ no · int · κ [ ϕ , G ] + Γ int [ ϕ , G ] Δ S κ ( ϕ ) , Γ no · int · κ [ ϕ , G ] = i 2 ϕ G no · int · κ 1 ϕ + i 2 Tr ln G 1 + i 2 Tr G no · int · κ 1 G , Γ int [ ϕ , G ] = i 4 ! λ ϕ 4 i 4 λ ϕ G ϕ + Γ 2 [ ϕ , G ; λ ] ,
where Γ 2 means the set of all 2PI graphs with two and more loops and we have defined i G no · int · κ 1 = i G no · int 1 R ^ κ = ( m 2 + R ^ k ) . We have subtracted the regulator term so that the effective action corresponds to the classical action at the ultraviolet scale μ . To simplify the notation we will write Γ = i Φ where both Γ and Φ have the same subscripts, and we define an imaginary regulator function R κ = i R ^ κ (the extra factor i will be removed when we change to Euclidean space variables).
The effective action is extremized by solving the variational equations of motion for the self consistent 1 and 2 point functions. These self consistent κ dependent solutions are denoted ϕ κ and G κ , but since we work with the symmetric theory we will set ϕ κ = 0 . We calculate n-point kernels by functionally differentiating the effective action
Λ ( m ) = 2 m δ m δ G m Φ int , Λ κ ( m ) = Λ ( m ) | G = G κ .
To simplify the notation we use special names for certain kernels: Λ ( 1 ) = Σ , Λ ( 2 ) = Λ , Λ ( 3 ) = Υ .

3. Flow Equations

Φ int and Λ ( m ) do not depend explicitly on κ and therefore we can use the chain rule to obtain
κ Λ κ ( m ) = 1 2 κ G κ Λ κ ( m + 1 ) .
In momentum space the equation becomes
κ Λ κ ( m ) ( P 1 , P 2 , P m ) = 1 2 d Q κ G κ ( Q ) Λ κ ( m + 1 ) ( P 1 , P 2 , P m + 1 , Q ) .
We will show below that this infinite hierarchy of coupled integral equations for the n-point kernels truncates at the level of the action. The flow equations can be rewritten in a more useful form using the stationary condition
δ Φ κ [ ϕ , G ] δ G | G = G κ = 0
which gives by a straightforward calculation
κ G κ = G κ ( κ G κ 1 ) G κ = G κ κ ( R κ + Σ κ ) G κ .
The first two equations in the hierarchy (8) now take the form
κ Σ κ ( P ) = 1 2 d Q κ Σ κ ( Q ) + R κ ( Q ) G κ 2 ( Q ) Λ κ ( P , Q ) ,
κ Λ κ ( P , K ) = 1 2 d Q κ R κ ( Q ) + Σ κ ( Q ) G κ 2 ( Q ) Υ κ ( P , K , Q ) .
By iterating Equation (11), we can reformulate the flow equation for the 2 point function Σ so that the kernel contains a Bethe-Salpeter (BS) vertex:
κ Σ κ ( P ) = 1 2 d Q κ R κ ( Q ) G κ 2 ( Q ) M κ ( P , Q ) ,
with
M κ ( P , K ) = Λ κ ( P , K ) + 1 2 d Q Λ κ ( P , Q ) G κ 2 ( Q ) M κ ( Q , K ) .
A different class of non-perturbative vertices can be defined by considering variations of the effective action with respect to the field. The 4 point function that is obtained in this way is related to the BS vertex as V = λ + 3 ( M Λ ) . The vertex V contains terms from all three (s, t and u) channels, and the shorthand notation which suppresses indices combines the three channels to give the factor (3) in Equation (3).
We rotate to Euclidean space for the numerical calculation, and to simplify the notation we do not introduce subscripts to denote Euclidean space quantities. The flow Equations (11) and (12) and the BS Equation (14) have the same form in Euclidean space. The Dyson equation has the form G 1 ( P ) = G no · int 1 ( P ) + Σ ( P ) , and the equation for the physical vertex in Euclidean space is V = λ + 3 ( M Λ ) . The regulator function becomes
R κ ( Q ) = Q 2 e Q 2 / κ 2 1 .
At the 4 loop level, the hierarchy of flow equations can be truncated at the level of the second equation (this is explained below). The n-point functions for the quantum theory can be obtained by starting from initial conditions defined at κ = μ and solving the integro-differential flow Equations (11) and (12). We choose the regulator function R κ so that the theory is described by the classical action at the ultraviolet scale κ = μ . The initial conditions are therefore obtained from the bare masses and couplings of the Lagrangian. The values of the bare parameters are unknown, but the values of the renormalized parameters are specified by the renormalization conditions
G 0 1 ( 0 ) = m 2 , M 0 ( 0 , 0 ) = λ
that are enforced by choice on the n-point functions that will be obtained at the quantum end of the flow. The method is to start from an initial guess for the bare parameters, solve the flow equations, extract the renormalized parameters, and then adjust the bare parameters (either up or down depending on the result). We then resolve the flow equations and repeat the procedure, continuing until the renormalization conditions are satisfied (to some numerically specified accuracy).
It can be shown [19] that consistency between the initial conditions and the renormalization conditions requires
Z = lim ( P 1 , P 2 ) 0 Λ ˜ 0 ( m ) ( P 1 , P 2 ) Λ ˜ 0 ( m ) ( 0 , 0 ) 0 .
If the hierarchy in (8) is truncated correctly, the condition (17) will be satisfied. This statement is proved by showing that if a given kernel obtained from functional differentiation satisfies the condition (17), it will also satisfy Λ 0 ( m ) ( 0 , 0 ) = λ and Λ μ ( m ) ( 0 , 0 ) = λ μ [19]. The result is that the flow equation for this kernel does not have to be solved. We therefore need to find the smallest value of m for which (17) is satisfied, and then solve self consistently the set of flow equations for the kernels with 2 × ( 1 , 2 , 3 , m 1 ) legs.
It is straightforward to show that any kernel that contains a diagram with a loop that is not forced by the structure of the diagram to carry one of the external momenta, will not satisfy (17), and the flow equation for this kernel must be solved [19]. If the effective action is truncated at the 3 loop level, the self energy will include the sunset diagram which will not satisfy (17). On the other hand, the kernel Λ has the tree graph and two 1 loop contributions that always carry external momenta, which means that Λ does not have to be flowed but can be simply substituted into the Σ flow equation. We have only to replace the tree vertex with the bare vertex ( λ μ ) to satisfy the initial condition. At the 4 loop level the kernel Λ does not satisfy (17), but the 6-leg kernel Υ does, and can be substituted directly into the Λ flow equation. There is no bare 6-vertex in the Lagrangian and therefore the integration constant is set to zero. The result is that at the 4 loop level we must solve the Σ and Λ flow equations self consistently.

4. Numerical Method

We start the flow of the 2 and 4 kernels from the initial conditions
Σ μ ( P ) = m μ 2 m 2 , Λ μ ( P , K ) = λ μ ,
and the propagator in the ultraviolet limit is G μ 1 ( P ) = P 2 + m μ 2 . We replace κ with the variable t = ln κ / μ so that we approach the quantum theory more slowly. We use κ max = μ = 100 , κ min = 10 2 and N κ = 50 and we have tested the insensitivity of our results to these choices. We have also used a generalized form of (15) to verify that our results are not dependent on the form of the regulator. The renormalized mass and coupling are obtained from the quantum functions
m found 2 = G 0 1 ( 0 ) = m 2 + Σ 0 ( 0 ) , λ found = M 0 ( 0 , 0 ) ,
and are then compared with the values specified in the renormalization conditions, adjusted, and tuned, by repeating the procedure until the renormalization conditions are satisfied to specified accuracy.
The 4-dimensional momentum integrals are written
d K f ( k 0 , k ) = n d 3 k ( 2 π ) 3 f ( m t n , k ) ,
with m t = 2 π T . There are N t terms in the summation with β = 1 T = N t a t and a t is the lattice spacing in the temporal direction. We use spherical coordinates and Gauss-Legendre integration to do the integrals over the 3-momenta.

5. Results and Discussion

We use N x = N ϕ = 8 points for the integrations over the cosine of the polar angle and the azimuthal angle, and we have checked that all results are stable when we increase the number of grid points in these dimensions. The momentum space grid spacing is Δ p 1 a s N s where a s is the spatial lattice spacing and N s is the number of lattice points for the momentum magnitude. The UV momentum cutoffs are ( p 0 ) max = π / a t and p max = π / a s . We use a t = a s = 1 / 8 so that ( p 0 ) max = p max = 8 π μ = 100 . The numerics are stable if results are unchanged when Δ p decreases while p max is held fixed, and we have checked that this is true if N s 14 . To test the renormalization we increase p max while holding Δ p 1 / L fixed. In Figure 1 we show V ( 0 ) versus p max . For purposes of comparison we also show a calculation that is done incorrectly, by working at 3 loop level and replacing one of the vertices in the 4 kernel with a bare vertex. We have checked that dependence on the renormalization scale is very small.
To evaluate the 2, 3 and 4 loop approximations in the context of a physical quantity, we calculate the pressure, which can be obtained from the effective action using P = T V Φ where V is the 3-volume. To 4 loop order, the contributions to the pressure are
P 0 = 1 2 d Q ln G no · int 1 ( Q ) π 2 T 4 90 P 1 = 1 2 d Q ln G 1 ( Q ) 1 Q 2 + m 2 P 2 = 1 2 d Q ( Q 2 + m b 2 ) G ( Q ) 1 P 3 = 1 8 λ b d Q G ( Q ) d L G ( L ) P 4 = 1 48 λ ( λ 2 λ b ) d P d K d Q G ( P ) G ( K ) G ( Q ) G ( P + K + Q ) P 5 = 1 48 λ 3 d Q d S G ( S ) G ( S + Q ) d L G ( L ) G ( L + Q ) d M G ( M ) G ( M + Q )
P sum = P 0 + P 1 + P 2 + P 3 + P 4 + P 5 .
There is a temperature independent divergence that can be subtracted off with a ‘cosmological constant’ renormalization, by setting the vacuum pressure to zero: Δ P = P sum P sum ( T = 0 ) . The arrow on the right side of (21) indicates that we have dropped a temperature independent constant which would have been removed by this shift. The term P 0 is the non-interacting ( λ = 0 ) pressure and since we want to compare Δ P to the non-interacting expression, we define P = Δ P P 0 . In Figure 2 we show our results for the pressure as a function of the coupling at the 2, 3 and 4 loop orders of approximation.

6. Conclusions

In this paper we present results from a 4 loop 2PI calculation in a symmetric ϕ 4 theory with the renormalization done using the RG method of [20]. No counterterms are introduced, and all divergences are absorbed into the bare parameters of the Lagrangian, the structure of which is fixed and independent of the order of the approximation. Our main goal is to use our method to do a calculation with the 4PI effective theory. The basic method is the same, since the form of the flow and Bethe-Salpeter equations are similar [22], but at the 4PI level we must introduce a flow equation for the variational 4 vertex. This calculation is currently in progress.

Author Contributions

Both authors have contributed equally to this work.

Funding

This work has been supported by the Natural Sciences and Engineering Research Council of Canada.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cornwall, J.M.; Jackiw, R.; Tomboulis, E. Effective action for composite operators. Phys. Rev. D 1974, 10, 2428–2445. [Google Scholar] [CrossRef]
  2. Norton, R.E.; Cornwall, J.M. On the formalism of relativistic many-body theory. Ann. Phys. 1975, 91, 106. [Google Scholar] [CrossRef]
  3. Arrizabalaga, A.; Smit, J. Gauge-fixing dependence of Φ-derivable approximations. Phys. Rev. D 2002, 66, 065014. [Google Scholar] [CrossRef]
  4. Carrington, M.E.; Kunstatter, G.; Zaraket, H. 2PI effective action and gauge invariance problems. Eur. Phys. J. C 2005, 42, 253–259. [Google Scholar] [CrossRef]
  5. Pilaftsis, A.; Teresi, D. Exact RG invariance and symmetry improved 2PI effective potential. Nucl. Phys. B 2017, 920, 298–318. [Google Scholar] [CrossRef]
  6. Pilaftsis, A.; Teresi, D. Symmetry-improved 2PI approach to the Goldstone-boson IR problem of the SM effective potential. Nucl. Phys. B 2016, 906, 381–407. [Google Scholar] [CrossRef] [Green Version]
  7. Markó, G.; Reinosa, U.; Szép, Z. Loss of solution in the symmetry improved Φ-derivable expansion scheme. Nucl. Phys. B 2016, 913, 405–424. [Google Scholar] [CrossRef] [Green Version]
  8. Markó, G.; Reinosa, U.; Szép, Z. O(N) model within the ϕ-derivable expansion to order λ2: On the existence and UV/IR sensitivity of the solutions to self-consistent equations. Phys. Rev. D 2015, 92, 125035. [Google Scholar] [CrossRef]
  9. Brown, M.J.; Whittingham, I.B. Soft symmetry improvement of two particle irreducible effective actions. Phys. Rev. D 2017, 95, 025018. [Google Scholar] [CrossRef] [Green Version]
  10. Van Hees, H.; Knoll, J. Renormalization of self-consistent approximation schemes at finite temperature. II. Applications to the sunset diagram. Phys. Rev. D 2002, 65, 105005. [Google Scholar] [CrossRef]
  11. Van Hees, H.; Knoll, J. Renormalization in self-consistent approximation schemes at finite temperature: Theory. Phys. Rev. D 2001, 65, 025010. [Google Scholar] [CrossRef]
  12. Blaizot, J.-P.; Iancu, E.; Reinosa, U. Renormalization of Φ-derivable approximations in scalar field theories. Nucl. Phys. A 2004, 736, 149–200. [Google Scholar] [CrossRef] [Green Version]
  13. Berges, J.; Borsányi, S.; Reinosa, U.; Serreau, J. Nonperturbative renormalization for 2PI effective action techniques. Ann. Phys. 2005, 320, 344–398. [Google Scholar] [CrossRef] [Green Version]
  14. Reinosa, U.; Serreau, J. 2PI functional techniques for gauge theories: QED. Ann. Phys. 2010, 325, 969–1017. [Google Scholar] [CrossRef] [Green Version]
  15. Carrington, M.E.; Kovalchuk, E. Leading order QED electrical conductivity from the three-particle irreducible effective action. Phys. Rev. D 2008, 77, 025015. [Google Scholar] [CrossRef]
  16. Carrington, M.E.; Kovalchuk, E. Leading order QCD shear viscosity from the three-particle irreducible effective action. Phys. Rev. D 2009, 80, 085013. [Google Scholar] [CrossRef]
  17. Carrington, M.E.; Fu, W.-J.; Mikula, P.; Pickering, D. Four-point vertices from the 2PI and 4PI effective actions. Phys. Rev. D 2014, 89, 025013. [Google Scholar] [CrossRef]
  18. Carrington, M.E.; Meggison, B.A.; Pickering, D. 2PI effective action at four loop order in ϕ4 theory. Phys. Rev. D 2016, 94, 025018. [Google Scholar] [CrossRef]
  19. Carrington, M.E.; Friesen, S.A.; Meggison, B.A.; Phillips, C.D.; Pickering, D.; Sohrabi, K. 2PI effective theory at next-to-leading order using the functional renormalization group. Phys. Rev. D 2018, 97, 036005. [Google Scholar] [CrossRef]
  20. Carrington, M.E.; Fu, W.-J.; Pickering, D.; Pulver, J.W. Renormalization group methods and the 2PI effective action. Phys. Rev. D 2015, 91, 025003. [Google Scholar] [CrossRef]
  21. Wetterich, C. Exact evolution equation for the effective potential. Phys. Lett. B 1993, 301, 90–94. [Google Scholar] [CrossRef] [Green Version]
  22. Carrington, M.E.; Fu, W.-J.; Fugleberg, T.; Pickering, D.; Russell, I. Bethe-Salpeter equations from the 4PI effective action. Phys. Rev. D 2013, 88, 085024. [Google Scholar] [CrossRef]
Figure 1. The physical vertex V ( 0 ) versus p max with λ = 2 , T = 2 and L = 4 at 4 loop level in the skeleton expansion. To set the scale we also show the results of an incorrect calculation (see text for more explanation).
Figure 1. The physical vertex V ( 0 ) versus p max with λ = 2 , T = 2 and L = 4 at 4 loop level in the skeleton expansion. To set the scale we also show the results of an incorrect calculation (see text for more explanation).
Universe 05 00009 g001
Figure 2. The pressure as a function of coupling. The right panel shows a close up of the large coupling region where the three approximations start to diverge from each other.
Figure 2. The pressure as a function of coupling. The right panel shows a close up of the large coupling region where the three approximations start to diverge from each other.
Universe 05 00009 g002

Share and Cite

MDPI and ACS Style

Carrington, M.E.; Phillips, C.D. Four Loop Scalar ϕ4 Theory Using the Functional Renormalization Group. Universe 2019, 5, 9. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5010009

AMA Style

Carrington ME, Phillips CD. Four Loop Scalar ϕ4 Theory Using the Functional Renormalization Group. Universe. 2019; 5(1):9. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5010009

Chicago/Turabian Style

Carrington, Margaret E., and Christopher D. Phillips. 2019. "Four Loop Scalar ϕ4 Theory Using the Functional Renormalization Group" Universe 5, no. 1: 9. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5010009

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop