Molecular dynamics for full QCD simulations with an improved action

Xiang-Qian Luo

DOI: https://doi.org/10.1016/0010-4655%2895%2900128-X

2004-03-27

Abstract:I derive the equation of motion in molecular dynamics for doing full lattice QCD simulations with clover quarks. The even-odd preconditioning technique, expected to significantly reduce the computational effort, is further developed for the simulations.

High Energy Physics - Lattice

What problem does this paper attempt to address?

Dynamical Simulations with Highly Improved Staggered Quarks

K. Y. Wong,R. M. Woloshyn

DOI: https://doi.org/10.48550/arXiv.0710.0737

2007-10-11

Abstract:It is well established that lattice artifacts can be suppressed substantially by the use of SU(3)-projected smeared links in the fermion action. An example is the Highly Improved Staggered Quark action where the ASQ-like effective links are constructed from reunitarized Fat7 links. A general procedure is presented for computing the derivative of the fermion action with respect to the base links (fermion force) - a key component in dynamical simulations using molecular dynamics evolution. The method is iterative and can be applied to actions with arbitrary levels of smearing and reunitarization. The cost of calculating the fermion force is determined for the ASQ action and the HISQ action. Test results show that calculating the HISQ force is about two times more expensive than the ASQ force.

High Energy Physics - Lattice
Quantum Simulation of Lattice QCD with Improved Hamiltonians

Anthony N. Ciavarella

2023-11-07

Abstract:Quantum simulations of lattice gauge theories are anticipated to directly probe the real time dynamics of QCD, but scale unfavorably with the required truncation of the gauge fields. Improved Hamiltonians are derived to correct for the effects of gauge field truncations on the SU(3) Kogut-Susskind Hamiltonian. It is shown in $1+1D$ that this enables low chromo-electric field truncations to quantitatively reproduce features of the untruncated theory over a range of couplings and quark masses. In $3+1D$, an improved Hamiltonian is derived for lattice QCD with staggered massless fermions. It is shown in the strong coupling limit that the spectrum qualitatively reproduces aspects of two flavor QCD and simulations of a small system are performed on IBM's {\tt Perth} quantum processor.

High Energy Physics - Lattice,Nuclear Theory,Quantum Physics
A Numerical Study of Improved Quark Actions on Anisotropic Lattices

SQ Su,LM Liu,X Li,CA Liu

DOI: https://doi.org/10.1142/s0217751x06024967

2007-01-01

Modern Physics Letters A

Abstract:Tadpole improved Wilson quark actions with clover terms on anisotropic lattices are studied numerically. Using asymmetric lattice volumes, the pseudo-scalar meson dispersion relations are measured for 8 lowest lattice momentum modes with quark mass values ranging from the strange to the charm quark with various values of the gauge coupling β and 3 different values of the bare speed of light parameter ν. These results can be utilized to extrapolate or interpolate to obtain the optimal value for the bare speed of light parameter ν_opt(m) at a given gauge coupling for all bare quark mass values m. In particular, the optimal values of ν at the physical strange and charm quark mass are given for various gauge couplings. The lattice action with these optimized parameters can then be used to study physical properties of hadrons involving either light or heavy quarks.
Droems: experimental mathematics, informatics and infinite dimensional geometry

Denis V. Juriev

DOI: https://doi.org/10.48550/arXiv.cs/9809119

1998-09-29

Abstract:The article is devoted to a problem of elaboration of the real-time interactive videosystems for accelerated nonverbal cognitive computer and telecommunications. The proposed approach is based on the using of droems (dynamically reconstructed objects of experimental mathematics) and interpretational figures as pointers to them. Four paragraphs of the article are devoted to (1) an exposition of basic notions of the interpretational geometry, (2) the operator methods in the theory of interactive dynamical videosystems, (3) the general concept of organization of the integrated interactive real-time videocognitive systems, (4) the droems and processes of their dynamical reconstruction, where the general notions are illustrated by a concrete example related to the infinite dimensional geometry. The exposition is presumably heuristic and conceptual (the first and the third paragraphs) though some particular aspects such as content of the second and the fourth paragraphs, which allow deeper formalization and detailing in present, are exposed on the mathematical level of rigor.

