BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Ion-Olimpiu Stamatescu 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
A formula expressing the fermionic determinant as an infinite product of smaller determinants is derived and discussed. These smaller determinants are of a fixed size, independent of the size of the lattice and are indexed by loops of increasing length.
A formula expressing the fermionic determinant (a large order polynomial) as an infinite product of smaller determinants is derived and discussed. These smaller determinants are of a fixed size, independent of the size of the lattice and are indexed by loops of increasing length.
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.
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.
In this paper we discuss one-dimensional models reproducing some features of quantum electrodynamics and quantum chromodynamics at nonzero density and temperature. Since a severe sign problem makes a numerical treatment of QED and QCD at high density difficult, such models help to explore various effects peculiar to the full theory. Studying them gives insights into the large density behavior of the Polyakov loop by taking both bosonic and fermionic degrees of freedom into account, although in one dimension only the implementation of a global gauge symmetry is possible. For these models we evaluate the respective partition functions and discuss several observables as well as the Silver Blaze phenomenon.
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.
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.
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.
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.
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.
We employ a new method, "gauge cooling", to stabilize complex Langevin simulations of QCD with heavy quarks. The results are checked against results obtained with reweigthing; we find agreement within the estimated errors. The method allows us to go to previously unaccessible high densities.
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.
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 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.
We calculate Landau gauge ghost and gluon propagators in pure SU(2) lattice gauge theory in two, three and four dimensions. The gauge fixing method we use, sc. stochastic quantisation, serves as a viable alternative approach to standard gauge fixing algorithms. We also investigate the spectrum of the Faddeev-Popov operator. At insufficiently accurate gauge fixing, we find evidence that stochastic quantisation samples configurations close to the Gribov horizon. Standard gauge fixing does so only at specific parameters; otherwise, there is a clear difference. However, this difference disappears if the gauge is fixed to sufficient accuracy. In this case, we confirm previous lattice results for the gluon and ghost propagator in two, three and four dimensions.
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.
We study the high density region of QCD within an effective model obtained in the frame of the hopping parameter expansion. The model still acknowledges the sign problem peculiar to non-zero chemical potential, but it permits the development of refined algorithms which ensure a good overlap of the Monte Carlo ensemble with the true one. We review the main features of the model, including the most explicit form of the resumed expansion, and present calculations concerning the dependence of various observables on the chemical potential and on the temperature, in particular of the charge density and the diquark susceptibility, which may be used to characterize the various phases expected at high baryonic density.
We discuss the contribution of OZI-suppressed diagrams to the hyperfine splitting of charmonium in lattice QCD. We study valence quark mass regions from strange to charm quark masses. No contribution of the disconnected diagram is seen in the vector meson channel. In the pseudo-scalar channel and for valence quark masses around the strange quark, the disconnected contribution induces a considerable increase of the meson mass. This contribution quickly decreases as the quark mass increases. For charmonium the effect is very small although a decrease of the pseudoscalar mass induced by the disconnected contribution cannot be ruled out.
We study the κ^2 corrections to the quenched limit of μ>0 QCD. We use an improved reweighting procedure.
We use a cooling algorithm based on an improved action with scale invariant instanton solutions, which needs no monitoring or calibration and has a inherent cut off for dislocations. We present results for SU(2) Yang-Mills theory where the method provides good susceptibility data and physical size distributions of instantons.