Understanding the Effect of Conserved Charges on the Coalescence Sum Rule of Directed Flow

Abstract

Recently, the rapidity-odd directed flow ($v_1$) of produced hadrons ($K^{-}$, $\varphi$, $\overline{p}$, $\overline{\Lambda }$, ${\overline{\Xi }}^{+}$, ${\Omega }^{-}$, and ${\overline{\Omega }}^{+}$) has been studied. Several combinations of these produced hadrons, with very small mass differences but differences in the net electric charge ($\Delta q$) and net strangeness ($\Delta S$) on the two sides, have been considered. A difference in $v_1$ between the two sides of these combinations ($\Delta v_1$) has been proposed as a consequence of the electromagnetic field produced in relativistic heavy-ion collisions, especially if $\Delta v_1$ increases with $\Delta \mathrm{q}$. Our study is performed to understand the effect of the coalescence sum rule (CSR) on $\Delta v_1$. We point out that the CSR leads to $\Delta v_1=c_q\Delta q+c_S\Delta S$, where the coefficients $c_q$ and $c_S$ reflect the $\Delta v_1$ of produced quarks. Equivalently, one can write $\Delta v_1=c_q\Delta q+c_B\Delta B$, involving the difference in the net baryon number $\Delta B$, where the CSR gives $c_B=-3c_S$. We then propose two methods to extract the coefficients for the $\Delta q$ and $\Delta S$ dependences of $\Delta v_1$.

1. Introduction

The properties of strongly interacting quark-gluon plasma, produced by relativistic heavy-ion collisions, can be studied using anisotropic flows, including the directed flow ($v_1$) [1,2,3]. $v_1$ has been found to be a sensitive probe of the equation of state of the produced medium [4,5]. The anisotropic flows of identified hadrons are expected to follow the coalescence sum rule (CSR) when the produced matter is initially in parton degrees of freedom and hadronizes via quark coalescence. The CSR states that the sum of constituent quarks’ flow equals the corresponding flow of the hadron; it leads to the number-of-constituent-quark scaling and provides information about the particle production mechanism and partonic collectivity [6,7]. This paper is based on our recent study [8], which was motivated by Ref. [9].

A new method of testing the coalescence sum rule using the difference in $v_1$ in different combinations of produced hadrons [$K^{-}\left(\overline{u}s\right)$, $\varphi \left(s\overline{s}\right)$, $\overline{p}\left(\overline{u}\overline{u}\overline{d}\right)$, $\overline{\Lambda }\left(\overline{u}\overline{d}\overline{s}\right)$, ${\overline{\Xi }}^{+}(\overline{d}\overline{s}\overline{s}$), ${\Omega }^{-}\left(sss\right)$ and ${\overline{\Omega }}^{+}\left(\overline{s}\overline{s}\overline{s}\right)$] has been previously proposed [9]. These selected hadrons are all produced, consisting of $\overline{u}$, $\overline{d}$, s, and $\overline{s}$ quarks. In contrast to the produced hadrons, transported hadrons such as ${\pi }^{\pm }$, p, and $\Lambda$ receive contributions from initial-state (incoming) u and d quarks, along with the produced quarks, which complicates the interpretations of the CSR. Table 1 shows a selection of five independent hadron sets [8]. Sets 1–3 are identical, with the same quark contents on both sides, whereas sets 4–5 are non-identical. Indeed, different choices of five independent hadron sets can be made [8]. One can also obtain all ten sets from Table 1 in Ref. [9]; however, they are not all independent. In a given set, the difference in $v_1$ between the left and right sides (after including the weighting factors of each hadron as listed in Table 1) is termed $\Delta v_1$. Similarly, the differences in the net electric charge ($\Delta q$ or $\Delta q_{ud}$) and net strangeness ($\Delta S$) are obtained.