Human-Computer Interaction,Graphics,Representation Theory
Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Ripeng Luo,Kuang Yu

DOI: https://doi.org/10.1063/5.0022013

2020-01-01

Abstract:Molecular dynamics (MD) is a powerful (and the most viable) tool to compute the thermal conductivities of solid disordered materials. However, conventional classical MD fails to describe the nuclear quantum effects (NQEs), so it may give inaccurate results for light materials at low temperatures. While the importance of NQE has been widely acknowledged, yet we do not have a fully reliable method to account for NQE in the MD thermal conductivity calculations. In this work, we will investigate and analyze the performances of a number of path-integral-based quantum MD methods, using ordered ice as a test case. To establish the validity of these methods, we will compare the MD results with the lattice dynamics results, in both classical and quantum limits. Through such a comparison, we will show that methods such as ring polymer MD stand as a good approach for a complex solid with short phonon lifetimes but could be problematic when describing long-living acoustic phonons. In addition, we will show that the rigid water model, which is the state-of-the-art model in the studies of ice/water systems, fails to capture most of the NQEs in ice thermal conductivity. Neglecting librational and translational NQEs leads to essential errors, which clearly demonstrates the importance of a true quantum simulation method that treats all modes at a consistent quantum level.
Simulation of quantum chromodynamics on the lattice with exactly chiral lattice fermions

sinya aoki,tingwai chiu,guido cossu,xu feng,hidenori fukaya,shoji hashimoto,tunghan hsieh,takashi kaneko,hideo matsufuru,j noaki,tetsuya onogi,eigo shintani,kouhei takeda

DOI: https://doi.org/10.1093/ptep/pts006

2012-01-01

Progress of Theoretical and Experimental Physics

Abstract:Numerical simulation of the low-energy dynamics of quarks and gluons is now feasible based on the fundamental theory of strong interaction, i.e. quantum chromodynamics (QCD). With QCD formulated on a 4D hypercubic lattice (called lattice QCD or LQCD), one can simulate the QCD vacuum and hadronic excitations on the vacuum using teraflop-scale supercomputers, which have become available in the last decade. A great deal of work has been done on this subject by many groups around the world; in this article we summarize the work done by the JLQCD and TWQCD collaborations since 2006. These collaborations employ Neuberger's overlap fermion formulation, which preserves the exact chiral and flavor symmetries on the lattice, unlike other lattice fermion formulations. Because of this beautiful property, numerical simulation of the formulation can address fundamental questions on the QCD vacuum, such as the microscopic structure of the quark-antiquark condensate in the chirally broken phase of QCD and its relation to SU(3) gauge field topology. Tests of the chiral effective theory, which is based on the assumption that the chiral symmetry is spontaneously broken in the QCD vacuum, can be performed, including the pion-loop effect test. For many other phenomenological applications, we adopt the all-to-all quark propagator technique, which allows us to compute various correlation functions without substantial extra cost. The benefit of this is not only that the statistical signal is improved but that disconnected quark-loop diagrams can be calculated. Using this method combined with the overlap fermion formulation, we study a wide range of physical quantities that are of both theoretical and phenomenological interest.
Generalized HMC using Nambu mechanics for lattice QCD

Erik Lundstrum

2024-10-18

Abstract:I describe a generalization of the Hybrid Monte Carlo (HMC) algorithm in which the molecular dynamics (MD) steps utilize Nambu generalized Hamiltonian dynamics. Characterized by multiple Hamiltonian functions, this formalism allows me to include forces from non-local objects in the MD evolution while preserving the target probability distribution. In this way, the changes proposed by the MD at one location can be made using instantaneous knowledge of the long-distance behavior of the gauge field to a degree beyond that usually provided by the fermion determinant. This represents a promising method for reducing critical slowing down in lattice QCD simulations.

High Energy Physics - Lattice
Coarsest-level improvements in multigrid for lattice QCD on large-scale computers

Jesus Espinoza-Valverde,Andreas Frommer,Gustavo Ramirez-Hidalgo,Matthias Rottmann

DOI: https://doi.org/10.48550/arXiv.2205.09104

2023-03-27

