Continuum limit of the Green function in scaled affine $\varphi^4_4$ quantum Euclidean covariant relativistic field theory

We prove through path integral Monte Carlo computer experiments that the affine quantization of the $\varphi_4^4$ scaled Euclidean covariant relativistic scalar field theory is a valid quantum field theory with a well defined continuum limit of the one- and two-point-function. Affine quantization leads to a completely satisfactory quantization of field theories using situations that involve scaled behavior leading to an unexpected, $\hbar^2/\varphi^2$ which arises only in the quantum aspects.

The aim of the present work is to show that as we approach the continuum on the computer, the one-and two-point-function converge to well defined results.In other words we prove the validity of the continuum limit for the field theory quantized through affine quantization.
Our result could become important in the physics of the standard model where the longstanding problem of the triviality of canonical quantum ϕ 4 theory is crucial for particle physics since it undermines the Higgs mechanism.It is also very important for progresses in quantum gravity where the role of the field is played by the metric tensor which must be positive definite [20].
We then introduce the classical Hamiltonian expressed in affine variables.This leads us to where r is a positive, even, integer and g ≥ 0 is the bare coupling constant such that for g → 0 we fall into the free field theory.With these variables we do not let ϕ(x) = ∞ otherwise ϕ(x) −2 = 0 which is not fair to κ(x) and, as we already observed, we must forbid also ϕ(x) = 0 which would admit ϕ(x) −2 = ∞ giving again an undetermined kinetic term.
Accordingly for the self adjoint kinetic term κ(x) φ(x (note that the factor 3/4 that holds for ϕ > 0 should be replaced by a factor 2 if |ϕ| > 0 [17]) and one finds for the quantum Hamiltonian operator The affine action is found adding time, x 0 = ct, where c is the speed of light constant and t is the Euclidean imaginary time, so that S = β 0 H dx 0 , with H the semi-classical Hamiltonian corresponding to the one of Eq. (2.2), will then read (see the appendix) where with an abuse of notation we here use x for (x 0 , x 1 , x 2 , . . ., x s ) and β = 1/k B T , with k B the Boltzmann's constant and T the absolute temperature.In this work we will set The vacuum expectation value of an observable O[ϕ] will then be given by the following expression where the functional integrals will be calculated on a lattice using the PIMC method as explained later on.
The theory considers a real scalar field ϕ taking the value ϕ(x) on each site x of a periodic n-dimensional lattice, with n = s + 1 space-time dimensions, of lattice spacing a, the ultraviolet cutoff, spacial periodicity L = Na and temporal periodicity β = Na.The field path is a closed loop on an n-dimensional closed surface of an (n + 1)-dimensional β-periodic cylinder of radius L: an (n + 1)-dimensional torus.We used a lattice formulation of the AQ field theory of Eq. (2.3) (also studied in Eq. ( 8) of [12]) using additionally the scaling ϕ → a −s/2 ϕ and g → a s(r−2)/2 g which is necessary to eliminate the Dirac delta factor δ(0) = a −1 divergent in the continuum limit a → 0. The affine action for the field (in the primitive approximation [23]) is then approximated by where e µ is a vector of length a in the +µ direction with µ = 0, 1, 2, . . ., s.We will have In this work we are interested in reaching the continuum limit by taking Na fixed and letting N → ∞ at fixed volume L s .
We performed a PIMC [23][24][25][26] calculation for the AQ field theory described by the action of Eq. (2.5) in natural Planck units c = h = k B = 1.Specifically we studied the s = 3 and r = 4 case.We calculated the renormalized coupling constant g R and mass m R defined in Eqs. ( 11) and ( 13) of [12] respectively, measuring them in the PIMC through vacuum expectation values like in Eq. (2.4).
In particular: and at zero momentum where φ(p) = d n x e ip•x ϕ(x) is the Fourier transform of the field and we choose the 4momentum p 0 with one spacial component equal to 2π/Na and all other components equal to zero.
We also calulated the one-, two-point-, and two-point-connected-function, respectively given by (2.10) By construction, these are periodic functions, D(z) = D(z + L), of period L.Moreover, since the action S contains only even powers of the field these functions must be symmetric

III. THE SCALING
As we have seen we decided to work with a scaled field ϕ ′ (x), related to the variable ϕ(x), used for example in [12], by In other words, we are renormalizing the bare field.This can be compared with the standard renormalization formula ϕ ren (x) is referred to as the renormalized field and Z is called the renormalization constant.
In this language, we are setting Z = a −3 .
At the same time, we are rescaling the coupling constant with In the Standard Model, the various coupling constants also need to be renormalized for the continuum limit to exist, but the renormalization is not simply given by a power of the lattice spacing.Instead, it needs to carefully be tuned to the cutoff and to the couplings.In perturbation theory of canonical ϕ 4 , the bare coupling constant can be expressed in terms of the renormalized one, order by order.The result consists of a series that starts with g ren : The standard renormalization procedure is based on the fact that the Fourier transform of the renormalized two-point-function contains a pole at p 2 = M 2 , where M is the physical mass of the particle.The renormalization constant Z is chosen such that the residue of this pole is equal to 1.This ensures, in particular, that ϕ ren (x) and ϕ(x) as well as g ren and g have the same dimension.Note that our rescaling (3.1) and (3.3) instead changes the dimension of these objects.
We will soon see, in our first case study below, that the expectation value ϕ ′ (x) tends to a constant when N becomes large.This means that the expectation value of the unscaled field, ϕ(x) , tends to infinity in proportion to N 3/2 [12,14].
As we are holding g ′ constant, the unscaled coupling constant g tends to zero in proportion to 1/N 3 .This suggests that, for the parameter values we consider, the connected Green's functions of the unscaled model tend to those of a free scalar field.

IV. NUMERICAL RESULTS
Our PIMC simulations use the Metropolis algorithm [24,25] to calculate the ensemble average of Eq. (2.4) which is a N n multidimensional integral.The simulation is started from the initial condition ϕ(x) = ǫ > 0 for all lattice points x, with ǫ a small positive number.
One PIMC step consisted in a random displacement of each one of the N n field values, ϕ(x), as follows where η is a uniform pseudo random number in [0, 1] and δ is the amplitude of the dis- Notice that the minima of the two symmetric potential wells in the semi-classical Hamiltonian density described by the function f (ϕ) = 1 2 ϕ 2 + ϕ 4 + 3 8 ϕ −2 , are at ϕ ± = ±2 −1/2 ≈ ±0.707107.From our Monte Carlo simulations (see Teble I), it seems that the vacuum expectation value of the field (one-point-value), V = x ϕ(x) /N n , tends to these values in the continuum limit, a = 1/N → 0. Note that in some of our previous works [16,18,19] where, instead of keeping the bare mass m constant, we tuned it so to have a constant renormalized mass m R we found V = 0 in all cases.This is due to the fact that as N increases so does the necessary bare mass which keeps constant m R .So that the two symmetric potential wells in the semi-classical Hamiltonian density has minima that tends to zero and one experiences tunneling of the potential barrier at ϕ = 0.
In Table II we show the values for the renormalized mass, m R , and coupling constant g R at increasing values of N = 1/a.We see that in the continuum limit lim a→0 m R = 0 and lim a→0 g R = 2, meaning that lim a→0 φ(0) 4 / φ(0) 2 2 = 1  In Figure 1 we show D(z) at increasing values of N = 1/a.From the plot of the simulation data we see that the function is symmetric respect to z = 1/2 as expected, since the action only contains even powers of the field.In Figure 2  Alternatively we could have adjusted, at each change of N, the value of the bare mass m so to have a fixed value for the renormalized mass m R .This would have resulted in a convergence towards a unique two-point-connected-function in the continuum limit N → ∞.
We did not choose did strategy because it is numerically problematic to tune the bare mass so to have at each N the same value for the renormalized mass.This was what we did in some of our previous papers [12,13,[16][17][18][19][20][21][22].As explained in those works keeping the renormalized mass constant is extremely cumbersome due to the unavoidable sistematic numerical error that one faces.It would have been then extremely difficult to obtain a reasonable comparison between the two point function at different N.  Notice that the minima of the two symmetric potential wells in the semi-classical Hamiltonian density described by the function f So with our choice of m = √ N we will find, in the continuum limit, lim N →∞ ϕ = lim N →∞ ϕ + ( √ N ) = 0 in agreement with the results in Refs.[12,13,16,[18][19][20][21][22].
From Fig. 3 we see the continuum limit N → ∞, of the scaled two-point-connectedfunction where, with an abuse of notation, we dropped from D c , the prime, adopted rigorously in Section III.Respect to the work [14] (see Fig. 3 there) which dealt with the unscaled free field case and D c (0) was found to increase with increasing N we see how now the scaling has the effect of letting the value of D ′ c (0) decrease with increasing N, as shown in Fig. 3, since ϕ ′ = ϕ/N 3/2 and D ′ c = D c /N 3 .It is only tuning the bare mass m so to have a constant renormalized mass m R for each N, that we would find true convergence.
Unfortunately this procedure is not easily accomplished numerically since for each N we would have to make several test runs with different values of m in order to find the value which keeps m R approximately constant.This procedure was nonetheless carried out in the following works [12,13,16,[18][19][20][21][22].

V. CONCLUSIONS
In this paper, we represent π(x) by k(x)/ϕ(x).To insure proper values for π(x) it is necessary to restrict 0 < ϕ(x) < ∞ as well as 0 ≤ |k(x)| < ∞.Indeed such symbol change is able to treat Hamiltonian densities with an interaction ϕ(x) 4 .This leads to a completely satisfactory quantization of field theories using situations that involve scaled behavior leading to an unexpected, h2 /ϕ(x) 2 which arises only in the quantum aspects.Indeed, it is fair to claim that this symbol change leads to valid field theory quantizations.
Respect to the work [14] which dealt with the free field case we here repeat that analysis but now for the ϕ 4 interacting case.
We prove through path integral Monte Carlo computer experiments that the affine quantization of the ϕ 4 4 scaled Euclidean covariant relativistic field theory is a well defined quantum field theory with a well defined continuum limit of the one-and two-point-function, the Green's function.
The simple pseudo-potential ∝ h2 /ϕ 2 stemming from the affine quantization procedure placement.Each one of these N n moves is accepted if exp(−∆S) > η where ∆S is the change in the action due to the move (it can be efficiently calculated considering how the kinetic part and the potential part change by the displacement of a single ϕ(x)) and rejected otherwise.The amplitude δ is chosen in such a way to have acceptance ratios as close as possible to 1/2 and is kept constant during the evolution of the simulation.One simulation consisted of M PIMC steps.The statistical error on the average O will then depend on the correlation time necessary to decorrelate the property O, τ O , and will be determined asτ O σ 2 O /(MN n ),where σ 2 O is the intrinsic variance for O.We used up to a lattice of N n = 25 4 = 390625 points (N = 25) and up to M = 2 × 10 6 corresponding to MN n PIMC displacement moves.A. First case study In our simulation we first chose the following study case m = g = L = β = 1 and ǫ = 10 −10 .

B.
Second case studyFor the parameter values we just used, the box plays a crucial role: the bare Compton wavelength (1/m) is equal to the size L of the box.In order for the box to be a purely technical device introduced to regularize the theory, it must be large compared to the correlation length of the model.At the same time, the lattice spacing must be small compared to it:1/L ≪ m ≪ 1/a.(4.2) Therefore, next we considered the study case with g = L = β = 1, m = √ N /L and ǫ = 10 −10 , which should be much less affected by the presence of the box than the previous choice m = 1/L.