Processing math: 100%

Signals of α clusters in 16O+16O collisions at the LHC from relativistic hydrodynamic simulations

  • In relativistic heavy ion collisions, the fluctuations of initial entropy density convert to the correlations of final state hadrons in momentum space through the collective expansion of strongly interacting QCD matter. Using a (3+1)D viscous hydrodynamic program, CLVisc, we consider whether the nuclear structure, which provides initial state fluctuations as well as correlations, can affect the final state of heavy ion collisions, and whether one can find signals of α cluster structures in oxygen using final state observables in 16O+16O collisions at the CERN Large Hadron Collider. For the initial nucleon distributions in oxygen nuclei, we compare three different configurations, a tetrahedral structure with four-α clusters, the deformed Woods-Saxon distribution, and a spherical symmetric Woods-Saxon distribution. Our results show that the charged multiplicity as a function of centrality and the elliptic flow at the most central collisions using the four-α structure differs from those with the Woods-Saxon and deformed Woods-Saxon distributions, which may help to identify α clustering structures in oxygen nuclei.
  • [1] Boris A. Gelman, Edward V. Shuryak, and Ismail Zahed, Phys. Rev. C 74, 044908 (2006), arXiv:nuclth/0601029 doi: 10.1103/PhysRevC.74.044908
    [2] Ulrich W. Heinz, J. Phys. A 42, 214003 (2009), arXiv:0810.5529[nucl-th doi: 10.1088/1751-8113/42/21/214003
    [3] H. Niemi, G. S. Denicol, H. Holopainen et al., Phys. Rev. C 87, 054901 (2013), arXiv:1212.1008[nucl-th doi: 10.1103/PhysRevC.87.054901
    [4] Fernando G. Gardim, Frederique Grassi, Matthew Luzum et al., Phys. Rev. C 85, 024908 (2012), arXiv:1111.6538[nucl-th doi: 10.1103/PhysRevC.85.024908
    [5] Long-Gang Pang, Kai Zhou, and Xin-Nian Wang, Interpretable deep learning for nuclear deformation in heavy ion collisions, (2019), arXiv: 1906.06429[nucl-th]
    [6] Lu-Meng Liu, Chun-Jian Zhang, Jia Zhou et al., Phys. Lett. B 834, 137441 (2022), arXiv:2203.09924[nucl-th doi: 10.1016/j.physletb.2022.137441
    [7] Jiangyong Jia, Phys. Rev. C 105, 014905 (2022), arXiv:2106.08768[nucl-th doi: 10.1103/PhysRevC.105.014905
    [8] Hao-jie Xu, Wenbin Zhao, Hanlin Li et al., Probing nuclear structure with mean transverse momentum in relativistic isobar collisions, (2021), arXiv: 2111.14812[nucl-th]
    [9] Jiangyong Jia, Phys. Rev. C 105, 044905 (2022), arXiv:2109.00604[nucl-th doi: 10.1103/PhysRevC.105.044905
    [10] Fei Li, Yu-Gang Ma, Song Zhang et al., Phys. Rev. C 106, 014906 (2022), arXiv:2201.10994[nucl-th doi: 10.1103/PhysRevC.106.014906
    [11] Giuliano Giacalone, Phys. Rev. Lett. 124, 202301 (2020), arXiv:1910.04673[nucl-th doi: 10.1103/PhysRevLett.124.202301
    [12] S. S. Wang, Y. G. Ma, D. Q. Fang et al., Phys. Rev. C 105, 034616 (2022), arXiv:2201.01116[nucl-th doi: 10.1103/PhysRevC.105.034616
    [13] Chun-Wang Ma, Yi-Pu Liu, Huiling Wei et al., Nucl. Sci. Tech. 33, 6 (2022), arXiv:2201.01442[nucl-th doi: 10.1007/s41365-022-00997-0
    [14] Alfredo Poves, Frederic Nowacki, and Yoram Alhassid, Phys. Rev. C 101, 054307 (2020), arXiv:1906.07542[nucl-th
    [15] Yi-An Li, Dong-Fang Wang, Song Zhang et al., Phys. Rev. C 104, 044906 (2021), arXiv:2110.04784[hep-ph doi: 10.1103/PhysRevC.104.044906
    [16] Yi-An Li, Song Zhang, and Yu-Gang Ma, Phys. Rev. C 102, 054907 (2020), arXiv:2010.10003[hep-ph doi: 10.1103/PhysRevC.102.054907
    [17] Jasmine Brewer, Aleksas Mazeliauskas, and Wilke van der Schee, Opportunities of OO and pO collisions at the LHC, (2021) arXiv: 2103.01939[hep-ph]
    [18] Nicholas Summerfield, Bing-Nan Lu, Christopher Plumberg et al., Phys. Rev. C 104, L041901 (2021), arXiv:2103.03345[nucl-th doi: 10.1103/PhysRevC.104.L041901
    [19] S. Zhang, Y. G. Ma, J. H. Chen et al., Phys. Rev. C 95, 064904 (2017), arXiv:1702.02507[nucl-th doi: 10.1103/PhysRevC.95.064904
    [20] Bo-Song Huang and Yu-Gang Ma, Phys. Rev. C 103, 054318 (2021), arXiv:2109.09430[nucl-th doi: 10.1103/PhysRevC.103.054318
    [21] Kun Wang and Bing-Nan Lu, Commun. Theor. Phys. 74, 015303 (2022), arXiv:2201.05953[nucl-th doi: 10.1088/1572-9494/ac3999
    [22] Roland Katz, Caio A. G. Prado, Jacquelyn Noronha-Hostler et al., Phys. Rev. C 102, 041901 (2020), arXiv:1907.03308[nucl-th doi: 10.1103/PhysRevC.102.041901
    [23] Govert Nijs and Wilke van der Schee, Predictions and postdictions for relativistic lead and oxygen collisions with Trajectum, (2021), arXiv: 2110.13153[nucl-th]
    [24] Debadatta Behera, Neelkamal Mallick, Sushanta Tripathy et al., Eur. Phys. J. A 58, 175 (2022), arXiv:2110.04016[hep-ph doi: 10.1140/epja/s10050-022-00823-6
    [25] G. F. Bertsch and W. Bertozzi, Nucl. Phys. A 165, 199-210 (1971) doi: 10.1016/0375-9474(71)90158-8
    [26] W. von Oertzen, Martin Freer, and Yoshiko Kanada-En’yo, Physics Reports 432, 43-113 (2006) doi: 10.1016/j.physrep.2006.07.001
    [27] Enrique Ruiz Arriola and Wojciech Broniowski, J. Phys. Conf. Ser. 630, 012060 (2015), arXiv:1411.5807[nucl-th doi: 10.1088/1742-6596/630/1/012060
    [28] Maciej Rybczyński, Milena Piotrowska, and Wojciech Broniowski, Phys. Rev. C 97, 034912 (2018), arXiv:1711.00438[nucl-th doi: 10.1103/PhysRevC.97.034912
    [29] Chen-Zhong Shi and Yu-Gang Ma, Nucl. Sci. Tech. 32, 66 (2021), arXiv:2109.09938[nucl-th doi: 10.1007/s41365-021-00897-9
    [30] G. Gamow and Ernest Rutherford, Mass defect curve and nuclear constitution, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 126, 632–644 (1930)
    [31] John Archibald Wheeler, Phys. Rev. 52, 1083-1106 (1937) doi: 10.1103/PhysRev.52.1083
    [32] L. R. Hafstad and E. Teller, Phys. Rev. 54, 681-692 (1938) doi: 10.1103/PhysRev.54.681
    [33] A. Tohsaki, H. Horiuchi, P. Schuck et al., Phys. Rev. Lett. 87, 192501 (2001), arXiv:nucl-th/0110014 doi: 10.1103/PhysRevLett.87.192501
    [34] B. S. Huang, Y. G. Ma, and W. B. He, Phys. Rev. C 95, 034606 (2017), arXiv:1803.07972[nuclth doi: 10.1103/PhysRevC.95.034606
    [35] P. I. Zarubin, “tomography” of the cluster structure of light nuclei via relativistic dissociation, in Clusters in Nuclei, Volume 3, edited by Christian Beck (Springer International Publishing, Cham, 2014) pp. 51–93
    [36] C. Beck, Acta Phys. Polon. B 42, 747-756 (2011), arXiv:1011.3423[nucl-ex doi: 10.5506/APhysPolB.42.747
    [37] W. B. He, Y. G. Ma, X. G. Cao et al., Phys. Rev. Lett. 113, 032506 (2014), arXiv:1407.5414[nucl-th doi: 10.1103/PhysRevLett.113.032506
    [38] Wojciech Broniowski and Enrique Ruiz Arriola, Phys. Rev. Lett. 112, 112501 (2014), arXiv:1312.0289[nucl-th doi: 10.1103/PhysRevLett.112.112501
    [39] Wojciech Broniowski, Piotr Bożek, Maciej Rybczyński et al., Nucl. Phys. A 1005, 121763 (2021), arXiv:2001.10260[nucl-th doi: 10.1016/j.nuclphysa.2020.121763
    [40] S. Zhang, Y. G. Ma, J. H. Chen et al., Eur. Phys. J. A 54, 161 (2018), arXiv:1808.10265[nucl-th doi: 10.1140/epja/i2018-12597-y
    [41] Pingal Dasgupta, Guo-Liang Ma, Rupa Chatterjee et al., Eur. Phys. J. A 57, 134 (2021), arXiv:2007.09543[nucl-th doi: 10.1140/epja/s10050-021-00441-8
    [42] L. Ma, Y. G. Ma, and S. Zhang, Phys. Rev. C 102, 014910 (2020), arXiv:2007.10063[hep-ph doi: 10.1103/PhysRevC.102.014910
    [43] Pingal Dasgupta, Rupa Chatterjee, and Guo-Liang Ma, Production and anisotropic flow of thermal photons in most-central α-clustered C+Au collisions at BNL Relativistic Heavy Ion Collider, (2022), arXiv: 2204.00235[nucl-th]
    [44] Serdar Elhatisari, Evgeny Epelbaum, Hermann Krebs et al., Phys. Rev. Lett. 119, 222505 (2017), arXiv:1702.05177[nucl-th doi: 10.1103/PhysRevLett.119.222505
    [45] Bing-Nan Lu, Ning Li, Serdar Elhatisari et al., Phys. Lett. B 797, 134863 (2019), arXiv:1812.10928[nucl-th doi: 10.1016/j.physletb.2019.134863
    [46] Chun Shen and Li Yan, Nucl. Sci. Tech. 31, 122 (2020), arXiv:2010.12377[nucl-th doi: 10.1007/s41365-020-00829-z
    [47] Long-Gang Pang, Hannah Petersen, and Xin-Nian Wang, Phys. Rev. C 97, 064918 (2018), arXiv:1802.04449[nucl-th doi: 10.1103/PhysRevC.97.064918
    [48] P. Mller, J. R. Nix, W. D. Myers et al., Atomic Data and Nuclear Data Tables 59, 1-204 (1993)
    [49] I. Angeli and K.P. Marinova, Atomic Data and Nuclear Data Tables 99, 69-95 (2013) doi: 10.1016/j.adt.2011.12.006
    [50] J. Scott Moreland, Jonah E. Bernhard, and Steffen A. Bass, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708[nucl-th doi: 10.1103/PhysRevC.92.011901
    [51] Xiang-Yu Wu, Guang-You Qin, Long-Gang Pang et al., Phys. Rev. C 105, 034909 (2022), arXiv:2107.04949[hep-ph doi: 10.1103/PhysRevC.105.034909
    [52] Pasi Huovinen and Pter Petreczky, Nucl. Phys. A 837, 26-53 (2010), arXiv:0912.2541[hep-ph doi: 10.1016/j.nuclphysa.2010.02.015
    [53] Meng Wang, Jun-Qi Tao, Hua Zheng et al., Nucl. Sci. Tech. 33, 37 (2022), arXiv:2203.10353[hep-ph doi: 10.1007/s41365-022-01019-9
    [54] Shao-Wei Lan and Shu-Su Shi, Nucl. Sci. Tech. 33, 21 (2022) doi: 10.1007/s41365-022-01006-0
    [55] Hai Wang and Jin-Hui Chen, Nucl. Sci. Tech. 33, 15 (2022) doi: 10.1007/s41365-022-00999-y
    [56] Jaroslav Adam et al. (ALICE), Phys. Lett. B 772, 567-577 (2017), arXiv:1612.08966[nucl-ex doi: 10.1016/j.physletb.2017.07.017
    [57] J. L. Klay et al. (E-0895), Phys. Rev. C 68, 054905 (2003), arXiv:nucl-ex/0306033 doi: 10.1103/PhysRevC.68.054905
    [58] L. Ahle et al. (E-802), Phys. Rev. C 57, R466-R470 (1998) doi: 10.1103/PhysRevC.57.R466
    [59] S. V. Afanasiev et al. (NA49), Phys. Rev. C 66, 054902 (2002), arXiv:nucl-ex/0205002 doi: 10.1103/PhysRevC.66.054902
    [60] M. C. Abreu et al. (NA50), Phys. Lett. B 530, 43-55 (2002) doi: 10.1016/S0370-2693(02)01353-9
    [61] B. Alver et al. (PHOBOS), Phys. Rev. C 83, 024913 (2011), arXiv:1011.1940[nucl-ex doi: 10.1103/PhysRevC.83.024913
    [62] I. G Bearden et al. (BRAHMS), Phys. Lett. B 523, 227-233 (2001), arXiv:nucl-ex/0108016 doi: 10.1016/S0370-2693(01)01333-8
    [63] I. G. Bearden et al. (BRAHMS), Phys. Rev. Lett. 88, 202301 (2002), arXiv:nucl-ex/0112001 doi: 10.1103/PhysRevLett.88.202301
    [64] Ehab Abbas et al. (ALICE), Phys. Lett. B 726, 610-622 (2013), arXiv:1304.0347[nucl-ex doi: 10.1016/j.physletb.2013.09.022
    [65] Iu. A. Karpenko, P. Huovinen, H. Petersen et al., Phys. Rev. C 91, 064901 (2015), arXiv:1502.01978[nucl-th doi: 10.1103/PhysRevC.91.064901
    [66] ALICE physics projections for a short oxygen-beam run at the LHC, (2021)
    [67] Maciej Rybczyński and Wojciech Broniowski, Phys. Rev. C 100, 064912 (2019), arXiv:1910.09489[hep-ph doi: 10.1103/PhysRevC.100.064912
    [68] S. H. Lim, J. Carlson, C. Loizides et al., Phys. Rev. C 99, 044904 (2019), arXiv:1812.08096[nucl-th doi: 10.1103/PhysRevC.99.044904
    [69] Chi Ding, Wei-Yao Ke, Long-Gang Pang et al., Chin. Phys. C 45, 074102 (2021), arXiv:2101.02356[nucl-th doi: 10.1088/1674-1137/abf645
    [70] Huichao Song, Causal Viscous Hydrodynamics for Relativistic Heavy Ion Collisions, Other thesis (2009), arXiv: 0908.3656[nucl-th]
  • [1] Boris A. Gelman, Edward V. Shuryak, and Ismail Zahed, Phys. Rev. C 74, 044908 (2006), arXiv:nuclth/0601029 doi: 10.1103/PhysRevC.74.044908
    [2] Ulrich W. Heinz, J. Phys. A 42, 214003 (2009), arXiv:0810.5529[nucl-th doi: 10.1088/1751-8113/42/21/214003
    [3] H. Niemi, G. S. Denicol, H. Holopainen et al., Phys. Rev. C 87, 054901 (2013), arXiv:1212.1008[nucl-th doi: 10.1103/PhysRevC.87.054901
    [4] Fernando G. Gardim, Frederique Grassi, Matthew Luzum et al., Phys. Rev. C 85, 024908 (2012), arXiv:1111.6538[nucl-th doi: 10.1103/PhysRevC.85.024908
    [5] Long-Gang Pang, Kai Zhou, and Xin-Nian Wang, Interpretable deep learning for nuclear deformation in heavy ion collisions, (2019), arXiv: 1906.06429[nucl-th]
    [6] Lu-Meng Liu, Chun-Jian Zhang, Jia Zhou et al., Phys. Lett. B 834, 137441 (2022), arXiv:2203.09924[nucl-th doi: 10.1016/j.physletb.2022.137441
    [7] Jiangyong Jia, Phys. Rev. C 105, 014905 (2022), arXiv:2106.08768[nucl-th doi: 10.1103/PhysRevC.105.014905
    [8] Hao-jie Xu, Wenbin Zhao, Hanlin Li et al., Probing nuclear structure with mean transverse momentum in relativistic isobar collisions, (2021), arXiv: 2111.14812[nucl-th]
    [9] Jiangyong Jia, Phys. Rev. C 105, 044905 (2022), arXiv:2109.00604[nucl-th doi: 10.1103/PhysRevC.105.044905
    [10] Fei Li, Yu-Gang Ma, Song Zhang et al., Phys. Rev. C 106, 014906 (2022), arXiv:2201.10994[nucl-th doi: 10.1103/PhysRevC.106.014906
    [11] Giuliano Giacalone, Phys. Rev. Lett. 124, 202301 (2020), arXiv:1910.04673[nucl-th doi: 10.1103/PhysRevLett.124.202301
    [12] S. S. Wang, Y. G. Ma, D. Q. Fang et al., Phys. Rev. C 105, 034616 (2022), arXiv:2201.01116[nucl-th doi: 10.1103/PhysRevC.105.034616
    [13] Chun-Wang Ma, Yi-Pu Liu, Huiling Wei et al., Nucl. Sci. Tech. 33, 6 (2022), arXiv:2201.01442[nucl-th doi: 10.1007/s41365-022-00997-0
    [14] Alfredo Poves, Frederic Nowacki, and Yoram Alhassid, Phys. Rev. C 101, 054307 (2020), arXiv:1906.07542[nucl-th
    [15] Yi-An Li, Dong-Fang Wang, Song Zhang et al., Phys. Rev. C 104, 044906 (2021), arXiv:2110.04784[hep-ph doi: 10.1103/PhysRevC.104.044906
    [16] Yi-An Li, Song Zhang, and Yu-Gang Ma, Phys. Rev. C 102, 054907 (2020), arXiv:2010.10003[hep-ph doi: 10.1103/PhysRevC.102.054907
    [17] Jasmine Brewer, Aleksas Mazeliauskas, and Wilke van der Schee, Opportunities of OO and pO collisions at the LHC, (2021) arXiv: 2103.01939[hep-ph]
    [18] Nicholas Summerfield, Bing-Nan Lu, Christopher Plumberg et al., Phys. Rev. C 104, L041901 (2021), arXiv:2103.03345[nucl-th doi: 10.1103/PhysRevC.104.L041901
    [19] S. Zhang, Y. G. Ma, J. H. Chen et al., Phys. Rev. C 95, 064904 (2017), arXiv:1702.02507[nucl-th doi: 10.1103/PhysRevC.95.064904
    [20] Bo-Song Huang and Yu-Gang Ma, Phys. Rev. C 103, 054318 (2021), arXiv:2109.09430[nucl-th doi: 10.1103/PhysRevC.103.054318
    [21] Kun Wang and Bing-Nan Lu, Commun. Theor. Phys. 74, 015303 (2022), arXiv:2201.05953[nucl-th doi: 10.1088/1572-9494/ac3999
    [22] Roland Katz, Caio A. G. Prado, Jacquelyn Noronha-Hostler et al., Phys. Rev. C 102, 041901 (2020), arXiv:1907.03308[nucl-th doi: 10.1103/PhysRevC.102.041901
    [23] Govert Nijs and Wilke van der Schee, Predictions and postdictions for relativistic lead and oxygen collisions with Trajectum, (2021), arXiv: 2110.13153[nucl-th]
    [24] Debadatta Behera, Neelkamal Mallick, Sushanta Tripathy et al., Eur. Phys. J. A 58, 175 (2022), arXiv:2110.04016[hep-ph doi: 10.1140/epja/s10050-022-00823-6
    [25] G. F. Bertsch and W. Bertozzi, Nucl. Phys. A 165, 199-210 (1971) doi: 10.1016/0375-9474(71)90158-8
    [26] W. von Oertzen, Martin Freer, and Yoshiko Kanada-En’yo, Physics Reports 432, 43-113 (2006) doi: 10.1016/j.physrep.2006.07.001
    [27] Enrique Ruiz Arriola and Wojciech Broniowski, J. Phys. Conf. Ser. 630, 012060 (2015), arXiv:1411.5807[nucl-th doi: 10.1088/1742-6596/630/1/012060
    [28] Maciej Rybczyński, Milena Piotrowska, and Wojciech Broniowski, Phys. Rev. C 97, 034912 (2018), arXiv:1711.00438[nucl-th doi: 10.1103/PhysRevC.97.034912
    [29] Chen-Zhong Shi and Yu-Gang Ma, Nucl. Sci. Tech. 32, 66 (2021), arXiv:2109.09938[nucl-th doi: 10.1007/s41365-021-00897-9
    [30] G. Gamow and Ernest Rutherford, Mass defect curve and nuclear constitution, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 126, 632–644 (1930)
    [31] John Archibald Wheeler, Phys. Rev. 52, 1083-1106 (1937) doi: 10.1103/PhysRev.52.1083
    [32] L. R. Hafstad and E. Teller, Phys. Rev. 54, 681-692 (1938) doi: 10.1103/PhysRev.54.681
    [33] A. Tohsaki, H. Horiuchi, P. Schuck et al., Phys. Rev. Lett. 87, 192501 (2001), arXiv:nucl-th/0110014 doi: 10.1103/PhysRevLett.87.192501
    [34] B. S. Huang, Y. G. Ma, and W. B. He, Phys. Rev. C 95, 034606 (2017), arXiv:1803.07972[nuclth doi: 10.1103/PhysRevC.95.034606
    [35] P. I. Zarubin, “tomography” of the cluster structure of light nuclei via relativistic dissociation, in Clusters in Nuclei, Volume 3, edited by Christian Beck (Springer International Publishing, Cham, 2014) pp. 51–93
    [36] C. Beck, Acta Phys. Polon. B 42, 747-756 (2011), arXiv:1011.3423[nucl-ex doi: 10.5506/APhysPolB.42.747
    [37] W. B. He, Y. G. Ma, X. G. Cao et al., Phys. Rev. Lett. 113, 032506 (2014), arXiv:1407.5414[nucl-th doi: 10.1103/PhysRevLett.113.032506
    [38] Wojciech Broniowski and Enrique Ruiz Arriola, Phys. Rev. Lett. 112, 112501 (2014), arXiv:1312.0289[nucl-th doi: 10.1103/PhysRevLett.112.112501
    [39] Wojciech Broniowski, Piotr Bożek, Maciej Rybczyński et al., Nucl. Phys. A 1005, 121763 (2021), arXiv:2001.10260[nucl-th doi: 10.1016/j.nuclphysa.2020.121763
    [40] S. Zhang, Y. G. Ma, J. H. Chen et al., Eur. Phys. J. A 54, 161 (2018), arXiv:1808.10265[nucl-th doi: 10.1140/epja/i2018-12597-y
    [41] Pingal Dasgupta, Guo-Liang Ma, Rupa Chatterjee et al., Eur. Phys. J. A 57, 134 (2021), arXiv:2007.09543[nucl-th doi: 10.1140/epja/s10050-021-00441-8
    [42] L. Ma, Y. G. Ma, and S. Zhang, Phys. Rev. C 102, 014910 (2020), arXiv:2007.10063[hep-ph doi: 10.1103/PhysRevC.102.014910
    [43] Pingal Dasgupta, Rupa Chatterjee, and Guo-Liang Ma, Production and anisotropic flow of thermal photons in most-central α-clustered C+Au collisions at BNL Relativistic Heavy Ion Collider, (2022), arXiv: 2204.00235[nucl-th]
    [44] Serdar Elhatisari, Evgeny Epelbaum, Hermann Krebs et al., Phys. Rev. Lett. 119, 222505 (2017), arXiv:1702.05177[nucl-th doi: 10.1103/PhysRevLett.119.222505
    [45] Bing-Nan Lu, Ning Li, Serdar Elhatisari et al., Phys. Lett. B 797, 134863 (2019), arXiv:1812.10928[nucl-th doi: 10.1016/j.physletb.2019.134863
    [46] Chun Shen and Li Yan, Nucl. Sci. Tech. 31, 122 (2020), arXiv:2010.12377[nucl-th doi: 10.1007/s41365-020-00829-z
    [47] Long-Gang Pang, Hannah Petersen, and Xin-Nian Wang, Phys. Rev. C 97, 064918 (2018), arXiv:1802.04449[nucl-th doi: 10.1103/PhysRevC.97.064918
    [48] P. Mller, J. R. Nix, W. D. Myers et al., Atomic Data and Nuclear Data Tables 59, 1-204 (1993)
    [49] I. Angeli and K.P. Marinova, Atomic Data and Nuclear Data Tables 99, 69-95 (2013) doi: 10.1016/j.adt.2011.12.006
    [50] J. Scott Moreland, Jonah E. Bernhard, and Steffen A. Bass, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708[nucl-th doi: 10.1103/PhysRevC.92.011901
    [51] Xiang-Yu Wu, Guang-You Qin, Long-Gang Pang et al., Phys. Rev. C 105, 034909 (2022), arXiv:2107.04949[hep-ph doi: 10.1103/PhysRevC.105.034909
    [52] Pasi Huovinen and Pter Petreczky, Nucl. Phys. A 837, 26-53 (2010), arXiv:0912.2541[hep-ph doi: 10.1016/j.nuclphysa.2010.02.015
    [53] Meng Wang, Jun-Qi Tao, Hua Zheng et al., Nucl. Sci. Tech. 33, 37 (2022), arXiv:2203.10353[hep-ph doi: 10.1007/s41365-022-01019-9
    [54] Shao-Wei Lan and Shu-Su Shi, Nucl. Sci. Tech. 33, 21 (2022) doi: 10.1007/s41365-022-01006-0
    [55] Hai Wang and Jin-Hui Chen, Nucl. Sci. Tech. 33, 15 (2022) doi: 10.1007/s41365-022-00999-y
    [56] Jaroslav Adam et al. (ALICE), Phys. Lett. B 772, 567-577 (2017), arXiv:1612.08966[nucl-ex doi: 10.1016/j.physletb.2017.07.017
    [57] J. L. Klay et al. (E-0895), Phys. Rev. C 68, 054905 (2003), arXiv:nucl-ex/0306033 doi: 10.1103/PhysRevC.68.054905
    [58] L. Ahle et al. (E-802), Phys. Rev. C 57, R466-R470 (1998) doi: 10.1103/PhysRevC.57.R466
    [59] S. V. Afanasiev et al. (NA49), Phys. Rev. C 66, 054902 (2002), arXiv:nucl-ex/0205002 doi: 10.1103/PhysRevC.66.054902
    [60] M. C. Abreu et al. (NA50), Phys. Lett. B 530, 43-55 (2002) doi: 10.1016/S0370-2693(02)01353-9
    [61] B. Alver et al. (PHOBOS), Phys. Rev. C 83, 024913 (2011), arXiv:1011.1940[nucl-ex doi: 10.1103/PhysRevC.83.024913
    [62] I. G Bearden et al. (BRAHMS), Phys. Lett. B 523, 227-233 (2001), arXiv:nucl-ex/0108016 doi: 10.1016/S0370-2693(01)01333-8
    [63] I. G. Bearden et al. (BRAHMS), Phys. Rev. Lett. 88, 202301 (2002), arXiv:nucl-ex/0112001 doi: 10.1103/PhysRevLett.88.202301
    [64] Ehab Abbas et al. (ALICE), Phys. Lett. B 726, 610-622 (2013), arXiv:1304.0347[nucl-ex doi: 10.1016/j.physletb.2013.09.022
    [65] Iu. A. Karpenko, P. Huovinen, H. Petersen et al., Phys. Rev. C 91, 064901 (2015), arXiv:1502.01978[nucl-th doi: 10.1103/PhysRevC.91.064901
    [66] ALICE physics projections for a short oxygen-beam run at the LHC, (2021)
    [67] Maciej Rybczyński and Wojciech Broniowski, Phys. Rev. C 100, 064912 (2019), arXiv:1910.09489[hep-ph doi: 10.1103/PhysRevC.100.064912
    [68] S. H. Lim, J. Carlson, C. Loizides et al., Phys. Rev. C 99, 044904 (2019), arXiv:1812.08096[nucl-th doi: 10.1103/PhysRevC.99.044904
    [69] Chi Ding, Wei-Yao Ke, Long-Gang Pang et al., Chin. Phys. C 45, 074102 (2021), arXiv:2101.02356[nucl-th doi: 10.1088/1674-1137/abf645
    [70] Huichao Song, Causal Viscous Hydrodynamics for Relativistic Heavy Ion Collisions, Other thesis (2009), arXiv: 0908.3656[nucl-th]
  • 加载中

Figures(6) / Tables(1)

Get Citation
Chi Ding, Long-Gang Pang, Song Zhang and Yu-Gang Ma. Signals of α clusters in 16O+16O collisions at LHC from relativistic hydrodynamic simulations[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac9fb8
Chi Ding, Long-Gang Pang, Song Zhang and Yu-Gang Ma. Signals of α clusters in 16O+16O collisions at LHC from relativistic hydrodynamic simulations[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac9fb8 shu
Milestone
Received: 2022-08-17
Article Metric

Article Views(1351)
PDF Downloads(57)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Signals of α clusters in 16O+16O collisions at the LHC from relativistic hydrodynamic simulations

  • 1. Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China
  • 2. Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Institute of Modern Physics, Fudan University, Shanghai 200433, China
  • 3. Shanghai Research Center for Theoretical Nuclear Physics NSFC and Fudan University, Shanghai 200438, China

Abstract: In relativistic heavy ion collisions, the fluctuations of initial entropy density convert to the correlations of final state hadrons in momentum space through the collective expansion of strongly interacting QCD matter. Using a (3+1)D viscous hydrodynamic program, CLVisc, we consider whether the nuclear structure, which provides initial state fluctuations as well as correlations, can affect the final state of heavy ion collisions, and whether one can find signals of α cluster structures in oxygen using final state observables in 16O+16O collisions at the CERN Large Hadron Collider. For the initial nucleon distributions in oxygen nuclei, we compare three different configurations, a tetrahedral structure with four-α clusters, the deformed Woods-Saxon distribution, and a spherical symmetric Woods-Saxon distribution. Our results show that the charged multiplicity as a function of centrality and the elliptic flow at the most central collisions using the four-α structure differs from those with the Woods-Saxon and deformed Woods-Saxon distributions, which may help to identify α clustering structures in oxygen nuclei.

    HTML

    I.   INTRODUCTION
    • The Relativistic Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory (BNL) and the Large Hadron Collider (LHC) at CERN have been used to create a near perfect fluid known as quark-gluon plasma (QGP) [1, 2]. QGP expands rapidly and collectively such that the collision geometry and initial state fluctuations are transformed into long-range multi-particle correlations at the final state [3, 4]. For example, the momentum anisotropy of final state hadrons is sensitive to spatial eccentricity at the initial state. A question to consider is whether the final state observables are also affected by the initial state nuclear structure, such as the nuclear size, nuclear shape deformation, neutron skin, and α clusters inside nucleus. Furthermore, will the nuclear structure signal survive the violent collision, especially since the colliding nuclei are highly relativistic and the initial state might be dominated by color glass condensates?

      Relativistic heavy ion collisions, especially isobar collisions, have been employed to study nuclear shape deformation [510], neutron skin thickness [1113], the nucleon-nucleon correlation [14], and the α clustering structure inside the nucleus [12, 1521]. Both the RHIC and LHC are planning to collect data for 16O+16O collisions, which is a medium-sized system [17, 22] that can serve as a control group for the peripheral collisions of large nuclei and high multiplicity p+p and p+Pb events in small systems. There have been many theoretical studies on 16O+16O collisions [23, 24].

      16O+16O collisions provide a unique opportunity to study the α clustering structure in the oxygen nucleus because 16O at the ground state is suggested to be a tetrahedral structure composed of four α clusters [2529]. The concept of the α cluster was proposed by Gamow in the 1930s [30], who noticed that light nuclei are more stable if their nucleons can form several αs. Increasing evidence indicates that the mean field cannot break the α cluster structure in light nuclei, such as 12C and 16O [3134]. The experimental evidence for clustering originates from fragmentation studies [35]. Nuclear lattice effective field theory (NLEFT) also suggests that nucleons in oxygen cluster together like α particles.

      α clusters inside a nucleus bring a significantly different nuclear shape deformation from the Woods-Saxon (WS) distribution [36, 37]. Wojciech proposed the study of α clusters in light nuclei using relativistic heavy ion collisions [38, 39]. Others have proposed the study of collision systems between light nuclei (12C) with α clustering structures and heavy nuclei without them, such as Pb and Au [4043]. Many studies have focused on observables such as collective flows with particle multiplicity and the eccentricity coefficient as functions of the number of wounded nucleons.

      Relativistic hydrodynamic simulations of 16O+16O collisions have shown signals of α-clustering in the oxygen nucleus [18]. In this study, the Trento Monte Carlo model is used to generate the initial conditions of 16O+16O collisions for iEBE-VISHNU hydrodynamic simulations. The nucleon distributions in 16O are provided by ab initio calculations using NLEFT [44, 45]. This ab initio calculation suggests the existence of α-clusters in the oxygen nucleus. It has been observed that α-clustering suppresses v3{2}/v2{2} and enhances v4{2}/v2{2}, whereas sub-nucleon fluctuations always have the opposite effect.

      Moreover, using a multi phase transport model (AMPT) to simulate 16O+16O collisions, it was discovered that the forward-backward (FB) multiplicity correlation is smaller for the four-α tetrahedron structure than the WS distribution. The FB correlation has thus been proposed as a probe for α-clustering inside the oxygen nucleus [15].

      In this study, we attempt to search for more signals of α-clusters in the oxygen nucleus using relativistic hydrodynamic simulations of 16O+16O collisions [46]. We calculate the multiplicity, transverse momentum spectra, and anisotropy of charged particles using relativistic hydrodynamic simulations with different nucleon distributions of 16O. We use the Trento Monte Carlo model to generate the initial entropy density distributions in the transverse plane for 16O+16O collisions with the WS distribution, deformed WS distribution, and four-α tetrahedral structure for the oxygen nucleus. The cross-section of nucleon-nucleon collisions at sNN=7 TeV (LHC) is set to σNN=7.1 fm2. The initial entropy density profile from Trento is fed into CLVisc model (a 3+1D hydrodynamic model) to simulate the expansion of QGP [47] and its transformation into hadrons. The remainder of this paper is organized as follows. In Section II, we introduce three different distributions of nucleons inside the oxygen nucleus and the relativistic hydrodynamic model CLVisc. Section III compares the charged particle multiplicities and their pT spectra and pT-differential vn using three different oxygen structures. Signals of α-clustering in the oxygen nucleus are collectively presented in the summary.

    II.   MODEL SETUP

      A.   Woods-Saxon and deformed Woods-Saxondistributions

    • The WS distribution is a good approximation to describe the coordinates of nucleons inside a double magic nucleus, e.g., 208Pb, whose shape is approximately a perfect sphere, i.e.,

      ρ(r)=ρ01+exp[rRa],

      (1)

      where r is the radial distance to the center of the nucleus, ρ0 is the nuclear density determined using A=2πρ(r,θ)r2sinθdrdθ, a is the surface diffuseness parameter, and R is the effective radius of the nucleus.

      For most nuclei whose number of protons or neutrons is not a magic number, the nuclear shape is deformed. The deformed WS distribution is usually employed to describe the nucleons inside non-spherical nuclei,

      ρ(r,θ)=ρ01+exp[rR(1+β2Y20(θ)+β4Y40(θ))a],

      (2)

      where θ is the polar angle with respect to the symmetry axis of the nucleus, and Y20(θ)=145π(3cos2θ1) and Y40(θ)=316π(35cos4θ30cos2θ+3) are spherical harmonic functions. The quadrupole deformation parameter β2 and hexadecapole parameter β4 control the magnitude of deformation from the spherical shape. As shown in the left window of Fig. 1, the shape change is caused by deformation on the transverse plane, and unlike the spherical Pb nucleus, the O nucleus is more like a square.

      Figure 1.  (color online) Left: Shape of Pb and oxygen. The Woods-Saxon and deformed Woods-Saxon distributions have been used for oxygen. Right: Four-α tetrahedral structure of oxygen, with four nucleons in each α sampled from the Woods-Saxon distribution.

      In this study, we use the WS and deformed WS distributions to sample the coordinates of 16 nucleons inside each oxygen nucleus. For the deformed WS distribution, the values of the parameters for oxygen are A=16, R=3.18 fm, a=0.662 fm, β2=0.021, and β4=0.108 [48]. For the WS distribution, we set β2=β4=0 and fix other parameters as in the deformed WS distribution. The resulting nuclear shapes from these two groups of parameters are shown on the left of Fig. 1, along with the 208Pb nucleus as a reference. For deformed nuclei, different orientations correspond to different collision geometries. To generate realistic collisions of deformed nuclei, we rotate the entire oxygen nucleus randomly after all nucleons inside the oxygen are sampled from the deformed WS distribution.

    • B.   α cluster and tetrahedral frame

    • The four-α tetrahedron structure of 16O is provided by the extended quantum molecular dynamics (EQMD) model [16, 37] with the effective Pauli potential. In this configuration, four αs are located at the vertices of a regular tetrahedron with a side length of 3.42 fm, as shown on the right of Fig. 1. The side length is fixed to produce a similar RMS-radius (2.699 fm) [16] as the experimental data (2.6991 fm) [49]. For each α cluster, the coordinates of two protons and two neutrons are sampled from a three-parameter distribution function,

      ρ(r)=ρ0(1+wr2R2)[1+exp(rRa)]1,

      (3)

      where ρ0 is the nuclear density at the center of the alpha, w=0.517, R=0.964 fm, and a=0.322 fm. Note that we apply three successive Eulerian rotations to each sampled oxygen nucleus before the collision, and the minimum distance between two nucleons inside the same oxygen is larger than dr=0.5 fm to avoid overlapping. This configuration has also been used in AMPT [16] simulations to search for signatures of α clustering in 16O. Using the same initial state nuclear configuration in relativistic hydrodynamic simulations helps to further investigate the effect of dynamical evolution on the signature of α clustering.

      We verify that the simulated radial distribution of charge density agrees well with the empirical result [45], as shown in Fig. 2.

      Figure 2.  (color online) Charge density distributions of 16O from empirical result (red line) and samples from the Trento model using a regular tetrahedron composed of four α clusters (black line with dots).

    • C.   Relativistic hydrodynamics

    • The initial entropy density in the transverse plane is locally proportional to the reduced thickness function TR,

      dSdy|τ=τ0TR(p;TA,TB),

      (4)

      where

      TR(p;TA,TB)={max(TA,TB),p+,TA+TB,p=+1,TATB,p=0,2TATB/(TA+TB),p=1,min(TA,TB),p.

      (5)

      where TA and TB indicate the participant nucleon thickness functions of the two colliding nuclei. In this study, we employ the default option, p=0, because of its excellent agreement with the IP-Glasma results and experimental data. Using this configuration, relativistic hydrodynamic simulations can describe not only the centrality dependence of charged multiplicity, but also the v3 to v2 ratio in Pb+Pb collisions [50].

      For the longitudinal direction, we use an envelope function that extends along the space-time rapidity,

      H(ηs)=exp[(|ηs|ηw)22σ2ηθ(|ηs|ηw)],

      (6)

      where ση=2.0 and ηw=1.7 are used for 16O+16O collisions at sNN=7 TeV.

      The CLVisc model is a (3+1)D viscous hydrodynamic model that is parallelized on GPU using OpenCL. This model simulates the hydrodynamic evolution of hot and dense QCD matter (QGP) and hadron resonance gas (HRG) before kinetic freeze-out by solving the following equations together with the Israel-Stewart equations for πμν:

      μTμν=0, with Tμν=(e+P)uμuνPgμν+πμν,

      (7)

      where in the energy-momentum tensor Tμν, e is the energy density, P is the pressure as a function of energy density given by the equation of state (EoS), uμ is the fluid four-velocity obeying uμuμ=1, and πμν is the shear stress tensor [51].

      We assume that hydrodynamic evolution starts at τ0=0.6 fm. In this study, we use the partial chemical equilibrium EoS with a chemical freeze-out temperature of 165 MeV and a smooth crossover between QGP at high temperatures and the HRG EoS at low temperatures [52], as inspired by the lattice QCD study in [52] (s95p-pce).

      The momentum distribution of hadron freeze-out from the hyper-surface follows the Cooper-Frye formula,

      dNidypTdpTdϕ=gi(2π)3pμdΣμf(pu)(1+δf),

      (8)

      where gi=2si+1 is the spin degeneracy of particle species i with spin si, pμ is the four-momenta of these particles in the lab frame, Σμ is their freeze-out hyper-surface, f(pu) is their Fermi-Dirac/Bose-Einstein distribution function,

      f(pu)=1exp[(puμi)/Tfrz]±1,

      (9)

      δf is responsible for the non-equilibrium correction and is given by considering the contribution of the shear stress tensor πμν,

      δf=(1feq)pμpνπμν2T2frz(ε+P).

      (10)

      We choose the freeze-out temperature Tfrz=137 MeV for light hadrons.

      The anisotropic collective flow vn is calculated from the Fourier expansion of the particle momentum spectrum in the transverse plane of the collisions [5355],

      d3NpTdpTdydϕ=d2N2πpTdpTdy[1+n=12vncos(n(ϕΨEP))],

      (11)

      where ΨEP denotes the azimuth angle of the event plane.

      In this study, we simulate event-by-event 16O+16O collisions at sNN=7 TeV for three centralities, with 1000 events in each centrality.

    III.   RESULTS
    • To calibrate the scale factor of the initial entropy density from the Trento model, we use the empirical formula

      Total  Nch/<Npart/2>=Asplogs+B

      (12)

      for the 0–5% most central collisions [56]. This expresses the total number of charged hadrons, Total Nch, produced in the most central heavy-ion collisions as a function of the squared collision energy s([GeV]2), which is the square of the collision energy of a pair of nucleons in their center of mass frame. The data for the most central collisions at the AGS (0–5% Au-Au) [57, 58], SPS (0–5% Pb-Pb) [59, 60], RHIC (0–5% and 0–6% Au-Au) [6163], and LHC (0–5% Pb-Pb) [64] are used for calibration. The parameters A=0.632, B=0.244, and p= 0.137 are used, and the value of the scale factor for 16O+16O collisions at sNN=7 TeV is 205.

      Shown in the left window of Fig. 3 is the charged multiplicity as a function of pseudo-rapidity from CLVisc relativistic hydrodynamic simulations of 16O+16O collisions at sNN=7 TeV, with nucleons in oxygen sampled from the WS distribution (blue-dashed line), deformed WS distribution (red-solid line), and four-α tetrahedron distribution (black-solid line with dots). The multiplicity as a function of pseudo-rapidity is shown to overlap for the WS and deformed WS distributions in 05%, 20%30%, and 40%60% collisions. The small deformation parameters β2 and β4 shown in the initial condition do not seem to make any difference in the distribution of charged particle multiplicity as a function of centrality. This is also confirmed by the probability distribution of initial total entropy shown in the right window of Fig. 3, which is calculated from 1 million events.

      Figure 3.  (color online) Left: Charged particle multiplicity as a function of pseudo-rapidity from CLVisc relativistic hydrodynamic simulations of 16O+16O collisions at sNN=7 TeV, with nucleons in oxygen sampled from the Woods-Saxon distribution (blue-dashed line), deformed Woods-Saxon distribution (red-solid line), and four-α tetrahedron distribution (black-solid line with dots). Right: Probability density distribution of total entropy calculated from Trento for these three different initial configurations of the oxygen nucleus.

      Using the same number of nucleons in oxygen and the same initial scaling factor, Trento produces more entropy using the four-α tetrahedron distribution than using WS and deformed WS for oxygen in 16O+16O collisions. It was studied in [65] that increasing the Gaussian smearing width of hot spots in the initial condition of heavy ion collisions leads to a higher total entropy at the initial state and more particles at the final state. The present study shows that different distributions of nucleons inside nuclei also bring different values of the total entropy and the final state charged multiplicity of the colliding system. As shown in the left window of Fig. 3, the charged multiplicity in the most central collisions using the four-α structure is 10% higher than that using the WS and deformed WS distributions.

      Note that the enhancement of charged multiplicity is not sufficient to tell whether α clusters exist in oxygen. It was verified that the oxygen structure calculated from the harmonic oscillator potential also increases the charged multiplicity in the most central 16O+16Ocollisions by 15% compared with the WS distribution [24]. However, using the harmonic oscillator potential to calculate the nucleon distribution inside the oxygen nucleus gives a different nucleon density profile. The enhancement of charged multiplicity can be used together with other observations to search for α clusters inside the oxygen nucleus.

      Figure 3 shows that the charged particle multiplicity as a function of centrality is sensitive to the nucleon distribution inside the oxygen nucleus. The ratios of charged particle multiplicity between different centralities are distinct for different nuclear structures, as shown in Table 1. The ratio of charged particle multiplicity between 05% and 20%30% for the four-α equilateral structure is approximately 2.29, which is smaller than that for the WS and deformed WS distributions by 6% percent. The ratio of charged multiplicity between 05% and 40%60% for the four-α equilateral structure is 5.38, which is approximately 9% smaller than that for the WS and deformed WS distributions. The experimental measurements of these two ratios can easily reveal whether the α-cluster structure plays a role in final state observables in heavy ion collisions. Usually, the total multiplicity in 05% collisions is used to determine the overall scaling factor. Owing to the lack of experimental data, the current calculation is merely a prediction for 16O+16O collisions at the LHC. There might be uncertainties in the scaling factor, but the ratio between different centrality bins seems to be a sensitive probe for the nuclear structure.

      Cent1/Cent2WSDeformed WSFour-α
      (0–5%)/(20%–30%)2.442.442.29
      (0–5%)/(40%–60%)5.995.915.38

      Table 1.  Ratio of charged multiplicity at middle pseudo-rapidity between different centralities.

      Figure 4 shows the pT spectrum of charged hadrons given by event-by-event hydrodynamic simulations of 16O+16O collisions for three different initial conditions, with nucleons in oxygen sampled from the WS distribution (black-dashed line), deformed WS distribution (red-solid line), and four-α tetrahedron distribution (black-solid line with dots). No significant differences are observed between the pT spectra with the WS and deformed WS distributions. This is consistent with the small differences between the shapes of oxygen and deformed oxygen, as demonstrated in Fig. 2. The quadrupole deformation of deformed oxygen can be neglected because the β2 parameter is too small to make a difference. Using the four-α tetrahedral structure, the pT spectrum is slightly enhanced at pT>1 GeV. This is also in agreement with the higher charged multiplicity using oxygen with α-clusters than WS and deformed WS distributions.

      Figure 4.  (color online) Transverse momentum spectra of charged hadrons given by event-by-event hydrodynamic simulations of 16O+16O collisions at sNN=7 TeV, with nucleons in oxygen sampled from the Woods-Saxon distribution (black-dashed line), deformed Woods-Saxon distribution (red-solid line), and four-α equilateral distribution (black-solid line with dots).

      The centrality dependence of the pT differential vn shown in Fig. 5 indicates that another difference may be used to look for α clusters inside the oxygen nucleus. Comparisons among the WS, deformed WS, and four-α tetrahedron distributions show that the α cluster structure decreases the pT differential v2 in 05% and 20%30% collisions and the pT differential v3 in 20%30% collisions, but slightly increases v4 in 20%30% collisions compared with the other two initial nuclear structures without α clusters. A previous study shows that the eccentricity in 16O+16O collisions is only weakly dependent on the α clustering structure [66]. The ordering of the pT differential vn seems to be a better signal than the absolute values of vn at different centralities. v3(pT) in 16O+16O collisions deceases as the centrality increases from central (0–5%) to peripheral (40%–60%) collisions. This feature is universal for the three different oxygen nuclear configurations and thus cannot be used to detect the α cluster in the oxygen nucleus.

      Figure 5.  (color online) Anisotropic flow of charged hadrons in the CLVisc simulation of 16O+16O collisions at sNN=7 TeV with and without the α-cluster structure using the event-plane method.

      It was observed in Ref. [67] that two-nucleon correlations in oxygen have only a moderate effect, and no significant differences should be expected in a comparison between 16O+16O and Pb+Pb collisions. In another study [22], a system size scan was conducted for the D meson v2{2}(pT) using initial conditions + hydrodynamics + the heavy flavor Langevin model. It is argued that the D meson v2{2}(pT) is independent of system size because of the competing effects between a larger geometric eccentricity, which will increase v2, and a smaller system size, which will suppress v2, for 16O+16O collisions compared with the case of Pb+Pb collisions. A realistic hydrodynamic simulation in the present study shows a considerably larger pT differential v2 for 16O+16O than Pb+Pb collisions in the 05% centrality because of the dominant large initial eccentricity. The effect of system size can be observed from the centrality dependence of v2(pT), which decreases with decreasing centrality in 16O+16O collisions but increases in Pb+Pb collisions. Our results are supported by the study in [68], which also observes a large v2(pT) at 05% 16O+16O collisions using the AMPT. If this is verified by experimental data in the near future, it may also suggest a considerably larger D meson v2{2}(pT) in 16O+16O collisions than Pb+Pb collisions if low pT D mesons approach local equilibrium [69].

      Because of larger event-by-event fluctuations in smaller collision systems, as shown in Fig. 6, the initial state eccentricity in 05% 16O+16O collisions is considerably larger than that in Pb+Pb collisions. In previous studies [70], the ratio vn/εn is typically compared between large and small collision systems, where εn is the n-th order geometric eccentricity, to eliminate the effect of initial state fluctuations as well as the differences between nuclear structures. This makes sense for the study of the properties of QGP. However, to study the effect of nuclear structure and initial state fluctuations, it is better to compare the pT differential vn between large and small systems directly. Indeed, vn has a strong dependence on the size of the colliding nuclei, as shown in Fig. 5.

      Figure 6.  (color online) Initial entropy density distributions of 16O+16O collisions and Pb+Pb collisions at a centrality of 0–5%. The entropy density distribution in the transverse plane of a single colliding event for (a) O-O collisions with the Woods-Saxon distribution, (b) O-O collisions with the deformed Woods-Saxon distribution, (c) O-O collisions with the four-α tetrahedral distribution, and (d) Pb-Pb collisions. Sub-figures (e), (f), (g), and (h) are the corresponding averaging entropy density distributions with each event rotated to their participant plane.

    IV.   SUMMARY
    • We study 16O+16Ocollisions using relativistic hydrodynamics with initial conditions given by the Trento Monte Carlo model for three different nucleon distributions inside the oxygen nucleus, i.e., the WS distribution, deformed WS distribution, and four-α tetrahedral structure. We observe that the four-α structure in the oxygen nucleus generates more entropy at the initial state of collisions than the WS and deformed WS distributions. As a result, more charged hadrons are produced at the 05% most central collisions using the four-α structure than with the other two configurations. The ratios between 05% and 20%30% and between 05% and 40%60% using the four-α structure are smaller than those for the WS and deformed WS distributions by 6%9%. The pT differential elliptic flow v2(pT) in 05% collisions using the four-α structure is also smaller than that for the WS and deformed WS distributions; however, v4(pT) in 20%30% collisions is larger than that for the WS and deformed WS distributions. These observations thus provide multiple signals of the α-cluster structure inside the oxygen nucleus.

    ACKNOWLEDGMENTS
    • LG Pang and C Ding also acknowledge the support provided by Huawei Technologies Co., Ltd.

Reference (70)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return