Prerpints.org logo
Preprint
Communication

Continuum Limit of the Green Function in Scaled Affine φ4 4 Quantum Euclidean Covariant Relativistic Field Theory

Altmetrics

Downloads

105

Views

26

Comments

1

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

28 December 2023

Posted:

28 December 2023

You are already at the latest version

Alerts
Abstract
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.
Keywords: 
Subject: Physical Sciences  -   Theoretical Physics

1. Introduction

It is well known that φ 4 4 quantum Euclidean covariant relativistic field theory when quantized through canonical (Dirac [1]) quantization (CQ) is non-renormalizable since its corresponding renormalized theory tends to a free theory in the continuum limit [2,3,4,5,6].
Recently J. R. Klauder [7,8,9,10,11] noticed that this difficulty can be overcome by using a different kind of quantization method, namely affine quantization (AQ).
In a sequel of recent papers [12,13,14,15,16,17,18,19,20,21,22] we proved, through path integral Monte Carlo (PIMC), that indeed affine quantization is able to make the φ 4 4 theory non-trivial. A crucial point left unanswered in these papers was the validity of the continuum limit at the level of the one- and two-point-functions.
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.

2. Field theory formulation

For a scalar field, φ , with spacial degrees of freedom x = ( x 1 , x 2 , … , x s ) and canonical momentum π ( x ) , the classical affine variables are κ ( x ) ≡ π ( x ) φ ( x ) and φ ( x ) ≠ 0 . The reason we insist that φ ( x ) ≠ 0 is because if φ ( x ) = 0 then x ) = 0 whatever is π ( x ) .
We then introduce the classical Hamiltonian expressed in affine variables. This leads us to
H ( κ , φ ) = ∫ { 1 2 [ κ ( x ) 2 φ ( x ) − 2 + ( ▿ φ ( x ) ) 2 + m 2 φ ( x ) 2 ] + g φ ( x ) r } d s x ,
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. Therefore the AQ bounds 0 < φ ( x ) < ∞  forbid any nonrenormalizability [12,13,14,15,16,17,18,19,20,21,22] which is otherwise possible for CQ [2,3,4,5].
The quantum affine operators are the scalar field φ ^ ( x ) = φ ( x ) and the dilation operator κ ^ ( x ) = [ φ ^ ( x ) π ^ ( x ) + π ^ ( x ) φ ^ ( x ) ] / 2 where the momentum operator is π ^ ( x ) = − i ℏ δ / δ φ ( x ) . Accordingly for the self adjoint kinetic term κ ^ ( x ) φ ^ ( x ) − 2 κ ^ ( x ) = π ^ ( x ) 2 + ( 3 / 4 ) ℏ δ ( 0 ) 2 s φ ( x ) − 2 (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
H ^ ( κ ^ , φ ^ ) = ∫ 1 2 [ π ^ ( x ) 2 + ( ∇ φ ( x ) ) 2 + m 2 φ ( x ) 2 ] + g φ ( x ) r + 3 8 ℏ 2 δ ( 0 ) 2 s φ ( x ) 2 d s x .
The affine action is found adding time, x 0 = c t , where c is the speed of light constant and t is the Euclidean imaginary time, so that S = ∫ 0 β H d x 0 , with H the semi-classical Hamiltonian corresponding to the one of Eq. (2), will then read
S [ φ ] = ∫ 0 β d x 0 ∫ L s d s x 1 2 ∑ μ = 0 s ∂ φ ( x ) ∂ x μ 2 + m 2 φ ( x ) 2 + g φ ( x ) r + 3 8 ℏ δ ( 0 ) 2 s φ ( x ) 2 ,
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 β = L .
The vacuum expectation value of an observable O [ φ ] will then be given by the following expression
〈 O 〉 = ∫ O [ φ ] exp ( − S [ φ ] ) D φ ( x ) ∫ exp ( − S [ φ ] ) D φ ( x ) ,
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 = N a and temporal periodicity β = N a . 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. (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
S [ φ ] a = 1 2 ∑ x , μ a − 2 [ φ ( x ) − φ ( x + e μ ) ] 2 + m 2 ∑ x φ ( x ) 2 + ∑ x g φ ( x ) r + 3 8 ℏ 2 φ ( x ) 2 ,
where e μ is a vector of length a in the + μ direction with μ = 0 , 1 , 2 , … , s . We will have S ≈ S .
In this work we are interested in reaching the continuum limit by taking N a 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. (5) in natural Planck units c = ℏ = 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. (4). In particular m R 2 = p 0 2 〈 | φ ˜ ( p 0 ) | 2 〉 / [ 〈 φ ˜ ( 0 ) 2 〉 − 〈 | φ ˜ ( p 0 ) | 2 〉 ] and g R = [ 3 〈 φ ˜ ( 0 ) 2 〉 2 − 〈 φ ˜ ( 0 ) 4 〉 ] / 〈 φ ˜ ( 0 ) 2 〉 2 , where φ ˜ ( p ) = ∫ d n x e i p · x [ φ ( x ) − 〈 φ ( x ) 〉 ] is the Fourier transform of the field and we choose the 4-momentum p 0 with one spacial component equal to 2 π / N a and all other components equal to zero. We also calulated the one-, two-point-, and two-point-connected-function, respectively given by
V = ∑ x 〈 φ ( x ) 〉 / N n ,
D ( z ) = ∑ x 〈 φ ( x ) φ ( x + z ) 〉 / N n ,
D c ( z ) = ∑ x ( 〈 φ ( x ) φ ( x + z ) 〉 − 〈 φ ( x ) 〉 2 ) / N n = D ( z ) − V 2 .
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 respect to z = L / 2 , namely D ( z ) = D ( L − z ) .

3. 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
φ ( x ) = a − 3 / 2 φ ′ ( x ) .
In other words, we are renormalizing the bare field. This can be compared with the standard renormalization formula
φ ( x ) = Z 1 / 2 φ ren ( x ) .
φ 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
g = a 3 g ′ .
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 r e n :
g = g ren + c 2 ( g ren ) 2 + c 3 ( g ren ) 3 + …
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 (9) and (11) 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.

4. Numerical results

Our PIMC simulations use the Metropolis algorithm [24,25] to calculate the ensemble average of Eq. (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
φ → φ + ( 2 η − 1 ) δ ,
where η is a uniform pseudo random number in [ 0 , 1 ] and δ is the amplitude of the displacement. 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 σ O 2 / ( M N n ) , where σ O 2 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 M N n PIMC displacement moves.

4.1. First case study

In our simulation we first chose the following study case m = g = L = β = 1 and ϵ = 10 − 10 .
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 1), 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 2 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 we show D c ( z ) at increasing values of N = 1 / a . From the plot of the simulation data we see that lim a → 0 D ( 1 / 2 ) = 0 . The width of the spike of D c ( z ) at z = 0 seems to be related to the value of the renormalized mass m R .
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].

4.2. Second case study

For 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 .
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 .
Table 2. Renormalized mass m R , renormalized coupling constant g R , and one-point-value (vacuum expectation value of the field) V = ∑ x 〈 φ ( x ) 〉 / N n for n = 3 + 1 , g = L = β = 1 , m = N / L and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Table 2. Renormalized mass m R , renormalized coupling constant g R , and one-point-value (vacuum expectation value of the field) V = ∑ x 〈 φ ( x ) 〉 / N n for n = 3 + 1 , g = L = β = 1 , m = N / L and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
N m R g R V
4 0.1461(2) 2.01199(2) 0.6672(4)
7 0.0627(1) 2.001154(4) 0.5992(3)
13 0.02463(5) 1.999844(2) 0.5169(2)
25 0.00867(5) 2.000069(7) 0.4359(3)
Notice that the minima of the two symmetric potential wells in the semi-classical Hamiltonian density described by the function f ( φ ) = 1 2 m 2 φ 2 + φ 4 + 3 8 φ − 2 , φ ± ( m ) , are such that
φ ± ( m ) = ± 3 1 / 4 ( 2 m ) − 1 / 2 + O ( m − 7 / 2 ) f o r m ≫ 1 .
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 Figure 3 we see the continuum limit N → ∞ , of the scaled two-point-connected-function where, with an abuse of notation, we dropped from D c , the prime, adopted rigorously in Section 3. 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 Figure 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].

