## Isospin splitting of pion elliptic flow in relativistic heavy-ion collisions

He Liu , ∗ 1, 2, 3 Feng-Tao Wang , † 1, 2 Kai-Jia Sun , ‡ 3 Jun Xu , § 4, 1 and Che Ming Ko ¶ 3

1 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China 2 University of Chinese Academy of Sciences, Beijing 100049, China 3 Cyclotron Institute and Department of Physics and Astronomy, Texas A &amp; M University, College Station, TX77843-3366, USA 4

Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China (Dated: October 4, 2019)

Based on the framework of an extended multiphase transport model with mean-field potentials in both the partonic phase and the hadronic phase, we explain the elliptic flow difference between π + and π -in the beam-energy scan program at the relativistic heavy-ion collider by incorporating the vector-isovector potential for quarks and antiquarks with different isospins. It is found that the isospin splitting of pion elliptic flow favors a strong vector-isovector interaction, and thus serves as a probe of the quark matter equation of state as well as the QCD phase structure at finite baryon and isospin chemical potentials.

Understanding the properties of the strongly interacting quark-gluon plasma (QGP) as well as the hadronquark phase transition is one of the main goals of relativistic heavy-ion collision experiments. The elliptic flow ( v 2 ) characterizing the different numbers of particles moving in the reaction plane and out of the reaction plane is a good probe of the QGP properties [1, 2]. The number-of-constituent-quark (NCQ) scaling of v 2 [3-8] in ultrarelativistic heavy-ion collisions serves as an evidence of the existence of the QGP, indicating that the elliptic flow is mostly developed in the partonic phase. Recently, 'low-energy' relativistic heavy-ion collisions were carried out in the beam-energy scan (BES) program at the relativistic heavy-ion collider (RHIC), in order to search for the signal of the critical point of the hadron-quark phase transition. Among various interesting phenomena different from those observed in ultrarelativistic heavy-ion collisions is the splitting of v 2 between protons and antiprotons, K + and K -, and π + and π -[9], showing the breakdown of the NCQ scaling law. The observed v 2 splitting can be due to the larger v 2 for transported quarks than that for produced quarks, or similarly, their different rapidity dependencies [10-12], hydrodynamic evolution of the QGP at finite baryon chemical potential [13, 14], and the smaller radial flow of particles than their antiparticles [15]. Although the v 2 splitting between π -and π + as well as its linear dependence on the charge asymmetry was proposed to be attributed to the electric quadrupole moment in the produced QGP induced by the chiral magnetic wave [16], it is disfavored by transport simulations based on the chiral kinetic equations of motion [17, 18]. In our previous study [19, 20], we have shown that the v 2 splitting between protons and antiprotons as well as between K + and K -can be attributed to the different mean-field potentials for particles and their antiparticles. We will demonstrated in the present study that the v 2 splitting between π + and π -can be explained by the different mean-field potentials for particles with different isospins.

∗ liuhe@sinap.ac.cn

† wangfengtao@sinap.ac.cn

‡ kjsun@tamu.edu

§ corresponding author: xujun@zjlab.org.cn

¶ ko@comp.tamu.edu

The elliptic flow in heavy-ion collisions at intermediate energies was reguarded as a probe of the mean-field potential or the nuclear matter equation of state [21]. In noncentral heavy-ion collisions at relativistic energies, the produced QGP is of an almond shape in the transverse plane. Particles with attractive potentials are more likely to be trapped in the system and move in the direction perpendicular to the reaction plane, while those with repulsive potentials are more likely to leave the system and move along the reaction plane, thus reducing and enhancing their respective elliptic flows. In the baryonrich matter produced in low-energy relativistic heavyion collisions, light quarks and baryons are affected by a more repulsive potential compared to light antiquarks and antibaryons as a result of the vector coupling, leading to the v 2 splitting between protons and antiprotons as well as between K + and K -[19, 20]. Since heavy-ions are neutron-rich nuclei, the produced matter is not only baryon-rich but also neutron-rich or d -quark-rich, where d (¯) quarks are expected to have a more repulsive potenu tial than u ( ¯ ) d in the presence of a vector-isovector interaction. This will enhance the v 2 of π -while reduce the v 2 of π + , qualitatively consistent with that observed experimentally. Since the stopping power becomes larger with the decreasing collision energy, the medium at midrapidities is expected to be not only more baryon-rich but also more neutron-rich, again qualitatively consistent with the larger v 2 splitting between π + and π -at lower collision energies. The isospin splitting of charged pion v 2 can thus be used to constrain the strength of the vector-isovector interaction and to help determine the quark matter equation of state as well as the QCD phase structure at finite baryon and isospin chemical potentials [22].

