
The nuclear structure parameterization includes half-density radius (\(R_0\) ), surface diffuseness (\(a_0\) ), and axial symmetric quadrupole deformation (\(\beta _2\) ) and octupole deformation (\(\beta _3\) ). Model studies of elliptic flow (\(v_2\) ) and triangular flow (\(v_3\) ) suggest \(\beta _2\sim 0.15\)  for \(^{96}\) Ru and \(\beta _3\sim 0.2\)  for \(^{96}\) Zr, respectively [1]}. In the mid-central collisions, the comparison of \(v_2\)  and the multiplicity distribution between the two isobaric systems supports the idea that \(^{96}\) Zr has a thicker neutron skin than \(^{96}\) Ru [2]}, [3]}. Predictions for other observables and their sensitivities to deformation and neutron skin have also been made, such as mean transverse momentum \(\left\langle p_{\mathrm {T}}\right\rangle \)  [4]}, \(\left\langle p_{\mathrm {T}}\right\rangle \)  fluctuations [5]}, and \(v_n\) –\(p_{\mathrm {T}}\)  correlations [6]}, [7]}, [8]}.