5. 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, ℏ 2 / φ ( 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 ∝ ℏ 2 / φ 2 stemming from the affine quantization procedure [9] not only does not disturb the continuum limit, as we proved here, but in addition is able to render renormalizable the φ 4 4 theory which is known [2] to be non-renormalizable when treated with the more commonly known [1] canonical quantization.

References

  1. Dirac, P.A.M. The Principles of Quantum Mechanics; Claredon Press: Oxford, 1958; in a footnote on page 114. [Google Scholar]
  2. Freedman, B.; Smolensky, P.; Weingarten, D. Monte Carlo Evaluation of the Continuum Limit of ϕ 4 4 and ϕ 3 4 . Physics Letters 1982, 113B, 481. [Google Scholar] [CrossRef]
  3. Aizenman, M. Proof of the Triviality of ϕ d 4 Field Theory and Some Mean-Field Features of Ising Models for d>4. Phys. Rev. Lett. 1981, 47, 886. [Google Scholar] [CrossRef]
  4. Fröhlich, J. On the Triviality of λ ϕ d 4 Theories and the Approach to the Critical Point in d≥4 Dimensions. Nuclear Physics B 1982, 200, 281. [Google Scholar] [CrossRef]
  5. Siefert, J.; Wolff, U. Triviality of φ4 theory in a finite volume scheme adapted to the broken phase. Physics Letters B 2014, 733, 11. [Google Scholar] [CrossRef]
  6. Wolff, U. Triviality of four dimensional ϕ4 theory on the lattice. Scholarpedia 2014, 9, 7367. [Google Scholar] [CrossRef]
  7. Klauder, J.R. A Modern Approach to Functional Integration; Springer, 2010.
  8. Klauder, J.R. Beyond Conventional Quantization; Cambridge University Press, 2000. Chap. 5.
  9. Klauder, J.R. The Benefits of Affine Quantization. Journal of High Energy Physics, Gravitation and Cosmology 2020, 6, 175. [Google Scholar] [CrossRef]
  10. Klauder, J.R. Using Affine Quantization to Analyze Non-renormalizable Scalar Fields and the Quantization of Einstein’s Gravity 2020. arXiv:2006.09156.
  11. Klauder, J.R. An Ultralocal Classical and Quantum Gravity Theory. Journal of High Energy Physics, Gravitation and Cosmology 2020, 6, 656. [Google Scholar] [CrossRef]
  12. Fantoni, R.; Klauder, J.R. Affine quantization of (φ4)4 succeeds while canonical quantization fails. Phys. Rev. D 2021, 103, 076013. [Google Scholar] [CrossRef]
  13. Fantoni, R. Monte Carlo evaluation of the continuum limit of (ϕ12)3. J. Stat. Mech. 2021; 083102. [Google Scholar] [CrossRef]
  14. Fantoni, R.; Klauder, J.R. Monte Carlo evaluation of the continuum limit of the two-point function of the Euclidean free real scalar field subject to affine quantization. J. Stat. Phys. 2021, 184, 28. [Google Scholar] [CrossRef]
  15. Fantoni, R.; Klauder, J.R. Monte Carlo evaluation of the continuum limit of the two-point function of two Euclidean Higgs real scalar fields subject to affine quantization. Phys. Rev. D 2021, 104, 054514. [Google Scholar] [CrossRef]
  16. Fantoni, R.; Klauder, J.R. Eliminating Nonrenormalizability Helps Prove Scaled Affine Quantization of φ 4 4 is Nontrivial. Int. J. Mod. Phys. A 2022, 37, 2250029. [Google Scholar] [CrossRef]
  17. Fantoni, R.; Klauder, J.R. Kinetic Factors in Affine Quantization and Their Role in Field Theory Monte Carlo. Int. J. Mod. Phys. A 2022, 37, 2250094. [Google Scholar] [CrossRef]
  18. Fantoni, R.; Klauder, J.R. Scaled Affine Quantization of φ 4 4 in the Low Temperature Limit. Eur. Phys. J. C 2022, 82, 843. [Google Scholar] [CrossRef]
  19. Fantoni, R.; Klauder, J.R. Scaled Affine Quantization of Ultralocal φ 2 4 a comparative Path Integral Monte Carlo study with scaled Canonical Quantization. Phys. Rev. D 2022, 106, 114508. [Google Scholar] [CrossRef]
  20. Klauder, J.R.; Fantoni, R. The Magnificent Realm of Affine Quantization: valid results for particles, fields, and gravity. Axioms 2023, 12, 911. [Google Scholar] [CrossRef]
  21. Fantoni, R. Scaled Affine Quantization of φ 3 12 is Nontrivial. Modern Physics Letters A (submitted) 2023, arXiv:2011.09862v4. [Google Scholar] [CrossRef]
  22. Klauder, J.R.; Fantoni, R. A modest redirection of quantum field theory solves all current problems. J. Stat. Phys. (submitted) 2023, arXiv:2308.13475. [Google Scholar]
  23. D. M. Ceperley. Rev. Mod. Phys. 1995, 67, 279.
  24. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.M.; Teller, E. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 1087, 21. [Google Scholar] [CrossRef]
  25. Kalos, M.H.; Whitlock, P.A. Monte Carlo Methods; Wiley-Vch Verlag GmbH & Co.: Germany, 2008. [Google Scholar]
  26. Fantoni, R. Localization of acoustic polarons at low temperatures: A path integral Monte Carlo approach. Phys. Rev. B 2012, 86, 144304. [Google Scholar] [CrossRef]
Preprints 94770 g001
Figure 1. Two-point-function D ( z ) = ∑ x 〈 φ ( x ) φ ( x + z ) 〉 / N n for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Figure 1. Two-point-function D ( z ) = ∑ x 〈 φ ( x ) φ ( x + z ) 〉 / N n for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Preprints 94770 g001
Preprints 94770 g002
Figure 2. Two-point-connected-function D c ( z ) = D ( z ) − V 2 for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Figure 2. Two-point-connected-function D c ( z ) = D ( z ) − V 2 for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Preprints 94770 g002
Preprints 94770 g003
Figure 3. Two-point-connected-function D c ( z ) = D ( z ) − V 2 for n = 3 + 1 , g = L = β = 1 , m = N / L and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Figure 3. Two-point-connected-function D c ( z ) = D ( z ) − V 2 for n = 3 + 1 , g = L = β = 1 , m = N / L and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Preprints 94770 g003
Table 1. Renormalized mass m R , renormalized coupling constant g R , and one-point-value (vacuum expectation value of the field) V = ∑ x 〈 φ ( x ) 〉 / N n for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
Table 1. Renormalized mass m R , renormalized coupling constant g R , and one-point-value (vacuum expectation value of the field) V = ∑ x 〈 φ ( x ) 〉 / N n for n = 3 + 1 , m = g = L = β = 1 and N = L / a = 4 , 7 , 13 , 25 . In our PIMC simulations we used Eqs. (4) and (5).
N m R g R V
4 0.1421(2) 2.01178(2) 0.7501(4)
7 0.0602(1) 2.001129(4) 0.7339(2)
13 0.0224(2) 1.999894(5) 0.7216(2)
25 0.0084(1) 2.000048(4) 0.7148(3)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated