Dynamical stabilisation of complex Langevin simulations of QCD

Abstract

The ability to describe strongly interacting matter at finite temperature and baryon density provides the means to determine, for instance, the equation of state of QCD at non-zero baryon chemical potential. From a theoretical point of view, direct lattice simulations are hindered by the numerical sign problem, which prevents the use of traditional methods based on importance sampling. Despite recent successes, simulations using the complex Langevin method have been shown to exhibit instabilities, which cause convergence to wrong results. We introduce and discuss the method of dynamic stabilisation (DS), a modification of the complex Langevin process aimed at solving these instabilities. We present results of DS being applied to the heavy-dense approximation of QCD, as well as QCD with staggered fermions at zero chemical potential and finite chemical potential at high temperature. Our findings show that DS can successfully deal with the aforementioned instabilities, opening the way for further progress.