The present study is based on the framework of an extended multiphase transport model (extended

AMPT) [19, 20, 23]. In this model the momentum distribution of initial partons is from melting hadrons generated by the heavy-ion jet interaction generator (HIJING) model [24] through the Lund string fragmentation model with its parameters fitted to particle multiplicities at RHIC-BES energies [25]. The coordinates of these partons in the transverse plane are set to be the same as those of the colliding nucleons from which they are produced, while their coordinates in the longitudinal direction are uniformly distributed within the finite thickness of the participant matter by considering the Lorentz contraction. The evolution of the partonic phase is then described by a relativistic Boltzmann transport approach, with the mean-field potentials of partons from a 3-flavor Nambu-Jona-Lasinio (NJL) model, and isotropic two-body scatterings with the cross section fitted by the final charged-particle elliptic flow. The meanfield potentials depend on the local phase-space distribution functions of partons, which are calculated numerically from averaging over parallel events using the test-particle method [26, 27]. The partonic phase ends when the central energy density is below 0.8 GeV/fm , 3 and then a coalescence model is used to describe the hadronization procedure. Different from the naive spatial coalescence algorithm for hadronization, we now use an improved coalescence algorithm [28] by selecting the combination of a quark-antiquark pair, three quarks, or three antiquarks, which has the largest Wigner function rather than the closest distance in coordinate space. In the improved coalescence scenario, the baryon, charge, and strangeness numbers for a single event are conserved, while mixed-event combinations are allowed for choosing partons close in phase space from different events in the spirit of the test-particle method, so that the parton combinations can have a large Wigner function. The improved hadronization algorithm prefers combination in not only coordinate space but also momentum space, thus preserves better the flow of partons during the hadronization procedure. The species of hadrons are determined by the flavor of its constituent quarks as well as their invariant mass [29]. The evolution of these hadrons is described by a relativistic transport (ART) model [30], with the mean-field potentials for pions, kaons, and baryons as well as antibaryons incorporated [31], in addition to various elastic collisions, inelastic collisions, and decay channels satisfying the detailed balance condition. The isovector hadronic mean-field potentials relevant for the present study are the pion s -wave potential and the symmetry potential for baryons with different isospins. The pion s -wave potential, which is repulsive for π -and attractive for π + in neutron-rich medium, is taken from Ref. [32]. The pion p -wave potential has shown to be less important compared to the s -wave potential [33-35]. The momentum-independent symmetry potential incorporated in the hadronic phase leads to a more repulsive potential for neutron-like baryons, e.g., ∆ -, and a more attractive potential for proton-like baryons, e.g., ∆ ++ , in the neutron-rich medium. The effects of these hadronic isovector potentials are, however, expected to be largely suppressed due to the low baryon density of the hadronic phase, originating from a more realistic finite thickness of initial partonic medium and the hadronization criterion employed in the present study compared with those used in Ref. [31].

The mean-field potentials of partons are given by the 3-flavor NJL model with isovector interactions, and its Lagrangian can be written as [47]

