Influence of α-clustering nuclear structure on the rotating collision system

Influence of α -clustering nuclear structure on the rotating collision system 1 2 first-author 3 4 corresp zhangsong@sinap.ac.cn zhangsong@sinap.ac.cn zhangsong@sinap.ac.cn 3 4 corresp ygma@sinap.ac.cn ygma@sinap.ac.cn ygma@sinap.ac.cn 3 4 3 4 Department of Physics, Fudan University, Shanghai 200433, China Department of Physics, Fudan University, Shanghai 200433, China Department of Physics and Astronomy, University of California, Los Angeles, CA 90095, USA Department of Physics and Astronomy, University of California, Los Angeles, CA 90095, USA Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Institute of Modern Physics, Fudan University, Shanghai 200433, China Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Institute of Modern Physics, Fudan University, Shanghai 200433, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China In recent years, the collective motion properties of global rotation of the symmetric colliding system in relativistic energies have been investigated. In addition, the initial geometrical shape effects on the collective flows have been explored using a hydrodynamical model, a transport model etc. In this work, we study the asymmetric 12 C + 197 Au collision at 200 GeV/ c and the effect of the exotic nuclear structure on the global rotation using a multi-phase transport model. The global angular momentum and averaged angular speed were calculated and discussed for the collision system at different evolution stages. Chiral magnetic effect Chiral vortical effect Initial geometrical effect Quark-gluon plasma Relativistic heavy-ion collisions 1 Introduction 1 In relativistic heavy-ion collisions, properties of the collision system are often investigated using the identified particle spectra, collective motion in transverse momentum plane (collective flow) [ 1 1 - 4 4 1-4 ], jet-quenching [ 5 5 , 6 6 ], Hanbury–Brown–Twiss (HBT) correlation [ 7 7 , 8 8 ], etc. Many experimental results suggest that hot-dense matter at the partonic level with large collective motion is created in the early stages of relativistic heavy-ion collisions; this is a new type of strong-coupling quark–gluon plasma (QGP) matter with almost the nearly lowest shear viscosity. On the other hand, global rotation of the collision system is important to understand the global and local polarization of the quark matter. By using a multi-phase transport (AMPT) model and hard-sphere model, evolution of angular momentum and vorticity fields have been calculated for Au + Au collisions [ 9 9 ]. Some theoretical works have predicted globally polarized QGP in relativistic heavy-ion collisions [ 10 10 , 11 11 ], and it was suggested that system vorticity and particle production mechanism should be studied by measuring global spin alignment of vector mesons in relativistic heavy-ion collisions [ 12 12 , 13 13 ]. Recently, RHIC-STAR collaboration [ 14 14 ] reported the first measurements for global Λ hyperon polarization in Au + Au collisions, indicating that the vorticity of the QGP might reach ω ~ (9±1) × 10 21 /s, which far surpasses the vorticity of all other known fluids. The initial geometry distribution could influence initial dynamical fluctuation and intrinsic structure in the collided system, which affects the collective flow, HBT correlation, etc. Hydrodynamical models [ 15 15 - 18 18 15-18 ] as well as transport models [ 19 19 - 23 23 19-23 ] were employed to investigate these effects. Some observables or physical quantities that are sensitive to the initial geometry distribution were proposed, such as collective flows [ 24 24 - 26 26 24-26 ], fluctuation of conserved quantities [ 27 27 , 28 28 ], density fluctuations [ 29 29 ], and charge separation [ 30 30 ]. In References [ 19 19 , 20 20 ] the carbon was considered with a 3- α structure and collided against a heavy nucleus at very high energies; the results implied that the final collective flow was sensitive to the intrinsic geometry distribution of carbon. The ratio of triangular flow to elliptic flow can be used to detect the intrinsic structure of α -clustering nuclei within an α -clustered 12 C colliding against heavy-ion by using AMPT model [ 31 31 ]. The α -cluster model, which was originally proposed by Gamow [ 32 32 ], considered some light nuclei made of N- α , such as 12 C with 3- α and 16 O with 4- α . It was suggested that α clustering configurations can be identified by giant dipole resonance [ 33 33 , 34 34 ] or photo-nuclear reaction in the quasi-deuteron region [ 35 35 , 36 36 ] by an extended quantum molecular dynamics (EQMD) model. In this work, the angular momentum, inertia, and average angular velocity for an asymmetric collision system, namely 12 C + 197 Au, was investigated for comparing the three different clustering configurations of 12 C. In the second section, an introduction to AMPT model and the calculation methods for angular momentum and other physical quantities in the many-body system is presented. In the third section, we first discuss the main calculation results of the angular momentum, inertia, and velocity for the three clustering configurations. Then, we explore the density distribution of the system in detail. The paper is concluded in the fourth section. 2 model and calculation methods 1 A step-by-step simulation of heavy-ion collisions in phase-space was obtained by AMPT model [ 37 37 ]. AMPT is a successful model used to study colliding energy ranges of relativistic heavy-ion collisions at RHIC [ 37 37 ] and LHC [ 38 38 ]. The model can be used to study pion-HBT correlations [ 39 39 ], collective flow [ 40 40 , 41 41 ], di-hadron azimuthal correlations [ 42 42 ], strangeness production [ 43 43 , 28 28 ], and chiral magnetic effects [ 44 44 , 45 45 ], etc. AMPT model was developed to simulate the collision system. It is applied to different transport theories to describe the many-body interaction from the initial to final states. The model is initialized by employing the HIJING model [ 46 46 , 47 47 ] that generated hard mini-jet partons and soft strings to form initial states. It was followed by a parton cascade model (ZPC) [ 48 48 ] that simulates the interaction of the melted partons. Then, a quark coalescence model was adopted to describe the formation of hadrons; these hadrons participate in hadronic rescattering using a relativistic transport (ART) model [ 49 49 ]. The initial nucleon distribution in 12 C would be a Woods–Saxon distribution, which originated from the HIJING model [ 46 46 , 47 47 ], or α -cluster configurations [ 33 33 , 34 34 , 50 50 ]. The latter configurations include two cases: the three α -clustering chain structure and the three α -clustering triangle structure. The parameters for the α -clustered 12 C were configured from the EQMD model [ 33 33 , 34 34 , 50 50 ] and are discussed in our previous work [ 31 31 ] in detail. Therefore, we can obtain the phase-space in 12 C + 197 Au collisions at <math id="M1"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV for the calculations using AMPT model. In this work, 12 C is configured as a chain structure, triangle structure with 3- α and the Woods–Saxon distribution of nucleons from the HIJING model [ 46 46 , 47 47 ] packaged in AMPT model. For 197 Au, its initial configuration was just the simple Woods–Saxon distribution. The distribution of the radial center of the α clusters in 12 C is assumed to be a Gaussian function, <math id="M2"><mrow><msup><mi>e</mi><mrow><mo>−</mo><mn>0.5</mn><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>r</mi><mo>−</mo><msub><mi>r</mi><mtext>c</mtext></msub></mrow><mrow><msub><mi>σ</mi><mrow><msub><mi>r</mi><mtext>c</mtext></msub></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow></msup></mrow></math> , where r c is average radial center of an α cluster and σ r c is the width of the distribution. In addition, the nucleon inside each α cluster will be given by the Woods–Saxon distribution. The parameters of r c and σ r c can be obtained using the EQMD calculations [ 33 33 - 36 36 33-36 ]. For the triangle structure, r c = 1.8 fm and σ r c = 0.1 fm. For the chain structure, r c = 2.5 fm and σ r c = 0.1 fm for the two outer α clusters, and the other cluster will be at the center in 12 C. Once the radial center of the α cluster is determined, the centers of the three clusters will be placed in an equilateral triangle for the triangle structure or in a line for the chain structure. Further details can be found in our previous work [ 31 31 , 51 51 ]. In heavy-ion collisions, impact parameter b is defined as the perpendicular distance between the path center of projectile nucleus and the target nucleus. The plane formed by b and the beam direction is called the reaction plane. Because the direction of b is random for every collision, the reaction plane direction will also be random. Hence, the observables that are sensitive to the reaction plane direction should be corrected to that. The participant plane is considered as the reasonable proximation to the reaction plane, and the participant plane angle Ψ n PP is defined by the following equation [ 52 52 - 54 54 52-54 ]: <math display="block" id="M3"><mrow><msub><mi>Ψ</mi><mi>n</mi></msub><mo>{</mo><mi>P</mi><mi>P</mi><mo>}</mo><mo>=</mo><mfrac><mrow><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>tan</mi></mrow></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mrow><mfrac><mrow><mo><</mo><msup><mi>r</mi><mn>2</mn></msup><mi>sin</mi><mrow><mo>(</mo><mrow><mi>n</mi><msub><mi>ϕ</mi><mrow><mtext>part</mtext></mrow></msub></mrow><mo>)</mo></mrow><mo>></mo></mrow><mrow><mo><</mo><msup><mi>r</mi><mn>2</mn></msup><mi>cos</mi><mrow><mo>(</mo><mrow><mi>n</mi><msub><mi>ϕ</mi><mrow><mtext>part</mtext></mrow></msub></mrow><mo>)</mo></mrow><mo>></mo></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mi>π</mi></mrow><mi>n</mi></mfrac><mn>,</mn></mrow></math> where Ψ n PP is the n th-order participant plane angle, r and ϕ part are the coordinate position and azimuthal angle of participants in the collision zone in the initial state, respectively, and the average 〈…〉 denotes the density weighting. For the calculation, the system will be rotated to the participant plane direction in the transverse plane. In AMPT model, we consider a many-body system with discrete particles, and the total angular momentum <math id="M4"><mover accent="true"><mi>J</mi><mo>→</mo></mover></math> can be calculated by summing each particle’s contribution as done for low-energy heavy-ion collisions [ 55 55 ]. <math display="block" id="M5"><mrow><mover accent="true"><mi>J</mi><mo>→</mo></mover><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo>×</mo><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub><mn>,</mn></mrow></math> where <math id="M6"><mrow><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> and <math id="M7"><mrow><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> are the coordinates and momentum of each particle, respectively. The relation between vorticity <math id="M8"><mrow><msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> and velocity <math id="M9"><mrow><msub><mover accent="true"><mi>v</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> is <math id="M10"><mrow><msub><mover accent="true"><mi>v</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> = <math id="M11"><mrow><msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mi>i</mi></msub><mo>×</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> for the particle and <math id="M12"><mrow><msub><mover accent="true"><mi>v</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> = <math id="M13"><mrow><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub><mo>/</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></math> , where E i is the particle’s energy. In this model, the system is symmetric around the system rotational axis <math id="M14"><mover accent="true"><mi>ω</mi><mo>^</mo></mover></math> with a constant vorticity <math id="M15"><mover accent="true"><mi>ω</mi><mo>→</mo></mover></math> . If the system is considered as an approximate rigid body in a fixed stage, each particle’s vorticity <math id="M16"><mrow><msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> has a component projected onto the axis <math id="M17"><mover accent="true"><mi>ω</mi><mo>^</mo></mover></math> , namely the <math id="M18"><mover accent="true"><mi>ω</mi><mo>→</mo></mover></math> and the vertical component will cancel each other out within the system. From the above discussion, the angular momentum can be written as <math display="block" id="M19"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mover accent="true"><mi>J</mi><mo>→</mo></mover><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo>×</mo><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mo>=</mo></mtd><mtd columnalign="left"><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mo stretchy="false">(</mo><mo>|</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><msup><mo>|</mo><mn>2</mn></msup><msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mi>i</mi></msub><mo>−</mo><mo stretchy="false">(</mo><msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mi>i</mi></msub><mo>⋅</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mo>=</mo></mtd><mtd columnalign="left"><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mo stretchy="false">(</mo><mo>|</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mover accent="true"><mi>ω</mi><mo>^</mo></mover><mo>⋅</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><msub><mi>E</mi><mi>i</mi></msub><mover accent="true"><mi>ω</mi><mo>→</mo></mover><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where <math id="M20"><mrow><msub><mi>I</mi><mn>0</mn></msub><mo>=</mo><mstyle displaystyle="true"><mo>∑</mo><mrow><msub><mrow/><mi>i</mi></msub></mrow></mstyle><mo stretchy="false">(</mo><mo>|</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mover accent="true"><mi>ω</mi><mo>^</mo></mover><mo>⋅</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></math> is the moment-of-inertia. Because the phase-space of the system is corrected to event plane angle, the event plane direction will be along the y -axis, and <math id="M21"><mover accent="true"><mi>ω</mi><mo>→</mo></mover></math> will also be along the y -axis. The above discussion and deducing are in the assumption of approximate rigid body system and for the expanding system with discrete particles, the averaged vorticity component along the y -axis can be defined using a similar method from [ 9 9 ] as <math display="block" id="M22"><mrow><mo>〈</mo><msub><mi>ω</mi><mi>y</mi></msub><mo>〉</mo><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mo stretchy="false">(</mo><mo>|</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mover accent="true"><mi>ω</mi><mo>^</mo></mover><mo>⋅</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><msub><mi>E</mi><mi>i</mi></msub><msub><mi>ω</mi><mrow><mi>i</mi><mn>,</mn><mi>y</mi></mrow></msub></mrow><mrow><msub><mi>I</mi><mn>0</mn></msub></mrow></mfrac><mn>.</mn></mrow></math> In the earlier stage, our observation focused on the u( <math id="M23"><mover accent="true"><mi>u</mi><mo>¯</mo></mover></math> ), d( <math id="M24"><mover accent="true"><mi>d</mi><mo>¯</mo></mover></math> ), and s( <math id="M25"><mover accent="true"><mi>s</mi><mo>¯</mo></mover></math> ) quarks, which are primarily composed of partons; therefore, implicated most physics in the collision system. After hardronization, the system experiences a phase transfer to later stages; then, we focused on (anti-)protons and other mesons, including π ± and K ± . The phases we considered for the transport simulation of the collision are the initial partons, freeze-out partons, hadrons without hadronic rescattering, and final-state hadrons. 3 Results and discussion 1 For non-central collisions, the total angular momentum of the collision system would be mainly along the y -axis [ 9 9 ], which in the case of 12 C + 197 Au collision would be of the order of magnitude of 10 5 . Most of the total angular momentum was carried by the spectators on the outer edge of the initial nucleus that passed on without being intervened, while the colliding parts, known as the participants, formed the QGP matter under high temperature and pressure [ 31 31 , 37 37 , 51 51 ].Therefore, because we only considered the participants using Eq. (3), the remnant angular momentum is reduced to the order of magnitude of 10 3 , as indicated in Fig. 3(a) Fig. 3(a) . From the hydrodynamic perspective, the colliding system was nearly isolated [ 56 56 , 57 57 ]. Therefore, the angular momentum would be conserved during the system evolvement after the collision. From Fig. 1(a) to (d) Fig. 1(a) to (d) , it can clearly be seen that J y did not change significantly during the four stages: initialed participants, partons at freeze-out stage, hadrons without hadron rescattering, and hadrons after hadron rescattering. However, the QGP matter soon cooled down and then froze out, scattering out the particles that had a high momentum [ 51 51 , 56 56 ]. Therefore, the system expanded, and the total angular inertia increased significantly from the order of magnitude of 10 5 to 10 7 , which is shown in Fig. 2 (a) to (d) Fig. 2 (a) to (d) . On the other hand, the decreasing average angular velocity ⟨ ω y ⟩ calculated using Eq. (4) could also manifest the system expansion in Fig. 3 Fig. 3 , corresponding to the increasing I 0 . In Fig. 1 (a) Fig. 1 (a) , the angular momentum of the participants first increased then decreased when the impact parameter b was increased gradually, which is in agreement with the results of other works [ 9 9 ]. This tendency can be understood from Eq. (2); in near-central collision, the angular momentum was dominated by the distance between each discrete particle and the center of mass, the increase of which contributed to the rising total angular momentum. However, the 197 Au nucleus employed the Woods–Saxon distribution [ 21 21 - 23 23 21-23 ], which possessed the Gaussian density drop against the radius; therefore, fewer initial nucleons (participants) could be contributed to the colliding process as the impact parameter increased. Because 12 C occupied a smaller region than 197 Au [ 19 19 , 20 20 ], the Woods–Saxon distribution of Au also induced the decreasing tendency of the total angular inertia at a larger impact parameter, as shown in Fig. 2 Fig. 2 . While for peripheral collision with larger b, the number of participants decreased rapidly, and at approximately 6 fm, which is approximately equal to the radius of the nucleus of 197 Au in the Woods–Saxon model [ 37 37 ], leading to a more rapid decay for a larger impact parameter. Studies on heavy-ion collision [ 9 9 , 58 58 ] have revealed a global rotation system with large angular momentum. In Fig.4 Fig.4 , we showed the distribution of the average angular velocity along the y-axis ω y in the rapidity-transverse momentum plane, which implied that most participants carried a small angular velocity of the order of magnitude of approximately 10 -2 . Correspondingly, in the second panel of Fig. 3 Fig. 3 , for the collective behavior, the total system at the parton freeze-out state was estimated to have an average angular velocity approximately equal to 0.021/fm, where we discussed the case of a middle range of impact parameter at b =4.5 fm. In most panel of Fig.4 Fig.4 , the range of ω y distribution in rapidity at relatively high p T was large compared to the low transverse momentum, although most of them remained around the mid-rapidity range. However, for the particles with relatively lower transverse momentum, rapidity range was narrower and at a larger magnitude, implying a strong collective structure, which could probably be seen as collective flow [ 31 31 , 51 51 ]. In addition, from the density distribution in Fig.5 Fig.5 , it can be seen that the center of collision region that composed of participants stayed with collective flow before being scattered out in the early stage. In the scenario [ 31 31 , 51 51 ], the 197 Au nucleus would move towards the negative direction along the rapidity axis, while 12 C took the opposite direction, which after collision transferred to a QGP matter and formed the global rotation [ 9 9 ]. In Fig. 5 Fig. 5 , for the impact parameter at 4.5 fm, most of the particle density was distributed in the range of low p T and mid-rapidity in Fig.5 Fig.5 . In the left column of Fig. 5 Fig. 5 , the two most dense regions at the parton freeze-out state represent the original 12 C and 197 Au nucleus [ 19 19 - 23 23 19-23 ]. With time evolution, the two nuclei in opposite directions, which is shown in the middle panels, and their transverse momentum expanded towards a larger range. After the particles scattered out to the final-state hadrons, the system had a lower density that was chiefly distributed in the range of low p T and mid-rapidity. Comparing the left and the middle column of Fig. 4 Fig. 4 , the angular velocity ω y moved towards the region of larger p T , which manifested the expansion of the system. At the same time, the magnitude of ω y decrease to 10 -3 , which corresponds to the collective ⟨ ω y ⟩ in Fig. 3 (c) and (d) Fig. 3 (c) and (d) . In the right column, some scattered hadrons at a higher p T and rapidity obtained higher angular momentum, which consequently led to a larger ω y . Few hadrons at large rapidity edge even carried large ω y that represent the high velocity scattered particles at the edge of rotation. In a previous work on the collective flow [ 51 51 ], we reported the configurations of Woods–Saxon distribution and two α -clustering distributions for 12 C to have very distinctive behaviors in the collective flow such as v 2 and v 3 . However, in this study, for angular velocity, we found that the differences among the three structures were not as evident as the previous cases. For instance, global behaviors of the system such as I 0 and 〈 ω y 〉 in Fig. 2 Fig. 2 - 3 3 perceived a rather small discrepancy among the three different cases of 12 C configuration. On the other hand, J y seemed to show a more significant difference, where Woods–Saxon distribution possessed the largest angular momentum J y , followed by the triangle clusters structure, and then the chain clusters. The discrepancy between Woods–Saxon distribution and chain could be as large as approximately 20-30% at the peak, with an impact parameter equal to 6 fm. The internal configurations of α -cluster inside 12 C [ 33 33 , 34 34 ] determined the initial geometry of the collision system. At and after the system phase transition, it is highly probable that such geometry remained in the system [ 51 51 ] and could be manifested in later evolvement. Even in the final state, shown in (d) of Fig. 1 Fig. 1 , the discrepancy of the angular momentum could be large as 20% at peak. However, because of the expansion for flow, the average speed of system was expected to gradually decrease, making the distinction difficult. It is much likely that more interesting mechanism could be suggested in physical behavior of collective behavior. 4 Summary 1 In summary, we simulated the collisions of 12 C against 197 Au at <math id="M28"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV using AMPT model, where three different structure patterns of C were considered separately and compared. We found the effects on the angular momentum of the chain and triangle α -clustered cases as well as the Woods–Saxon distribution case to have a discrepancy. Other collective behaviors, including angular inertia and average angular velocity, were also examined and were found to have much less distinction. This is consistent with the density distribution and angular velocity distribution for the three patterns, suggesting possible analysis for experimental study. References L. Adamczyk , J. K. Adkins , G. Agakishiev et al . (STAR Collaboration) , Beam-energy dependence of the directed flow of protons, antiprotons, and pions in Au+Au collisions . Phys. Rev. Lett . 112 , 162301 ( 2014 ). doi: 10.1103/PhysRevLett.112.162301 http://doi.org/10.1103/PhysRevLett.112.162301 Beam-energy dependence of the directed flow of protons, antiprotons, and pions in Au+Au collisions L. Adamczyk , J. K. Adkins , G. Agakishiev et al . (STAR Collaboration) , Elliptic flow of identified hadrons in Au + Au collisions at <math id="M29"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 7.7-62.4 GeV . Phys. Rev. C 88 , 014902 ( 2013 ). doi: 10.1103/PhysRevC.88.014902 http://doi.org/10.1103/PhysRevC.88.014902 Elliptic flow of identified hadrons in Au + Au collisions at sNN = 7.7-62.4 GeV L. Adamczyk , J. K. Adkins , G. Agakishiev , et al . (STAR Collaboration) , Third harmonic flow of charged particles in Au + Au collisions at <math id="M30"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. C 88 , 014904 ( 2013 ). doi: 10.1103/PhysRevC.88.014904 http://doi.org/10.1103/PhysRevC.88.014904 Third harmonic flow of charged particles in Au + Au collisions at sNN = 200 GeV L. Adamczyk , J.K. Adkins , G. Agakishiev , et al . (STAR Collaboration) , Observation of charge asymmetry dependence of pion elliptic flow and the possible chiral magnetic wave in heavy-ion collisions . Phys. Rev. Lett . 114 , 252302 ( 2015 ). doi: 10.1103/PhysRevLett.114.252302 http://doi.org/10.1103/PhysRevLett.114.252302 Observation of charge asymmetry dependence of pion elliptic flow and the possible chiral magnetic wave in heavy-ion collisions J. Adams , C. Adler , M. M. Aggarwal , et al . (STAR Collaboration) , Evidence from d+ Au measurements for final-state suppression of high-pT hadrons in Au+Au collisions at RHIC . Phys. Rev. Lett . 91 , 072304 ( 2003 ). doi: 10.1103/PhysRevLett.91.072304 http://doi.org/10.1103/PhysRevLett.91.072304 Evidence from d+ Au measurements for final-state suppression of high-pT hadrons in Au+Au collisions at RHIC J. Adams , C. Adler , M. M. Aggarwal , et al . (STAR Collaboration) , Transverse-momentum and collision-energy dependence of high-pT hadron suppression in Au+Au collisions at ultrarelativistic energies . Phys. Rev. Lett . 91 , 172302 ( 2003 ). doi: 10.1103/PhysRevLett.91.172302 http://doi.org/10.1103/PhysRevLett.91.172302 Transverse-momentum and collision-energy dependence of high-pT hadron suppression in Au+Au collisions at ultrarelativistic energies L. Adamczyk , J. K. Adkins , G. Agakishiev et al . (STAR Collaboration) , Beam-energy-dependent two-pion interferometry and the freeze-out eccentricity of pions measured in heavy ion collisions at the STAR detector . Phys. Rev. C 92 , 014904 ( 2015 ). doi: 10.1103/PhysRevC.92.014904 http://doi.org/10.1103/PhysRevC.92.014904 Beam-energy-dependent two-pion interferometry and the freeze-out eccentricity of pions measured in heavy ion collisions at the STAR detector A. Adare , et al . (PHENIX Collaboration) , Beam-energy and system-size dependence of the space-time extent of the pion emission source produced in heavy ion collisions . arXiv:1410.2559 ( 2014 ). Beam-energy and system-size dependence of the space-time extent of the pion emission source produced in heavy ion collisions Y. Jiang , Z.-W. Lin , J. Liao , Rotating quark-gluon plasma in relativistic heavy-ion collisions . Phys. Rev. C 94 , 044910 ( 2016 ) doi: 10.1103/PhysRevC.94.044910 http://doi.org/10.1103/PhysRevC.94.044910 Rotating quark-gluon plasma in relativistic heavy-ion collisions Z.-T. Liang , X.-N. Wang , Globally polarized quark-gluon plasma in noncentral A+A collisions . Phys. Rev. Lett . 94 , 102301 ( 2005 ). doi: 10.1103/PhysRevLett.94.102301 http://doi.org/10.1103/PhysRevLett.94.102301 Globally polarized quark-gluon plasma in noncentral A+A collisions X.-G. Huang , P. Huovinen , X.-N. Wang , Quark polarization in a viscous quark-gluon plasma . Phys. Rev. C 84 , 054910 ( 2011 ). doi: 10.1103/PhysRevC.84.054910 http://doi.org/10.1103/PhysRevC.84.054910 Quark polarization in a viscous quark-gluon plasma A. H. Tang , B. Tu , C. S. Zhou , Practical considerations for measuring global spin alignment of vector mesons in relativistic heavy ion collisions . Phys. Rev. C 98 , 044907 ( 2018 ). doi: 10.1103/PhysRevC.98.044907 http://doi.org/10.1103/PhysRevC.98.044907 Practical considerations for measuring global spin alignment of vector mesons in relativistic heavy ion collisions Z.-T. Liang , X.-N. Wang , Spin alignment of vector mesons in non-central A+A collisions . Phys. Lett. B 629 , 20 - 26 ( 2005 ). doi: 10.1016/j.physletb.2005.09.060 http://doi.org/10.1016/j.physletb.2005.09.060 Spin alignment of vector mesons in non-central A+A collisions L. Adamczyk and et al . (STAR Collaboration) , Global Λ hyperon polarization in nuclear collisions . Nature , 541 , 62 - 65 ( 2017 ). doi: 10.1038/nature23004 http://doi.org/10.1038/nature23004 Global Λ hyperon polarization in nuclear collisions R. D. de Souza , J. Takahashi , T. Kodama , et al ., Effects of initial state fluctuations in the final state elliptic flow measurements using the NeXSPheRIO model . Phys. Rev. C 85 , 054909 ( 2012 ). doi: 10.1103/PhysRevC.85.054909 http://doi.org/10.1103/PhysRevC.85.054909 Effects of initial state fluctuations in the final state elliptic flow measurements using the NeXSPheRIO model A. Chaudhuri , Fluctuating initial conditions and fluctuations in elliptic and triangular flow . Phys. Lett. B 710 , 339 - 342 ( 2012 ). doi: 10.1016/j.physletb.2012.02.067 http://doi.org/10.1016/j.physletb.2012.02.067 Fluctuating initial conditions and fluctuations in elliptic and triangular flow H. Song , S. A. Bass , U. Heinzi , Elliptic flow in <math id="M31"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV Au + Au collisions and <math id="M32"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 2.76 TeV Pb + Pb collisions: Insights from viscous hydrodynamics + hadron cascade hybrid model . Phys. Rev. C 83 , 054912 ( 2011 ). doi: 10.1103/PhysRevC.83.054912 http://doi.org/10.1103/PhysRevC.83.054912 Elliptic flow in s = 200 GeV Au + Au collisions and s = 2.76 TeV Pb + Pb collisions: Insights from viscous hydrodynamics + hadron cascade hybrid model S. Floerchinger , U. A. Wiedemann , Mode-by-mode fluid dynamics for relativistic heavy ion collisions . Phys. Lett. B , 728 , 407 - 411 ( 2014 ). doi: 10.1016/j.physletb.2013.12.025 http://doi.org/10.1016/j.physletb.2013.12.025 Mode-by-mode fluid dynamics for relativistic heavy ion collisions P. Bozek , W. Broniowski , E. R. Arriola , et al ., α clusters and collective flow in ultrarelativistic carbon-heavy-nucleus collisions . Phys. Rev. C 90 , 064902 ( 2014 ). doi: 10.1103/PhysRevC.90.064902 http://doi.org/10.1103/PhysRevC.90.064902 α clusters and collective flow in ultrarelativistic carbon-heavy-nucleus collisions W. Broniowski and E. R. Arriola , Signatures of α clustering in light nuclei from relativistic nuclear collisions . Phys. Rev. Lett . 112 , 112501 ( 2014 ). doi: 10.1103/PhysRevLett.112.112501 http://doi.org/10.1103/PhysRevLett.112.112501 Signatures of α clustering in light nuclei from relativistic nuclear collisions L. Ma , G. L. Ma , Y. G. Ma , Initial partonic eccentricity fluctuations in a multiphase transport model . Phys. Rev. C 94 , 044915 ( 2016 ). doi: 10.1103/PhysRevC.94.044915 http://doi.org/10.1103/PhysRevC.94.044915 Initial partonic eccentricity fluctuations in a multiphase transport model L. Ma , G. L. Ma , Y. G. Ma , Anisotropic flow and flow fluctuations for Au + Au at <math id="M33"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV in a multiphase transport model . Phys. Rev. C 89 , 044907 ( 2014 ). doi: 10.1103/PhysRevC.89.044907 http://doi.org/10.1103/PhysRevC.89.044907 Anisotropic flow and flow fluctuations for Au + Au at sNN=200 GeV in a multiphase transport model L.X. Han , G.L. Ma , Y.G. Ma , et al ., Initial fluctuation effect on harmonic flows in high-energy heavy-ion collisions . Phys. Rev. C 84 , 064907 ( 2011 ). doi: 10.1103/PhysRevC.84.064907 http://doi.org/10.1103/PhysRevC.84.064907 Initial fluctuation effect on harmonic flows in high-energy heavy-ion collisions H.-C. Song , Y. Zhou , K. Gajdosova , Collective flow and hydrodynamics in large and small systems at the LHC . Nucl. Sci. Tech . 28 , 99 ( 2017 ). doi: 10.1007/s41365-017-0245-4 http://doi.org/10.1007/s41365-017-0245-4 Collective flow and hydrodynamics in large and small systems at the LHC J. Wang , Y. G. Ma , G. Q. Zhang , et al ., Initial fluctuation effect on elliptic flow in Au+Au collision at 1 GeV/A . Nucl. Sci. Tech . 24 , 030501 ( 2013 ). doi: 10.13538/j.1001-8042/nst.2013.03.004 http://doi.org/10.13538/j.1001-8042/nst.2013.03.004 Initial fluctuation effect on elliptic flow in Au+Au collision at 1 GeV/A C. C. Guo , W. B. He , Y. G. Ma , Collective flows of 16O+16O collisions with α-clustering configurations . Chinese Phys. Lett . 34 , 092101 ( 2017 ). doi: 10.1088/0256-307X/34/9/092101 http://doi.org/10.1088/0256-307X/34/9/092101 Collective flows of 16O+16O collisions with α-clustering configurations X.-F. Luo , N. Xu , Search for the QCD critical point with fluctuations of conserved quantities in relativistic heavy-ion collisions at RHIC: An overvie . Nucl. Sci. Tech . 28 , 112 ( 2017 ). doi: 10.1007/s41365-017-0257-0 http://doi.org/10.1007/s41365-017-0257-0 Search for the QCD critical point with fluctuations of conserved quantities in relativistic heavy-ion collisions at RHIC: An overvie X. Jin , J. Chen , Z. Lin , et al ., Explore the QCD phase transition phenomena from a multiphase transport model . Sci. China Phys. Mech . 62 , 11012 ( 2018 ). doi: 10.1007/s11433-018-9272-4 http://doi.org/10.1007/s11433-018-9272-4 Explore the QCD phase transition phenomena from a multiphase transport model C. M. Ko , F. Li , Density fluctuations in baryon-rich quark matter . Nucl. Sci. Tech . 27 , 140 ( 2016 ). doi: 10.1007/s41365-016-0141-3 http://doi.org/10.1007/s41365-016-0141-3 Density fluctuations in baryon-rich quark matter Q.Y. Shou , G.L. Ma , Y.G. Ma , Charge separation with fluctuating domains in relativistic heavy-ion collisions . Phys. Rev. C 90 , 047901 ( 2014 ). doi: 10.1103/PhysRevC.90.047901 http://doi.org/10.1103/PhysRevC.90.047901 Charge separation with fluctuating domains in relativistic heavy-ion collisions S. Zhang , Y. G. Ma , J. H. Chen , et al ., Nuclear cluster structure effect on elliptic and triangular flows in heavy-ion collisions . Phys. Rev. C 95 , 064904 ( 2017 ). doi: 10.1103/PhysRevC.95.064904 http://doi.org/10.1103/PhysRevC.95.064904 Nuclear cluster structure effect on elliptic and triangular flows in heavy-ion collisions F. D. Murnaghan , Review: G. Gamow, Constitution of Atomic Nuclei and Radioactivity . Bull. Amer. Math. Soc . 39 , 487 ( 1993 ) Review: G. Gamow, Constitution of Atomic Nuclei and Radioactivity W.B. He , Y.G. Ma , X.G. Cao , et al ., Giant dipole resonance as a fingerprint of α clustering configurations in 12C and 16O . Phys. Rev. Lett. bf 113 , 032506 ( 2014 ). doi: 10.1103/PhysRevLett.113.032506 http://doi.org/10.1103/PhysRevLett.113.032506 Giant dipole resonance as a fingerprint of α clustering configurations in 12C and 16O W.B. He , Y.G. Ma , X.G. Cao , et al ., Dipole oscillation modes in light α-clustering nuclei . Phys. Rev. C , 94 , 014301 ( 2016 ). doi: 10.1103/PhysRevC.94.014301 http://doi.org/10.1103/PhysRevC.94.014301 Dipole oscillation modes in light α-clustering nuclei B. S. Huang , Y. G. Ma , W. B. He , Photonuclear reaction as a probe for α-clustering nuclei in the quasi-deuteron region . Phys.Rev. C , 95 , 034606 ( 2017 ). doi: 10.1103/PhysRevC.95.034606 http://doi.org/10.1103/PhysRevC.95.034606 Photonuclear reaction as a probe for α-clustering nuclei in the quasi-deuteron region B. S. Huang , Y. G. Ma , W. B. He , Alpha-clustering effects on 16O(4γ, np)14N in quasi-deuteron region . Eur. Phys. J. A 53 , 119 ( 2017 ). Alpha-clustering effects on 16O(4γ, np)14N in quasi-deuteron region Z.-W. Lin , C. M. Ko , B.-A. Li , et al ., Multiphase transport model for relativistic heavy ion collisions . Phys. Rev. C 72 , 064901 ( 2005 ). doi: 10.1103/PhysRevC.72.064901 http://doi.org/10.1103/PhysRevC.72.064901 Multiphase transport model for relativistic heavy ion collisions G.-L. Ma , Z.-W. Lin , Predictions for <math id="M34"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt><mo>=</mo><mn>5.02</mn></mrow></math> TeV Pb + Pb collisions from a multiphase transport model . Phys. Rev. C 93 , 054911 ( 2016 ). doi: 10.1103/PhysRevC.93.054911 http://doi.org/10.1103/PhysRevC.93.054911 Predictions for sNN=5.02 TeV Pb + Pb collisions from a multiphase transport model Z.-W. Lin , C. M. Ko , S. Pal , Partonic effects on pion interferometry at the relativistic heavy-ion collider . Phys. Rev. Lett . 89 , 152301 ( 2002 ). doi: 10.1103/PhysRevLett.89.152301 http://doi.org/10.1103/PhysRevLett.89.152301 Partonic effects on pion interferometry at the relativistic heavy-ion collider B. I. Abelev , M. M. Aggarwal , Z. Ahammed , et al . (STAR Collaboration) , System-size independence of directed flow measured at the BNL relativistic heavy-ion collider . Phys. Rev. Lett . 101 , 252301 ( 2008 ). doi: 10.1103/PhysRevLett.101.252301 http://doi.org/10.1103/PhysRevLett.101.252301 System-size independence of directed flow measured at the BNL relativistic heavy-ion collider A. Bzdak and G.-L. Ma , Elliptic and triangular flow in p-Pb and peripheral Pb-Pb collisions from parton scatterings . Phys. Rev. Lett . 113 , 252301 ( 2014 ). doi: 10.1103/PhysRevLett.113.252301 http://doi.org/10.1103/PhysRevLett.113.252301 Elliptic and triangular flow in p-Pb and peripheral Pb-Pb collisions from parton scatterings G.-L. Ma , S. Zhang , Y.-G. Ma , et al ., Di-hadron azimuthal correlation and Mach-like cone structure in a parton/hadron transport model . Phys. Lett. B 641 , 362 - 367 ( 2006 ). Di-hadron azimuthal correlation and Mach-like cone structure in a parton/hadron transport model X.-H. Jin , J.-H. Chen , Y.-G. Ma , et al ., Ω and ϕ production in Au+Au collisions at =11.5 and 7.7 GeV in a dynamical quark coalescence model . Nucl. Sci. Tech . 29 , 54 ( 2018 ). doi: 10.1007/s41365-018-0393-1 http://doi.org/10.1007/s41365-018-0393-1 Ω and ϕ production in Au+Au collisions at =11.5 and 7.7 GeV in a dynamical quark coalescence model K. Hattori , X.-G. Huang , Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions . Nucl. Sci. Tech . 28 , 26 ( 2017 ). Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions X.-L. Zhao , Y.-G. Ma , G.-L. Ma , Electromagnetic fields in small systems from a multiphase transport model . Phys. Rev. C 97 , 024910 ( 2018 ). doi: 10.1103/PhysRevC.97.024910 http://doi.org/10.1103/PhysRevC.97.024910 Electromagnetic fields in small systems from a multiphase transport model X.-N. Wang , M. Gyulassy , hijing: A Monte Carlo model for multiple jet production in pp, pA, and AA collisions . Phys. Rev. D 44 , 3501 - 3516 .( 1991 ). doi: 10.1103/PhysRevD.44.3501 http://doi.org/10.1103/PhysRevD.44.3501 hijing: A Monte Carlo model for multiple jet production in pp, pA, and AA collisions M. Gyulassy , X.-N. Wang , HIJING 1.0: A Monte Carlo program for parton and particle production in high energy hadronic and nuclear collisions . Comput. Phys. Commun . 83 , 307 ( 1994 ). doi: 10.1016/0010-4655(94)90057-4 http://doi.org/10.1016/0010-4655(94)90057-4 HIJING 1.0: A Monte Carlo program for parton and particle production in high energy hadronic and nuclear collisions B. Zhang , ZPC 1.0.1: a parton cascade for ultrarelativistic heavy ion collisions . Computer Physics Communications 109 , 193 - 206 ( 1998 ). ZPC 1.0.1: a parton cascade for ultrarelativistic heavy ion collisions B.-A. Li , C. M. Ko , Formation of superdense hadronic matter in high energy heavy-ion collisions . Phys. Rev. C 52 , 2037 - 2063 ( 1995 ). doi: 10.1103/PhysRevC.52.2037 http://doi.org/10.1103/PhysRevC.52.2037 Formation of superdense hadronic matter in high energy heavy-ion collisions T. Maruyama , K. Niita , A. Iwamoto , Extension of quantum molecular dynamics and its application to heavy-ion collisions . Phys. Rev. C 53 , 297 - 304 ( 1996 ). doi: 10.1103/PhysRevC.53.297 http://doi.org/10.1103/PhysRevC.53.297 Extension of quantum molecular dynamics and its application to heavy-ion collisions S. Zhang , Y. G. Ma , J. H. Chen , et al ., Collective flows of α-clustering 12C + 197Au by using different flow analysis methods . Euro. Phys. J. A 54 , 161 ( 2018 ). doi: 10.1140/epja/i2018-12597-y http://doi.org/10.1140/epja/i2018-12597-y Collective flows of α-clustering 12C + 197Au by using different flow analysis methods S. A. Voloshin , A. M. Poskanzer , A. Tang , et al ., Elliptic flow in the Gaussian model of eccentricity fluctuations . Phys. Lett. B 659 , 537 - 541 ( 2008 ). doi: 10.1016/j.physletb.2007.11.043 http://doi.org/10.1016/j.physletb.2007.11.043 Elliptic flow in the Gaussian model of eccentricity fluctuations B. Alver , G. Roland , Collision-geometry fluctuations and triangular flow in heavy-ion collisions . Phys. Rev. C 81 , 054905 ( 2010 ). doi: 10.1103/PhysRevC.81.054905 http://doi.org/10.1103/PhysRevC.81.054905 Collision-geometry fluctuations and triangular flow in heavy-ion collisions R. A. Lacey , R. Wei , J. Jia , et al ., Initial eccentricity fluctuations and their relation to higher-order flow harmonics . Phys. Rev. C 83 , 044902 ( 2011 ). doi: 10.1103/PhysRevC.83.044902 http://doi.org/10.1103/PhysRevC.83.044902 Initial eccentricity fluctuations and their relation to higher-order flow harmonics Y. G. Ma , W. Q. Shen , Z. Y. Zhu , Collective motions of reverse reaction system in the intermediate energy domain via the quantum molecular dynamics arrpoach . Phys. Rev. C 51 , 1029 ( 1995 ). doi: 10.1103/PhysRevC.51.1029 http://doi.org/10.1103/PhysRevC.51.1029 Collective motions of reverse reaction system in the intermediate energy domain via the quantum molecular dynamics arrpoach L. He , T. Edmonds , Z.-W. Lin , et al ., Phys. Lett. B 753 , 506 - 510 ( 2016 ). doi: 10.1016/j.physletb.2015.12.051 http://doi.org/10.1016/j.physletb.2015.12.051 Phys. Lett. B G.-L. Ma , A. Bzdak , Flow in small systems from parton scatterings . Nucl. Phys. A 956 , 745 - 748 ( 2016 ). doi: 10.1016/j.nuclphysa.2016.01.057 http://doi.org/10.1016/j.nuclphysa.2016.01.057 Flow in small systems from parton scatterings W.-T. Deng , X.-G. Huang , Vorticity in heavy-ion collisions . Phys. Rev. C 93 , 064907 ( 2016 ). doi: 10.1103/PhysRevC.93.064907 http://doi.org/10.1103/PhysRevC.93.064907 Vorticity in heavy-ion collisions 10.1007/s41365-018-0523-9 This work was supported in part by National Key R&D Program of China (No. 2016YFE0100900), the National Natural Science Foundation of China (Nos. 11421505, 11220101005, 11775288 and U1232206), the Major State Basic Research Development Program in China (No. 2014CB845400), the Key Research Program of Frontier Sciences of the CAS (No. QYZDJ-SSW-SLH002), and the Key Research Program of the Chinese Academy of Sciences (NO. XDPB09). 2018-11-20 2018-11-26 2018-11-26 2018-12-03 2018-12-01 2021-12-20 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0523-9 http://dx.doi.org/10.1007/s41365-018-0523-9 . 2018 Nuclear Science and Techniques 12 29