Abstract:Numerical simulations of quantum chromodynamics (QCD) on a lattice require the frequent solution of linear systems of equations with large, sparse and typically ill-conditioned matrices. Algebraic multigrid methods are meanwhile the standard for these difficult solves. Although the linear systems at the coarsest level of the multigrid hierarchy are much smaller than the ones at the finest level, they can be severely ill-conditioned, thus affecting the scalability of the whole solver. In this paper, we investigate different novel ways to enhance the coarsest-level solver and demonstrate their potential using DD-$\alpha$AMG, one of the publicly available algebraic multigrid solvers for lattice QCD. We do this for two lattice discretizations, namely clover-improved Wilson and twisted mass. For both the combination of two of the investigated enhancements, deflation and polynomial preconditioning, yield significant improvements in the regime of small mass parameters. In the clover-improved Wilson case we observe a significantly improved insensitivity of the solver to conditioning, and for twisted mass we are able to get rid of a somewhat artificial increase of the twisted mass parameter on the coarsest level used so far to make the coarsest level solves converge more rapidly.

Numerical Analysis,Distributed, Parallel, and Cluster Computing,High Energy Physics - Lattice
Controlling atom-atom interaction at ultralow temperatures by dc electric fields

M. Marinescu,L. You

DOI: https://doi.org/10.1103/PhysRevLett.81.4596

1998-10-15

Abstract:We propose a physical mechanism for tuning the atom-atom interaction strength at ultra-low temperatures. In the presence of a dc electric field the interatomic potential is changed due to the effective dipole-dipole interaction between the polarized atoms. Detailed multi-channel scattering calculations reveal features never before discussed for ultra-cold atomic collisions. We demonstrate that optimal control of the effective atom-atom interactions can be achieved under reasonable laboratory conditions. Implications of this research on the physics of atomic Bose-Einstein condensation (BEC) and on the pursuit for atomic degenerate fermion gases will be discussed.

Quantum Physics,Condensed Matter
Calculating the I=2 pion scattering length using tadpole improved clover Wilson action on coarse anisotropic lattices

Chuan Liu,Junhua Zhang,Ying Chen,JianPing Ma

DOI: https://doi.org/10.1016/S0550-3213(01)00662-9

IF: 3.045

2002-01-01

Nuclear Physics B

Abstract:In an exploratory study, using the tadpole improved clover Wilson quark action on small, coarse and anisotropic lattices, the ππ scattering length in the I=2 channel is calculated within quenched approximation. A new method is proposed which enables us to make chiral extrapolation of our lattice results without calculating the decay constant on the lattice. Finite volume and finite lattice spacing errors are analyzed and the results are extrapolated towards the infinite volume and continuum limit. Comparisons of our lattice results with the new experiment and the results from Chiral Perturbation Theory are made. Good agreements are found.
Reducing Discretization Errors in Lattice QCD Spectroscopy

Sara Collins,R. G. Edwards,U. M. Heller,John Sloan,R.G. Edwards,U.M. Heller

DOI: https://doi.org/10.48550/arXiv.hep-lat/9611022

1996-11-22

High Energy Physics - Lattice

Abstract:The improved Wilson quark action - the clover action - is constructed to have smaller discretization errors than the normal Wilson quark action. We test this in a quenched spectroscopy computation on 6 lattice ensembles with spacings from 0.15 to 0.43 fm. To ensure that the dominant scaling violations come from the fermions we use an $O(a^2)$ improved 6-link $SU(3)$ pure gauge action. We find evidence that fermionic scaling violations are consistent with $O(a^2)$ for clover fermions and $O(a)$ with a nonnegligible $O(a^2)$ term for standard Wilson fermions. This latter mixed ansatz makes a reliable continuum extrapolation problematic for Wilson fermions. For clover fermions, on the other hand, we obtain accurate predictions for hadron masses in quenched continuum QCD. We find that the slopes of the scaling violations are roughly 200 MeV for both Wilson and clover fermions.
Full QCD with milder topological freezing

Claudio Bonanno,Giuseppe Clemente,Massimo D'Elia,Lorenzo Maio,Luca Parente

2024-09-02