$$\text{ac} n _ { \text{$s$-navor invalid mooel with isovector interactions, and its} } \\ \text{long} _ { \text{$L$_NIL} } \, = \, \bar { q } ( i \hat { \emptyset} - \hat { m } ) q + \frac { G _ { S } } { 2 } \sum _ { a = 0 } ^ { 8 } [ ( \bar { q } \lambda _ { a } q ) ^ { 2 } + ( \bar { q } i \gamma _ { 5 } \lambda _ { a } q ) ^ { 2 } ] \\ \text{ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$d$-n} } \\ \text{$f$-n} \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { \text{$f$-n} } \\ \text{$ans} _ { i } } \, = \, ( u, d, s ) ^ { T } \text{ and } \hat { m } = \text{diag} ( m _ { u }, m _ { d }, m _ { s } ) \\ \text{use} _ { \text{$e the quark fields and the current quark mass matrix} }$$

In the above, q = ( u, d, s ) T and ˆ m = diag( m ,m ,m u d s ) are the quark fields and the current quark mass matrix for u d , , and s quarks, respectively; λ a are the Gell-Mann matrices with λ 0 = √ 2 / I 3 in the 3-flavor space with the SU(3) symmetry; G S and G V are respectively the scalarisoscalar and the vector-isoscalar coupling constant; and the K term represents the six-point Kobayashi-Maskawat'Hooft interaction that breaks the axial U (1) A symmetry [36]. The additional G IS and G IV terms represent the scalar-isovector and the vector-isovector interactions, with G IS and G IV the corresponding coupling constants, respectively. Since the Gell-Mann matrices with a = 1 2 3 are identical to the Pauli matrices in , , u and d space, the isovector couplings break the SU(3) symmetry while keeping the isospin symmetry. In the present study, we employ the parameters m u = m d = 3 6 MeV, . m s = 87 MeV, G S Λ 2 = 3 6, . K Λ 5 = 8 9, and the cutoff value in the . momentum integral Λ = 750 MeV given in Refs. [37, 38]. As is known, the position of the critical point for the chiral phase transition is sensitive to G V [38-40], which was later constrained within 0 5 . G S &lt; G V &lt; 1 1 . G S from the relative v 2 splitting between protons and antiprotons as well as between K + and K -in relativistic heavy-ion collisions [19]. In the present study, we choose G V = 1 1 . G S throughout the calculation.

In the mean-field approximation and considering only the flavor-singlet state for the vector-isoscalar term, the single-particle Hamiltonian for partons with flavor i ( i = u, d, s ) can be written as

$$\dot { \vec { H } } _ { i } \, = \, & \, \sqrt { M _ { i } ^ { 2 } + p _ { i } ^ { * 2 } } \pm \frac { 2 } { 3 } G _ { V } ( \rho _ { u } ^ { 0 } + \rho _ { d } ^ { 0 } + \rho _ { s } ^ { 0 } ) \\ & \pm G _ { I V } \tau _ { 3 i } ( \rho _ { u } ^ { 0 } - \rho _ { d } ^ { 0 } ). \quad \, ( 2 ) \\ \text{Here we take the convention that the inner (lower)} \sin$$

Here we take the convention that the upper (lower) sign is for quarks (antiquarks). M i is the constituent quark

mass determined by the gap equation

$$M _ { i } = m _ { i } - 2 G _ { s } \sigma _ { i } + 2 K \sigma _ { j } \sigma _ { k } - 2 G _ { 1 S } \tau _ { 3 i } ( \sigma _ { u } - \sigma _ { d } ), \ ( 3 )$$

where σ i = 〈 q q i ¯ i 〉 is the quark condensate, ( , i j , k ) is any permutation of ( u , d , s ), and τ 3 i is the isospin quantum number of quarks, i.e., τ 3 u = 1, τ 3 d = -1, and τ 3 s = 0. /vector p ∗ i = /vector p i ∓ 2 3 G V ( /vector ρ u + /vector ρ d + /vector ρ s ) ∓ G IV τ 3 i ( /vector ρ u -/vector ρ d ) is the real momentum of partons, with /vector p i being the canonical momentum and /vector ρ i ≡ 〈 ¯ q /vector q i γ i 〉 being the space components of the net quark density. ρ 0 i ≡ 〈 ¯ q γ i 0 q i 〉 in Eq. (2) is the time component of the net quark density. The quark condensate and the 4-component net quark density can be calculated from the phase-space distributions of quarks ( f i ) and antiquarks ( ¯ ) through the relations f i

