Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMonte Carlo simulation of lattice ${rm CP}^{N-1}$ models at large N
71
0
0.0
(
0
)
Authors
Ettore Vicari
Ask ChatGPT about the research
No Arabic abstract
In order to check the validity and the range of applicability of the 1/N expansion, we performed numerical simulations of the two-dimensional lattice CP(N-1) models at large N, in particular we considered the CP(20) and the CP(40) models. Quantitative agreement with the large-N predictions is found for the correlation length defined by the second moment of the correlation function, the topological susceptibility and the string tension. On the other hand, quantities involving the mass gap are still far from the large-$N$ results showing a very slow approach to the asymptotic regime. To overcome the problems coming from the severe form of critical slowing down observed at large N in the measurement of the topological susceptibility by using standard local algorithms, we performed our simulations implementing the Simulated Tempering method.
rate research
Read More
We develop numerical tools for Diagrammatic Monte-Carlo simulations of non-Abelian lattice field theories in the tHooft large-N limit based on the weak-coupling expansion. First we note that the path integral measure of such theories contributes a bare mass term in the effective action which is proportional to the bare coupling constant. This mass term renders the perturbative expansion infrared-finite and allows to study it directly in the large-N and infinite-volume limits using the Diagrammatic Monte-Carlo approach. On the exactly solvable example of a large-N O(N) sigma model in D=2 dimensions we show that this infrared-finite weak-coupling expansion contains, in addition to powers of bare coupling, also powers of its logarithm, reminiscent of re-summed perturbation theory in thermal field theory and resurgent trans-series without exponential terms. We numerically demonstrate the convergence of these double series to the manifestly non-perturbative dynamical mass gap. We then develop a Diagrammatic Monte-Carlo algorithm for sampling planar diagrams in the large-N matrix field theory, and apply it to study this infrared-finite weak-coupling expansion for large-N U(N)xU(N) nonlinear sigma model (principal chiral model) in D=2. We sample up to 12 leading orders of the weak-coupling expansion, which is the practical limit set by the increasingly strong sign problem at high orders. Comparing Diagrammatic Monte-Carlo with conventional Monte-Carlo simulations extrapolated to infinite N, we find a good agreement for the energy density as well as for the critical temperature of the deconfinement transition. Finally, we comment on the applicability of our approach to planar QCD at zero and finite density.
We consider properties of the inhomogeneous solution found recently for mbox{$mathbb{CP}^{,N-1}$} model. The solution was interpreted as a soliton. We reevaluate its energy in three different ways and find that it is negative contrary to the previous claims. Hence, instead of the solitonic interpretation it calls for reconsideration of the issue of the true ground state. While complete resolution is still absent we show that the energy density of the periodic elliptic solution is lower than the energy density of the homogeneous ground state. We also discuss similar solutions for the ${mathbb{O}}(N)$ model and for SUSY extensions.
The lattice provides a powerful tool to non-perturbatively investigate strongly coupled supersymmetric Yang-Mills (SYM) theories. The pure SU(2) SYM theory with one supercharge is simulated on large lattices with small Majorana gluino masses down to about $am_{tilde g}=0.068$ with lattice spacing $asimeq 0.125$ fm. The gluino dynamics is simulated by the Two-Step Multi-Boson (TSMB) and the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) algorithms. Supersymmetry (SUSY) is broken explicitly by the lattice and the Wilson term and softly by the presence of a non-vanishing gluino mass. However, the recovery of SUSY is expected in the infinite volume continuum limit by tuning the bare parameters to the SUSY point in the parameter space. This scenario is studied by the determination of the low-energy mass spectrum and by means of lattice SUSY Ward-Identities (WIs).
The topological charge distribution P(Q) is calculated for lattice ${rm CP}^{N-1}$ models. In order to suppress lattice cut-off effects we employ a fixed point (FP) action. Through transformation of P(Q) we calculate the free energy $F(theta)$ as a function of the $theta$ parameter. For N=4, scaling behavior is observed for P(Q), $F(theta)$ as well as the correlation lengths $xi(Q)$. For N=2, however, scaling behavior is not observed as expected. For comparison, we also make a calculation for the ${rm CP}^{3}$ model with standard action. We furthermore pay special attention to the behavior of P(Q) in order to investigate the dynamics of instantons. For that purpose, we carefully look at behavior of $gamma_{it eff}$, which is an effective power of P(Q)($sim exp(-CQ^{gamma_{it eff}})$), and reflects the local behavior of P(Q) as a function of Q. We study $gamma_{it eff}$ for two cases, one of which is the dilute gas approximation based on the Poisson distribution of instantons and the other is the Debye-Huckel approximation of instanton quarks. In both cases we find similar behavior to the one observed in numerical simulations.
We present a method for direct hybrid Monte Carlo simulation of graphene on the hexagonal lattice. We compare the results of the simulation with exact results for a unit hexagonal cell system, where the Hamiltonian can be solved analytically.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
