An Exact Integration of A\Phi^ 4 Quantum Field Theory
Total Page:16
File Type:pdf, Size:1020Kb
An exact integration of a φ4 quantum field theory Timothy D. Andersen (Dated: Received: date / Accepted: date) Abstract Most quantum field theories are not exactly solvable. In this paper show the statistical equiva- lence of the standard exponential path integral to products of Heaviside functions, i.e. a product of specially tuned uniform distributions. This allows exact integrations of certain quantum field theories. I apply the equivalence to calculate the exact, non-perturbative path integral for a 3+1-D scalar (real) phi-4 field theory. PACS numbers: 03.70.+k, 05.20.Gg arXiv:1306.6836v1 [physics.gen-ph] 27 Jun 2013 1 The standard or “exponential” path-integral formulation of quantum field theory is the functional integral over a Hilbert space, i = Dφe ~ S[φ], Z Z where S[φ] = d4x , the action functional in the quantum field φ, is only exactly com- L putable for a smallR number of Lagrangians such as the free field, [φ]= ∂µφ∂ φ m2φ2 + Jφ, Lfree µ − for a scalar field φ and source J both in L2 Hilbert space over Minkowski spacetime [13]. The massive, real scalar φ4 theory is one of the simplest theories in quantum field theory without an exact integration. Its Lagrangian is given by, [φ]= ∂µφ∂ φ m2φ2 λφ4 + Jφ, L µ − − where λ > 0. There are no known mathematical techniques for calculating for the φ4 Z theory. The main analytical approaches are Feynman diagrams, which give perturbations of the free particle ensemble for small coupling constant, λ, or perturbation series for strong coupling and small kinetic term[7]. The former perturbation series, however, do not converge for the quartic interaction theory and are asymptotic at best, while the latter have limited applicability in high energy physics. What we need for strong couplings and high energy is an exact solution. Under a Wick rotation, the 3 + 1-D φ4 quantum field theory becomes a 4-D statistical mechanical theory [11], − 1 (J)= Dφe ~ S[φ], Z Z where S[φ]= d4x and [φ]= ∂µφ∂ φ+m2φ2 +λφ4 Jφ. The functional is now a 4-D H H µ − H energy functional,R the amplitude is a partition function, and Planck’s reduced constant ~ Z is equivalent to temperature. The path integral is equivalent to the “canonical” ensemble. For experiments, we are only interested in how the quantization changes relative to the vacuum. Therefore, the quantity of interest is the ratio: = Z(J) . Most expectations of A Z(0) observables can be found from this ratio or a scaling of it. Perturbation theory, specifically Feynman diagrams, is a time-tested and useful technique for evaluating this ratio at high energy, but it is also tedious and tends to blow up if the coupling is too strong. Therefore, there is a strong motivation in quantum field theory 2 to calculate the ratio exactly for as many theories as possible, even ones as simple as the φ4, because of the possibility that it will lead to the development of new forms of non- perturbative calculus for making quantum predictions in QED and QCD, useful at high energies and densities where couplings become too strong. With this motivation in mind, I apply the principle of ensemble equivalence to derive an exact solution. The principle of ensemble equivalence is a common tool for simplifying computations in both statistical mechanics and quantum field theory. Let Ω (J)= dφ dφ δ(A σ [φ, J]), N 1 ··· N N − N Z be the N-dimensional microcanonical ensemble and − ~ Z (J)= dφ dφ e σN [φ]/ , N 1 ··· N Z the N-dimensional canonical ensemble, with σN the N dimensional action. By definition, two statistical ensembles are equivalent if, in the “thermodynamic” limit of infinite degrees of freedom, N , they generate identical expectations of observables. Equivalence is → ∞ typically established by an asymptotic relation (e.g., with the method of steepest-descent or saddle point method [5]) that becomes exact in the limit. The equivalence between microcanonical and canonical ensembles has been known since at least the 1930’s [10] and criteria for non-equivalence established in the 1970’s [9][12][8]. For efficiency, lattice gauge simulations frequently implement the microcanonical quantum field theory either in “demon” Monte Carlo as in the work of Creutz et al. [3][4] or the Hamiltonian flow method of Callaway et al. [1][2]. Perturbation theories for microcanonical quantum field theory have also been established [11][6]. What is relatively unknown is that the microcanonical ensemble is not the only one for which an equivalence to the canonical can be shown. As I show below, a specially tuned product of uniform distributions, i.e. an ensemble density given by a product of step or Heaviside functions, can also be shown to be equivalent to the standard quantization by invoking the multidimensional method of steepest descent. In this paper, I show the following exact equivalence: Z(J) d4k C = exp F [J] (1) Z(0) (2π)4 Z where 3 √ 2/3 2 3( h + g ) −1 1/3 1/3 F [J] = J − + λp λp (4α + h/g + g ) wg1/3 λ−1( 2α + h/g1/3 + g1/3) − × p − h q q g − 108 + w 1√3( h + g2/3)(h +4αg1/3 + g2/3) − h 2αg1/3 + g2/3 − × − i −1 1/3 2/3 −2 −1 1/3 1/3 −1 1/3 1/3 (h +4αg + g ) 8 λp ( 2α + h/g + g )+ λp (4α + h/g + g ) (2), ! × − − ! hq q i given that g α3 + 36Aαλ +6√3w, h α2 12Aλ , ≡ p ≡ − p 1 w Aλ (α2 +4Aλ ), α (k2 + m2). ≡ p p ≡ 2 The value for A is given by thep implicit equation, 1 = √3α4 + 72√3Aα2λ + 144√3A2λ2 + 24α3 Aλ + 288α(Aλ )3/2 (√3h + 12α√Aλ)g2/3 ~ p p p p − × h i λ−1( 2α + h/g1/3 + g1/3) λ−1(4αp+ h/g1/3 + g1/3) p − − − p × hq q i √ 4/3 −1 1/3 1/3 −1 1/3 1/3 2 Aλg λp ( 2α + h/g + g ) λp (4α + h/g + g ) " − − × q q −1 −1 1/3 1/3 −1 1/3 1/3 λp ( 2α + h/g + g )+ λp (4α + h/g + g ) , (3) − − # q q and λp is the renormalized or “physical” coupling constant. The proportionality constant C is a trivial scaling. N Let S[φ, J] = limN→∞ σN where σN = ǫ j=1 Lj[φj,Jj] and, working in momentum ∆k 4 space, ǫ = 1/N is the hypercubic volumeP element of the 4-D momentum space 2π ∼ which, along with N, is a regularization. Lemma 1. Given the limit exists, the canonical ensemble and an ensemble defined by the product of Heaviside functions are equivalent for particular Aj, i.e., Z[J] N dφ θ(A ǫL [φ ,J ]) C = lim j j − j j j , Z[0] N→∞ dφ θ(A ǫL [φ , 0]) j=1 R j j j j Y − R where C is a trivial proportionality constant. The constant is trivial if and only if, given a 4 Hilbert space functional O = lim →∞ O such that O O(1), N N N ∼ dφ dφ O ρ (φ ,...,φ ) O = lim 1 ··· N N θ 1 N h i N→∞ dφ dφ ρ (φ ,...,φ ) R 1 ··· N θ 1 N dφ dφ O ρ (φ ,...,φ ) = lim R 1 ··· N N exp 1 N , (4) N→∞ dφ dφ ρ (φ ,...,φ ) R 1 ··· N exp 1 N − ~ for ρ = θ(A ǫL ) and ρ = R e ǫLj / for some choice of A . θ j j − j exp j j Proof. TheQ Heaviside has numerousQ representations, but the best one to use here is, x x γ+i∞ dη −η(x0−t) θ(x x0)= dt δ(t x0)= dt e , − − − ∞ 2πi Z0 Z0 Zγ i where we have used the inverse Laplacian representation for the Dirac delta function [5]. Now, given that dφ θ(A ǫL [φ ,J ]) = j j − j j j Z Aj γ+i∞ dηj ηj (aj −ǫLj [φj ,Jj ]) dφj daj e , − ∞ 2πi Z Z0 Zγ i we can move the integrals over φj inside: Aj γ+i∞ dηj ηj aj −ηj ǫLj [φj ,Jj ]) daj e dφje = − ∞ 2πi j 0 γ i Y Z Z Z Aj γ+i∞ dηj ηj aj −ηj Fj daj e e , (5) − ∞ 2πi j 0 γ i Y Z Z −ηj Fj −ηj ǫLj [φj ,Jj] where e = dφje . Note that Fj is a function of ηj. The methodR of steepest descent applies here because the “entropy” expression increases with N, S = η (a F ) O(N). Thus, the exponent increases with N, and the integral N j j j − j ∼ can be expandedP asymptotically in N. The N-vector saddle-point η0, where ∂/∂ηSN = 0 is the maximal point of the entropy. As N , we obtain the first order asymptotic relation, →∞ γ+i∞ γ+i∞ dη1 dηN SN [η] SN [η0] e CN e , (6) − ∞ 2πi ··· − ∞ 2πi ≈ Zγ i Zγ i where CN is related to the Hessian matrix of SN /N, − 2π N/2 ∂2S 1/2 C = det N −1 N . N N − ∂η η i j Since S [η ]= ηT a ηT F , normalizing, N 0 0 − 0 A ηT a−ηT F [J] −ηT F da da C [J]e 0 0 C [J]e 0 [J] 0 1 ··· N N = N . A ηT a−ηT F [0] −ηT F [0] da da C [0]e 0 0 CN [0]e 0 R 0 1 ··· N N R 5 With all ensemble equivalence proofs, trivial scaling factors vanish in normalization (see pg.