As $v_1$ develops in the early stage of collisions, it is sensitive to the strong electromagnetic field produced by incoming protons in the two colliding nuclei. The correlation between the difference in $v_1$ on $\Delta q$ has been considered as breaking the coalescence sum rule and being a consequence of the electromagnetic fields [9,10], especially when the difference in $v_1$ increases with $\Delta q$. This motivated us to critically examine the consequence of the coalescence sum rule on the difference in $v_1$ [8]. This paper is organized as follows. The derivation of the CSR relation is given in Section 2. Two methods for extracting the dependence of $\Delta v_1$ on $\Delta q$ and $\Delta S$ are presented in Section 3. A summary is finally given in Section 4.

2. The CSR Relation for the Difference in ${\mathit{v}}_1$ of a Hadron Set

The relation between the $v_1$ of a hadron (H) and those of its constituent quarks from the CSR can be written as

$$v_{1,\mathrm{H}}\left(p_{\mathrm{T},\mathrm{H}}\right)=\sum _j{v}_{1,j}\left(p_{T,j}\right),$$

(1)

which is a sum over each constituent quark j of the hadron. The simplest case for the coalescence of comoving equal-mass quarks gives the usual relation [7] $v_{1,\mathrm{H}}\left(N_{cq}{p}_{\mathrm{T},\mathrm{q}}\right)=N_{cq}{v}_{1,q}\left(p_{T,q}\right)$, where $N_{cq}$ is the number of constituent quarks of the hadron.

For each of the five hadron sets in

Table 1. List of five independent hadron sets, including identical ($\Delta q=0$ and $\Delta S=0$) and non-identical ($\Delta q\ne 0$ and/or $\Delta S\ne 0$) constituent quark combinations with similar mass.

Set	$\mathbf{\Delta }\mathit{q}$	$\mathbf{\Delta }{\mathit{q}}_{\mathit{ud}}$	$\mathbf{\Delta }\mathit{S}$	$\mathbf{\Delta }\mathit{B}$	Left Side (L)	Right Side (R)
1	0	0	0	0	$v_1\left[K^{-}\left(\overline{u}s\right)\right]+v_1\left[\overline{\Lambda }\left(\overline{u}\overline{d}\overline{s}\right)\right]$	$v_1\left[\overline{p}\left(\overline{u}\overline{u}\overline{d}\right)\right]+v_1\left[\varphi \left(s\overline{s}\right)\right]$
2	0	0	0	0	$v_1\left[\overline{\Lambda }\left(\overline{u}\overline{d}\overline{s}\right)\right]$	$\frac{v_1}{2}\left[{\overline{\Xi }}^{+}\left(\overline{d}\overline{s}\overline{s}\right)\right]+\frac{v_1}{2}\left[\overline{p}\left(\overline{u}\overline{u}\overline{d}\right)\right]$
3	0	0	0	0	$\frac{v_1}{3}\left[{\Omega }^{-}\left(sss\right)\right]+\frac{v_1}{3}\left[{\overline{\Omega }}^{+}\left(\overline{s}\overline{s}\overline{s}\right)\right]$	$v_1\left[\varphi \left(s\overline{s}\right)\right]$
4	1/3	0	1	−1/3	$\frac{v_1}{2}\left[\varphi \left(s\overline{s}\right)\right]$	$\frac{v_1}{3}\left[{\Omega }^{-}\left(sss\right)\right]$
5	2/3	1/3	1	−1/3	$\frac{v_1}{2}\left[\varphi \left(s\overline{s}\right)\right]+\frac{v_1}{3}\left[\overline{p}\left(\overline{u}\overline{u}\overline{d}\right)\right]$	$v_1\left[K^{-}\left(\overline{u}s\right)\right]$

, the number of quarks of flavor i ($N_i$) for each side (L: left; R: right) is calculated as the sum of the number of quarks of flavor i in each hadron, multiplied by the weighting factor. With $\Delta N_i=N_i^L-N_i^R$, each set satisfies the following relations:

$$\Delta N_{\overline{u}}+\Delta N_{\overline{d}}=0,\Delta N_s+\Delta N_{\overline{s}}=0.$$

(2)

In other words, the two sides have the same number of $\overline{u}+\overline{d}$ quarks and the same number of $s+\overline{s}$ quarks. However, they may have different $N_{\overline{u}}$ and/or $N_s$. Therefore, they can have a different total net electric charge (q), total net strangeness (S), or total net electric charge in light quarks ($q_{ud}$), defined as $\Delta q=q^L-q^R,\Delta S=S^L-S^R$, and $\Delta q_{ud}=q_{ud}^L-q_{ud}^R$, respectively. For the hadron sets in Table 1, Equation (2) leads to

$$\Delta q_{ud}=\Delta N_{\overline{d}},\Delta S=2\Delta N_{\overline{s}},\Delta q=\Delta q_{ud}+\Delta S/3.$$

(3)

The total $v_1$ from each side of a hadron set can be written as $v_1^{L,R}={\sum }_i{N}_i^{L,R}{v}_{1,\mathrm{i}}$. It is important to note that we neglect the effect of different constituent quark masses in this study. Therefore, the above relation applies to light and strange (anti)quarks at the same $p_{\mathrm{T}}$ value or $p_{\mathrm{T}}$ range, and consequently, it also applies to baryons at a $p_{\mathrm{T}}$ value or range that is 3/2 times that of the meson’s. The difference in $v_1$ between the two sides is then given by

$$\begin{array}{ccc}\hfill \Delta v_1& =& v_1^L-v_1^R=\sum _i\Delta N_i{v}_{1,\mathrm{i}}\hfill \end{array}$$

$$\begin{array}{ccc}& =& (v_{1,\overline{d}}-v_{1,\overline{u}})\Delta q_{ud}+\left(\frac{v_{1,\overline{s}}-v_{1,s}}{2}\right)\Delta S\hfill \end{array}$$

(4)

$$\begin{array}{ccc}& =& \left(v_{1,\overline{d}}-v_{1,\overline{u}}\right)\Delta q+\left(\frac{v_{1,\overline{s}}-v_{1,s}}{2}-\frac{v_{1,\overline{d}}-v_{1,\overline{u}}}{3}\right)\Delta S.\hfill \end{array}$$

(5)

Equations (4) and (5) show the linear dependence of $\Delta v_1$ of a hadron set on both $\Delta q$ and $\Delta S$, where the coefficients are simply determined by the quark-level differences in $v_1$. We can also see that the coefficient for the $\Delta S$ dependence in Equation (5) is not as clean as the corresponding coefficient in Equation (4).

The difference in $v_1$ in the hadron sets also depends on $\Delta B$, which represents the difference in the net baryon number between the two sides of a hadron set. Under the condition from Equation (2), we have

$$\Delta B=-2\Delta N_{\overline{s}}/3=-\Delta S/3.$$

(6)

Therefore, $\Delta v_1$ depends linearly on both $\Delta q$ and $\Delta B$ as follows:

$$\begin{array}{ccc}\hfill \Delta v_1& =& (v_{1,\overline{d}}-v_{1,\overline{u}})\Delta q_{ud}-3\left(\frac{v_{1,\overline{s}}-v_{1,s}}{2}\right)\Delta B\hfill \end{array}$$

(7)

$$\begin{array}{ccc}& =& \left(v_{1,\overline{d}}-v_{1,\overline{u}}\right)\Delta q-3\left(\frac{v_{1,\overline{s}}-v_{1,s}}{2}-\frac{v_{1,\overline{d}}-v_{1,\overline{u}}}{3}\right)\Delta B.\hfill \end{array}$$

(8)