Abstract:We simulate $N_f=2+1$ QCD at the physical point combining open and periodic boundary conditions in a parallel tempering framework, following the original proposal by M. Hasenbusch for $2d$ $\mathrm{CP}^{N-1}$ models, which has been recently implemented and widely employed in $4d$ $\mathrm{SU}(N)$ pure Yang-Mills theories too. We show that using this algorithm it is possible to achieve a sizable reduction of the auto-correlation time of the topological charge in dynamical fermions simulations both at zero and finite temperature, allowing to avoid topology freezing down to lattice spacings as fine as $a \sim 0.02$ fm. Therefore, this implementation of the Parallel Tempering on Boundary Conditions algorithm has the potential to substantially push forward the investigation of the QCD vacuum properties by means of lattice simulations.

High Energy Physics - Lattice
Lattice QCD with exponentially small chirality breaking

A.A. Slavnov

DOI: https://doi.org/10.1016/s0920-5632(00)91763-8

2000-04-01

Abstract:A new multifield formulation of lattice QCD is proposed. The model is free of spectrum doubling and preserves all nonanomalous chiral symmetries up to exponentially small corrections. It is argued that a small number of fields may provide a good approximation making computer simulations feasible.
BQCD - Berlin quantum chromodynamics program

Yoshifumi Nakamura,Hinnerk Stüben

DOI: https://doi.org/10.48550/arXiv.1011.0199

2010-11-01

Abstract:We publish BQCD as free software under the GNU General Public License. BQCD is a Hybrid Monte-Carlo program that simulates lattice QCD with dynamical Wilson fermions. It is one of the main production programs of the QCDSF collaboration. The program can simulate 2 and 2 + 1 fermion flavours with pure, clover improved, and stout smeared fat link Wilson fermions as well as standard plaquette, and an improved (rectangle) gauge action. The single flavour is simulated with the Rational Hybrid Monte-Carlo algorithm.

High Energy Physics - Lattice
Spontaneous chiral-symmetry breaking of lattice QCD with massless dynamical quarks

XiangQian Luo

DOI: https://doi.org/10.1007/s11433-007-2015-5

2007-02-01

Science in China Series G: Physics, Mechanics and Astronomy

Abstract:One of the most challenging issues in QCD is the investigation of spontaneous chiral-symmetry breaking, which is characterized by the non-vanishing chiral condensate when the bare fermion mass is zero. In standard methods of the lattice gauge theory, one has to perform expensive simulations at multiple bare quark masses, and employ some modeled functions to extrapolate the data to the chiral limit. This paper applies the probability distribution function method to computing the chiral condensate in lattice QCD with massless dynamical quarks, without any ambiguous mass extrapolation. The results for staggered quarks indicate that this might be a promising and efficient method for investigating the spontaneous chiral-symmetry breaking in lattice QCD, which deserves further investigation.
Improved Dirichlet boundary conditions for lattice gauge-fermion theories

Tobias Rindlisbacher,Kari Rummukainen,Ahmed Salami

2024-01-11

Abstract:Hybrid Monte Carlo (HMC) simulations of lattice gauge theories with fermionic matter rely on the invertibility of the lattice Dirac operator. Near-zero modes of the latter can therefore significantly slow down the update algorithm and cause instabilities. This is in particular a problem when dealing with massless fermions. Homogeneous temporal Dirichlet boundary conditions can be used to remove zero modes from massless lattice Dirac operators, but the standard implementation of these boundary conditions can cause severe finite-volume cutoff effects in regions of parameter space where the physics at the ultraviolet (UV) cutoff scale is dominated by the fermionic instead of the gauge action. In lattice quantum chromodynamics (QCD) this is usually not an issue, as the gauge action dominates the UV physics and the problem does not show up. In studies of beyond standard model (BSM) theories, on the other hand, the finite-volume artifacts can be severe. We have identified the origin of these IR cutoff effects and propose a simple improvement on the homogeneous temporal Dirichlet boundary conditions to prevent them. We demonstrate the benefits of using our improved boundary conditions at the example of SU(2) lattice gauge theory with $N_f=24$ massless Wilson-clover flavors. Due to the large number of fermions in this theory, the boundary-related finite volume artifacts are particularly strong, and the effect from switching from the normal to our improved homogeneous Dirichlet boundary conditions is therefore distinct.

