# QCD Phase Transition with Chiral Quarks and Physical Quark Masses

**INSPIRE:** [1282138](https://inspirehep.net/literature/1282138)
**arXiv:** [1402\.5175](https://arxiv.org/abs/1402.5175)
**DOI:** [10\.1103/PhysRevLett\.113\.082001](https://doi.org/10.1103/PhysRevLett.113.082001)

**Authors:** Bhattacharya, Tanmoy, Buchoff, Michael I\., Christ, Norman H\., Ding, H\.\-T\., Gupta, Rajan, Jung, Chulwoo, Karsch, F\., Lin, Zhongjie, Mawhinney, R\.D\., McGlynn, Greg, Mukherjee, Swagato, Murphy, David, Petreczky, P\., Schroeder, Chris, Soltz, R A\., Vranas, P\.M\., Yin, Hantao

**Submitted:** 20 February 2014

**Subjects:**
- hep\-lat
- Lattice
- Phenomenology\-HEP
- Other

**Journal reference:** Phys\.Rev\.Lett\. 113 \(8\) 082001 \(2014\)

**Citations:** 504

## Abstract

\<p\>We report on the first lattice calculation of the QCD phase transition using chiral fermions with physical quark masses\. This calculation uses \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mn\>2\</mml:mn\>\<mml:mo\>\+\</mml:mo\>\<mml:mn\>1\</mml:mn\>\</mml:math\>\</inline\-formula\> quark flavors, spatial volumes between \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mrow\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>4\</mml:mn\>\<mml:mtext\> \</mml:mtext\>\<mml:mtext\> \</mml:mtext\>\<mml:mi\>fm\</mml:mi\>\<mml:msup\>\<mml:mrow\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\</mml:mrow\>\<mml:mrow\>\<mml:mn\>3\</mml:mn\>\</mml:mrow\>\</mml:msup\>\</mml:mrow\>\</mml:math\>\</inline\-formula\> and \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mrow\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>11\</mml:mn\>\<mml:mtext\>  \</mml:mtext\>\<mml:mi\>fm\</mml:mi\>\<mml:msup\>\<mml:mrow\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\</mml:mrow\>\<mml:mrow\>\<mml:mn\>3\</mml:mn\>\</mml:mrow\>\</mml:msup\>\</mml:mrow\>\</mml:math\>\</inline\-formula\> and temperatures between 139 and 196 MeV\. Each temperature is calculated at a single lattice spacing corresponding to a temporal Euclidean extent of \<inline\-formula\>\<mml:math display\="inline"\>\<mml:msub\>\<mml:mi\>N\</mml:mi\>\<mml:mi\>t\</mml:mi\>\</mml:msub\>\<mml:mo\>\=\</mml:mo\>\<mml:mn\>8\</mml:mn\>\</mml:math\>\</inline\-formula\>\. The disconnected chiral susceptibility, \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mrow\>\<mml:msub\>\<mml:mrow\>\<mml:mi\>χ\</mml:mi\>\</mml:mrow\>\<mml:mrow\>\<mml:mtext\>disc\</mml:mtext\>\</mml:mrow\>\</mml:msub\>\</mml:mrow\>\</mml:math\>\</inline\-formula\> shows a pronounced peak whose position and height depend sensitively on the quark mass\. We find no metastability near the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm\. Each result is strong evidence that the QCD “phase transition” is not first order but a continuous crossover for \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mrow\>\<mml:msub\>\<mml:mrow\>\<mml:mi\>m\</mml:mi\>\</mml:mrow\>\<mml:mrow\>\<mml:mi\>π\</mml:mi\>\</mml:mrow\>\</mml:msub\>\<mml:mo\>\=\</mml:mo\>\<mml:mn\>135\</mml:mn\>\<mml:mtext\>  \</mml:mtext\>\<mml:mi\>MeV\</mml:mi\>\</mml:mrow\>\</mml:math\>\</inline\-formula\>\. The peak location determines a pseudocritical temperature \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mrow\>\<mml:msub\>\<mml:mrow\>\<mml:mi\>T\</mml:mi\>\</mml:mrow\>\<mml:mrow\>\<mml:mi\>c\</mml:mi\>\</mml:mrow\>\</mml:msub\>\<mml:mo\>\=\</mml:mo\>\<mml:mn\>155\</mml:mn\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>1\</mml:mn\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>8\</mml:mn\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\<mml:mtext\>  \</mml:mtext\>\<mml:mi\>MeV\</mml:mi\>\</mml:mrow\>\</mml:math\>\</inline\-formula\>, in agreement with earlier staggered fermion results\. However, the peak height is 50% greater than that suggested by previous staggered results\. Chiral \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mi\>S\</mml:mi\>\<mml:mi\>U\</mml:mi\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>2\</mml:mn\>\<mml:msub\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\<mml:mi\>L\</mml:mi\>\</mml:msub\>\<mml:mo\>×\</mml:mo\>\<mml:mi\>S\</mml:mi\>\<mml:mi\>U\</mml:mi\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>2\</mml:mn\>\<mml:msub\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\<mml:mi\>R\</mml:mi\>\</mml:msub\>\</mml:math\>\</inline\-formula\> symmetry is fully restored above 164 MeV, but anomalous \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mi\>U\</mml:mi\>\<mml:mo stretchy\="false"\>\(\</mml:mo\>\<mml:mn\>1\</mml:mn\>\<mml:msub\>\<mml:mo stretchy\="false"\>\)\</mml:mo\>\<mml:mi\>A\</mml:mi\>\</mml:msub\>\</mml:math\>\</inline\-formula\> symmetry breaking is nonzero above \<inline\-formula\>\<mml:math display\="inline"\>\<mml:msub\>\<mml:mi\>T\</mml:mi\>\<mml:mi\>c\</mml:mi\>\</mml:msub\>\</mml:math\>\</inline\-formula\> and vanishes as \<inline\-formula\>\<mml:math display\="inline"\>\<mml:mi\>T\</mml:mi\>\</mml:math\>\</inline\-formula\> is increased to 196 MeV\.\</p\>