Arxiv:2012.02926V2 [Hep-Ph] 26 Feb 2021 of This Critical Point Was Established at M 0.33G [60– 62]; the Resulting Phase Diagram Is Shown In≈ Fig.1

Real-time dynamics of Chern-Simons fluctuations near a critical point Kazuki Ikeda,1, ∗ Dmitri E. Kharzeev,2, 3, 4, y and Yuta Kikuchi3, z 1Research Institute for Advanced Materials and Devices, Kyocera Corporation, Soraku, Kyoto 619-0237, Japan 2Center for Nuclear Theory, Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA 3Department of Physics, Brookhaven National Laboratory, Upton, New York 11973-5000 4RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000 The real-time topological susceptibility is studied in (1+1)-dimensional massive Schwinger model with a θ-term. We evaluate the real-time correlation function of electric field that represents the topological Chern-Pontryagin number density in (1+1) dimensions. Near the parity-breaking critical point located at θ = π and fermion mass m to coupling g ratio of m=g ≈ 0:33, we observe a sharp maximum in the topological susceptibility. We interpret this maximum in terms of the growth of critical fluctuations near the critical point, and draw analogies between the massive Schwinger model, QCD near the critical point, and ferroelectrics near the Curie point. I. INTRODUCTION (1 + 1)-dimensional Minkowski space is Z 2 1 µν gθ µν ¯ Schwinger model [1] is quantum electrodynamics in S = d x F Fµν + Fµν + (iD= m) ; −4 4π − (1 + 1) space-time dimensions. For massless fermions, (1) Schwinger model is analytically solvable and equivalent µ to the theory of a free massive boson field [1{8]; however with D= = γ (@µ igAµ). Note that the gauge field Aµ the model with massive fermions presents a challenge for and the coupling− constant g have mass dimensions 0 and analytical methods and has a rich dynamics. 1, respectively. Upon a chiral transformation, eiγ5θ Recently, quantum algorithms have emerged as an ef- and ¯ ¯eiγ5θ, the action is transformed to,1 ! ficient (and potentially superior) way to explore the dy- ! Z 1 2 µν ¯ µ iγ5θ namics of quantum field theories, including the Schwinger S = d x F Fµν + (iγ Dµ me ) : (2) model [9{56]. Previously, we have addressed the −4 − real-time dynamics of vector current [47] (a (1 + 1)- It is clear from (2) that the massive theory with a positive dimensional analog of the chiral magnetic current [57, mass m > 0 at θ = π is equivalent to the theory at θ = 0 58]) induced by the \chiral quench" { an abrupt change but with a negative mass m. of the θ-angle, or the chiral chemical potential. In this pa- − per, we will explore the connection between the real-time topological fluctuations and criticality using Schwinger model as a testing ground. The massive Schwinger model possesses a quantum phase transition at θ = π between the phases with opposite orientations of the electric field, see Fig.1. For m g, this phase transition is first order. However, it was shown by Coleman [59] that the line of the first order phase transition terminates at some critical value m∗, where the phase transition is second order. The position ∗ arXiv:2012.02926v2 [hep-ph] 26 Feb 2021 of this critical point was established at m 0:33g [60{ 62]; the resulting phase diagram is shown in≈ Fig.1. The phase diagram of the theory in the (m=g; θ) plane is thus FIG. 1: Phase diagram of the massive Schwinger model in reminiscent of the phase diagram of QCD in the (T; µ) the (θ; m=g) plane. At θ = π and large masses m > m∗, the plane of temperature T and baryon chemical potential µ ferroelectric phases with opposite orientations of electric field [63, 64]. are separated by the line of the first order phase transition. ∗ To understand the physics behind this phase diagram, This line terminates at m ≈ 0:33g at the critical point, where the phase transition is second order. For small masses m we need to recall the role of θ-angle in the model. The ∗ action of the massive Schwinger model with θ term in m , the electric field is screened by the production of light fermion-antifermion pairs. ∗kazuki.ikeda.gtATkyocera.jp ydmitri.kharzeevATstonybrook.edu 1 The action is invariant under this transformation only up to the zykikuchiATbnl.gov boundary term. 2 II. TOPOLOGICAL FLUCTUATIONS NEAR phase. At some critical value m m∗, the effective THE CRITICAL POINT potential becomes flat { this corresponds≈ to the critical point with a second order phase transition. Let us now discuss topological fluctuations in massive Schwinger model. For that purpose, it will be conve- For m > m∗, the minima are separated by a potential nient to use the Hamiltonian formalism with the tempo- barrier, and the transition between them (corresponding ral gauge, A0 = 0. From the action (2) the canonical to the change in the direction of the electric dipole mo- momentum conjugate to A1 can be read off as Π = A_ 1. ment of the system) is first order { for example, having a The corresponding Hamiltonian is then given by domain with an opposite orientation of the electric dipole moment would cost an additional energy due to the \sur- Z 1 2 ¯ 1 iγ5θ face tension". Due to this extra energy, the fluctuations H = dx Π (iγ D1 me ) ; (3) 2 − − of the electric dipole moment at m m∗ are suppressed. At the critical point, where m m∗, the potential bar- ≈ with commutation relations [A1(x); Π(y)] = iδ(x y), rier between the two minima disappears, and the \surface ¯ − and (x); (y) = γ0δ(x y). tension" of the domains with opposite orientations of the f g − electric dipole moments vanishes. Because of this, the In the bosonized description of the theory, the Hamil- fluctuations of the electric dipole moment near the crit- tonian is given by (see AppendixA for details) ical point are strongly enhanced. Here one can draw a useful analogy to the physics of ferroelectrics, where the Z 2 2 h1 2 1 2 µ θ electric susceptibility exhibits critical behavior near the H = dx '_ + (@1') + ' + 2 2 2 2pπ (4) Curie point [65]. i cmµ cos(2pπ') : When the ratio m=g becomes very small, the elec- − tric field is easily screened by the production of fermion- The potential of the model is thus given by antifermion pairs (or kink-antikink pairs in the bosonized description), and the fluctuations of electric dipole mo- µ2 θ 2 ment again become small. Basing on this qualitative U(') = ' + cmµ cos (2pπ'); (5) 2 2pπ − picture (that we will confirm below with a more formal treatment), we expect to see a maximum in the electric where µ = g=pπ is the mass of the scalar boson, and the susceptibility near the critical point. We will show that dimensionless coefficient c is given by c = eγ =2π with the this is indeed the case. Euler constant γ = 0:5774. The electric field is related to the boson field as In order to characterize the topological fluctuation we θ compute the static topological susceptibility. It is the E = µ ' + ; (6) zero frequency and wavelength limit of the real-time two- − 2pπ point correlation function of the topological charge: where the first term is the quantum, dynamical contri- bution and the second one is the classical background χCS induced by the θ angle, Ecl = gθ=2π. − Z The physics is periodic as a function of θ in the absence = lim lim d2x ei!t−ikx Q(x)Q(0) Q(x) Q(0) of boundaries { as θ increases, the electric field becomes !!0 k!0 h i − h ih i capable of producing a fermion-antifermion pair (or kink- Z = 2Re d2xΘ(t) Q(x)Q(0) Q(x) Q(0) ; (8) antikink pair, in bosonic description) that screens it. At h i − h ih i θ = π the potential takes the form 2 µ2 pπ where Θ(t) is the Heaviside's step function. Q(x) = U(') = ' + cmµ cos (2pπ'): (7) µ 2 2 − @µK is the density of the Chern-Pontryagin number given by the divergence of the Chern-Simons current µ µ µν When m µ, the potential has two well separated K . In (1 + 1) dimensions, K = (g=2π) Aν , and µν minima at ' 0; pπ, associated with spontaneous Q = (g=4π) Fµν , where Fµν is the electromagnetic field symmetry breaking.≈ − As follows from (6), the minima strength tensor. The corresponding change of Chern- R 2 0 R 2 ' = 0; pπ correspond to the electric fields E = g=2 Simons number is ∆NCS = d xK = (g=2π) d xA1; and E =− g=2, respectively. − the Chern-Pontryagin number density is given by the When m µ, there is a single minimum at ' = electric field: Q(x) = (g=2π)E(x). We call χCS (8) pπ=2, which according to (6) corresponds to the phase the real-time topological susceptibility because it is com- with− no electric field. This is because at small m, the elec- puted from the real-time two-point function unlike the tric field is easily screened by the production of fermion- conventional topological susceptibility which is computed antifermion pairs - so we are dealing with the screened in the Euclidean space-time. 3 III. LATTICE SCHWINGER MODEL at the Hamiltonian, 2 2 N−1 " n i # ag X X y 1 ( 1) A. Lattice Hamiltonian H = χ χi − − 2 i − 2 n=1 i=1 N−1 For the purpose of numerical simulation, we place the i X χy χ χy χ (14) theory (3) on a spatial lattice. We introduce the stag- 2a n+1 n n n+1 y − n=1 − gered fermion χn and χn [66, 67], and lattice gauge field N operators Ln and Un with an integer n labeling a lat- X n y tice site; the lattice spacing is a.

Load more