High Energy Physics - Lattice
Quantum benefit of the quantum equation of motion for the strongly coupled many-body problem

Manqoba Q. Hlatshwayo,John Novak,Elena Litvinova

2023-09-19

Abstract:We investigate the quantum equation of motion (qEOM), a hybrid quantum-classical algorithm for computing excitation properties of a fermionic many-body system, with a particular emphasis on the strong-coupling regime. The method is designed as a stepping stone towards building more accurate solutions for strongly coupled fermionic systems, such as medium-heavy nuclei, using quantum algorithms to surpass the current barrier in classical computation. Approximations of increasing accuracy to the exact solution of the Lipkin-Meshkov-Glick Hamiltonian with $N=8$ particles are studied on digital simulators and IBM quantum devices. Improved accuracy is achieved by applying operators of growing complexity to generate excitations above the correlated ground state, which is determined by the variational quantum eigensolver (VQE). We demonstrate explicitly that the qEOM exhibits a quantum benefit due to the independence of the number of required quantum measurements from the configuration complexity. Post-processing examination shows that quantum device errors are amplified by increasing configuration complexity and coupling strength. A detailed error analysis is presented, and error mitigation based on zero noise extrapolation is implemented.

Nuclear Theory,Other Condensed Matter,Quantum Physics
Improved lattice actions for supersymmetric quantum mechanics

Sebastian Schierenberg,Falk Bruckmann

DOI: https://doi.org/10.1103/PhysRevD.89.014511

2012-10-19

Abstract:We analyze the Euclidean version of supersymmetric quantum mechanics on the lattice by means of a numerical path integral. We consider two different lattice derivatives and improve the actions containing them with respect to supersymmetry by systematically adding interaction terms with non-zero extent. To quantize this improvement, we measure boson and fermion masses and Ward identities for the naive as well as the improved models. The masses are degenerate in all models, but the magnitude of the Ward identities decreases significantly for both derivative operators using the improved actions. This is a clear sign that the breaking of supersymmetry due to lattice artifacts is reduced.

High Energy Physics - Lattice,High Energy Physics - Theory
Improved actions and asymptotic scaling in lattice Yang-Mills theory

Kurt Langfeld

DOI: https://doi.org/10.1103/PhysRevD.76.094502

2007-09-05

Abstract:Improved actions in SU(2) and SU(3) lattice gauge theories are investigated with an emphasis on asymptotic scaling. A new scheme for tadpole improvement is proposed. The standard but heuristic tadpole improvement emerges from a mean field approximation from the new approach. Scaling is investigated by means of the large distance static quark potential. Both, the generic and the new tadpole scheme yield significant improvements on asymptotic scaling when compared with loop improved actions. A study of the rotational symmetry breaking terms, however, reveals that only the new improvement scheme efficiently eliminates the leading irrelevant term from the action.

High Energy Physics - Lattice,High Energy Physics - Theory

Molecular dynamics for full QCD simulations with an improved action

Dynamical Simulations with Highly Improved Staggered Quarks

Quantum Simulation of Lattice QCD with Improved Hamiltonians

A Numerical Study of Improved Quark Actions on Anisotropic Lattices

Droems: experimental mathematics, informatics and infinite dimensional geometry

Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Simulation of quantum chromodynamics on the lattice with exactly chiral lattice fermions

Generalized HMC using Nambu mechanics for lattice QCD

Coarsest-level improvements in multigrid for lattice QCD on large-scale computers

Controlling atom-atom interaction at ultralow temperatures by dc electric fields

Calculating the I=2 pion scattering length using tadpole improved clover Wilson action on coarse anisotropic lattices

Reducing Discretization Errors in Lattice QCD Spectroscopy

Full QCD with milder topological freezing

Lattice QCD with exponentially small chirality breaking

BQCD - Berlin quantum chromodynamics program

Spontaneous chiral-symmetry breaking of lattice QCD with massless dynamical quarks

Improved Dirichlet boundary conditions for lattice gauge-fermion theories

Quantum benefit of the quantum equation of motion for the strongly coupled many-body problem

Improved lattice actions for supersymmetric quantum mechanics

Improved actions and asymptotic scaling in lattice Yang-Mills theory