No Arabic abstract
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 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.
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.
We explore the QCD phase diagram at finite density with four-flavor staggered fermions using the complex Langevin method, which is a promising approach to overcome the sign problem. In our previous work on an $8^3 times 16$ lattice at $beta = 5.7$ with the quark mass $m = 0.01$, we have found that the baryon number density has a clear plateau as a function of the chemical potential. In this study, we use a $16^3 times 32$ lattice to reduce finite volume effects and find that the plateau structure survives. Moreover, the number of quarks in the plateau region turns out to be 24, which is exactly the same as the one obtained previously on the $8^3 times 16$ lattice. We provide a simple interpretation of this number, which suggests that the Fermi sphere is starting to form.
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.