Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β₂ and triaxiality γ . The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ɛ₂ and the density gradient d_⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances «ɛ₂²» or «(δ^{d_⊥/d_⊥) 2}» and skewnesses «ɛ₂²δ d_⊥/d_⊥» or «(δ^{d_⊥/d_⊥) 3}» have a simple analytical form of a^'+b^'β₂² and a^'+[b^'+c^'cos(3 γ ) ] β₂³ , respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β₂ and show that the correlations between a normalized skewness and a variance can constrain simultaneously β₂ and γ . Assuming a linear relation with elliptic flow v₂ and mean-transverse momentum [p_T] of final-state particles, v₂∝ɛ₂ and δ [p_T] /[p_T] ∝δ d_⊥/d_⊥ , similar conclusions are also expected for the variances and skewnesses of v₂ and [p_T] , i.e., a +b β₂² for «v₂²» and «(δ^{[p_T] /[p_T] ) 2}» and a +[b +c cos(3 γ ) ] β₂³ for «v₂²δ [p_T] /[p_T] » or «(δ^{[p_T] /[p_T] ) 3}» . My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos(3 γ )≈1 , and near oblate nuclei, cos(3 γ )≈−1 , with known β₂ values, followed by collisions of other species of interest with similar mass number. The (β₂,γ ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. The results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.