Assuming that the rapidity of a hadron formed by quark coalescence is the same as that of the coalescing quarks, the $v_1$ slope at $y=0$, denoted as $v_1^{\prime }=d{v}_1/dy(y=0)$, satisfies exactly the same relations as Equations (4), (5), (7) and (8), where one just needs to replace $v_1$ with $v_1^{\prime }$ [8].

3. Extracting the $\Delta \mathit{q}$ and $\Delta \mathit{S}$ Dependences of the Coefficients

The coalescence sum rule may not be satisfied for certain collision energies and/or systems, e.g., in cases where $v_1$ is not dominated by parton dynamics. Since Equations (4) and (5) from the coalescence sum rule predict $\Delta v_1$ = 0 for hadron sets with $\Delta S=\Delta q_{ud}$ = 0 or $\Delta S=\Delta q$ = 0, we propose that the following modified equations are used to fit the $\Delta v_1$ data from the five independent sets:

$$\begin{array}{ccc}\hfill \Delta v_1& =& c_0+c_q\Delta q_{ud}+c_S\Delta S,\hfill \end{array}$$

(9)

$$\begin{array}{ccc}\hfill \Delta v_1& =& c_0^{✶}+c_q^{✶}\Delta q+c_S^{✶}\Delta S.\hfill \end{array}$$

(10)

This way, a non-zero value of the intercept parameter, denoted as $c_0$ or $c_0^{✶}$, would mean breaking the coalescence sum rule. Note that even for a given collision system, the coalescence sum rule may be satisfied around midrapidity but not satisfied at large rapidities; therefore, these coefficients are rapidity-dependent.

To extract the coefficients of the $\Delta q$ and $\Delta S$ dependences of $\Delta v_1$ for a given collision system, one can utilize the five-set method by simply fitting the 5 data points using Equation (9) or Equation (10). It should be noted that the fit function represents a two-dimensional plane over the $\Delta q-\Delta S$ space. Therefore, one should not simply fit the $\Delta v_1$ data with a one-dimensional function of $\Delta q$ without taking into account the different $\Delta S$ values of different hadron sets.

Alternatively, we can use the three-set method [8], where we take the average of sets 1–3, because they have the same $\Delta q$ and $\Delta S$ values, and denote the average as set A. Let us also denote set 4 and set 5 as set B and set C, respectively. Equation (9) then leads to ($\Delta v_1{)}_A$ = $c_0$, ($\Delta v_1{)}_B$ = $c_0$ + $c_S$, and ($\Delta v_1{)}_C$ = $c_0$ + $c_q/3$ + $c_S$. Thus, we obtain

$$\begin{array}{c}\hfill c_0={(\Delta v_1)}_A;c_q=3[{(\Delta v_1)}_C-{(\Delta v_1)}_B];c_S={(\Delta v_1)}_B-{(\Delta v_1)}_A.\end{array}$$

(11)

Similarly, the three-set method yields the coefficients in Equation (10) as

$$\begin{array}{c}\hfill \begin{array}{c}\hfill c_0^{★}={(\Delta v_1)}_A=c_0;c_q^{★}=3[{(\Delta v_1)}_C-{(\Delta v_1)}_B]=c_q;\\ \hfill c_S^{★}=2{(\Delta v_1)}_B-{(\Delta v_1)}_A-{(\Delta v_1)}_C=c_S-c_q/3.\end{array}\end{array}$$

(12)

We see that the three-set method has an advantage over the five-set method in that we can extract the $\Delta q$ and $\Delta S$ coefficients analytically instead of performing a two-dimensional fit. It should be noted that the CSR predicts the following coefficients [8]:

$$\begin{array}{ccc}& & c_0=c_0^{✶}=0;c_q=c_q^{★}=v_{1,\overline{d}}-v_{1,\overline{u}};\hfill \end{array}$$

(13)

$$\begin{array}{ccc}& & c_S=(v_{1,\overline{s}}-v_{1,s})/2,c_S^{★}=(v_{1,\overline{s}}-v_{1,s})/2-(v_{1,\overline{d}}-v_{1,\overline{u}})/3.\hfill \end{array}$$

