-
In this study, we use the string melting AMPT model, which can give a better description of the correlation function in HBT research [23, 24]. The string melting AMPT model contains four main parts. The first part uses the HIJING model to produce the partons and strings, where the strings will fragment into partons. The second part uses the Zhang’s parton-cascade (ZPC) model to describe the interactions among these partons. Then, in the third part, a quark coalescence model is used to combine these partons into hadrons. In the last part, the interaction of these hadrons is described by a relativistic transport (ART) model till the hadrons freeze out. With the models in these four main parts, the string melting AMPT model can be used to describe heavy-ion collisions.
The HBT radii are important in HBT research, and they can be extracted by the HBT three-dimensional correlation function, whereas the normal form can be written as [7]
$ \begin{array}{l} C({\boldsymbol{q}},{\boldsymbol{K}})=1 + \lambda{{\rm e}}^ {-q_{\rm{ o}}^2 R_{\rm{ o}}^2({{\boldsymbol{K}}})-q_{\rm {s}}^2 R_{\rm {s}}^2({{\boldsymbol{K}}})- q_{\rm{ l}}^2 R_{\rm{ l}}^2({{\boldsymbol{K}}})-2q_{\rm{ o}}q_{\rm{ l}} R_{\rm{ ol}}^2({{\boldsymbol{K}}})}. \end{array} $
(1) In addition, in this study, we use the form with Coulomb interaction as [25]
$ \begin{array}{l} C({\boldsymbol{q}},{\boldsymbol{K}}) = 1 - \lambda + \lambda K_{{\rm coul}}\big[ 1 + {{\rm e}}^ {-q_{\rm{ o}}^2 R_{\rm{ o}}^2({{\boldsymbol{K}}})-q_{\rm {s}}^2 R_{\rm {s}}^2({{\boldsymbol{K}}})- q_{\rm{ l}}^2 R_{\rm{ l}}^2({{\boldsymbol{K}}})-2q_{\rm{ o}}q_{\rm{ l}} R_{\rm{ ol}}^2({{\boldsymbol{K}}})}\big], \end{array} $
(2) where
$ {{{\boldsymbol{q}}}}={{{\boldsymbol{p}}}}_1-{{{\boldsymbol{p}}}}_2 $ ,$ {{{\boldsymbol{K}}}}=({{{\boldsymbol{p}}}}_1+{{{\boldsymbol{p}}}}_2)/2 $ , C is the two pion correlation function,$ K_{{\rm coul}} $ is the squared Coulomb wave function, and λ is the coherence parameter. R is for the HBT radius. It is usually used in the 'out-side-long' coordinate system, and this system is used for the pair particles, as shown in Fig. 2. l represents the longitudinal direction and beam direction. o and s are outward and sideward directions that are both defined on the transverse plane. The pair particles momentum direction is defined as the outward direction, and the side direction is the direction perpendicular to the outward direction.$ R_{\rm{ ol}} $ is the cross term, which will vanish at mid-rapidity in a symmetric system. In this study, we set the biggest impact parameter as 1.4 fm and chose the mid-rapidity range$ -0.5<\eta<0.5 $ . We use Eq. (1) to fit the correlation function without Coulomb interaction and Eq. (2) to fit the correlation function with Coulomb interaction.We use the Correlation After Burner (CRAB) code to generate the correlation function, which can read the phase space information from the string melting AMPT model [26]. It generates the correlation function by the formula
$ \begin{array}{l} C({\boldsymbol{q}},{\boldsymbol{K}})=1+\frac{\displaystyle\int {\rm d}^4x_1 {\rm d}^4x_2 S_1(x_1,{{{\boldsymbol{p}}}}_2) S_2(x_2,{{{\boldsymbol{p}}}}_2) {\left|\psi_{\rm{rel}} \right|}^{2}} {\displaystyle\int {\rm d}^{4}x_{1}{\rm d}^{4}x_{2}S_{1}(x_{1},{{{\boldsymbol{p}}}}_{2})S_{2}(x_{2},{{{\boldsymbol{p}}}}_{2})}, \end{array} $
(3) where
$ \psi_{\rm{rel}} $ is the relative two particles wave function, which includes the interaction between two particles.$ S(x,{{{\boldsymbol{p}}}}) $ is the single particle emission function. -
We produced the particles of Au+Au collisions at
$ \sqrt{S_{NN}}= $ 14.5, 19.6, 27, 39, 62.4, 200 GeV by the string melting AMPT model. These are the BES energies, each one of them has more than sixty thousand collision events. We separate the pair transverse momentum 125–625 MeV/c into 9 bins. While the last bin has fewer particles, we set it for 525–625 MeV/c, and the width is twice that of the others. Then we use the CRAB code to calculate the correlation functions for the$ \pi^+ $ . The Coulomb interaction can change the correlation functions, which are shown in Fig. 3. If the Coulomb interaction is neglected, because of the collective flow created by the collisions, two$ \pi^+ $ are easily frozen out along the same direction, and the correlation function will be a Gaussian form. While considering the Coulomb interaction, the two$ \pi^+ $ have the Coulomb potential between them, and it will change the two$ \pi^+ $ wave function$ \psi_{\rm{rel}} $ . The closer are the two$ \pi^+ $ , the bigger the Coulomb potential between them will be. When two$ \pi^+ $ are close to each other, the probability density of the two$ \pi^+ $ will decrease, and will cause the decrease of the Correlation function in the direction along the momentum difference direction.$ q_{\rm {s}} $ is related to the momentum difference on the transverse plane, so there is a little gap in the correlation function in the$ q_{\rm {s}} $ direction.Figure 3. (color online) Correlation functions of
$ \pi^+ $ in the$ q_{\rm{ o}} $ and$ q_{\rm {s}} $ directions, generated by the string melting AMPT model for the Au+Au collisions at$ \sqrt{S_{NN}}= $ 200 GeV. The$ K_{\rm {T}} $ range is 325–375 MeV/c, and the$ q_{\rm{ l}} $ range is -3–3 MeV/c. Figure (a) neglects the Coulomb interaction, and figure (b) is the situation with the Coulomb interaction.We only focus on the HBT radius
$ R_{\rm {s}} $ ; it is least affected by other physical factors and related to the transverse size of the source [27]. In addition, the transverse momentum dependence of$ R_{\rm {s}} $ with and without Coulomb interaction is shown in Fig. 4. The obvious difference are the values of$ R_{\rm {s}} $ calculated with the Coulomb interaction, which are smaller than those without the Coulomb interaction at small pair momentum. Besides, the strength of$ K_{\rm {T}} $ dependence of$ R_{\rm {s}} $ , especially for the lower collision energies, can barely be distinguished. And it is even worse for the situation with the Coulomb interaction. At low collision energies, there are similar collective flow velocities in each$ K_{\rm {T}} $ bin, which leads to similar single-particle momentum-space angle distributions [19]. Moreover, the single-particle momentum-space angle distributions can affect the HBT radii [18], so the values of$ R_{\rm {s}} $ are similar. To distinguish the strength of$ K_{\rm {T}} $ dependence of$ R_{\rm {s}} $ , we introduced a parameter b in the fit function$ \begin{array}{l} R_{\rm {s}} = a K_{\rm {T}}^{b}, \end{array} $
(4) where parameter a is just a common constant. The parameter b can be used to describe the strength for the
$ K_{\rm {T}} $ dependence of$ R_{\rm {s}} $ , as shown in Fig. 5.In Fig. 5, All the
$ |b| $ values increase with increasing collision energies. It indicates that, at high collision energies, the changes in$ R_{s} $ values are more intense with the transverse pair momenta. The black dots for all pions are taken from our last work, which treated all pions as one kind, generated with the 0 fm impact parameter. Comparing to the 0–1.4 fm for all pions, the$ |b| $ values decrease with a higher impact parameter. The latter can decrease the HBT radii [28], which indicates that the impact parameter has different influences on$ R_{\rm {s}} $ in different$ K_{\rm {T}} $ sections. Comparing$ \pi^+ $ without Coulomb interaction and all pions, shows that the kind of pions can also affect the strength of the$ K_{\rm {T}} $ dependence of$ R_{\rm {s}} $ . The$ |b| $ values are smaller with the Coulomb interaction, which means the Coulomb interaction inhibits the strength of the transverse momentum dependence of HBT radii. Furthermore, with increasing the collision energies, the differences of the$ |b| $ values for$ \pi^+ $ between the two situations are also increased, which indicates this inhibition is stronger at high collision energies for the Coulomb interaction. -
In HBT research, the correlation function is often approximated by the on-shell momenta [29, 30]. So the correlation can also be written as [31]
$ \begin{array}{l} C({\boldsymbol{q}},{\boldsymbol{K}})\approx 1+\dfrac{\bigg|\displaystyle\int {\rm d}^4x S(x,K)) {\rm e}^{{\rm i} {{\boldsymbol{q}}} \cdot {{\boldsymbol{x}}}}\bigg|^2} {\bigg|\displaystyle\int {\rm d}^{4}xS(x,K)\bigg|^2} \equiv 1+|<{\rm e}^{{\rm i} {{\boldsymbol{q}}} \cdot {{\boldsymbol{x}}}}>|^2, \end{array} $
(5) where
$ K_0=E_K=\sqrt{(m^2+|{{\boldsymbol{K}}}|^2)} $ ,$ {{\boldsymbol{q}}} = {{\boldsymbol{p}}}_1-{{\boldsymbol{p}}}_2 $ , then we get$ \begin{array}{l} {\rm e}^{{\rm i}{{\boldsymbol{q}}} \cdot {{\boldsymbol{x}}}}=\exp[{\rm i}({\boldsymbol{p}}_1 \cdot {\boldsymbol{x}}-{\boldsymbol{p}}_2 \cdot {\boldsymbol{x}})], \end{array} $
(6) for the single-particle, we have already defined the angle
$ \Delta\varphi $ , which is between the momentum direction and space direction, then we define two more angles,$ \cos(\Delta\alpha) \equiv \frac{{{\boldsymbol{p}}}_2 \cdot {{{\boldsymbol{x}}}}}{|{{\boldsymbol{p}}}_2||{{{\boldsymbol{x}}}}|}, $
(7) $ \cos(\Delta\beta) \equiv \frac{{{\boldsymbol{p}}}_1 \cdot {{\boldsymbol{p}}}_2}{|{{\boldsymbol{p}}}_1||{{\boldsymbol{p}}}_2|}, $
(8) where
$ \Delta\alpha $ is the angle between the momentum direction and radius direction for two particles, and$ \Delta\beta $ is the angle between the momentum directions for two particles. The smoothness assumption is$ |{{\boldsymbol{p}}}_1| \approx |{{\boldsymbol{p}}}_2| \approx\frac{1}{2} |{{\boldsymbol{K}}}| $ [32], then that$ \begin{array}{l} {\rm e}^{{\rm i}{{\boldsymbol{q}}} \cdot {{\boldsymbol{x}}}}\approx\exp\left[\dfrac{{\rm i}}{2}\left(|{{\boldsymbol{K}}}||x|\cos(\Delta\varphi) - |{{\boldsymbol{K}}}||x|\cos(\Delta\alpha)\right)\right], \end{array} $
(9) the
$ \Delta\alpha $ angle is not independent, its values are related to the$ \Delta\beta $ and the$ \Delta\varphi $ ,$ \Delta\alpha=f(\Delta\varphi,\Delta\beta) $ , so the correlation function can be written as$ \begin{aligned}[b] C({{\boldsymbol{K}}}, \Delta\varphi, \Delta\beta) \approx & 1+\Big|\bigg<\exp\bigg[ \frac{{\rm i}}{2}|{{\boldsymbol{K}}}||{{\boldsymbol{x}}}|[\cos(\Delta\varphi)\\&-\cos(f(\Delta\varphi,\Delta\beta))] \bigg]\bigg >\Big|^2. \end{aligned} $
(10) When we calculate the HBT correlation function, we need to limit the values of momentum difference,
$ |{{\boldsymbol{q}}}|<q_{\rm {max}} $ . If we choose the pair particles in one transverse momentum section, and because we focus on the mid-rapidity range, the pair momenta are also screened,$ k_{\rm {min}}<|{{\boldsymbol{K}}}|<k_{\rm {max}} $ . At the same time, the angle$ \Delta\beta $ is also limited,$ \cos(\Delta\beta)>\dfrac{2k_{\rm {min}}^2}{k_{\rm {min}}^2+q_{\rm {max}}^2}-1 $ . So the$ \Delta\beta $ values are related to the$ |{{\boldsymbol{K}}}| $ values and the interaction betweenthe particles. Therefore, the single-particle space-momentum angle distribution can directly affect the correlation function. Our research is focused on the transverse plane, so we discuss the connection between the$ \Delta\theta $ distribution and$ R_{\rm {s}} $ .For a source, the collective flow is created by the expansion, and the transverse flow velocities are different at different layers of the source, which leads to different normalized
$ \Delta\theta $ distributions, as shown in Fig. 6.Figure 6. (color online) The normalized
$ \Delta\theta $ distributions in different transverse radii,$ \pi^{+} $ are generated by the string melting AMPT model for the Au+Au collisions at$ \sqrt{S_{NN}}= $ 62.4 GeV/c.In Fig. 6, the average transverse flow velocity can be calculated by
$ <v_{\rm {T}}> = \left< \frac{{{\boldsymbol{p}}}_{\rm {T}} \cdot {{\boldsymbol{r}}}_{\rm {T}}}{E|{{\boldsymbol{r}}}_{\rm {T}}|} \right>. $
(11) $N$ and$N_{\rm R}$ are the particle numbers in each bin, whereas$N_{\rm R}$ is obtained with the random p and r particles. The normalization process is using the$ \cos(\Delta\theta) $ distribution divided by the random$ \cos(\Delta\theta) $ distribution, and we let$\sum N =\sum N_{\rm R}$ . The particles with lower transverse momenta are closer to$ N/N_{\rm R} = 1 $ . This phenomenon indicates that the source is approaching a random freeze-out source. While with higher flow velocities, the distributions are closer to$ \cos(\Delta\theta)=1 $ , it means the particles tend to freeze out along the radius direction. For the particles located at 8–9 and 11–12 fm, their$ <\nu_{\rm {T}}> $ values are closer, so they have similar normalized$ \Delta\theta $ distributions.Furthermore, for the particles with different transverse momenta
$ p_{\rm {T}} $ , they correspond to the parts of the whole source. For the particles with higher$ p_{\rm {T}} $ , the part sources have bigger collective flow, as shown in Fig. 7. The$ <\nu_{\rm {T}}> $ increase with$ r_{\rm {T}}=\sqrt{x^2+y^2} $ inside of the source, and at the outside of the source,$ <\nu_{\rm {T}}> $ almost becomes a constant. So the sizes of part sources are changing with the$ p_{\rm {T}} $ .Figure 7. (color online) The transverse flow velocities for the transverse radii in different transverse momenta,
$ \pi^{+} $ are generated by the string melting AMPT model for the Au+Au collisions at$ \sqrt{S_{NN}}= $ 62.4 GeV/c.Before particle freezing out from the source, the system goes through a series of processes, containing parton production and interaction, hadronization, and hadron cascade. The expansion of the source will be reflected in the collective flow of the freeze-out particles. The collective flow is different at different locations in the source, meanwhile, the
$ \Delta\theta $ is related to the flow. Thus, different$ \Delta\theta $ have different freeze-out sources, as shown in Fig. 8.Figure 8. (color online) The freeze-out
$ \pi^+ $ sources in different$ \cos(\Delta\theta) $ sections. The particles are generated by the string melting AMPT model for the Au+Au collisions at$ \sqrt{S_{NN}}= $ 62.4 GeV, and the$ p_{\rm {T}} $ range is 250–450 MeV/c.For the
$ \pi^+ $ that have large$ \Delta\theta $ angles, their freeze-out directions are reversed to the space directions. They have a small freeze-out source, which indicates they are almost freezing out from the center of the whole$ \pi^+ $ source. Those$ \pi^+ $ freeze-out directions are perpendicular to the space directions and have a bigger source, as shown in Fig. 8(b). While in Fig. 8(c), for the$ \pi^+ $ whose momentum directions are along to the space directions, their freeze-out source has a ring shape, and has the biggest number density. So in this$ p_{\rm {T}} $ section, most of the$ \pi^+ $ freeze out from the shell. The closer to the center of the source, the higher the probability of producing large$ \Delta\theta $ angle particles. For we limited the differences of the pair particles$ q_{\rm max} $ , the normalized$ \cos(\Delta\theta) $ distribution in$ K_{\rm {Tmin}}<K_{\rm {T}}<K_{\rm {Tmax}} $ section only contains$ K_{\rm Tmin}- q_{\rm max}/2 <p_{\rm {T}}< K_{\rm Tmax}+q_{\rm max}/2 $ particles. Besides, different$ \cos(\Delta\theta) $ in the$ K_{\rm {T}} $ section has different partial sources, similar to Fig. 8. The normalized$ \cos(\Delta\theta) $ distribution corresponds to the superposition of these partial sources. Therefore, the distribution is related to the HBT radius$ R_{\rm {s}} $ .If there is no correlation between the space and momentum, the relation between
$ \Delta\theta $ and the source will break, and the differences in the normalized$ \cos(\Delta\theta) $ distributions will disappear. Then the$ \Delta\theta $ angle will be completely random, as shown in Fig. 9. The normalized "no x-p correlation" distribution becomes a line and is the same in other$ K_{\rm {T}} $ sections.Figure 9. (color online) The original distribution and the no x-p correlation distribution at
$ \sqrt{S_{NN}}= $ 62.4 GeV, and the$ K_{\rm {T}} $ range is 325–375 MeV/c.In Fig. 9, we disrupt the space and momentum correlation in each collision event. For example, original data are
$ x_{1}p_{1} $ and$ x_{2}p_{2} $ , and after the disruption, they are$ x_{1}p_{2} $ and$ x_{2}p_{1} $ . This is a rough but effective method, and it can change the correlation between the space and the momentum, which means the particles can freeze out from the source in any direction. Then we calculate the HBT radius$ R_{\rm {s}} $ , as shown in Fig. 10. The results show that the phenomenon of transverse dependence of$ R_{\rm {s}} $ for "no x-p correlation" has almost disappeared. The$ R_{\rm {s}} $ values for "no x-p correlation" are closer to each other, because they have the same normalized$ \cos(\Delta\theta) $ distribution. This method destroys the physical process, there will be fluctuations in calculating the correlation function, which leads to the emergence of the big error bars.Figure 10. (color online) The original
$ R_{\rm {s}} $ and no x-p correlation$ R_{\rm {s}} $ at$ \sqrt{S_{NN}}= $ 62.4 GeV.Moreover, we can build a numerical connection between the normalized
$ \cos(\Delta\theta) $ distribution and HBT radius$ R_{\rm {s}} $ . The AMPT model does not contain the Coulomb interaction at the stage of the hadron cascade [23], so we discuss the Coulomb interaction only for the HBT analysis. Moreover, the two conditions are using the same string melting AMPT data, leading to the same normalized$ \cos(\Delta\theta) $ distribution for different$ K_{\rm {T}} $ regions. In order to describe the normalized$ \cos(\Delta\theta) $ distributions, we introduce the fit function$ \begin{array}{l} f=0.002\exp\bigg\{ c_{1}\exp\Big[c_{2}\cos(\Delta\theta) \Big] \bigg\}. \end{array} $
(12) $ c_1 $ and$ c_2 $ are the fit parameters.$ c_1 $ is influenced by the proportion of particles whose$ \cos(\Delta\theta)=0 $ , and$ c_2 $ is influenced by the strength of the distribution approaching$ \cos(\Delta\theta)=1 $ . The value of 0.002 is settled by us to get good fitting results, and this value is different from that in our last work because the normalized$ \cos(\Delta\theta) $ distributions have changed, as shown in Fig. 11. We can see in two$ K_{\rm {T}} $ sections,$ \pi^+ $ and all pions with 0–1.4 fm impact parameter have similar distributions, while the impact parameter has a big influence on the distributions. Thus, the impact parameter plays an important role in the single-particle space-momentum angle correlation. For we change this settled value, it leads to the changing of the fitting results in the subsequent analysis. The parameters$ c_1 $ and$ c_2 $ are changing with the$ K_{\rm {T}} $ , so we can use the fitting functions to describe their changing pattern. The fit functions areFigure 11. (color online) The normalized
$ \Delta\theta $ distributions, generated by the string melting AMPT model for the Au+Au collisions at$ \sqrt{S_{NN}}= $ 39 GeV/c.$ c_{1}=k_{1} \exp\Big[-6\times\left(\frac{K_{\rm {T}}}{1000}\right)^{2} \Big]+j_{1}, $
(13) $ c_{2}=k_{2} \exp\Big[-4.5\times\left(\frac{K_{\rm {T}}}{1000}\right)^{2} \Big]+j_{2}, $
(14) where –6 and –4.5 were chosen by us to obtain good fitting results, k and j are fit parameters.
$ k_1 $ and$ j_1 $ are related to$ c_1 $ .$ k_1 $ is influenced by the changing range of the proportion of pions with$ \cos(\Delta\theta)=0 $ , and$ j_1 $ is influenced by this lowest proportion. Besides,$ k_2 $ and$ j_2 $ are related to$ c_2 $ ,$ k_2 $ is influenced by the changing range of the strength of the distribution approaching$ \cos(\Delta\theta)=1 $ , and$ j_2 $ is influenced by its highest strength.The collision energies not only have influences on the
$ K_{{\rm T}} $ dependence of$ R_{{\rm s}} $ , but also on the$ K_{{\rm T}} $ dependence of normalized$ \cos(\Delta\theta) $ distributions. Then we can plot the parameter b, as a function of the parameters k and j, as shown in Fig. 12. The changing patterns indicate parameter b has an extremum; we set$ b_1=-0.2 $ for the situation without the Coulomb interaction, and$ b_2=-0.3 $ for considering the Coulomb interaction. We use the following functions to fit these patterns,Figure 12. (color online) The numerical connections between the strength of the transverse momentum dependence of
$ R_{\rm {s}} $ and the parameters related to the normalized$ \cos(\Delta \theta) $ distribution, the lines are fit lines.$ b(k_m) = \mu_{m1} |k_m| ^{\mu_{m2}}+b_n, $
(15) $ b(j_m) = \nu_{m1} |j_m| ^{\nu_{m2}}+b_n, $
(16) where
$ m= $ 1 or 2, for distinguishing the parameters of$ c_1 $ or$ c_2 $ . Also$ n= $ 1 or 2, then 1 is for the situation without the Coulomb interaction, and 2 refers to the Coulomb interaction. μ and ν are fit parameters, and the fitting results are shown in Table 1. The results show that the Coulomb interaction has an influence on these fit parameters, it can decrease$ |\mu_{m1}| $ and$ |\nu_{m1}| $ , and increase$ |\mu_{m2}| $ and$ |\nu_{m2}| $ . It is caused by the changing of the differences of b values. The differences in b values between the two situations are increasing with increasing collision energies. Compared to our previous work, the changing patterns are different, so these numerical connections are also related to the kinds of particles.$ n=1 $ $ n=2 $ m=1 $ \mu_{11} $ $ -0.0004\pm0.0003 $ $ -0.00004\pm0.00003 $ $ \mu_{12} $ $ 3.3\pm0.5 $ $ 4.7\pm0.4 $ $ \nu_{11} $ $ -0.086\pm0.004 $ $ -0.084\pm0.003 $ $ \nu_{12} $ $ -0.68\pm0.09 $ $ -0.94\pm0.08 $ m=2 $ \mu_{21} $ $ -0.024\pm0.003 $ $ -0.014\pm0.003 $ $ \mu_{22} $ $ 2.3\pm0.3 $ $ 3.2\pm0.3 $ $ \nu_{21} $ $ -0.021\pm0.004 $ $ -0.012\pm0.002 $ $ \nu_{22} $ $ 2.4\pm0.3 $ $ 3.3\pm0.3 $ Table 1. Fit results for parameters μ and ν.
With these fit functions and parameters, the numerical connections between the normalized
$ \cos(\Delta\theta) $ distribution and the strength of the transverse momentum dependence of$ R_{\rm {s}} $ have been built. Comparing to our previous work, we found that the fit parameters are related to the impact parameter, the kinds of pions, and the Coulomb interaction. We improve the fit parameters, which makes our numerical correlation applied to the experiment closer. With further research, the experiment data can be directly used to study the single-particle space-momentum correlation.
Single-particle space-momentum angle distribution effect on two-pion HBT correlation with Coulomb interaction
- Received Date: 2022-06-03
- Available Online: 2022-12-15
Abstract: We calculate the HBT radius