Inhomogeneous confinement and chiral symmetry breaking induced by imaginary angular velocity
Abstract
We investigate detailed properties of imaginary rotating matter with gluons and quarks at high temperature. Previously, we showed that imaginary rotation induces perturbative confinement of gluons at the rotation center. We perturbatively calculate the Polyakov loop potential and find inhomogeneous confinement above a certain threshold of imaginary angular velocity. We also evaluate the quark contribution to the Polyakov loop potential and confirm that spontaneous chiral symmetry breaking occurs in the perturbatively confined phase.
keywords:
Confinement, Rotation, Chiral Symmetry Breaking, Inhomogeneous States[1]organization=School of Physics and Astronomy, University of Minnesota, addressline=Minneapolis, city=MN, postcode=55455, country=U.S.A.
[2]organization=Department of Physics, The University of Tokyo, addressline=7-3-1 Hongo, Bunkyo-ku, city=Tokyo, postcode=113-0033, country=Japan
1 Introduction
The physical mechanism and the interplay of confinement and chiral symmetry breaking are long-standing puzzles in Quantum Chromodynamics (QCD) that is a fundamental theory in terms of quarks and gluons. At finite temperature, if the quark mass is infinitely large to quench dynamical quarks, confinement and deconfinement of gluons are classified by center symmetry. Dynamical quarks explicitly break center symmetry [1], and the order parameter for center symmetry, i.e., the Polyakov loop is only an approximate measure of confinement. In the limit of zero quark mass, chiral symmetry is exact and the chiral condensate characterizes the state of matter. Various phases of QCD have been considered as functions of external environmental parameters such as the temperature , the quark chemical potential , the isospin chemical potential [2, 3], the electric and magnetic fields [4, 5, 6], and the rotation angular velocity [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29], and their mixture [30, 31, 32]. In this way, the QCD phase diagram has been intensively studied in theoretical and experimental contexts [33].
Among various parameters, the angular velocity has attracted much attention in recent years. An extraordinary value of is estimated based on the data in the heavy-ion collision [34]. In nonrelativistic theories rotation effects are similar to those induced by magnetic fields, while relativistic rotation of quark matter has features analogous to finite density [30, 8, 9, 35]. In the same way as chiral symmetry restoration at high density, chiral effective models such as the (Polyakov–)Nambu–Jona-Lasino model exhibit the chiral phase transition induced by rotation [7, 8, 9, 10, 11]. Interestingly, rotation directly affects gluons, which makes a contrast to effects of finite density and magnetic fields, and the angular velocity turns out to be a useful probe to confinement and deconfinement. For the purpose to investigate a chance of (de)confinement caused by rotation, the hadron resonance gas model and the holographic QCD model have been adopted [11, 12, 13, 14, 15, 16], which predicted that real rotation favors deconfinement. Thus, it is expected that the deconfinement critical temperature should decrease for larger .
The numerical results from lattice-QCD simulations [18, 19, 20] have invoked controversies reporting a trend opposite to the model predictions. That is, the Polyakov loop decreases as a result of real rotation, so that the deconfinement critical temperature should increase. A technical subtlety lies in the treatment of analytical continuation; the sign problem can be evaded with imaginary angular velocity, , and then the critical temperature is inferred from if the Wick rotation does not hit singularities.
The subtle connection between real and imaginary rotation has been studied recently [10, 20, 23, 24, 28, 29] but the validity condition for analytical continuation is only partially clarified for the moment. There are various scenarios; some support the lattice results [21, 22], some propose a modified model which can explain the lattice results [23], some suggest a non-monotonic function of [24, 25].
Previously in Ref. [36], we exploited the perturbative Polyakov loop potential to investigate rotation effects based on the first-principles approach. Thanks to the asymptotic freedom, QCD or the pure gluonic theory can be perturbatively analyzed if the temperature is sufficiently high. Without rotation, the perturbative Polyakov loop potential, known as the GPY-Weiss potential [37, 38, 39, 40, 41], spontaneously breaks center symmetry, which is consistent with deconfinement expected at high temperature; for a review, see Ref. [42]. We obtained the Polyakov loop potential in the rotating pure gluonic system and found that the Wick rotation, , is hindered by singularity which corresponds to the causality violation. To cure this obstacle, the system must have a finite boundary, and then a closed analytical expression is no longer available.
Interestingly enough, the theory with imaginary angular velocity contains rich physics. In our previous work [36], we discovered exotic realization of confinement even at arbitrarily high temperature, i.e., the perturbatively confined phase. Thus, is a novel theoretical device to tackle the confinement mechanism. In the present work, we will discuss two nontrivial extensions from our previous study. One is exploring inhomogeneous structures away from the imaginary rotation center. In Ref. [36], we focused on only (where is the radial distance) and argued that confinement is homogeneously realized for at . There are some related works [27, 28, 29] in favor of inhomogeneous confinement/deconfinement with real and imaginary rotation. We will show that the perturbative Polyakov loop potential exhibits an inhomogeneous distribution of the Polyakov loop at . Another extension is the inclusion of dynamical quarks, with which we can consider chiral symmetry breaking in the perturbatively confined phase. The relation between confinement and chiral symmetry breaking is not fully understood, but under reasonable assumptions, one can see that confinement leads to chiral symmetry breaking or spontaneous generation of quark mass. We will verify that quark mass is indeed nonzero in the perturbatively confined phase.
2 Inhomogeneous confinement and deconfinement
Real rotation is a real-time effect, while imaginary rotation is interpreted as a geometrical effect. The latter has a theoretical advantage as explicated below. In the presence of an angular velocity vector, , a thermal system is described by the following partition function:
(1) |
where is the inverse temperature and denotes the total angular momentum operator. We can regard the above expression as a thermal expectation value of the topological operator, . For pure imaginary , this topological operator,
(2) |
is a unitary operator that generates a rotation by angle along the rotation axis . In the Euclidean path integral, we can take account of in the partition function (1) by imposing an aperiodic thermal boundary condition which in cylindrical coordinates takes the form,
(3) |
Alternatively, we can perform the coordinate transformation to change the condition (3) to the ordinary periodic one, but then the metric gains nontrivial components in a rotating frame, i.e.,
(4) |
Both approaches lead us to the same partition function.
We shall discuss inhomogeneous structures in the perturbatively confined phase in pure gluonic matter with and . Although we have already given a derivation in Ref. [36], we quickly review the formulation to make this paper self-contained. In our calculation, we take the cylindrical coordinates, , and assume rigidly rotating matter along the -axis.
We take a constant background and then is replaced by the covariant derivative . The inhomogeneity will be studied within the local density approximation (in which the spatial derivatives on are neglected). Using the Cartan subalgebra of Lie algebra of gauge group , the covariant derivative is
(5) |
where is a vector of bases of and is normalized . Using this background field, we fix the gauge by setting .
In the actual calculations, the scalar Laplacian, , is the basic building block. We solve the equation of motion, , to find the spectrum as
(6) |
where quantum numbers are , , , and with denoting the union of the root system and the zero roots. The ghost contribution needs the determinant of , which is given in terms of .
The contribution from the gauge field requires the vector Laplacian, . For the explicit form of , see Ref. [36]. The eigenvalue spectrum is the same as the scalar one. Each eigenmode has four polarization degrees of freedom, and two out of four are canceled by the ghost contribution. After some calculations, we arrive at the Polyakov loop potential resulting from the two physical (transverse) modes as
(7) |
For notational brevity, we introduced a dimensionless imaginary angular velocity; . By expanding the logarithm, we can perform the momentum integration to simplify the above form into
(8) |
The potential is dependent on the dimensionless radial distance; . The potential is minimized at the optimal value of , and the Polyakov loop expectation value, , is evaluated accordingly. The denominator has singularities at even for (i.e., free theory), which is consistent with Ref. [43].
Specifically, for the Yang-Mills theory, the background gauge field is , where and constitute the Cartan subalgebra of . The traced Polyakov loop in the fundamental representation of is
(9) |
In our previous work [36], we showed that confinement, , is realized for at .
It is a straightforward exercise to find the global minima of the potential (8) for . Figure 1 shows the results from such extensive analyses of Eq. (8) for (left) and (right). It is notable that both cases generally develop -dependent spatial structures. For the case as shown in the left of Fig. 1, the Polyakov loop changes to zero, indicating confinement, with second-order phase transition as grows up. In the previous paper [36], we focused on two edges of and only. Our present results imply that, for for example, there should appear a spatial interface separating the confined phase for and the deconfined phase for . We can confirm a qualitatively similar trend for the case, but the detailed shape looks different as shown in the right of Fig. 1. In this case of , the homogeneously confined phase is realized around , and the interface between confinement and deconfinement emerges only when becomes larger.
We should emphasize that this dependence of inhomogeneity is qualitatively consistent with the lattice-QCD results [29]. As closely discussed in our previous work [36], our dependence of the Polyakov loop is opposite to the lattice-QCD case; goes smaller with larger imaginary rotation in our perturbative calculation, while goes smaller with larger real rotation in the lattice-QCD simulation. However, in the both cases of ours and the lattice-QCD calculations, becomes larger (smaller) with larger for a given imaginary (real) angular velocity. This makes a sharp contrast to preceding works [30, 14, 15, 27, 28] and might be a key observation to resolve the controversy of rotating QCD matter.
Now, let us turn our attention to the Polyakov loop potential and the symmetry properties. Since we sum up all the roots in , different pairs can have the same potential and the Polyakov loop value according to Weyl symmetry of the root lattice. We can see the characteristic patterns of the potential minima in Fig. 2. The repetition of the minima reflects Weyl symmetry. The red triangle region with three edges, , , , is the fundamental domain and we can identify the state of matter from the minimum inside this triangular domain. The triangle has geometric symmetry, including three-fold rotational symmetry, which manifests center symmetry, and two-fold reflective symmetry, which manifests charge conjugation. The potential shape visualized by the pattern in Fig. 2 does not change by rotation but the phase of the Polyakov loop, , does. The center of the triangle at corresponds to , that is, a center symmetric vacuum. Although the potential minima may break center symmetry, the charge conjugation symmetry is never broken.
We point out a nontrivial observation in Fig. 2; an emergent symmetry is realized at . We observe a reflective mirror on the line of . This emergent symmetry comes from the vanishing of odd- terms in the one-loop potential (8) at and exists for not only but any radius. It could be either a one-loop artifact or a genuine symmetry. In the latter case, it has to be a non-invertible symmetry like that in the 2D critical Ising model since it exchanges the unbroken and broken vacua.
In this work, we also quantify the spatial inhomogeneity. In Fig. 3, we plot the Polyakov loop potential for (confined phase) and (deconfined phase) at in the case. It is intriguing that the potential minima are located in a different way from the ordinary perturbative vacuum where three vertices of the triangle minimize the potential. In this sense, this deconfined phase discovered in the upper right region in the right panel of Fig. 1 may have exotic properties different from the ordinary one.
We mention the difference from the previous arguments [27] based on the Tolman-Ehrenfest (TE) effect. In our calculation, we treat as a Lagrange multiplier to conserve the total energy. By construction, is a global variable without dependence. Therefore, in our study, we do not need to introduce the apparent temperature as a result of the TE effect. Nevertheless, the calculations naturally lead to dependent structures.
3 Chiral symmetry breaking in the perturbatively confined phase
We can repeat similar calculations including dynamical quark contributions that break center symmetry explicitly. We can also address a relation between confinement and chiral symmetry breaking from two (approximate) order parameters, namely, the Polyakov loop and the dynamical quark mass, , which is rooted in the chiral condensate. In the perturbative treatment, the pressure (the free energy) is maximized (minimized) for , and the dynamical mass generation is energetically disfavored. It is quite interesting what would happen in the perturbatively confined phase with imaginary rotation.
The theoretical treatments of fermions in the rotating frame are found in Refs. [30, 44, 7]. We should be careful about the fact that the Polyakov loop coupling with quarks is given by the fundamental representation. We can calculate the fermionic partition function by imposing an aperiodic thermal boundary condition or equivalently considering the ordinary anti-periodic boundary condition in the rotating frame. In this paper, for fermions, we choose the latter and the fermionic partition function is the determinant of the Dirac operator in the rotating frame, i.e.,
(10) |
Here, is the covariant derivative including the background field with , where and . We denote the gamma matrices of the flat space-time by and . After all, we arrive at the expression for the Polyakov loop potential per one fermion (particle or anti-particle but without flavor degrees of freedom) as
(11) |
where denotes the set of the fundamental weights of . The limit of renders the above expression to the free quark and anti-quark gas energy with imaginary rotation derived from statistical mechanics.
First, we shall consider the interplay of and at by setting . Then, the total effective potential is given by (where comes from the sum over particle and anti-particle). In our numerical calculation, we consider and . We note that the former region for is sufficient thanks to Weyl symmetry. Because we do not model quark interactions to keep our analysis as model independent as possible, blows up in the phase where chiral symmetry is fully broken. We thus set a numerical cutoff on in the practical calculation.
Figure 4 shows the potential contour as functions of two order parameters in the case. The minimum of the potential is given by the darkest spot in the figures. In the left of Fig. 4 for , the minimum is still located at and (i.e., ). This is different from the pure gluonic case in the left of Fig. 1; fermions generally favor deconfinement. Interestingly in the right of Fig. 4 for , however, the minimum moves to and blows up to the cutoff value, which signifies confinement and spontaneous chiral symmetry breaking.
We make the plots for the behavior of the Polyakov loop and the dynamical mass in Fig. 5. The blue solid line represents , while the green dashed line represents . With increasing , the Polyakov loop goes smaller and the dynamical mass grows up. Then, a confined and chiral symmetry broken state is favored near . For the case as shown in the left of Fig. 5, the behavior is periodic because the gauge group covers the fermionic parity . Also we see that the behavior is symmetric around . For the case in the right of Fig. 5, in contrast, the order parameters depend on in a complicated way and the confined phase seems to be less favored. Also the periodicity becomes because the gauge group does not contain .
In the same way as the pure gluonic theory, we have quantified the inhomogeneous structures in Fig. 6. The left and the right panels show the Polyakov loop and the dynamical mass, respectively, as functions of and . We have performed this calculation for the case only. As compared to the left of Fig. 1 in the pure gluonic case, the perturbatively confined region is shrunk to the upper left corner. Interestingly, in view of the right of Fig. 6, confinement and chiral symmetry breaking are locked together and the phase transitions occur simultaneously even at large imaginary rotation. In addition, there appear some substructures, but the numerical convergence seems to be insufficient. We have tried, but the potential shape is quite complicated and it is difficult to obtain smooth results in the whole region at . This complication might be a precursory of possible fractals [45]. More investigations about these substructures as well as the full calculations are left as future problems.
4 Conclusions
We considered the inhomogeneous structures and the chiral symmetry breaking in the perturbatively confined phase which is realized by rotation with imaginary angular velocity, . In both color and cases, the perturbatively confined phase is induced at the rotation center, , for above a certain threshold. We found that the perturbative confinement may break down in the region away from the rotation center. This indicates the existence of phase interface between confinement and deconfinement and the critical temperature is also dependent accordingly. Our results predict that the critical temperature should decrease with increasing in the presence of imaginary rotation. This trend agrees with the latest results from the lattice simulation. Although the analytical continuation to real rotation has some subtle points, the inhomogeneity we have confirmed from the perturbative effective potential can presumably persist in real rotating systems.
We then added the quark contribution to the Polyakov loop effective potential, which is also a function of the dynamical quark mass, . Our analysis of the potential minimum search showed that confinement and chiral symmetry breaking are tightly correlated in such unconventional systems with imaginary rotation and even without fermionic interaction. In the presence of quarks that explicitly break center symmetry, the Polyakov loop is only an approximate order parameter. Still, in the region where the Polyakov loop is vanishingly small, departs from zero signifying spontaneous breaking of chiral symmetry. In our treatment, in fact, blows up when chiral symmetry is broken, so that large suppresses center symmetry breaking. In this way, we established simultaneous realization of confinement and chiral symmetry breaking. We also found highly complicated substructures in the deconfined regions. The imaginary-rotating inhomogeneous states with dynamical quarks should deserve further investigations in the future.
Acknowledgement
The authors thank Maxim Chernodub and Xu-Guang Huang for useful discussions. This work was supported by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Nos. 22H01216 (K.F.) and 22H05118 (K.F.) and JST SPRING, Grant Number JPMJSP2108 (Y.S.).
References
- Svetitsky [1986] B. Svetitsky, Phys. Rept. 132, 1 (1986).
- Son and Stephanov [2001a] D. T. Son and M. A. Stephanov, Phys. Atom. Nucl. 64, 834 (2001a), arXiv:hep-ph/0011365 .
- Son and Stephanov [2001b] D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001b), arXiv:hep-ph/0005225 .
- Yamamoto [2013] A. Yamamoto, Phys. Rev. Lett. 110, 112001 (2013), arXiv:1210.8250 [hep-lat] .
- Andersen et al. [2016] J. O. Andersen, W. R. Naylor, and A. Tranberg, Rev. Mod. Phys. 88, 025001 (2016), arXiv:1411.7176 [hep-ph] .
- Chen et al. [2017] H.-L. Chen, K. Fukushima, X.-G. Huang, and K. Mameda, Phys. Rev. D 96, 054032 (2017), arXiv:1707.09130 [hep-ph] .
- Jiang and Liao [2016] Y. Jiang and J. Liao, Phys. Rev. Lett. 117, 192302 (2016), arXiv:1606.03808 [hep-ph] .
- Wang et al. [2019] X. Wang, M. Wei, Z. Li, and M. Huang, Phys. Rev. D 99, 016018 (2019), arXiv:1808.01931 [hep-ph] .
- Chernodub and Gongyo [2017] M. N. Chernodub and S. Gongyo, JHEP 01, 136, arXiv:1611.02598 [hep-th] .
- Chen et al. [2023] H.-L. Chen, Z.-B. Zhu, and X.-G. Huang, Phys. Rev. D 108, 054006 (2023), arXiv:2306.08362 [hep-ph] .
- Sun et al. [2023] F. Sun, K. Xu, and M. Huang, Phys. Rev. D 108, 096007 (2023), arXiv:2307.14402 [hep-ph] .
- Fujimoto et al. [2021] Y. Fujimoto, K. Fukushima, and Y. Hidaka, Phys. Lett. B 816, 136184 (2021), arXiv:2101.09173 [hep-ph] .
- Chen et al. [2021] X. Chen, L. Zhang, D. Li, D. Hou, and M. Huang, JHEP 07, 132, arXiv:2010.14478 [hep-ph] .
- Braga et al. [2022] N. R. F. Braga, L. F. Faulhaber, and O. C. Junqueira, Phys. Rev. D 105, 106003 (2022), arXiv:2201.05581 [hep-th] .
- Yadav [2023] G. Yadav, Phys. Lett. B 841, 137925 (2023), arXiv:2203.11959 [hep-th] .
- Wang and Feng [2024] J.-H. Wang and S.-Q. Feng, (2024), arXiv:2403.01814 [hep-ph] .
- Yamamoto and Hirono [2013] A. Yamamoto and Y. Hirono, Phys. Rev. Lett. 111, 081601 (2013), arXiv:1303.6292 [hep-lat] .
- Braguta et al. [2020] V. V. Braguta, A. Y. Kotov, D. D. Kuznedelev, and A. A. Roenko, Pisma Zh. Eksp. Teor. Fiz. 112, 9 (2020).
- Braguta et al. [2021] V. V. Braguta, A. Y. Kotov, D. D. Kuznedelev, and A. A. Roenko, Phys. Rev. D 103, 094515 (2021), arXiv:2102.05084 [hep-lat] .
- Yang and Huang [2023] J.-C. Yang and X.-G. Huang, (2023), arXiv:2307.05755 [hep-lat] .
- Mameda and Takizawa [2023] K. Mameda and K. Takizawa, Phys. Lett. B 847, 138317 (2023), arXiv:2308.07310 [hep-ph] .
- Sun et al. [2024] F. Sun, J. Shao, R. Wen, K. Xu, and M. Huang, (2024), arXiv:2402.16595 [hep-ph] .
- Cao [2024] G. Cao, Phys. Rev. D 109, 014001 (2024), arXiv:2310.03310 [nucl-th] .
- Jiang [2023] Y. Jiang, (2023), arXiv:2312.06166 [hep-th] .
- Gaspar et al. [2023] I. I. Gaspar, L. A. Hernández, and R. Zamora, Phys. Rev. D 108, 094020 (2023), arXiv:2305.00101 [hep-ph] .
- Tabatabaee Mehr [2023] S. M. A. Tabatabaee Mehr, Phys. Rev. D 108, 094042 (2023), arXiv:2306.11753 [nucl-th] .
- Chernodub [2021] M. N. Chernodub, Phys. Rev. D 103, 054027 (2021), arXiv:2012.04924 [hep-ph] .
- Chernodub et al. [2023] M. N. Chernodub, V. A. Goy, and A. V. Molochkov, Phys. Rev. D 107, 114502 (2023), arXiv:2209.15534 [hep-lat] .
- Braguta et al. [2023] V. V. Braguta, M. N. Chernodub, and A. A. Roenko, (2023), arXiv:2312.13994 [hep-lat] .
- Chen et al. [2016] H.-L. Chen, K. Fukushima, X.-G. Huang, and K. Mameda, Phys. Rev. D 93, 104052 (2016), arXiv:1512.08974 [hep-ph] .
- Flachi and Fukushima [2018] A. Flachi and K. Fukushima, Phys. Rev. D 98, 096011 (2018), arXiv:1702.04753 [hep-th] .
- Zhang et al. [2020] H. Zhang, D. Hou, and J. Liao, Chin. Phys. C 44, 111001 (2020), arXiv:1812.11787 [hep-ph] .
- Fukushima and Hatsuda [2011] K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74, 014001 (2011), arXiv:1005.4814 [hep-ph] .
- Adamczyk et al. [2017] L. Adamczyk et al. (STAR), Nature 548, 62 (2017), arXiv:1701.06657 [nucl-ex] .
- Fukushima et al. [2020] K. Fukushima, T. Shimazaki, and L. Wang, Phys. Rev. D 102, 014045 (2020), arXiv:2004.05852 [hep-ph] .
- Chen et al. [2022] S. Chen, K. Fukushima, and Y. Shimada, Phys. Rev. Lett. 129, 242002 (2022), arXiv:2207.12665 [hep-ph] .
- Gross et al. [1981] D. J. Gross, R. D. Pisarski, and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981).
- Weiss [1981] N. Weiss, Phys. Rev. D 24, 475 (1981).
- Weiss [1982] N. Weiss, Phys. Rev. D 25, 2667 (1982).
- Korthals Altes [1994] C. P. Korthals Altes, Nucl. Phys. B 420, 637 (1994), arXiv:hep-th/9310195 .
- Gocksch and Pisarski [1993] A. Gocksch and R. D. Pisarski, Nucl. Phys. B 402, 657 (1993), arXiv:hep-ph/9302233 .
- Fukushima and Skokov [2017] K. Fukushima and V. Skokov, Prog. Part. Nucl. Phys. 96, 154 (2017), arXiv:1705.00718 [hep-ph] .
- Chernodub [2022] M. N. Chernodub, (2022), arXiv:2210.05651 [quant-ph] .
- Ebihara et al. [2017] S. Ebihara, K. Fukushima, and K. Mameda, Phys. Lett. B 764, 94 (2017), arXiv:1608.00336 [hep-ph] .
- Ambruş and Chernodub [2023] V. E. Ambruş and M. N. Chernodub, Phys. Rev. D 108, 085016 (2023), arXiv:2304.05998 [hep-th] .