Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFourier acceleration, the HMC algorithm and renormalizability
281
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The analysis developed by Luscher and Schaefer of the Hybrid Monte Carlo (HMC) algorithm is extended to include Fourier acceleration. We show for the $phi^4$ theory that Fourier acceleration substantially changes the structure of the theory for both the Langevin and HMC algorithms. When expanded in perturbation theory, each five-dimensional auto-correlation function of the fields $phi(x_i, t_i)$, $1le i le N $, corresponding to a specific 4-dimensional Feynman graph separates into two factors: one depending on the Monte-Carlo evolution times $t_i$ and the second depending on the space-time positions $x_i$. This separation implies that only auto-correlation times at the lattice scale appear, eliminating critical slowing down in perturbation theory.
rate research
Read More
Mass preconditioned HMC and DD-HMC are among the most popular algorithms to simulate Wilson fermions. We present a comparison of the performance of the two algorithms for realistic quark masses and lattice sizes. In particular, we use the locally deflated solver of the DD-HMC environment also for the mass preconditioned simulations.
We present results of a hybrid Monte-Carlo algorithm for dynamical, $n_f=2$, four-dimensional QCD with overlap fermions. The fermionic force requires careful treatment, when changing topological sectors. The pion mass dependence of the topological susceptibility is studied on $6^4$ and $12cdot 6^3$ lattices. The results are transformed into physical units.
For an asymptotically free theory, a promising strategy for eliminating Critical Slowing Down (CSD) is naive Fourier acceleration. This requires the introduction of gauge-fixing into the action, in order to isolate the asymptotically decoupled Fourier modes. In this article, we present our approach and results from a gauge-fixed Fourier-accelerated hybrid Monte Carlo algorithm, using an action that softly fixes the gauge links to Landau gauge. We compare the autocorrelation times with those of the pure hybrid Monte Carlo algorithm. We work on a small-volume lattice at weak coupling. We present preliminary results and obstacles from working with periodic boundary conditions, and then we present results from using fixed, equilibrated boundary links to avoid $mathbb{Z}_3$ and other topological barriers and to anticipate applying a similar acceleration to many small cells in a large, physically-relevant lattice volume.
Polynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in HMC simulations. The use of a multiple time-scale integration scheme allows the filtered pseudofermions to be evolved using a coarse step size. We introduce a novel generalisation of the nested leapfrog which allows for far greater flexibility in the choice of time scales. We observe a reduction in the computational expense of the molecular dynamics integration of between 3--5 which improves as the quark mass decreases.
We report on a study of 2+1 flavor lattice QCD with the $O(a)$-improved Wilson quarks on a $16^3times 32$ lattice at the lattice spacing $1/aapprox 2$GeV employing Lueschers domain-decomposed HMC(LDDHMC) algorithm. This is dedicated to a preliminary study for the PACS-CS project which plans to complete the Wilson-clover $N_f=2+1$ program lowering the up-down quark masses close to the physical values as much as possible. We focus on three issues: (i) how light quark masses we can reach with LDDHMC, (ii) efficiency of the algorithm compared with the conventional HMC, (iii) parameter choice for the production runs on PACS-CS.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
