Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userComplex Langevin simulations of a finite density matrix model for QCD
555
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study a random matrix model for QCD at finite density via complex Langevin dynamics. This model has a phase transition to a phase with nonzero baryon density. We study the convergence of the algorithm as a function of the quark mass and the chemical potential and focus on two main observables: the baryon density and the chiral condensate. For simulations close to the chiral limit, the algorithm has wrong convergence properties when the quark mass is in the spectral domain of the Dirac operator. A possible solution of this problem is discussed.
rate research
Read More
We study the Stephanov model, which is an RMT model for QCD at finite density, using the Complex Langevin algorithm. Naive implementation of the algorithm shows convergence towards the phase quenched or quenched theory rather than to intended theory with dynamical quarks. A detailed analysis of this issue and a potential resolution of the failure of this algorithm are discussed. We study the effect of gauge cooling on the Dirac eigenvalue distribution and time evolution of the norm for various cooling norms, which were specifically designed to remove the pathologies of the complex Langevin evolution. The cooling is further supplemented with a shifted representation for the random matrices. Unfortunately, none of these modifications generate a substantial improvement on the complex Langevin evolution and the final results still do not agree with the analytical predictions.
We demonstrate that the complex Langevin method (CLM) enables calculations in QCD at finite density in a parameter regime in which conventional methods, such as the density of states method and the Taylor expansion method, are not applicable due to the severe sign problem. Here we use the plaquette gauge action with $beta = 5.7$ and four-flavor staggered fermions with degenerate quark mass $m a = 0.01$ and nonzero quark chemical potential $mu$. We confirm that a sufficient condition for correct convergence is satisfied for $mu /T = 5.2 - 7.2$ on a $8^3 times 16$ lattice and $mu /T = 1.6 - 9.6$ on a $16^3 times 32$ lattice. In particular, the expectation value of the quark number is found to have a plateau with respect to $mu$ with the height of 24 for both lattices. This plateau can be understood from the Fermi distribution of quarks, and its height coincides with the degrees of freedom of a single quark with zero momentum, which is 3 (color) $times$ 4 (flavor) $times$ 2 (spin) $=24$. Our results may be viewed as the first step towards the formation of the Fermi sphere, which plays a crucial role in color superconductivity conjectured from effective theories.
Monte Carlo studies of QCD at finite density suffer from the sign problem, which becomes easily uncontrollable as the chemical potential $mu$ is increased even for a moderate lattice size. In this work we make an attempt to approach the high density low temperature region by the complex Langevin method (CLM) using four-flavor staggered fermions with reasonably small quark mass on a $8^3 times 16$ lattice. Unlike the previous work on a $4^3 times 8$ lattice, the criterion for correct convergence is satisfied within a wide range of $mu$ without using the deformation technique. In particular, the baryon number density exhibits a plateau behavior consistent with the formation of eight baryons, and it starts to grow gradually at some $mu$.
In this paper we test the complex Langevin algorithm for numerical simulations of a random matrix model of QCD with a first order phase transition to a phase of finite baryon density. We observe that a naive implementation of the algorithm leads to phase quenched results, which were also derived analytically in this article. We test several fixes for the convergence issues of the algorithm, in particular the method of gauge cooling, the shifted representation, the deformation technique and reweighted complex Langevin, but only the latter method reproduces the correct analytical results in the region where the quark mass is inside the domain of the eigenvalues. In order to shed more light on the issues of the methods we also apply them to a similar random matrix model with a milder sign problem and no phase transition, and in that case gauge cooling cooling solves the convergence problems as was shown before in the literature.
Statistical sampling with the complex Langevin (CL) equation is applied to (0+1)-dimensional Thirring model, and its uniform-field variant, at finite fermion chemical potential $mu$. The CL simulation reproduces a crossover behavior which is similar to but actually deviating from the exact solution in the transition region, where we confirm that the CL simulation becomes susceptible to the drift singularities, i.e., zeros of the fermion determinant. In order to simulate the transition region with the CL method correctly, we examine two approaches, a reweighting method and a model deformation, in both of which a single thimble with an attractive fixed point practically covers the integration domain and the CL sampling avoids the determinant zeros. It turns out that these methods can reproduce the correct crossover behavior of the original model with using reference ensembles in the complexified space. However, they need evaluation of the reweighting factor, which scales with the system size exponentially. We discuss feasibility of applying these methods to the Thirring model and to more realistic theories.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
