Do you want to publish a course? Click here

A New Exact Method for Dynamical Fermion Computations with Non-Local Actions

58   0   0.0 ( 0 )
 Added by Anthony D. Kennedy
 Publication date 1998
  fields
and research's language is English




Ask ChatGPT about the research

We introduce a new algorithm which we call the {Rational Hybrid Monte Carlo} Algorithm (RHMC). This method uses a rational approximation to the fermionic kernel together with a noisy Kennedy-Kuti acceptance step to give an efficient algorithm with no molecular dynamics integration step-size errors.



rate research

Read More

The fermion bag is a powerful idea that helps to solve fermion lattice field theories using Monte Carlo methods. Some sign problems that had remained unsolved earlier can be solved within this framework. In this work we argue that the fermion bag also gives insight into a new mechanism of fermion mass generation, especially at strong couplings where fermion masses are related to the fermion bag size. On the other hand, chiral condensates arise due to zero modes in the Dirac operator within a fermion bag. Although in traditional four-fermion models the two quantities seem to be related, we show that they can be decoupled. While fermion bags become small at strong couplings, the ability of zero modes of the Dirac operator within fermion bags to produce a chiral condensate, can be suppressed by the presence of additional zero modes from other fermions. Thus, fermions can become massive even without a chiral condensate. This new mechanism of mass generation was discovered long ago in lattice field theory, but has remained unappreciated. Recent work suggests that it may be of interest even in continuum quantum field theory.
We consider a field theoretical model where a SU(2) fermion doublet, subjected to non-Abelian gauge interactions, is also coupled to a complex scalar field doublet via a Yukawa and an irrelevant Wilson-like term. Despite the presence of these two chiral breaking operators in the Lagrangian, an exact symmetry acting on fermions and scalars prevents perturbative mass corrections. In the phase where fermions are massless (Wigner phase) the Yukawa coupling can be tuned to a critical value at which chiral transformations acting on fermions only become a symmetry of the theory (up to cutoff effects). In the Nambu-Goldstone phase of the critical theory a fermion mass term of dynamical origin is expected to arise in the Ward identities of the purely fermionic chiral transformations. Such a non-perturbative mechanism of dynamical mass generation can provide a natural (`a la t Hooft) alternative to the Higgs mechanism adopted in the Standard Model. Here we lay down the theoretical framework necessary to demonstrate the existence of this mechanism by means of lattice simulations.
We report on a study of the light quark spectrum using an improved gauge action and both Kogut-Susskind and Naik quark actions. We have studied six different lattice spacings, corresponding to plaquette couplings ranging from 6.8 to 7.9, with five to six quark masses per coupling. We compare the two quark actions in terms of the spectrum and restoration of flavor symmetry. We also compare these results with those from the conventional action.
We report on our first experiences in simulating Neuberger valence fermions on CLS $N_f=2$ configurations with light sea quark masses and small lattice spacings. Valence quark masses are considered that allow to explore the matching to (partially quenched) chiral perturbation theory both in the $epsilon$- and $p$-regimes. The setup is discussed, and first results are presented for spectral observables.
Lattice QCD calculations including the effects of one or more non-degenerate sea quark flavors are conventionally performed using the Rational Hybrid Monte Carlo (RHMC) algorithm, which computes the square root of the determinant of $mathscr{D}^{dagger} mathscr{D}$, where $mathscr{D}$ is the Dirac operator. The special case of two degenerate quark flavors with the same mass is described directly by the determinant of $mathscr{D}^{dagger} mathscr{D}$ --- in particular, no square root is necessary --- enabling a variety of algorithmic developments, which have driven down the cost of simulating the light (up and down) quarks in the isospin-symmetric limit of equal masses. As a result, the relative cost of single quark flavors --- such as the strange or charm --- computed with RHMC has become more expensive. This problem is even more severe in the context of our measurements of the $Delta I = 1/2$ $K rightarrow pi pi$ matrix elements on lattice ensembles with $G$-parity boundary conditions, since $G$-parity is associated with a doubling of the number of quark flavors described by $mathscr{D}$, and thus RHMC is needed for the isospin-symmetric light quarks as well. In this paper we report on our implementation of the exact one flavor algorithm (EOFA) introduced by the TWQCD collaboration for simulations including single flavors of domain wall quarks. We have developed a new preconditioner for the EOFA Dirac equation, which both reduces the cost of solving the Dirac equation and allows us to re-use the bulk of our existing high-performance code. Coupling these improvements with careful tuning of our integrator, the time per accepted trajectory in the production of our 2+1 flavor $G$-parity ensembles with physical pion and kaon masses has been decreased by a factor of 4.2.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا