BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Denes Sexty contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
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.
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.
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.
Recent progress of the complex Langevin method and the Lefschetz thimble in connection with the sign problem is reviewed. These methods rely on the complexification of the original field manifold and they allow direct simulations of theories with non-real measures. Similarities and differences of the two approaches are pointed out. Results using the complex Langevin method, which allows simulations to evade the sign problem in full QCD, are presented. Promising results of the thimble approach for non-gauge theories are also discussed.
Progress in the application of the complex Langevin method to full QCD at non-zero chemical potential is reported. The method evades the sign problem which makes naive simulations at nonzero density impossible. The procedure 'gauge cooling' is used to stabilize the simulations at small enough lattice spacings. The method allows simulations also at high densities, all the way up to saturation. Simulations in a systematic hopping parameter expansion are also performed and good convergence is observed, validating the full as well as the expanded 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 complex Langevin method is extended to full QCD at non-zero chemical potential. The use of gauge cooling stabilizes the simulations at small enough lattice spacings. At large fermion mass the results are compared to the HQCD approach, in which the spatial hoppings of fermionic variables are neglected, and good agreement is found. The method allows simulations also at high densities, all the way up to saturation.
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.
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 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.
Simulations of full QCD at nonzero baryon density using light quark masses are presented. The sign problem is evaded by the usage of the complex Langevin equation. The simulations are stabilized by the gauge cooling procedure for small lattice spacings. The method allows simulations at high densities, up to the saturation. The sign average is measured in the full as well as the phasequenched theory. Results are compared to the HQCD approach, in which the spatial hopping terms of fermionic variables are dropped, and good agreement is found at large masses.
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.
New aspects of parametrically resonant heating of a relativistic scalar O(2)-symmetric self-interacting field are presented. This process is a candidate for reheating at the end of the early-universe epoch of inflation. Although a model with a fully O(2)-symmetric ground state is used, transient, metastable spontaneous symmetry breaking can be observed. This manifests itself in the form of persistent regimes of opposite and, inside these, uniform charge overdensities separated by thin lines and walls similar to topological defects, in two and three spatial dimensions, respectively. The configuration is found to correspond to an attractive non-equilibrium fixed point of the underlying dynamic equations which prevents thermalisation over an extended period of time.
The non-equilibrium time evolution of an Anderson quantum dot is investigated. The quantum dot is coupled between two leads forming a chemical-potential gradient. We use Kadanoff-Baym dynamic equations within a non-perturbative resummation of the s-channel bubble chains. The effect of the resummation leads to the introduction of a frequency-dependent 4-point vertex. The tunneling to the leads is taken into account exactly. The method allows the determination of the transient as well as stationary transport through the quantum dot, and results are compared with different schemes discussed in the literature (fRG, ISPI, tDMRG and QMC).