# Quartic interaction

This article refers to a type of self-interaction in scalar field theory, a topic in quantum field theory. Other types of quartic interactions may be found under the topic of four-fermion interactions. A classical free scalar field φ satisfies the Klein–Gordon equation. If a scalar field is denoted φ, a quartic interaction is represented by adding a potential term $\frac{\lambda}{4!} \phi^4$. The coupling constant λ is dimensionless in 4-dimensional space-time.

## The Lagrangian

The Lagrangian for a real scalar field takes the form

$\mathcal{L}=\frac{1}{2}\partial^\mu \phi \partial_\mu \phi -\frac{m^2}{2}\phi^2 -\frac{\lambda}{4}\phi^4.$

This Lagrangian has a global Z2 symmetry mapping φ to −φ. For a complex scalar field the Lagrangian is,

$\mathcal{L}=\partial^\mu \phi^* \partial_\mu \phi -m^2 \phi^* \phi -\lambda (\phi^* \phi)^2.$

With N real scalar fields, we can have a φ4 model with a global SO(N) symmetry

$\mathcal{L}=\frac{1}{2}\partial^\mu \phi_a \partial_\mu \phi_a - \frac{m^2}{2}\phi_a \phi_a -\frac{\lambda}{4}(\phi_a \phi_a)^2.$

Expanding the complex field in real and imaginary parts shows that it is equivalent to the SO(2) model of real scalar fields.

In all of the models above, the coupling constant λ must be positive, since, otherwise, the potential would be unbounded below, and there would be no stable vacuum. Also, the Feynman path integral discussed below would be ill-defined. In 4 dimensions, φ4 theories have a Landau pole. This means that without a cut-off on the high-energy scale, renormalization would render the theory trivial.

## Feynman Integral Quantization

The Feynman diagram expansion may be obtained also from the Feynman path integral formulation.[1] The time ordered vacuum expectation values of polynomials in φ, known as the n-particle Green's functions, are constructed by integrating over all possible fields, normalized by the vacuum expectation value with no external fields,

$\langle\Omega|\mathcal{T}\{{\phi}(x_1)\cdots {\phi}(x_n)\}|\Omega\rangle=\frac{\int \mathcal{D}\phi \phi(x_1)\cdots \phi(x_n) e^{i\int d^4x \left({1\over 2}\partial^\mu \phi \partial_\mu \phi -{m^2 \over 2}\phi^2-{\lambda\over 4!}\phi^4\right)}}{\int \mathcal{D}\phi e^{i\int d^4x \left({1\over 2}\partial^\mu \phi \partial_\mu \phi -{m^2 \over 2}\phi^2-{\lambda\over 4!}\phi^4\right)}}.$

All of these Green's functions may be obtained by expanding the exponential in J(x)φ(x) in the generating function

$Z[J] =\int \mathcal{D}\phi e^{i\int d^4x \left({1\over 2}\partial^\mu \phi \partial_\mu \phi -{m^2 \over 2}\phi^2-{\lambda\over 4!}\phi^4+J\phi\right)} = Z[0] \sum_{n=0}^{\infty} \frac{1}{n!} \langle\Omega|\mathcal{T}\{{\phi}(x_1)\cdots {\phi}(x_n)\}|\Omega\rangle.$

A Wick rotation may be applied to make time imaginary. Changing the signature to (++++) then gives a φ4 statistical mechanics integral over a 4-dimensional Euclidean space,

$Z[J]=\int \mathcal{D}\phi e^{-\int d^4x \left({1\over 2}(\nabla\phi)^2+{m^2 \over 2}\phi^2+{\lambda\over 4!}\phi^4+J\phi\right)}.$

Normally, this is applied to the scattering of particles with fixed momenta, in which case, a Fourier transform is useful, giving instead

$\tilde{Z}[\tilde{J}]=\int \mathcal{D}\tilde\phi e^{-\int d^4p \left({1\over 2}(p^2+m^2)\tilde\phi^2-\tilde{J}\tilde\phi+{\lambda\over 4!}{\int d^4p_1d^4p_2d^4p_3\delta(p-p_1-p_2-p_3)\tilde\phi(p)\tilde\phi(p_1)\tilde\phi(p_2)\tilde\phi(p_3)}\right)}.$

where $\delta(x)$ is the Dirac delta function.

The standard trick to evaluate this functional integral is to write it as a product of exponential factors, schematically,

$\tilde{Z}[\tilde{J}]=\int \mathcal{D}\tilde\phi \prod_p \left[e^{-(p^2+m^2)\tilde\phi^2/2} e^{-\lambda/4!\int d^4p_1d^4p_2d^4p_3\delta(p-p_1-p_2-p_3)\tilde\phi(p)\tilde\phi(p_1)\tilde\phi(p_2)\tilde\phi(p_3)} e^{\tilde{J}\tilde\phi}\right].$

The second two exponential factors can be expanded as power series, and the combinatorics of this expansion can be represented graphically. The integral with λ = 0 can be treated as a product of infinitely many elementary Gaussian integrals, and the result may be expressed as a sum of Feynman diagrams, calculated using the following Feynman rules:

• Each field $\tilde{\phi}(p)$ in the n-point Euclidean Green's function is represented by an external line (half-edge) in the graph, and associated with momentum p.
• Each vertex is represented by a factor .
• At a given order λk, all diagrams with n external lines and k vertices are constructed such that the momenta flowing into each vertex is zero. Each internal line is represented by a factor 1/(q2 + m2), where q is the momentum flowing through that line.
• Any unconstrained momenta are integrated over all values.
• The result is divided by a symmetry factor, which is the number of ways the lines and vertices of the graph can be rearranged without changing its connectivity.
• Do not include graphs containing "vacuum bubbles", connected subgraphs with no external lines.

The last rule takes into account the effect of dividing by $\tilde{Z}[0]$. The Minkowski-space Feynman rules are similar, except that each vertex is represented by $-i\lambda$, while each internal line is represented by a factor i/(q2-m2 + i ε), where the ε term represents the small Wick rotation needed to make the Minkowski-space Gaussian integral converge.

## Renormalization

Main article: renormalization

The integrals over unconstrained momenta, called "loop integrals", in the Feynman graphs typically diverge. This is normally handled by renormalization, which is a procedure of adding divergent counter-terms to the Lagrangian in such a way that the diagrams constructed from the original Lagrangian and counterterms are finite.[2] A renormalization scale must be introduced in the process, and the coupling constant and mass become dependent upon it. It is this dependence that leads to the Landau pole mentioned earlier, and requires that the cutoff be kept finite. Alternatively, if the cutoff is allowed to go to infinity, the Landau pole can be avoided only if the renormalized coupling runs to zero, rendering the theory trivial.

## Spontaneous symmetry breaking

An interesting feature can occur if m2 happens to be negative, but with λ positive. In this case, the vacuum consists of two lowest-energy states, each of which spontaneously breaks the Z2 global symmetry of the original theory. This leads to the appearance of interesting collective states like domain walls. In the O(2) theory, the vacua would lie on a circle, and the choice of one would spontaneously break the O(2) theory. A continuous broken symmetry leads to a Goldstone boson. This type of spontaneous symmetry breaking is the essential component of the Higgs mechanism.[3]