$$\langle \bar { q } _ { i } q _ { i } \rangle \ & = \ - 2 N _ { c } \int _ { 0 } ^ { \Lambda } \frac { d ^ { 3 } p } { ( 2 \pi ) ^ { 3 } } \frac { M _ { i } } { \sqrt { M _ { i } ^ { 2 } + p _ { i } ^ { 2 } } } ( 1 - f _ { i } - \bar { f } _ { i } ) ( 4 ) \quad \text{FIG.} \\ \int _ { \bar { \sigma } \sim \mu \sigma \cdot \vee } \ & = \ \varepsilon \, N \, \int ^ { \Lambda } \frac { d ^ { 3 } p } { \left ( \frac { p ^ { \mu } } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gam { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gamma } { \left ( \frac { \gam { \left ( \frac { \gamma } { \left ( \frac { \gam { \left ( \left ( \frac { \gamma } { \left ( \frac { \gam { \left ( \left ( \frac { \gam { \left ( \left ( \left \bar { \left ( \left \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \bar { \left } } { \left ( \sigma } \sum \delta } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }, } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } | } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ( } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }. } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } & } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \\ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where 2 N c is the quark spin and color degeneracy. In the transport approach, the whole system is divided into grids with f i and ¯ at each cell calculated from averaging f i over parallel events using the test-particle method. In this way, the spatial distributions of the above quantities can be obtained.

$$\langle \bar { q } _ { i } \gamma ^ { \mu } q _ { i } \rangle \ = \ 2 N _ { c } \int _ { 0 } ^ { \Lambda } \frac { d ^ { 3 } p } { ( 2 \pi ) ^ { 3 } } \frac { p ^ { \mu } } { \sqrt { M _ { i } ^ { 2 } + p _ { i } ^ { 2 } } } ( f _ { i } - \bar { f } _ { i } ), \quad ( 5 ) \quad \text{$\text{$\text{$\text{$c}$} \\ \text{$\text{$\text{$model}$} } } \\ \text{where $2N_{\text{$is the mark sign and color degeneracy.\ In}$}$$

Given the single-particle Hamiltonian as Eq. (2), partons with flavor i evolve according to the following canonical equations

$$\frac { d \vec { r } _ { i } } { d t } \, = \, \frac { \partial H _ { i } } { \partial \vec { p } _ { i } } = \frac { \vec { p } _ { i } ^ { * } } { \sqrt { M _ { i } ^ { 2 } + p _ { i } ^ { * } 2 } }, \quad \quad \, ( 6 ) \, \begin{array} { c } \text{eq.} \\ \text{as a} \\ \text{as d} \\ \text{part} \end{array} \,.$$

$$\frac { d \vec { p } _ { i } ^ { * } } { d t } \, = \, - \frac { \partial H _ { i } } { \partial \vec { r } _ { i } } \mp \frac { d } { d t } [ \frac { 2 } { 3 } G _ { V } ( \vec { \rho } _ { u } + \vec { \rho } _ { d } + \vec { \rho } _ { s } ) \, & \quad \text{part} \\ + & G _ { I V } \tau _ { 3 i } ( \vec { \rho } _ { u } - \vec { \rho } _ { d } ) ]. \quad ( 7 ) \, & \quad \text{com} ] \\ \tau _ { u } \, = \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \delta \, \text{AMI} }$$

The time component of the vector-isovector potential G IV τ 3 i ( ρ 0 u -ρ 0 d ) is expected to give a more repulsive (attractive) potential for d ( u ) quarks in the baryon-rich and d -quark-rich partonic phase. The space component of the vector-isovector potential -G IV τ 3 i ( /vector ρ u -/vector ρ d ) may have an opposite effect, but its effect is small at the early stage when the net quark flux has not been developed. A similar argument for the vector-isoscalar potential can be found in Refs. [41, 42]. For the ease of discussion, we define R IV = G IV /G S as the reduced strength of the vector-isovector coupling and will discuss its effects on the isospin splitting of the elliptic flow. The scalarisovector interaction leads to the constituent mass splitting of u and d quarks, and the effect is only pronounced around the boundary of the chiral phase transition [22] at the later stage of the partonic phase when the density is rather low, and it has been found to have negligible effects on the isospin splitting of the elliptic flow.

We first fit the parton scattering cross section in order to reproduce the charged particle v 2 at RHIC-BES

/s118

/s84

FIG. 1: (color online) Comparison of the transverse momentum dependence of charged particle elliptic flows at midpseudorapidities ( | η | &lt; 1) in mid-central (20 -30%) Au+Au collisions at RHIC-BES energies from the extended AMPT model with those measured by the STAR Collaboration [43].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-5b14dbda-0d95-4882-9586-5df776bc7e14.jpg)

energies. As shown in Fig. 1, using the same eventplane method [1, 2] as in the experimental analysis, the extended AMPT model with an isotropic parton scattering cross section of 10 mb can reproduce reasonably well the transverse-momentum ( p T ) dependence of midpseudorapidity charged particle v 2 at all RHIC-BES energies. The cross section is larger than that in Ref. [20], as a result of using an improved coalescence approach. As discussed in Ref. [28], a larger cross section is needed to reproduce the same hadron v 2 once combinations of partons close in momentum space is favored, compared with the naive spatial coalescence scenario. We have also compared the p T spectra of mid-rapidity charged pions in central Au+Au collisions at RHIC-BES energies from AMPT calculations with the STAR data in Fig. 2. It is seen that these p T spectra are reproduced reasonably well, except that the AMPT results give slightly stiffer p T spectra at lower energies while smaller multiplicities at higher energies. The agreement of the p T dependence of v 2 as well as the p T spectra serves as a baseline for the following study on the v 2 splitting of charged pions.

We choose minibias Au+Au collisions at √ s NN = 7 7 . GeV as a representative system to discuss the effect of the vector-isovector interaction on the isospin splitting of the elliptic flow. The transverse-momentum dependence of v 2 for light quarks and antiquarks with different isospins are displayed in Fig. 3. The v 2 splitting between quarks and antiquarks is due to the vector-isoscalar potential. Without the vector-isovector interaction ( R IV = 0), u and d quarks as well as ¯ and u ¯ have the same elliptic d flow. With the vector-isovector interaction ( R IV = 2), d quarks have a slightly larger v 2 than u quarks. This shows that the time component of the vector-isovector potential has the dominate effect over the space compo-

/s32

/s32

/s84

FIG. 2: (color online) Comparison of transverse momentum spectra of mid-rapidity ( | y | &lt; 0 1) . π + and π -in central (0 -5%) Au+Au collisions at RHIC-BES energies from the extended AMPT model with those measured by the STAR Collaboration [44].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-5d5ee2da-dcc3-4b15-b58b-dcc9e039100b.jpg)

/s32

/s118

/s32

/s84

FIG. 3: (color online) Transverse momentum dependence of the elliptic flows for mid-pseudorapidity ( | η | &lt; 1) light quarks and antiquarks in minibias Au+Au collisions at √ s NN = 7 7 . GeV from the extended AMPT model with (right) and without (left) the vector-isovector interaction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-b6d5fec9-651d-4e46-8743-1a75cda66dc2.jpg)

nent, leading to a more repulsive (attractive) potential for d ( u ) quarks, as discussed above. Also, ¯ ( ¯ ) are afu d fected by a more repulsive (attractive) potential according to Eq. (2), and have thus a slightly larger (smaller) v 2 in the dominating lowp T region. We note that the increasing parton elliptic flows with transverse momentum are due to parton scatterings as the mean-field potentials mainly lead to a splitting of their elliptic flows.

Since the isospin splitting of light quark v 2 is preserved during hadronization, a larger splitting of charged pion v 2 is obtained as π + and π -are composed of ud ¯ and ¯ ud , respectively. For resonance states with different isospins, e.g., ρ +( -) or ∆ ++( -) , formed from the quark coalescence, they eventually decay into charged pions and carry information of the isospin-dependent potentials on their constituent light quarks. The final elliptic flows of mid-

/s32

/s118

/s32

/s84

FIG. 4: (color online) Transverse momentum dependence of the elliptic flows for mid-pseudorapidity ( | η | &lt; 1) π + and π -in minibias Au+Au collisions at √ s NN = 7 7 GeV from the . extended AMPT model with (right) and without (left) the vector-isovector interaction. The experimental results measured by the STAR Collaboration [45] are also plotted for reference.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-7be96700-e0ad-40cf-8610-2c243066bbc0.jpg)

pseudorapidity π + and π -are displayed in Fig. 4 for the same collision system as in Fig. 3. It is seen that the resulting v 2 as a function of transverse momentum from the extended AMPT model has a similar overall magnitude compared with those measured by the STAR Collaboration. A slightly larger v 2 for π -than π + is observed with the vector-isovector interaction, while the v 2 for π -and π + are similar without the vector-isovector interaction.

To illustrate more clearly the effect of the vectorisovector interaction, we display in the upper panel of Fig. 5 the p T -integrated elliptic flow difference between mid-pseudorapidity π + and π -obtained with different vector-isovector coupling strengths. Without the vectorisovector interaction, although the mean value of the v 2 difference is negative due to the weak hadronic potential, it is consistent with 0 within the statistical error. The v 2 splitting increases with the increasing coupling strength of the vector-isovector interaction until it is about to saturate for R IV &gt; 2. To understand the saturation effect from the vector-isovector interaction, we further show in the lower panel of Fig. 5 the evolution of the time component of the vector-isoscalar potential U 0 V = 2 3 G V ( ρ 0 u + ρ 0 d + ρ 0 s ) as well as the vector-isovector potential U 0 IV = G IV ( ρ 0 u -ρ 0 d ) for different coupling strengths in the central region of the partonic phase. Although the peak value of U 0 IV increases with increasing R IV , it decreases more rapidly for R IV = 5 compared with that for R IV = 2. This is due to the quicker reduction of the isospin asymmetry in the high-density phase as a result of the stronger vector-isovector potential. For the vector-isoscalar potential, it has a much larger magnitude compared to the vector-isovector potential. The saturation effect of the vector-isovector potential is also observed at higher collision energies, where the lifetimes of both potentials become even shorter as a result of a

/s45

/s32

FIG. 5: (color online) Upper: Elliptic flow difference between mid-pseudorapidity ( | η | &lt; 1) π + and π -with different vector-isovector coupling strengths; Lower: Evolution of the time component of the vector-isoscalar potential as well as the vector-isovector potential with different coupling strengths in the central region of the partonic phase. Results are for minibias Au+Au collisions at √ s NN = 7.7 GeV from the extended AMPT model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-dad066d3-0792-46db-8445-05c299c55776.jpg)

thinner thickness and more rapid expansion of the produced partonic matter.

In Fig. 6, the collision energy dependence of the p T -integrated v 2 difference between mid-pseudorapidity π + and π -from the extended AMPT model is compared with those measured by the STAR Collaboration. At each collision energy, tens of thousands events were generated from transport simulations. The decreasing trend of v 2 splitting with increasing collision energy is qualitatively reproduced by the extended AMPT model. Although the mean values of the v 2 difference between π + and π -have a smaller magnitude compared with the experimental data, they have overlaps within statistical errors. Our calculations show that the v 2 difference between π + and π -at RHIC-BES energies favors a strong vector-isovector interaction of R IV &gt; 2. It is worthy to point out that the strong vector-isovector coupling obtained in the present study does not affect the v 2 splitting between nucleons and antinucleons as well as between K + and K -obtained in Ref. [19] since they are due to the vector-isoscalar coupling.

/s45

FIG. 6: (color online) Comparison of the elliptic flow difference between mid-pseudorapidity ( | η | &lt; 1) π + and π -at different RHIC-BES energies from the extended AMPT model with those measured by the STAR Collaboration [9].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01156v2-figure-cdfb6be2-73a0-40d7-958c-02c532711e5d.jpg)

To summarize, we have studied the effect of the vectorisovector interaction on the elliptic flow difference between π + and π -based on the framework of an extended multiphase transport model, in which the partonic meanfield potentials are obtained from a 3-flavor NambuJona-Lasinio model including isovector couplings, and the hadronization algorithm is improved by considering the coalescence of partons close in phase space. We have found that the time component of the vector-isovector potential leads to a more repulsive (attractive) potential for d ( u ) quarks in the baryon-rich and d -quark-rich medium, thus enhancing the v 2 of d quarks and π -while reducing that of u quarks and π + . The space component of the vector-isovector potential and the hadronic potentials are, however, found to have less important effect on the isospin splitting of v 2 . With increasing strength of the vector-isovector coupling, the v 2 difference between π + and π -tends to saturate. Results from our transport approach reproduce the decreasing v 2 splitting with increasing collision energy seen in experiments, and the experimentally observed v 2 difference between π + and π -favors a strong vector-isovector coupling, with its coupling strength larger than twice the scalar-isoscalar coupling strength. In the d -quark-rich matter, a stronger vector-isovector coupling disfavors the splitting of u and d quark chiral phase transition boundary, leads to a critical point at higher temperatures, and results in a larger quark matter symmetry energy [22].

Although the total isospin splitting of pion v 2 can be explained by a strong vector-isovector interaction, its dependence on the charge asymmetry [46] needs further investigations. Also, the isospin-independent polyakovloop potential, which is essential for describing the QCD deconfinement phase transition and the thermodynamic properties of QGP [22], needs to be incorporated in transport simulations. These studies are in progress.

We thank Chen Zhong for maintaining the high-quality performance of the computer facility. J.X. acknowledges support from the Major State Basic Research Development Program (973 Program) of China under Contract No. 2015CB856904 and the National Natural Science

- [1] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58 , 1671 (1998).
- [2] S. A. Voloshin, A. M. Poskanzer, and R. Snellings, arXiv: 0809.2949 [nucl-ex].
- [3] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92 , 052302 (2004).
- [4] J. Adams et al. [STAR Collaboration], Phys. Rev. C 72 , 014904 (2005).
- [5] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 95 , 122301 (2005).
- [6] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98 , 162301 (2007).
- [7] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 99 , 052301 (2007).
- [8] J. Adams et al. [STAR Collaboration], Phys. Rev. C 77 , 054901 (2008).
- [9] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 110 , 142301 (2013).
- [10] J. C. Dunlop, M. A. Lisa, and P. Sorensen, Phys. Rev.C 84 , 044914 (2011).
- [11] V. Greco, M. Mitrovski, and G. Torrieri, Phys. Rev. C 86 , 044905 (2012).
- [12] Yu. B. Ivanov, Phys. Lett. B 723 , 475 (2013).
- [13] J. Steinheimer, V. Koch, and M. Bleicher, Phys. Rev. C 86 , 044903 (2012).
- [14] Y. Hatta, A. Monnai, and B. W. Xiao, Phys. Rev. D 92 , 114010 (2015).
- [15] X. Sun, H. Masui, A. M. Poskanzer, and A. Schmah, Phys. Rev. C 91 , 024903 (2015).
- [16] Y. Burnier, D. E. Kharzeev, J. Liao, and H. U. Yee, Phys. Rev. Lett. 107 , 052303 (2011).
- [17] Y. F. Sun, C. M. Ko, and F. Li, Phys. Rev. C 94 , 045204 (2016).
- [18] Z. Z. Han and J. Xu, Phys. Rev. C 99 , 044919 (2019).
- [19] J. Xu, T. Song, C. M. Ko, and F. Li, Phys. Rev. Lett. 112 , 012301 (2014).
- [20] J. Xu and C. M. Ko, Phys. Rev. C 94 , 054909 (2016).
- [21] P. Danielewicz, R. Lacey, and W. G. Lynch, Science 298 , 1592 (2002).
- [22] H. Liu, J. Xu, L. W. Chen, and K. J. Sun, Phys. Rev. D 94 , 065032 (2016).
- [23] C. Q. Guo, H. Liu, and J. Xu, Phys. Rev. C 98 , 024914 (2018).
- [24] X. N. Wang and M. Gyulassy, Phys. Rev. D 44 , 3501

Foundation of China under Grant Nos. 11922514 and 11421505. C.M.K. acknowledges support from the US Department of Energy under Contract No. de-sc0015266 and the Welch Foundation under Grant No. A-1358.

- (1991).
- [25] C. J. Zhang and J. Xu, Phys. Rev. C 93 , 024906 (2016).
- [26] C. Y. Wong, Phys. Rev. C 25 , 1460 (1982).
- [27] G. F. Bertsch and S. Das Gupta, Phys. Rep. 160 , 189 (1988).
- [28] F. T. Wang and J. Xu, arXiv: 1908.04956 [nucl-th].
- [29] Z. W. Lin, C. M. Ko, B. A. Li, B. Zhang, and S. Pal, Phys. Rev. C 72 , 064901 (2005).
- [30] B. A. Li and C. M. Ko, Phys. Rev. C 52 , 2037 (1995).
- [31] J. Xu, L. W. Chen, C. M. Ko, and Z. W. Lin, Phys. Rev. C 85 , 041901 (2012).
- [32] N. Kaiser and W. Weise, Phys. Lett. B 512 , 283 (2001).
- [33] J. Xu, C. M. Ko, and Y. Oh, Phys. Rev. C 81 , 024910 (2010).
- [34] J. Xu, L. W. Chen, C. M. Ko, B. A. Li, and Y. G. Ma, Phys. Rev. C 87 , 067601 (2013).
- [35] Z. Zhang and C. M. Ko, Phys. Rev. C 95 , 064604 (2017).
- [36] G. 't Hooft, Phys. Rev. D 14 , 3432 (1976); ibid. 18 , 2199(E) (1978).
- [37] M. Lutz, S. Klimt, and W. Weise, Nucl. Phys. A 542 , 521 (1992).
- [38] N. M. Bratovic, T. Hatsuda, and W. Weise, Phys. Lett. B 719 , 131 (2013).
- [39] M. Asakawa and K. Yazaki, Nucl. Phys. A 504 , 668 (1989).
- [40] K. Fukushima, Phys. Rev. D 77 , 114028 (2008); ibid. 78 , 039902 (E) (2008).
- [41] T. Song, S. Plumari, V. Greco, C. M. Ko, and F. Li, arXiv: 1211.5511 [nucl-th].
- [42] C. M. Ko, T. Song, F. Li, V. Greco, and S. Plumari, Nucl. Phys. A 928 , 234 (2014).
- [43] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 86 , 054908 (2012).
- [44] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 96 , 044904 (2017).
- [45] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 88 , 014902 (2013).
- [46] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 114 , 252302 (2015).
- [47] The signs for vector interaction terms are set to be negative to be consistent with the relativistic mean-field model, different from those in Ref. [22].
He Liu Feng-Tao Wang Kai-Jia Sun Jun Xu Che Ming Ko
Isospin splitting of pion elliptic flow in relativistic heavy-ion collisions
Based on the framework of an extended multiphase transport model with
mean-field potentials in both the partonic phase and the hadronic phase, we
explain the elliptic flow difference between $\pi^+$ and $\pi^-$ in the
beam-energy scan program at the relativistic heavy-ion collider by
incorporating the vector-isovector potential for quarks and antiquarks with
different isospins. It is found that the isospin splitting of pion elliptic
flow favors a strong vector-isovector interaction, and thus serves as a probe
of the quark matter equation of state as well as the QCD phase structure at
finite baryon and isospin chemical potentials.