(14)

When the difference in $v_1$ of the hadron sets is expressed in terms of $\Delta B$, we can use the following modified equations, which are similar to Equations (9) and (10):

$$\begin{array}{ccc}\hfill \Delta v_1& =& c_0+c_q\Delta q_{ud}+c_B\Delta B,\hfill \end{array}$$

(15)

$$\begin{array}{ccc}\hfill \Delta v_1& =& c_0^{✶}+c_q^{✶}\Delta q+c_B^{✶}\Delta B.\hfill \end{array}$$

(16)

Here, the CSR predicts the following $\Delta B$ coefficients:

$$c_B=-3(v_{1,\overline{s}}-v_{1,s})/2=-3c_S,c_B^{★}=-3(v_{1,\overline{s}}-v_{1,s})/2+(v_{1,\overline{d}}-v_{1,\overline{u}})=-3c_S^{★}.$$

(17)

Obviously, the CSR relates the $\Delta q$ and $\Delta S$ (or $\Delta B$) coefficients to the quark-level differences in $v_1$. It should be noted that the quark $v_1$ in Equations (13), (14) and (17) refers to its value after the partonic evolution just before the quark coalescence. Therefore, these coefficients could be non-zero due to the flavor dependence of the strong interaction or the effect of the electromagnetic field on the produced quarks. Additionally, the extraction methods outlined in this section apply to the $v_1^{\prime }$ data in exactly the same way. Furthermore, in the full study [8], we demonstrated these extraction methods for the $\Delta q$ and $\Delta S$ coefficients using results from the AMPT model for Au+Au collisions at different energies.

4. Summary

In this study, we derived the relations for $\Delta v_1$, representing the difference in $v_1$ between the two sides of a hadron set, using the coalescence sum rule. In earlier studies, seven produced hadron species ($K^{-}$, $\varphi$, $\overline{p}$, $\overline{\Lambda }$, ${\overline{\Xi }}^{+}$, ${\Omega }^{-}$, and ${\overline{\Omega }}^{+}$) were considered, and a non-zero $\Delta v_1$ dependence on the difference in the net electric charge ($\Delta q$) was considered as breaking the coalescence sum rule and being a consequence of the electromagnetic fields. Our study showed that the coalescence sum rule only leads to a zero $\Delta v_1$ for a hadron set if its two sides have identical constituent quark contents (i.e., $\Delta q$ = $\Delta S$ = 0). In general, $\Delta v_1$ depends linearly on both $\Delta q$ and $\Delta S$, or on both $\Delta q_{ud}$ and $\Delta S$, and the same applies to the difference in the $v_1$ slopes at midrapidity ($\Delta v_1^{\prime }$). Since there are only five such independent hadron sets, there will be five independent data points from the measurement of a given collision system. We propose using a two-dimensional function such as $c_0+c_q\Delta q_{ud}+c_S\Delta S$ to extract the coefficients of the $\Delta q_{ud}$ and $\Delta S$ dependences, where a non-zero value of the coefficient $c_0$ indicates the breaking of the coalescence sum rule. In the five-set method, one fits the five data points using this function. On the other hand, in the three-set method, the three data points with the same $\Delta q$ and $\Delta S$ values are combined into one set, allowing us to obtain the coefficients analytically. The coalescence sum rule relates these coefficients to the quark-level differences in $v_1$ just before quark coalescence. Therefore, the coefficients can be affected by both the strong interaction and the electromagnetic fields. Equivalently, one can also express $\Delta v_1$ or $\Delta v_1^{\prime }$ in terms of the net baryon difference $\Delta B$, where the function $c_0+c_q\Delta q_{ud}+c_B\Delta B$ can be used to extract the coefficients, and the coalescence sum rule yields $c_B=-3c_S$.