-
SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations
Authors:
Lukas Mazur,
Dennis Bollweg,
David A. Clarke,
Luis Altenkort,
Olaf Kaczmarek,
Rasmus Larsen,
Hai-Tao Shu,
Jishnu Goswami,
Philipp Scior,
Hauke Sandmeyer,
Marius Neumann,
Henrik Dick,
Sajid Ali,
Jangho Kim,
Christian Schmidt,
Peter Petreczky,
Swagato Mukherjee
Abstract:
The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HIS…
▽ More
The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ fermions in QCD studies. This situation encourages the development of software featuring a HISQ action that is flexible, high-performing, open source, easy to use, and easy to adapt. In this technical paper, we explain the design strategy, provide implementation details, list available algorithms and modules, and show key performance indicators for SIMULATeQCD, a simple multi-GPU lattice code for large-scale QCD calculations, mainly developed and used by the HotQCD collaboration. The code is publicly available on GitHub.
△ Less
Submitted 3 April, 2024; v1 submitted 1 June, 2023;
originally announced June 2023.
-
Scaling functions of the three-dimensional $Z(2)$, $O(2)$ and $O(4)$ models and their finite size dependence in an external field
Authors:
Frithjof Karsch,
Marius Neumann,
Mugdha Sarkar
Abstract:
We analyze scaling functions in the $3$-$d$, $Z(2)$, $O(2)$ and $O(4)$ universality classes and their finite size dependence using Monte Carlo simulations of improved $φ^4$ models. Results for the scaling functions are fitted to the Widom-Griffiths form, using a parametrization also used in analytic calculations. We find good agreement on the level of scaling functions and the location of maxima i…
▽ More
We analyze scaling functions in the $3$-$d$, $Z(2)$, $O(2)$ and $O(4)$ universality classes and their finite size dependence using Monte Carlo simulations of improved $φ^4$ models. Results for the scaling functions are fitted to the Widom-Griffiths form, using a parametrization also used in analytic calculations. We find good agreement on the level of scaling functions and the location of maxima in the universal part of susceptibilities. We also find that an earlier parametrization of the $O(4)$ scaling function, using 14 parameters, is well reproduced when using the Widom-Griffiths form with only three parameters. We furthermore show that finite size corrections to the scaling functions are distinctively different in the $Z(2)$ and $O(N)$ universality classes and determine the volume dependence of the peak locations in order parameter and mixed susceptibilities.
△ Less
Submitted 17 July, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
A machine learning approach to the classification of phase transitions in many flavor QCD
Authors:
Frithjof Karsch,
Anirban Lahiri,
Marius Neumann,
Christian Schmidt
Abstract:
Normalizing flows are generative machine learning models which can efficiently approximate probability distributions, using only given samples of a distribution. This architecture is used to interpolate the chiral condensate obtained from QCD simulations with five degenerate quark flavors in the HISQ action. From this a model for the probability distribution of the chiral condensate as function of…
▽ More
Normalizing flows are generative machine learning models which can efficiently approximate probability distributions, using only given samples of a distribution. This architecture is used to interpolate the chiral condensate obtained from QCD simulations with five degenerate quark flavors in the HISQ action. From this a model for the probability distribution of the chiral condensate as function of lattice volume, quark mass and gauge coupling is obtained. Using the model, first order and crossover regions can be classified and the boundary between these regions can be marked by a critical mass. An extension of this model to studies of phase transitions in QCD with variable number of flavors is expected to be possible.
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
Towards the chiral phase transition in the Roberge-Weiss plane
Authors:
F. Cuteri,
J. Goswami,
F. Karsch,
Anirban Lahiri,
M. Neumann,
O. Philipsen,
Christian Schmidt,
A. Sciarra
Abstract:
We discuss the interplay between chiral and center sector phase transitions that occur in QCD with an imaginary quark chemical potential $μ=i(2n+1) πT/3$. Based on a finite size scaling analysis in (2+1)-flavor QCD using HISQ fermions with a physical strange quark mass and a range of light quark masses, we show that the endpoint of the line of first-order Roberge-Weiss (RW) transitions between cen…
▽ More
We discuss the interplay between chiral and center sector phase transitions that occur in QCD with an imaginary quark chemical potential $μ=i(2n+1) πT/3$. Based on a finite size scaling analysis in (2+1)-flavor QCD using HISQ fermions with a physical strange quark mass and a range of light quark masses, we show that the endpoint of the line of first-order Roberge-Weiss (RW) transitions between center sectors is second order for light quark masses $m_l\ge m_s/320$, and that it belongs to the $3$-d, $Z(2)$ universality class. The operator for the chiral condensate behaves like an energy-like operator in an effective spin model for the RW phase transition. As a consequence, for any non-zero value of the quark mass, the chiral condensate will have an infinite slope at the RW phase transition temperature, $T_{RW}$. Its fluctuation, the disconnected chiral susceptibility, behaves like the specific heat in $Z(2)$ symmetric models and diverges in the infinite volume limit at the RW phase transition temperature $T_{RW}$ for any non-zero value of the light quark masses. Our analysis suggests the critical temperatures for the RW phase transition and the chiral phase transition coincide in the RW plane. On lattices with temporal extent $N_τ=4$, we find in the chiral limit $T_χ=T_{RW}=195(1)~$MeV.
△ Less
Submitted 25 May, 2022;
originally announced May 2022.
-
Lattice QCD at Imaginary Chemical Potential in the Chiral Limit
Authors:
D. A. Clarke,
Jishnu Goswami,
F. Karsch,
Anirban Lahiri,
M. Neumann,
C. Schmidt
Abstract:
We report on an ongoing study on the interplay between Roberge-Weiss (RW) and chiral transitions in simulations of (2+1)-flavor QCD with an imaginary chemical potential. We established that the RW endpoint belongs to the 3-$d$, $Z_2$ universality class when calculations are done with the Highly Improved Staggered Quark (HISQ) action in the RW plane with physical quark masses. We also have explored…
▽ More
We report on an ongoing study on the interplay between Roberge-Weiss (RW) and chiral transitions in simulations of (2+1)-flavor QCD with an imaginary chemical potential. We established that the RW endpoint belongs to the 3-$d$, $Z_2$ universality class when calculations are done with the Highly Improved Staggered Quark (HISQ) action in the RW plane with physical quark masses. We also have explored a range of quark masses corresponding to pion mass values, $m_π\geq40$~MeV and found that the transition is consistent with $Z_2$ universality class. We argue that observables that were usually used to determine the chiral phase transition temperature, e.g. the chiral condensate and chiral susceptibility, are sensitive to the RW transition and are energy-like observables for the $Z_2$ transition, contrary to the magnetic-like (order parameter) behavior at vanishing chemical potential. Moreover the calculations performed at $m_π\sim40$~MeV also put a stringent constraint for a critical pion mass at zero chemical potential for a possible first-order chiral phase transition.
△ Less
Submitted 30 November, 2021;
originally announced November 2021.
-
Critical end points in (2+1)-flavor QCD with imaginary chemical potential
Authors:
Jishnu Goswami,
Frithjof Karsch,
Anirban Lahiri,
Marius Neumann,
Christian Schmidt
Abstract:
We present here the results from an ongoing determination of the critical quark mass in simulations of (2+1)-flavor QCD with an imaginary chemical potential. Studies with unimproved actions found the existence of a critical quark mass value at which the crossover transition ends on a second order phase transition and becomes first order for smaller values of the quark mass for the case of both van…
▽ More
We present here the results from an ongoing determination of the critical quark mass in simulations of (2+1)-flavor QCD with an imaginary chemical potential. Studies with unimproved actions found the existence of a critical quark mass value at which the crossover transition ends on a second order phase transition and becomes first order for smaller values of the quark mass for the case of both vanishing and imaginary chemical potential. We use the Highly Improved Staggered Quark (HISQ) action and perform calculations in the Roberge-Weiss (RW) plane, where the value of the critical mass is expected to be largest. The lowest quark mass value used in our simulation corresponds to the pion mass $m_π$, down to $40$ MeV. Contrary to calculations performed with unimproved actions we find no evidence for the occurrence of first order transitions at the smallest quark mass values explored so far. Moreover we also show that the chiral observables are sensitive to the RW transition. Our results also indicate that the RW transition and chiral transition could coincide in the chiral limit.
△ Less
Submitted 8 May, 2019;
originally announced May 2019.