BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Gert Aarts appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment of this approach in the future.
We present the most recent results from the FASTSUM collaboration for hadron properties at high temperature. This includes the temperature dependence of the light and charmed meson and baryon spectrum, as well as properties of heavy quarkonia. The results are obtained using anisotropic lattices with a fixed scale approach. We also present the status of our next generation gauge ensembles.
We propose inverse renormalization group transformations within the context of quantum field theory that produce the appropriate critical fixed point structure, give rise to inverse flows in parameter space, and evade the critical slowing down effect in calculations pertinent to criticality. Given configurations of the two-dimensional $ϕ^{4}$ scalar field theory on sizes as small as $V=8^{2}$, we apply the inverse transformations to produce rescaled systems of size up to $V'=512^{2}$ which we utilize to extract two critical exponents. We conclude by discussing how the approach is generally applicable to any method that successfully produces configurations from a statistical ensemble and how it can give novel insights into the structure of the renormalization group.
The transition to Euclidean space and the discretization of quantum field theories on spatial or space-time lattices opens up the opportunity to investigate probabilistic machine learning within quantum field theory. Here, we will discuss how discretized Euclidean field theories, such as the $ϕ^{4}$ lattice field theory on a square lattice, are mathematically equivalent to Markov fields, a notable class of probabilistic graphical models with applications in a variety of research areas, including machine learning. The results are established based on the Hammersley-Clifford theorem. We will then derive neural networks from quantum field theories and discuss applications pertinent to the minimization of the Kullback-Leibler divergence for the probability distribution of the $ϕ^{4}$ machine learning algorithms and other probability distributions.
eXplainable Artificial Intelligence (XAI) is a sub-field of Artificial Intelligence (AI) that is at the forefront of AI research. In XAI, feature attribution methods produce explanations in the form of feature importance. People often use feature importance as guidance for intervention. However, a limitation of existing feature attribution methods is that there is a lack of explanation towards the consequence of intervention. In other words, although contribution towards a certain prediction is highlighted by feature attribution methods, the relation between features and the consequence of intervention is not studied. The aim of this paper is to introduce a new framework, called a peer influence framework to look deeper into explanations using graph representation for feature-to-feature interactions to improve the interpretability of black-box Machine Learning models and inform intervention.
We present results from the fastsum collaboration's programme to determine the spectrum of the bottomonium system as a function of temperature. Three different methods of extracting spectral information are discussed: a Maximum Likelihood approach using a Gaussian spectral function for the ground state, the Backus Gilbert method, and the Kernel Ridge Regression machine learning procedure. We employ the fastsum anisotropic lattices with 2+1 dynamical quark flavours, with temperatures ranging from 47 to 375 MeV.
We present a physical interpretation of machine learning functions, opening up the possibility to control properties of statistical systems via the inclusion of these functions in Hamiltonians. In particular, we include the predictive function of a neural network, designed for phase classification, as a conjugate variable coupled to an external field within the Hamiltonian of a system. Results in the two-dimensional Ising model evidence that the field can induce an order-disorder phase transition by breaking or restoring the symmetry, in contrast with the field of the conventional order parameter which causes explicit symmetry breaking. The critical behavior is then studied by proposing a Hamiltonian-agnostic reweighting approach and forming a renormalization group mapping on quantities derived from the neural network. Accurate estimates of the critical point and of the critical exponents related to the operators that govern the divergence of the correlation length are provided. We conclude by discussing how the method provides an essential step toward bridging machine learning and physics.
The precise equivalence between discretized Euclidean field theories and a certain class of probabilistic graphical models, namely the mathematical framework of Markov random fields, opens up the opportunity to investigate machine learning from the perspective of quantum field theory. In this contribution we will demonstrate, through the Hammersley-Clifford theorem, that the $ϕ^{4}$ scalar field theory on a square lattice satisfies the local Markov property and can therefore be recast as a Markov random field. We will then derive from the $ϕ^{4}$ theory machine learning algorithms and neural networks which can be viewed as generalizations of conventional neural network architectures. Finally, we will conclude by presenting applications based on the minimization of an asymmetric distance between the probability distribution of the $ϕ^{4}$ machine learning algorithms and target probability distributions.
New results on bottomonium at nonzero temperature are presented, using the FASTSUM Generation 2L ensembles. Preliminary results for spectral function reconstruction using Kernel Ridge Regression, a machine learning technique, are shown as well and compared to results from the Maximum Entropy Method.
Complex Langevin simulations provide an alternative to sample path integrals with complex weights and therefore are suited to determine the phase diagram of QCD from first principles. We use our proposed method of Dynamic Stabilisation (DS) to ensure improved convergence to the right limit and present new systematic tests of this technique. We also show results on QCD in the limit of heavy quarks and an analysis of DS compared to known results from reweighting.
Baryonic correlation functions provide an ideal tool to study parity doubling and chiral symmetry using lattice simulations. We present a study using $2+1$ flavors of anisotropic Wilson clover fermions on the FASTSUM ensembles and find clear evidence that parity doubling emerges in the quark-gluon plasma. This result is confirmed on the level of spectral functions, which are obtained using a MEM reconstruction. We further highlight the importance of Gaussian smearing in this study.
In the complex Langevin approach to lattice simulations at nonzero density, zeroes of the fermion determinant lead to a meromorphic drift and hence a need to revisit the theoretical derivation. We discuss how poles in the drift affect the formal justification of the approach and then explore the various potential issues in simple models, in a manner that is applicable to heavy dense and full QCD.
At zero temperature the negative-parity ground states of the nucleon and delta baryons are non-degenerate with the positive-parity partners due to spontaneous breaking of chiral symmetry. However, chiral symmetry is expected to be restored at sufficiently high temperature, in particular when going from the hadronic to the quark-gluon plasma (QGP) phase. This would imply that channels with opposite parity become degenerate. We study the nucleon (spin $1/2$) and $Δ$ (spin $3/2$) baryons in both parity sectors using lattice QCD. The range of temperatures spans both the hadronic and QGP phases. Using the FASTSUM anisotropic $N_f = 2 + 1$ ensembles, we analyze the correlation functions and the spectral functions using respectively exponential fits and the Maximum Entropy Method. We find clear evidence of in-medium effects in the hadronic phase, especially for the negative-parity ground state, and of parity doubling in the QGP phase.
Complex Langevin simulations have been able to successfully reproduce results from Monte Carlo methods in the region where the sign problem is mild and make predictions when it is exponentially hard. We present here our study of the QCD phase diagram and the boundary between the confined and deconfined phases in the limit of heavy and dense quarks (HDQCD) for 3 different lattice volumes. We also briefly discuss instabilities encountered in our simulations.
Complex Langevin simulations allow numerical studies of theories that exhibit a sign problem, such as QCD, and are thereby potentially suitable to determine the QCD phase diagram from first principles. Here we study QCD in the limit of heavy quarks for a wide range of temperatures and chemical potentials. Our results include an analysis of the adaptive gauge cooling technique, which prevents large excursions into the non-compact directions of the SL($3, \mathbb{C}$) manifold. We find that such excursions may appear spontaneously and change the statistical distribution of physical observables, which leads to disagreement with known results. Results whose excursions are sufficiently small are used to map the boundary line between confined and deconfined quark phases.
Following a recent lattice study of nucleon parity doubling at finite temperature from the computation of the two-point nucleon correlators, we study the spectral functions of free nucleons at finite temperature. Spectral densities in the continuum are presented along with a comparison to (free) results on the lattice. Particular attention is given to lattice artefacts at higher energies.
We present a lattice QCD calculation of the charge diffusion coefficient, the electrical conductivity and various susceptibilities of conserved charges, for a range of temperatures below and above the deconfinement crossover. The calculations include the contributions from up, down and strange quarks. We find that the diffusion coefficient is of the order of 1/(2πT) and has a dip around the crossover temperature. Our results are obtained with lattice simulations containing 2+1 dynamical flavours on anisotropic lattices. The Maximum Entropy Method is used to construct spectral functions from correlators of the conserved vector current.
Complex Langevin simulations provide an alternative to sample path integrals with complex weights and therefore are suited to determine the phase diagram of QCD from first principles. Adaptive step-size scaling and gauge cooling are used to improve the convergence of our simulations. We present results for the phase diagram of QCD in the limit of heavy quarks and discuss the order of the phase transitions, which are studied by varying the spatial simulation volume.
These lecture notes contain an elementary introduction to lattice QCD at nonzero chemical potential. Topics discussed include chemical potential in the continuum and on the lattice; the sign, overlap and Silver Blaze problems; the phase boundary at small chemical potential; imaginary chemical potential; and complex Langevin dynamics. An incomplete overview of other approaches is presented as well. These lectures are meant for postgraduate students and postdocs with an interest in extreme QCD. A basic knowledge of lattice QCD is assumed but not essential. Some exercises are included at the end.
It is expected that nucleons and their parity partners become degenerate when chiral symmetry is restored. We investigate this question in the context of the thermal transition from the hadronic phase to the quark-gluon plasma, using lattice QCD simulations with N_f=2+1 flavours. We observe a clear sign of parity doubling in the quark-gluon plasma. Besides, we find that the nucleon ground state is, within the uncertainty, largely independent of the temperature, whereas temperature effects are substantial in the negative-parity (N^*) channel, already in the confined phase.
The spectrum of nucleons and their parity partners is studied as a function of temperature spanning the deconfinement transition. We analyse our results using the correlation functions directly, exponential fits in the hadronic phase, and the Maximum Entropy Method. These techniques all indicate that there is degeneracy in the parity partners' channels in the deconfined phase. This is in accordance with the expectation that there is parity doubling and chiral symmetry in the deconfined phase. In the hadronic phase, we also find that the nucleon ground state is largely independent of temperature, whereas there are substantial temperature effects in the negative parity channel. All results are obtained using our FASTSUM 2+1 flavour ensembles.
Some recent developments to handle the numerical sign problem in QCD and related theories at nonzero density are reviewed. In this contribution I focus on changing the integration order to soften the severity of the sign problem, the density of states, and the extension into the complex plane (complex Langevin dynamics and Lefshetz thimbles).
The sign problem of QCD prevents standard lattice simulations to determine the phase diagram of strong interactions with a finite chemical potential directly. Complex Langevin simulations provide an alternative method to sample path integrals with complex weights. We report on our ongoing project to determine the phase diagram of QCD in the limit of heavy quarks (HDQCD) using Complex Langevin simulations.
Monte Carlo methods cannot probe far into the QCD phase diagram with a real chemical potential, due to the famous sign problem. Complex Langevin simulations, using adaptive step-size scaling and gauge cooling, are suited for sampling path integrals with complex weights. We report here on tests of the deconfinement transition in pure Yang-Mills SU(3) simulations and present an update on the QCD phase diagram in the limit of heavy and dense quarks.
The FASTSUM collaboration has initiated a detailed study of thermal QCD using 2+1 flavours of improved Wilson quarks on anisotropic lattices. Spatial volumes of (3fm)^3 and (4fm)^3 are used at fixed cut-off with temperatures ranging from 40 to 350 MeV (corresponding to temporal lattice extents of 128 to 16 lattice units). Results presented here include the deconfinement temperature and a study of the restoration of chiral symmetry, together with a brief summary of our collaboration's other results from these ensembles.
One of the yet unsolved questions of QCD in the context of the Standard Model is to explain the strong CP problem. A way to look for a better understanding of it is to investigate the theory in the presence of a non-zero topological theta term. On the lattice such a term is complex: hence it introduces a sign problem which, in general, limits the applicability of standard Monte Carlo methods. Here we will discuss the approach of complex Langevin dynamics and show results for both real and imaginary values of theta. We also report on our experience with the gradient flow for real and imaginary theta.
A brief overview of the QCD phase diagram at nonzero temperature and density is provided. It is explained why standard lattice QCD techniques are not immediately applicable for its determination, due to the sign problem. We then discuss a selection of recent lattice approaches that attempt to evade the sign problem and classify them according to the underlying principle: constrained simulations (density of states, histograms), holomorphicity (complex Langevin, Lefschetz thimbles), partial summations (clusters, subsets, bags) and change in integration order (strong coupling, dual formulations).
Simulations of QCD with a finite chemical potential typically lead to a severe sign problem, prohibiting any standard Monte Carlo approach. Complex Langevin simulations provide an alternative to sample path integrals with oscillating weight factors and therefore potentially enable the determination of the phase diagram of QCD. Here we present results for QCD in the limit of heavy quarks and show evidence that the phase diagram can be mapped out by direct simulation. We apply adaptive step-size scaling and adaptive gauge cooling to ensure the convergence of these simulations.
We summarise recent progress in simulating QCD at nonzero baryon density using complex Langevin dynamics. After a brief outline of the main idea, we discuss gauge cooling as a means to control the evolution. Subsequently we present a status report for heavy dense QCD and its phase structure, full QCD with staggered quarks, and full QCD with Wilson quarks, both directly and using the hopping parameter expansion to all orders.
The FASTSUM collaboration has been carrying out simulations of N_f=2+1 QCD at nonzero temperature in the fixed-scale approach using anisotropic lattices. Here we present the status of these studies, including recent results for electrical conductivity and charge diffusion, and heavy quarkonium (charm and beauty) physics.
Progress in simulating QCD at nonzero baryon density requires, amongst others, substantial numerical effort. Here we propose two different expansions to all orders in the hopping parameter, preserving the full Yang-Mills action, which are much cheaper to simulate. We carry out simulations using complex Langevin dynamics, both in the hopping expansions and in the full theory, for two flavours of Wilson fermions, and agreement is seen at sufficiently high order in the expansion. These results provide support for the use of complex Langevin dynamics to study QCD at nonzero density, both in the full and the expanded theory, and for the convergence of the latter.
Lefschetz thimbles and complex Langevin dynamics both provide a means to tackle the numerical sign problem prevalent in theories with a complex weight in the partition function, e.g. due to nonzero chemical potential. Here we collect some findings for the quartic model, and for U(1) and SU(2) models in the presence of a determinant, which have some features not discussed before, due to a singular drift. We find evidence for a relation between classical runaways and stable thimbles, and give an example of a degenerate fixed point. We typically find that the distributions sampled in complex Langevin dynamics are related to the thimble(s), but with some important caveats, for instance due to the presence of unstable fixed points in the Langevin dynamics.
We use the heavy dense formulation of QCD (HD-QCD) as the basis for an analytic expansion as systematic approximation to QCD at non-zero density, keeping the full Yang-Mills action. We analyse the structure of the baryonic density and other quantities and present data from the complex Langevin equation (CLE) and reweighting (RW) calculations for 2 flavours of Wilson fermions.
We present results on the bottomonium spectrum at temperatures above and below the deconfinement crossover temperature, $T_c$, from dynamical lattice QCD simulations. The heavy quark is treated with a non-relativistic effective field theory on the lattice and serves as a probe of the hot medium. Ensembles with a finer spatial lattice spacing and a greater range of temperatures below $T_c$ than those previously employed by this collaboration are used. In addition, there are $N_f=2+1$ flavours of Wilson clover quark in the sea with $M_π\approx400$ MeV and we perform a more careful tuning of the bottom quark mass in this work. We calculate the spectral functions of S and P wave bottomonium states using the maximum entropy method and confirm earlier findings on the survival of the ground state S wave states up to at least $2T_c$ and the immediate dissociation of the P wave states above $T_c$.
In the case of nonabelian gauge theories with a complex weight, a controlled exploration of the complexified configuration space during a complex Langevin process requires the use of SL(N,C) gauge cooling, in order to minimize the distance from SU(N). Here we show that adaptive gauge cooling can lead to an efficient implementation of this idea. First results for SU(3) Yang-Mills theory in the presence of a nonzero theta-term are presented as well.
We investigate the modification of S and P wave states in the bottomonium spectrum above and below the deconfinement crossover temperature through their spectral functions obtained from the maximum entropy method. Anisotropic ensembles with $N_f=2+1$ Wilson clover fermions with tadpole-improvement are used while the bottom quark is treated with an improved non-relativistic action. We observe significant modifications of the P wave $χ_{b1}$ ground state directly above the crossover temperature, $T_c$, while the S wave $Υ$ ground state remains relatively unchanged up to temperatures of almost $2T_c$. This evidence supports earlier conclusions from our FASTSUM collaboration of the immediate dissociation of the P wave states above $T_c$ and the survival of the S wave ground states up to $2T_c$.
I review the presence of the sign problem in lattice QCD at nonzero baryon density and its relation with the overlap and Silver Blaze problems. I then discuss progress in some cases where the sign problem can be handled, either because the sign problem is absent or because it is milder than in full QCD. Some time is spent on effective three-dimensional models, which can be treated with a variety of methods. I conclude with a discussion of the applicability of complex Langevin dynamics at nonzero density.
We first test the Complex Langevin method (CLE) on various simple models. We then introduce the method of Gauge Cooling to control the dynamics of the process and ensure thin distributions in the imaginary direction. We finally apply CLE with gauge cooling to a QCD-related lattice model (HQCD) and compare the results by CLE and by a refined Reweighting method (rRW). Very good agreement is found in all regions of physical interest.
At nonzero chemical potential the numerical sign problem in lattice field theory limits the use of standard algorithms based on importance sampling. Complex Langevin dynamics provides a possible solution, but it has to be applied with care. In this review, we first summarise our current understanding of the approach, combining analytical and numerical insight. In the second part we study SL(N,C) gauge cooling, which was introduced recently as a tool to control complex Langevin dynamics in nonabelian gauge theories. We present new results in Polyakov chain models and in QCD with heavy quarks and compare various adaptive cooling implementations.
The FASTSUM Collaboration presents its first results for the electric charge susceptibility in QCD using 2+1 dynamical flavours of Wilson quark on anisotropic lattices. Spatial volumes of (2.94 fm)^3 are used at fixed cut-off with temperatures ranging from below T_c to ~2 T_c.
A lattice calculation is presented for the electrical conductivity of the QCD plasma with 2+1 dynamical flavours at nonzero temperature. We employ the conserved lattice current on anisotropic lattices using a tadpole-improved clover action and study the behaviour of the conductivity over a wide range of temperatures, both below and above the deconfining transition. The conductivity is extracted from a spectral-function analysis using the Maximal Entropy Method and a discussion of its systematics is provided.
Lattice field theories with a complex action can be studied numerically by allowing a complexified configuration space to be explored. Here we compare the recently introduced formulation on a Lefschetz thimble with the result from stochastic quantisation (or complex Langevin dynamics) in the case of a simple model and contrast the distributions being sampled. We also study the role of the residual phase on the Lefschetz thimble.
Complex Langevin dynamics can solve the sign problem appearing in numerical simulations of theories with a complex action. In order to justify the procedure, it is important to understand the properties of the real and positive distribution, which is effectively sampled during the stochastic process. In the context of a simple model, we study this distribution by solving the Fokker-Planck equation as well as by brute force and relate the results to the recently derived criteria for correctness. We demonstrate analytically that it is possible that the distribution has support in a strip in the complexified configuration space only, in which case correct results are expected.
Complex Langevin dynamics can be used to perform numerical simulations of theories with a complex action. In order to justify the procedure, it is important to understand the properties of the real and positive distribution, which is effectively sampled during the stochastic process. In the context of a simple model, we study this distribution and relate the results to the recently derived criteria for correctness. We demonstrate analytically that if the distribution has support only on a strip in the complexified configuration space, correct results are expected.
We present an overview of bottomonium spectral functions in the quark-gluon plasma, obtained by the FASTSUM collaboration, using lattice QCD simulations with two light quark flavours on anisotropic lattices. The bottom quark is treated nonrelativistically. While we find that the S wave ground states survive up to the highest temperature we consider, we have strong indications that P wave states melt immediately above Tc.
The FASTSUM Collaboration has calculated several quantities relevant for QCD studies at non-zero temperature using the lattice technique. We report here our results for the (i) interquark potential in charmonium; (ii) bottomonium spectral functions; and (iii) electrical conductivity. All results were obtained with 2+1 flavours of dynamical fermions on an anisotropic lattice which allows greater resolution in the temporal direction.
We extend our study of bottomonium spectral functions in the quark-gluon plasma to nonzero momentum. We use lattice QCD simulations with two flavours of light quark on highly anisotropic lattices and treat the bottom quark with nonrelativistic QCD (NRQCD). We focus on S wave (Upsilon and eta_b) channels and consider nonrelativistic velocities, v/c < 0.2. A comparison with predictions from effective field theory is made.
The sign problem at nonzero chemical potential prohibits the use of importance sampling in lattice simulations. Since complex Langevin dynamics does not rely on importance sampling, it provides a potential solution. Recently it was shown that complex Langevin dynamics fails in the disordered phase in the case of the three-dimensional XY model, while it appears to work in the entire phase diagram in the case of the three-dimensional SU(3) spin model. Here we analyse this difference and argue that it is due to the presence of the nontrivial Haar measure in the SU(3) case, which has a stabilizing effect on the complexified dynamics. The freedom to modify and stabilize the complex Langevin process is discussed in some detail.
The FASTSUM collaboration presents a study on the temperature dependence of the electrical conductivity $σ$ in the quark-gluon plasma, using the methods of lattice QCD. Correlators of the exactly conserved vector current are measured at different temperatures across the deconfinement transition, using ensembles of $2+1$ flavours of dynamical fermions on anisotropic lattices. We then employ bayesian methods (MEM) to extract the relevant spectral functions, which are found to be consistent with $σ/T$ rising as a function of $T$. The robustness of the results is verified by a detailed analysis of the systematics involved in the bayesian reconstruction of the spectral functions.
The three-dimensional SU(3) spin model is an effective Polyakov loop model for QCD at nonzero temperature and density. It suffers from a sign problem at nonzero chemical potential. We revisit this model using complex Langevin dynamics and assess in particular the justification of this approach, using analyticity at small mu^2 and the criteria for correctness developed recently. Finite-stepsize effects are discussed in some detail and a higher-order algorithm is employed to eliminate leading stepsize corrections. Our results strongly indicate that complex Langevin dynamics is reliable in this theory in both phases, including the critical region. This is in sharp contrast to the case of the XY model, where correct results were obtained in only part of the phase diagram.
Using Non-Relativistic QCD (NRQCD), we study heavy quark mass dependence of S-wave and P-wave bottomonium correlators for 0.42Tc <= T <= 2.09Tc and study spectral functions of S-wave bottomonium states moving in a thermal bath at these temperatures using Maximum Entropy Method with NRQCD kernel. For the studied momentum range, the energy of moving states shows quadratic momentum-dependence and the width of moving states does not show significant changes as the momentum of bottomonium is increased. Also, we find that in correlator ratios, the temperature effect is larger than the effect caused by 20% change in the bottom quark mass.
The temperature dependence of bottomonium states at temperatures above and below $T_c$ is presented, using non-relativistic dynamics for the bottom quark and full relativistic lattice QCD simulations for two light flavors on a highly anisotropic lattice. We find that the S-waves ($Υ$ and $η_b$) show little temperature dependence in this range while the P wave propagators show a crossover from the exponential decay characterizing the hadronic phase to a power-law behavior consistent with nearly-free dynamics at approximately twice the critical temperature.
The complex Langevin method is a leading candidate for solving the sign problem occurring in various physical situations, notably QCD at finite chemical potential. Its most vexing problem is `convergence to the wrong limit', where the simulation gives a well defined, but incorrect, result. Here, we first outline a formal justification of the method and identify points at which it might fail. From these we derive a condition that must be satisfied in order for correct results to be obtained. We then apply these ideas to the three-imensional SU(3) spin model at finite chemical potential and show strong indications that complex Langevin dynamics yields correct results in this theory.
The complex Langevin method is a leading candidate for solving the so-called sign problem occurring in various physical situations. Its most vexing problem is that in some cases it produces `convergence to the wrong limit'. In the first part of the paper we go through the formal justification of the method, identify points at which it may fail and identify a necessary and sufficient criterion for correctness. This criterion would, however, require checking infinitely many identities, and therefore is somewhat academic. We propose instead a truncation to the check of a few identities; this still gives a necessary criterion, but a priori it is not clear whether it remains sufficient. In the second part we carry out a detailed study of two toy models: first we identify the reasons why in some cases the method fails, second we test the efficiency of the truncated criterion and find that it works perfectly at least in the toy models studied.
We study bottomonium spectral functions in the quark-gluon plasma in the Upsilon and eta_b channels, using lattice QCD simulations with two flavours of light quark on highly anisotropic lattices. The bottom quark is treated with nonrelativistic QCD (NRQCD). In the temperature range we consider, 0.42 < T/T_c < 2.09, we find that the ground states survive, whereas the excited states are suppressed as the temperature is increased. The position and width of the ground states are compared to analytical effective field theory (EFT) predictions. Systematic uncertainties of the maximum entropy method (MEM), used to construct the spectral functions, are discussed in some detail.
Stochastic quantization offers the opportunity to simulate field theories with a complex action. In some theories unstable trajectories are prevalent when a constant stepsize is employed. We construct algorithms for generating an adaptive stepsize in complex Langevin simulations and find that unstable trajectories are completely eliminated. To illustrate the generality of the approach, we apply it to the three-dimensional XY model at nonzero chemical potential and the heavy dense limit of QCD.
We demonstrate analytically that complex Langevin dynamics can solve the sign problem in one-dimensional QCD in the thermodynamic limit. In particular, it is shown that the contributions from the complex and highly oscillating spectral density of the Dirac operator to the chiral condensate are taken into account correctly. We find an infinite number of classical fixed points of the Langevin flow in the thermodynamic limit. The correct solution originates from a continuum of degenerate distributions in the complexified space.
We investigate N=4 SYM coupled to fundamental flavours at nonzero imaginary quark chemical potential in the strong coupling and large N limit, using gauge/gravity duality applied to the D3-D7 system, treating flavours in the probe approximation. The interplay between Z(N) symmetry and the imaginary chemical potential yields a series of first-order Roberge-Weiss transitions. An additional thermal transition separates phases where quarks are bound/unbound into mesons. This results in a set of Roberge-Weiss endpoints: we establish that these are triple points, determine the Roberge-Weiss temperature, give the curvature of the phase boundaries and confirm that the theory is analytic in mu^2 when mu^2~0.
It has long been known that QCD undergoes a deconfining phase transition at high temperature. One of the consequent features of this new, quark-gluon phase is that hadrons become unbounded. In this talk meson correlation functions at non-zero momentum are studied in the deconfined phase using the Maximum Entropy Method.
The three-dimensional XY model is studied at finite chemical potential using complex Langevin dynamics. The validity of the approach is probed at small chemical potential using imaginary chemical potential and continuity arguments, and at larger chemical potential by comparison with the world line method. While complex Langevin works for larger beta, we find that it fails for smaller beta, in the region of the phase diagram corresponding to the disordered phase. Diagnostic tests are developed to identify symptoms correlated with incorrect convergence. We argue that the erroneous behaviour at smaller beta is not due to the sign problem, but rather resembles dynamics observed in complex Langevin simulations of simple models with complex noise.
We analyze to what extent the complex Langevin method, which is in principle capable of solving the so-called sign problems, can be considered as reliable. We give a formal derivation of the correctness and then point out various mathematical loopholes. The detailed study of some simple examples leads to practical suggestions about the application of the method.
The three-dimensional XY model is studied at finite chemical potential using complex Langevin dynamics. An adaptive stepsize algorithm is implemented to cure the problem of runaway solutions that appears when using a constant stepsize. The validity of complex Langevin dynamics is tested against calculations using imaginary chemical potential and the world line formalism. While complex Langevin dynamics is found to work correctly at larger beta, it fails for smaller beta, in the region of the phase diagram corresponding to the disordered phase. Diagnostic tests are developed to identify behaviour symptomatic of incorrect convergence. These indicate that the failure does not depend on the severeness of the sign problem, but has a different origin.
I answer the question in the title for the relativistic Bose gas at finite chemical potential using numerical lattice simulations, complemented with analytical understanding.
Stochastic quantization can potentially be used to simulate theories with a complex action due to a nonzero chemical potential. We study complex Langevin dynamics in the relativistic Bose gas analytically, using a mean field approximation. We concentrate on the region with a Silver Blaze problem and discuss convergence, stability, fixed points, and the severeness of the sign problem. The real distribution satisfying the extended Fokker-Planck equation is constructed and its nonlocal form is explained. Finally, we compare the mean field results in finite volume with the numerical data presented in Ref. [1].
Lattice QCD at finite chemical potential is difficult due to the sign problem. We use stochastic quantization and complex Langevin dynamics to study this issue. First results for QCD in the hopping expansion are encouraging. U(1) and SU(3) one link models are used to gain further insight into why the method appears to be successful.
In order to enhance our understanding of spectral functions in lattice QCD obtained with the help of the Maximum Entropy Method, we study meson spectral functions for lattice fermions with chiral symmetry. In particular we analyse lattice artefacts for standard overlap, overlap hypercube and domain wall fermions in the free field limit. We also present first results for pseudoscalar spectral functions in dynamical QCD with 2+1 flavours of domain wall fermions, using data generated by the UKQCD and RBC collaborations on QCDOC machines.
We report first results for spectral functions of charmonium in 2-flavour QCD. The spectral functions are determined from vector and pseudoscalar correlators on a dynamical, anisotropic lattice. J/psi and eta_c are found to survive well into the deconfined phase before melting away at T<~2T_c. Current systematic uncertainties prevent us from drawing any definite conclusions at this stage.
In the high-temperature phase of QCD meson spectral functions at nonzero momentum are expected to have a nontrivial and interesting structure. In order to provide a reference point for lattice studies employing e.g. the Maximal Entropy Method, we discuss several characteristics of meson spectral functions in the infinite-temperature limit. We report on ongoing work in quenched QCD with staggered fermions.
We study the κ^2 corrections to the quenched limit of μ>0 QCD. We use an improved reweighting procedure.
We consider the time evolution of nonequilibrium quantum scalar fields in the O(N) model, using the next-to-leading order 1/N expansion of the 2PI effective action. A comparison with exact numerical simulations in 1+1 dimensions in the classical limit shows that the 1/N expansion gives quantitatively precise results already for moderate values of N. For sufficiently high initial occupation numbers the time evolution of quantum fields is shown to be accurately described by classical physics. Eventually the correspondence breaks down due to the difference between classical and quantum thermal equilibrium.
For a weakly coupled quantum field at high temperature the classical approximation offers a possibility to gain insight into nonperturbative real-time dynamics. I use this to present a nonperturbative approach to the computation of spectral functions in real time. Results are shown for a scalar field in 2+1 dimensions.