Nuclear cluster structure effect in $^{16}$O+$^{16}$O collisions at the top RHIC energy

Using the improved string-melting version of a Multi-Phase Transport model, we investigated the impact of nuclear geometry of $^{16}$O on anisotropic flows in O+O collisions at $\sqrt{s_{\rm NN}} = 200$ GeV. To evaluate the influence of nuclear structure and potential alpha clustering, we implemented four candidate configurations: Woods-Saxon, tetrahedron, square, and Nuclear Lattice Effective Field Theory. Initial-state geometry is quantified via the eccentricity cumulant ratio $\varepsilon_{2}\{4\}/\varepsilon_{2}\{2\}$, which provides a robust and evolution-independent measure sensitive to configuration differences. The model reproduces $v_{2}(p_{\rm T})$ at low $p_{\rm T}$ and $v_{3}(p_{\rm T})$ across the full $p_{\rm T}$ range, with integrated $v_{2}\{2\}$ and $v_{3}\{2\}$ matching the STAR data, demonstrating that transport dynamics captures the essential collectivity in this intermediate-size system. These findings establish a baseline for extending nuclear-structure studies in O+O collisions to other energies and differential observables within a unified transport model framework.