Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGoldstone mode singularities in O(n) models
90
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Monte Carlo (MC) analysis of the Goldstone mode singularities for the transverse and the longitudinal correlation functions, behaving as G_{perp}(k) simeq ak^{-lambda_{perp}} and G_{parallel}(k) simeq bk^{-lambda_{parallel}} in the ordered phase at k -> 0, is performed in the three-dimensional O(n) models with n=2, 4, 10. Our aim is to test some challenging theoretical predictions, according to which the exponents lambda_{perp} and lambda_{parallel} are non-trivial (3/2<lambda_{perp}<2 and 0<lambda_{parallel}<1 in three dimensions) and the ratio bM^2/a^2 (where M is a spontaneous magnetization) is universal. The trivial standard-theoretical values are lambda_{perp}=2 and lambda_{parallel}=1. Our earlier MC analysis gives lambda_{perp}=1.955 pm 0.020 and lambda_{parallel} about 0.9 for the O(4) model. A recent MC estimation of lambda_{parallel}, assuming corrections to scaling of the standard theory, yields lambda_{parallel} = 0.69 pm 0.10 for the O(2) model. Currently, we have performed a similar MC estimation for the O(10) model, yielding lambda_{perp} = 1.9723(90). We have observed that the plot of the effective transverse exponent for the O(4) model is systematically shifted down with respect to the same plot for the O(10) model by Delta lambda_{perp} = 0.0121(52). It is consistent with the idea that 2-lambda_{perp} decreases for large $n$ and tends to zero at n -> infty. We have also verified and confirmed the expected universality of bM^2/a^2 for the O(4) model, where simulations at two different temperatures (couplings) have been performed.
rate research
Read More
Correlation functions in the O(n) models below the critical temperature are considered. Based on Monte Carlo (MC) data, we confirm the fact stated earlier by Engels and Vogt, that the transverse two-plane correlation function of the O(4) model for lattice sizes about L=120 and small external fields h is very well described by a Gaussian approximation. However, we show that fits of not lower quality are provided by certain non-Gaussian approximation. We have also tested larger lattice sizes, up to L=512. The Fourier-transformed transverse and longitudinal two-point correlation functions have Goldstone mode singularities in the thermodynamic limit at k --> 0 and h=+0, i.e., G_perp(k) = a k^{-lambda_perp} and G_parallel(k) = b k^{-lambda_parallel}, respectively. Here a and b are the amplitudes, k is the magnitude of the wave vector. The exponents lambda_perp, lambda_parallel and the ratio b M^2/a^2, where M is the spontaneous magnetization, are universal according to the GFD (grouping of Feynman diagrams) approach. Here we find that the universality follows also from the standard (Gaussian) theory, yielding b M^2/a^2 = (n-1)/16. Our MC estimates of this ratio are 0.06 +/- 0.01 for n=2, 0.17 +/- 0.01 for n=4 and 0.498 +/- 0.010 for n=10. According to these and our earlier MC results, the asymptotic behavior and Goldstone mode singularities are not exactly described by the standard theory. This is expected from the GFD theory. We have found appropriate analytic approximations for G_perp(k) and G_parallel(k), well fitting the simulation data for small k. We have used them to test the Patashinski--Pokrovski relation and have found that it holds approximately.
Power-law singularities and critical exponents in n-vector models are considered from different theoretical points of view. It includes a theoretical approach called the GFD (grouping of Feynman diagrams) theory, as well as the perturbative renormalization group (RG) treatment. A non-perturbative proof concerning corrections to scaling in the two-point correlation function of the phi^4 model is provided, showing that predictions of the GFD theory rather than those of the perturbative RG theory can be correct. Critical exponents determined from highly accurate experimental data very close to the lambda-transition point in liquid helium, as well as the Goldstone mode singularities in n-vector spin models, evaluated from Monte Carlo simulation results, are discussed with an aim to test the theoretical predictions. Our analysis shows that in both cases the data can be well interpreted within the GFD theory.
A relation between O$(n)$ models and Ising models has been recently conjectured [L. Casetti, C. Nardini, and R. Nerattini, Phys. Rev. Lett. 106, 057208 (2011)]. Such a relation, inspired by an energy landscape analysis, implies that the microcanonical density of states of an O$(n)$ spin model on a lattice can be effectively approximated in terms of the density of states of an Ising model defined on the same lattice and with the same interactions. Were this relation exact, it would imply that the critical energy densities of all the O$(n)$ models (i.e., the average values per spin of the O$(n)$ Hamiltonians at their respective critical temperatures) should be equal to that of the corresponding Ising model; it is therefore worth investigating how different the critical energies are and how this difference depends on $n$. We compare the critical energy densities of O$(n)$ models in three dimensions in some specific cases: the O$(1)$ or Ising model, the O$(2)$ or $XY$ model, the O$(3)$ or Heisenberg model, the O$(4)$ model and the O$(infty)$ or spherical model, all defined on regular cubic lattices and with ferromagnetic nearest-neighbor interactions. The values of the critical energy density in the $n=2$, $n=3$, and $n=4$ cases are derived through a finite-size scaling analysis of data produced by means of Monte Carlo simulations on lattices with up to $128^3$ sites. For $n=2$ and $n=3$ the accuracy of previously known results has been improved. We also derive an interpolation formula showing that the difference between the critical energy densities of O$(n)$ models and that of the Ising model is smaller than $1%$ if $n<8$ and never exceeds $3%$ for any $n$.
On the example of the spherical model we study, as a function of the temperature $T$, the behavior of the Casimir force in O(n) systems with a diffuse interface and slab geometry $infty^{d-1}times L$, where $2<d<4$ is the dimensionality of the system. We consider a system with nearest-neighbor anisotropic interaction constants $J_parallel$ parallel to the film and $J_perp$ across it. The model represents the $ntoinfty$ limit of O(n) models with antiperiodic boundary conditions applied across the finite dimension $L$ of the film. We observe that the Casimir amplitude $Delta_{rm Casimir}(d|J_perp,J_parallel)$ of the anisotropic $d$-dimensional system is related to that one of the isotropic system $Delta_{rm Casimir}(d)$ via $Delta_{rm Casimir}(d|J_perp,J_parallel)=(J_perp/J_parallel)^{(d-1)/2} Delta_{rm Casimir}(d)$. For $d=3$ we find the exact Casimir amplitude $ Delta_{rm Casimir}= [ {rm Cl}_2 (pi/3)/3-zeta (3)/(6 pi)](J_perp/J_parallel)$, as well as the exact scaling functions of the Casimir force and of the helicity modulus $Upsilon(T,L)$. We obtain that $beta_cUpsilon(T_c,L)=(2/pi^{2}) [{rm Cl}_2(pi/3)/3+7zeta(3)/(30pi)] (J_perp/J_parallel)L^{-1}$, where $T_c$ is the critical temperature of the bulk system. We find that the effect of the helicity is thus strong that the Casimir force is repulsive in the whole temperature region.
In previous work with Scullard, we defined a graph polynomial P_B(q,T) that gives access to the critical temperature T_c of the q-state Potts model on a general two-dimensional lattice L. It depends on a basis B, containing n x m unit cells of L, and the relevant root of P_B(q,T) was observed to converge quickly to T_c in the limit n,m to infinity. Moreover, in exactly solvable cases there is no finite-size dependence at all. We reformulate this method as an eigenvalue problem within the periodic Temperley-Lieb algebra. This corresponds to taking m to infinity first, so the bases B are semi-infinite cylinders of circumference n. The limit implies faster convergence in n, while maintaining the n-independence in exactly solvable cases. In this setup, T_c(n) is determined by equating the largest eigenvalues of two topologically distinct sectors of the transfer matrix. Crucially, these two sectors determine the same critical exponent in the continuum limit, and the observed fast convergence is thus corroborated by results of conformal field theory. We obtain similar results for the dense and dilute phases of the O(N) loop model, using now a transfer matrix within the dilute periodic Temperley-Lieb algebra. The eigenvalue formulation allows us to double the size n for which T_c(n) can be obtained, using the same computational effort. We study in details three significant cases: (i) bond percolation on the kagome lattice, up to n = 14; (ii) site percolation on the square lattice, to n = 21; and (iii) self-avoiding polygons on the square lattice, to n = 19. Convergence properties of T_c(n) and extrapolation schemes are studied in details for the first two cases. This leads to rather accurate values for the percolation thresholds: p_c = 0.524404999167439(4) for bond percolation on the kagome lattice, and p_c = 0.59274605079210(2) for site percolation on the square lattice.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
