Refine
Year of publication
Document Type
- Conference Proceeding (35)
- Article (17)
Language
- English (52)
Has Fulltext
- yes (52)
Is part of the Bibliography
- no (52)
Keywords
- Lattice QCD (3)
- Lattice Quantum Field Theory (2)
- Effective Field Theories (1)
- FOS: Physical sciences (1)
- Finite baryon density (1)
- High Energy Physics - Lattice (hep-lat) (1)
- High Energy Physics - Phenomenology (hep-ph) (1)
- High Energy Physics - Theory (hep-th) (1)
- Hybrid Monte Carlo algorithm (1)
- Lattice field theory (1)
- Lattice gauge field theories (1)
- Lattice gauge theory (1)
- Nuclear Theory (nucl-th) (1)
- Phase Diagram of QCD (1)
- QCD equation of state (1)
- QCD phase diagram (1)
- Quark-gluon Plasmap (1)
- Strong coupling expansion (1)
- chiral symmetry restoration (1)
- cluster expansion model (1)
- finite baryon density (1)
- lattice QCD (1)
Institute
- Physik (52)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- ELEMENTS (2)
- Informatik (1)
The SU(3) spin model with chemical potential corresponds to a simplified version of QCD with static quarks in the strong coupling regime. It has been studied previously as a testing ground for new methods aiming to overcome the sign problem of lattice QCD. In this work we show that the equation of state and the phase structure of the model can be fully determined to reasonable accuracy by a linked cluster expansion. In particular, we compute the free energy to 14-th order in the nearest neighbour coupling. The resulting predictions for the equation of state and the location of the critical end points agree with numerical determinations to O(1%) and O(10%), respectively. While the accuracy for the critical couplings is still limited at the current series depth, the approach is equally applicable at zero and non-zero imaginary or real chemical potential, as well as to effective QCD Hamiltonians obtained by strong coupling and hopping expansions.
Lattice QCD with heavy quarks reduces to a three-dimensional effective theory of Polyakov loops, which is amenable to series expansion methods. We analyse the effective theory in the cold and dense regime for a general number of colours, Nc. In particular, we investigate the transition from a hadron gas to baryon condensation. For any finite lattice spacing, we find the transition to become stronger, i.e. ultimately first-order, as Nc is made large. Moreover, in the baryon condensed regime, we find the pressure to scale as p ∼ Nc through three orders in the hopping expansion. Such a phase differs from a hadron gas with p ∼ N0c, or a quark gluon plasma, p ∼ N2c, and was termed quarkyonic in the literature, since it shows both baryon-like and quark-like aspects. A lattice filling with baryon number shows a rapid and smooth transition from condensing baryons to a crystal of saturated quark matter, due to the Pauli principle, and is consistent with this picture. For continuum physics, the continuum limit needs to be taken before the large Nc limit, which is not yet possible in practice. However, in the controlled range of lattice spacings and Nc-values, our results are stable when the limits are approached in this order. We discuss possible implications for physical QCD.
Lattice constraints on the QCD chiral phase transition at finite temperature and baryon density
(2021)
The thermal restoration of chiral symmetry in QCD is known to proceed by an analytic crossover, which is widely expected to turn into a phase transition with a critical endpoint as the baryon density is increased. In the absence of a genuine solution to the sign problem of lattice QCD, simulations at zero and imaginary baryon chemical potential in a parameter space enlarged by a variable number of quark flavours and quark masses constitute a viable way to constrain the location of a possible non-analytic phase transition and its critical endpoint. In this article I review recent progress towards an understanding of the nature of the transition in the massless limit, and its critical temperature at zero density. Combined with increasingly detailed studies of the physical crossover region, current data bound a possible critical point to μB ≳ 3T.
The nature of the QCD chiral phase transition in the limit of vanishing quark masses has remained elusive for a long time, since it cannot be simulated directly on the lattice and is strongly cutoff-dependent. We report on a comprehensive ongoing study using unimproved staggered fermions with Nf ∈ [2, 8] mass-degenerate flavours on Nτ ∈ {4, 6, 8} lattices, in which we locate the chiral critical surface separating regions with first-order transitions from crossover regions in the bare parameter space of the lattice theory. Employing the fact that it terminates in a tricritical line, this surface can be extrapolated to the chiral limit using tricritical scaling with known exponents. Knowing the order of the transitions in the lattice parameter space, conclusions for approaching the continuum chiral limit in the proper order can be drawn. While a narrow first-order region cannot be ruled out, we find initial evidence consistent with a second-order chiral transition in all massless theories with Nf ≤ 6, and possibly up to the onset of the conformal window at 9 ≲ N∗f ≲ 12. A reanalysis of already published O(a)-improved Nf = 3 Wilson data on Nτ ∈ [4, 12] is also consistent with tricritical scaling, and the associated change from first to second-order on the way to the continuum chiral limit. We discuss a modified Columbia plot and a phase diagram for many-flavour QCD that reflect these possible features.
Recently, an approximate SU(4) chiral spin-flavour symmetry was observed in multiplet patterns of QCD meson correlation functions, in a temperature range above the chiral crossover. This symmetry is larger than the chiral symmetry of massless QCD, and can only arise effectively when colour-electric quark-gluon interactions dynamically dominate the quantum effective action. At temperatures about three times the crossover temperature, these patterns disappear again, indicating the screening of colour-electric interactions, and the expected chiral symmetry is recovered. In this contribution we collect independent evidence for such an intermediate temperature range, based on screening masses and the pion spectral function. Both kinds of observables behave non-perturbatively in this window, with resonance-like peaks for the pion and its first excitation disappearing gradually with temperature. Using symmetry arguments and the known behaviour of screening masses at small densities, we discuss how this chiral spin symmetric band continues into the QCD phase diagram.
In the strong coupling and heavy quark mass regime, lattice QCD dimensionally reduces to effective theories of Polyakov loops depending on the parameters of the original Wilson action β,κ and Nτ. We apply coarse graining techniques to such theories in 1d and 2d, corresponding to lattice QCD at finite temperature and non-zero chemical potential in 1+1d and 2+1d, respectively. In 1d the method is applied to the effective theories up to O(κ4). Using the transfer matrix, the recursion relations are solved analytically. The thermodynamic limit is taken for some observables. Afterwards, continuum extrapolation is performed numerically and results are discussed. In 2d the coarse graining method is applied in the pure gauge and static quark limit. Running couplings are obtained and the fixed points of the transformations are discussed. Finally, the critical coupling of the deconfinement transition is determined in both limits. Agreement to about 12% with Monte Carlo results of 2+1d Yang-Mills theory from the literature is observed.
The order of the chiral phase transition of lattice QCD with unimproved staggered fermions is known to depend on the number of quark flavours, their masses and the lattice spacing. Previous studies in the literature for Nf∈{3,4} show first-order transitions, which weaken with decreasing lattice spacing. Here we investigate what happens when lattices are made coarser to establish contact to the strong coupling region. For Nf∈{4,8} we find a drastic weakening of the transition when going from Nτ=4 to Nτ=2, which is consistent with a second-order chiral transition reported in the literature for Nf=4 in the strong coupling limit. This implies a non-monotonic behaviour of the critical quark or pseudo-scalar meson mass, which separates first-order transitions from crossover behaviour, as a function of lattice spacing.
A lot of effort in lattice simulations over the last years has been devoted to studies of the QCD deconfinement transition. Most state-of-the-art simulations use rooted staggered fermions, while Wilson fermions are affected by large systematic uncertainties, such as coarse lattices or heavy sea quarks. Here we report on an ongoing study of the transition, using two degenerate flavours of nonperturbatively O(a) improved Wilson fermions. We start with Nt = 12 and 16 lattices and pion masses of 600 to 450 MeV, aiming at chiral and continuum limits with light quarks.
We explore the phase diagram of two flavour QCD at vanishing chemical potential using dynamical O(a)-improved Wilson quarks. In the approach to the chiral limit we use lattices with a temporal extent of Nt = 16 and spatial extent L = 32;48 and 64 to enable the extrapolation to the thermodynamic limit with small discretisation effects. In addition to an update on the scans at constant k, reported earlier, we present first results from scans along lines of constant physics at a pion mass of 290 MeV.We probe the transition using the Polyakov loop and the chiral condensate, as well as spectroscopic observables such as screening masses.
The SU(3) pure gauge theory exhibits a first-order thermal deconfinement transition due to spontaneous breaking of its global Z3 center symmetry. When heavy dynamical quarks are added, this symmetry is broken explicitly and the transition weakens with decreasing quark mass until it disappears at a critical point. We compute the critical hopping parameter and the associated pion mass for lattice QCD with Nf=2 degenerate standard Wilson fermions on Nτ∈{6,8,10} lattices, corresponding to lattice spacings a=0.12 fm, a=0.09 fm, a=0.07 fm, respectively. Significant cutoff effects are observed, with the first-order region growing as the lattice gets finer. While current lattices are still too coarse for a continuum extrapolation, we estimate mcπ≈4 GeV with a remaining systematic error of ∼20%. Our results allow us to assess the accuracy of the leading-order and next-to-leading-order hopping expanded fermion determinant used in the literature for various purposes. We also provide a detailed investigation of the statistics required for this type of calculation, which is useful for similar investigations of the chiral transition.
The QCD equation of state at finite baryon density is studied in the framework of a Cluster Expansion Model (CEM), which is based on the fugacity expansion of the net baryon density. The CEM uses the two leading Fourier coefficients, obtained from lattice simulations at imaginary μB, as the only model input and permits a closed analytic form. Excellent description of the available lattice data at both μB = 0 and at imaginary μB is obtained. We also demonstrate how the Fourier coefficients can be reconstructed from baryon number susceptibilities.
Approaching the continuum limit of the deconfinement critical point for Nf=2 staggered fermions
(2022)
Quenched QCD at zero baryonic chemical potential undergoes a first-order deconfinement phase transition at a critical temperature Tc, which is related to the spontaneous breaking of the global center symmetry. The center symmetry is broken explicitly by including dynamical quarks, which weaken the first-order phase transition for decreasing quark masses. At a certain critical quark mass, which corresponds to the Z(2)-critical point, the first-order phase transition turns into a smooth crossover. We investigate the Z(2)-critical quark mass for Nf=2 staggered fermions on Nτ=8,10 lattices, where larger Nτ correspond to finer lattices. Monte-Carlo simulations are performed for several quark mass values and aspect ratios in order to extrapolate to the thermodynamic limit. We present final results for Nτ=8 and preliminary results for Nτ=10 for the critical mass, which are obtained from fitting to a kurtosis finite size scaling formula of the absolute value of the Polyakov loop.
For the exploration of the phase diagram of QCD, effective Polyakov loop theories derived from lattice QCD provide a valuable tool in the heavy quark mass regime. In practice, the evaluation of these theories is complicated by the appearance of long-range and multipoint interaction terms. On the other hand, it is well known that for theories with such kind of interactions mean field approximations can be expected to yield reliable results. Here, we apply this framework to the critical endpoint of the deconfinement transition and results are compared to the literature. This treatment can also be used to investigate the phase diagram at non-zero baryon and isospin chemical potential.
Spectral functions encode a wealth of information about the dynamics of any given system, and the determination of their non-perturbative characteristics is a long-standing problem in quantum field theory. Whilst numerical simulations of lattice QCD provide ample data for various Euclidean correlation functions, the inversion required to extract spectral functions is an ill-posed problem. In this work, we pursue previously established constraints imposed by field locality at finite temperature T, namely that spectral functions possess a non-perturbative representation which generalises the well-known Källén-Lehmann spectral form to T > 0. Using this representation, we analyse lattice QCD data of the spatial pseudo-scalar correlator in the temperature range 220–960 MeV, and obtain an analytic expression for the corresponding spectral function, with parameters fixed by the data. From the structure of this spectral function we find evidence for the existence of a distinct pion state above the chiral pseudo-critical temperature Tpc, and contributions from its first excitation, which gradually melt as the temperature increases. As a non-trivial test, we find that the extracted spectral function reproduces the corresponding temporal lattice correlator data for T = 220 MeV.
In quantum field theories at finite temperature spectral functions describe how particle systems behave in the presence of a thermal medium. Although data from lattice simulations can in principle be used to determine spectral function characteristics, existing methods rely on the extraction of these quantities from temporal correlators, which requires one to circumvent an illposed inverse problem. In these proceedings we report on a recent approach that instead utilises the non-perturbative constraints imposed by field locality to extract spectral function information directly from spatial correlators. In particular, we focus on the application of this approach to lattice QCD data of the spatial pseudo-scalar meson correlator in the temperature range 220−960 MeV, and outline why this data supports the conclusion that there exists a distinct pion state above the chiral pseudo-critical temperature Tpc.
Quenched QCD at zero baryonic chemical potential undergoes a first-order deconfinement phase transition at a critical temperature Tc, which is related to the spontaneous breaking of the global center symmetry. Including heavy, dynamical quarks breaks the center symmetry explicitly and weakens the first-order phase transition. For decreasing quark masses the first-order phase transition turns into a smooth crossover at a Z2-critical point. The critical quark mass corresponding to this point has been examined with Nf=2 Wilson fermions for several Nτ in a recent study within our group. For comparison, we also locate the critical point with Nf=2 staggered fermions on Nτ=8 lattices. For this purpose we perform Monte Carlo simulations for several quark mass values and various aspect ratios in order to extrapolate to the thermodynamic limit. The critical mass is obtained by fitting to a finite size scaling formula of the kurtosis of the Polyakov loop. Our results indicate large discretization effects, requiring simulations on lattices with Nτ>8.
The so-called Columbia plot summarises the order of the QCD thermal transition as a function of the number of quark flavours and their masses. Recently, it was demonstrated that the first-order chiral transition region, as seen for Nf∈[3,6] on coarse lattices, exhibits tricritical scaling while extrapolating to zero on sufficiently fine lattices. Here we extend these studies to imaginary baryon chemical potential. A similar shrinking of the first-order region is observed with decreasing lattice spacing, which again appears compatible with a tricritical extrapolation to zero.
Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interactions in an extended parameter space, and the possibility of model building for Dark Matter and Electro-Weak Symmetry Breaking. Methodological challenges are addressed as well, including new developments in Artificial Intelligence geared towards the identification of different phases and transitions.
The global center symmetry of quenched QCD at zero baryonic chemical potential is broken spontaneously at a critical temperature Tc leading to a first-order phase transition. Including heavy dynamical quarks breaks the center symmetry explicitly and weakens the first-order phase transition for decreasing quark masses until it turns into a smooth crossover at a Z(2)-critical point. We investigate the Z(2)-critical quark mass value towards the continuum limit for Nf=2 flavors using lattice QCD in the staggered formulation. As part of a continued study, we present results from Monte-Carlo simulations on Nτ=8,10 lattices. Several aspect ratios and quark mass values were simulated in order to obtain the critical mass from a fit of the Polyakov loop to a kurtosis finite size scaling formula. Moreover, the possibility to develop a Ginzburg-Landau effective theory around the Z(2)-critical point is explored.