Search

Magnetic fields in heavy-ion collisions are pivotal and subject to diverse factors. In this study, we quantitatively investigate the impact of fluid vorticity on the evolution of magnetic fields in the 20-50\% centrality class in Au+Au collisions, with collision energies of $\sqrt{s_{NN}}=(7.7, 14.5, 19.6, 27, 39, 62.4, 200)$ GeV. Our results indicate that fluid vorticity leads to a delay in the evolution of the magnetic field, in which this effect becomes more pronounced as the collision energy decreases. Additionally, we have calculated the mean magnetic field values on the freeze-out hypersurface for various collision energies. Our simulation results align with the values inferred from experimental data of $\bar{\Lambda}-\Lambda$, within the error margins.

We have implemented the 3+1 dimensional CLVisc hydrodynamics model with TRENTO-3D initial conditions to investigate the spin polarization of $\Lambda$ hyperons along the beam direction in p+Pb collisions at $\sqrt{s_{NN}} = 8.16$ TeV. Following our previous theoretical framework based on quantum kinetic theory, we consider three different scenarios: $\Lambda$ equilibrium, $s$ quark equilibrium, and iso-thermal equilibrium scenarios. We have computed the second Fourier sine coefficients of spin polarization along the beam direction, denoted as $\left\langle P_{z} \sin 2(\phi_{p} - \Psi_{2}) \right\rangle$, with $\phi_{p} - \Psi_{2}$ being the azimuthal angle relative to the second-order event plane $\Psi_{2}$, as functions of multiplicity, transverse momentum and pseudo-rapidity in the three scenarios. Additionally, we have also computed the spin polarization along the beam direction, $P_{z}$, as a function of the azimuthal angle. We find that the spin polarization induced by thermal vorticity always provides an opposite contribution compared to the shear-induced polarization in p+Pb collisions. The total spin polarization computed by the current hydrodynamic model disagrees with the data measured by LHC-CMS experiments. Our findings imply that other non-flow effects may play a crucial role in p+Pb collisions., Comment: 8 pages, 5 figures and 1 table. A new figure for local polarization as a function of pseudo-rapidity is added. submitted to PRC

A systematic study of intermediate invariant mass dilepton production in Pb+Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV is performed, using next-to-leading-order (NLO) thermal QCD dilepton emission rates with a multistage dynamical approach which includes event-by-event IP-Glasma initial conditions, relativistic viscous fluid dynamics, and a hadronic afterburner. Considering dilepton yield and anisotropic flow, special attention is paid to the out-of-equilibrium aspects, both thermal and chemical, and to the contribution of the Drell-Yan process. The relative contribution of each of those different channels to dilepton observables is calculated and discussed., Comment: 13 pages, 10 figures

In the SMASH-CLVisc-hybrid framework, including SMASH for the initial conditions and the hadronic rescattering stage, and CLVisc for the quark gluon plasma (QGP) evolution, we investigate net baryon number fluctuations via considering the equation of state (EoS) with and without a critical end point (CEP) in the QCD phase transition. Specifically, two distinct QCD EoS are utilized: one with smooth crossover derived from NEOS and another with a critical end point sourced from the rPNJL model. Our results show that non-monotonic behavior of $\kappa\sigma^2$ is not observed above the collision energy $7.7 {\rm GeV}$, nor are there explicit differences between the EoS characterized by crossover and that by CEP. This could be attributed to the significant deviation of the freeze-out line from the location of CEP. It is also found that the pure SMASH result of $\kappa \sigma^2$ is positive and close to zero at 3 GeV, which is different from the negative value observed from STAR.

We study the charmed meson and baryon production and elliptic flow in ultra-relativistic nucleus-nucleus collisions at the LHC energies. The space-time evolution of quark-gluon plasma (QGP) produced in these energetic collisions is obtained via the (3+1)-dimensional CLVisc hydrodynamics model, the heavy quark dynamics inside the QGP is simulated using an improved Langevin model that incorporates both elastic and inelastic parton energy loss processes, and the heavy quark hadronization is simulated utilizing a comprehensive coalescence-fragmentation model. Using our combined approach, we first calculate charmed hadron ratios, $\Lambda_c/D^0$ and $D_s/D^0$, as well as their elliptic flow ($v_2$) as a function of transverse momentum ($p_T$) for different centralities in Pb+Pb collisions at $\sqrt{s_{NN}}=5.02$~TeV. Due to strangeness enhancement and parton coalescence effects, $D_s/D^0$ and $\Lambda_c/D^0$ ratios increase from peripheral to central collisions, and such centrality dependence for $\Lambda_c/D^0$ is stronger than $D_s/D^0$. We further predict the $p_T$ and centrality dependences of charmed hadron chemistry and $v_2$ in smaller Xe+Xe, Ar+Ar and O+O collisions at the LHC energies. Strong centrality and system size dependences for $\Lambda_c/D^0$ and $D_s/D^0$ ratios are observed across four collision systems. As for charmed hadron flow, both system size and collision geometry are important to understand the centrality dependence of $v_2$ in different collision systems. Our study provides a significant reference for studying heavy quark evolution and hadronizaiton in large and small systems in relativistic nuclear collisions., Comment: 12 pages, 10 figures

Quantum fluctuations plays an essential role in forming the collective flow of hadrons observed in relativistic heavy-ion collisions. Event-by-event fluctuations of the collective flow can arise from various sources, such as the fluctuations in the initial geometry, hydrodynamic expansion, hadronization, and hadronic evolution of the nuclear matter, while the exact contribution from each source is still an open question. Using a (3+1)-dimensional relativistic hydrodynamic model coupled to a Monte-Carlo Glauber initial condition, Cooper-Frye particlization and a hadronic transport model, we explore the system size and shape dependences of the collective flow fluctuations in Au+Au, Cu+Au, and O+O collisions at $\sqrt{s_\mathrm{NN}}=200$~GeV. The particle yields, mean transverse momenta, 2-particle and 4-particle cumulant elliptic flows ($v_2\{2\}$ and $v_2\{4\}$) from our calculation agree with the currently existing data from RHIC. Different centrality dependences of the flow fluctuations, quantified by the $v_2\{4\}/v_2\{2\}$ ratio, are found for different collision systems due to their different sizes and shapes. By comparing $v_2\{4\}/v_2\{2\}$ between different hadron species, and comparing $v_2\{4\}/v_2\{2\}$ to the initial state geometric fluctuations quantified by the cumulant eccentricity ratio $\varepsilon_2\{4\}/\varepsilon_2\{2\}$, we find that while the initial state fluctuations are the main source of the $v_2$ fluctuations in large collision systems, other sources like nonlinear hydrodynamic response, hadronization, and hadronic afterburner can significantly affect the $v_2$ fluctuations in small systems., Comment: 12 pages, 6 figures

We study the transverse momentum ($p_T$) differential anisotropic flow and flow fluctuation in Pb+Pb collisions at $\sqrt{s_{NN}}$=5.02 TeV at the LHC. A (3+1)-dimensional CLVisc hydrodynamics framework with fluctuating TRENTO (or AMPT) initial conditions is utilized to simulate the space-time evolution of the quark-gluon plasma (QGP) medium. The effects of shear viscosity and the sub-nucleon structure on anisotropic flow and flow fluctuation are analyzed. Our result shows that shear viscosity tends to suppress both flow coefficients (${v_2\{2}\}$, ${v_2\{4\}}$, ${\langle v_2\rangle}$) and flow fluctuation (${\sigma_{v_2}}$) due to its smearing effect on local density fluctuation. The flow coefficients appear to be insensitive to the sub-nucleon structure, whereas for flow fluctuation ${\sigma_{v_2}}$, it tends to be suppressed by the sub-nucleon structure in central collisions but enhanced in peripheral collisions. After taking into account the sub-nucleon structure effect, our numerical result can quantitatively describe the relative flow fluctuations (${v_2\{4\}/v_2\{2\}}$, $F({v_2})$) measured by the ALICE Collaboration at the LHC. We further investigate the effects of shear viscosity, sub-nucleon structure and initial condition model on the flow angle and flow magnitude decorrelations (${A_2^f}$, ${M_2^f}$) using the four-particle correlation method. We find that the flow decorrelation effect is typically stronger in central collisions than in peripheral collisions. The flow angle decorrelation is found to be insensitive to the shear viscosity and sub-nucleon structure, whereas the flow magnitude decorrelation shows quite different behavior when using TRENTO or AMPT initial condition model. Our study sheds light on the anisotropic flow, transport properties and initial structure of the QGP created in high-energy nuclear collisions., Comment: 11 pages, 6 figures

We report our recent study on the global and local polarization of $\Lambda$ hyperons in Au+Au collisions at RHIC-BES energies within the (3+1)-dimensional CLVisc hydrodynamics framework. We present our numerical results for the global polarization as the function of collision energies and the local polarization along the beam direction as functions of azimuthal angle in $20-50$% centrality at $\sqrt{s_{NN}}$=7.7 GeV Au+Au collision energy. We have discussed the effects of initial conditions, Spin Hall effect and baryon diffusion., Comment: 4 pages, 3 figures. Contribution to the proceedings of Quark Matter 2023 (Houston, TX, 3-9 Sep. 2023)

We investigate the hyperon polarization and its relation with the directed flow of the quark-gluon plasma (QGP) in non-central Au+Au collisions at $\sqrt{s_{\textrm{NN}}} = 27$ GeV. A modified 3-dimensional (3D) Glauber model is developed and coupled to a (3+1)-D viscous hydrodynamic evolution of the QGP. Within this framework, we obtain a satisfactory simultaneous description of the directed flow of identified particles and $\Lambda$ polarization, and show sensitivity of polarization to both the tilted geometry and the longitudinal flow profile of the QGP. A non-monotonic transverse momentum dependence of the $\Lambda$ polarization is found in our calculation, which is absent from hydrodynamic simulation using other initialization methods and can be tested by future experimental data with higher precision. The relation between the global polarization and directed flow of $\Lambda$ is explored as the longitudinal flow field or the medium deformation varies. Due to the common origin of these two observables, their combination may provide a more stringent constraint on the initial condition of the QGP., Comment: 15 pages, 12 figures, published in Phy. Rev. C

In a quasi-particle model of QCD matter at finite temperature with thermal masses for quarks and gluons from hard thermal loops, the equation of state (EOS) can be described by an effective temperature dependence of the strong coupling $g(T)$. Assuming the same effective coupling between the exchanged gluon and thermal partons, the EOS can also be related to parton energy loss.} Based on the quasi-particle linear Boltzmann transport (QLBT) model coupled to a (3+1)-dimensional viscous hydrodynamic model of the quark-gluon plasma (QGP) evolution and a hybrid fragmentation-coalescence model for heavy quark hadronization, we perform a Bayesian analysis of the experimental data on $D$ meson suppression $R_{\rm AA}$ and anisotropy $v_2$ at RHIC and the LHC. We achieve a simultaneous constraint on the QGP EOS and the heavy quark transport coefficient, both consistent with the lattice QCD results., Comment: 8 pages, 5 figures

We investigate the hydrodynamic helicity polarization of $\Lambda$ hyperons, defined as the projection of the spin polarization vector along the directions of particle momenta, at RHIC-BES energies by utilizing the relativistic (3+1)D CLVisc hydrodynamics framework with SMASH initial conditions. As opposed to local spin polarization at high energy collisions, our hydrodynamic simulations demonstrate that the helicity polarization induced by the kinetic vorticity dominates over other contributions at intermediate and low collision energies. Our findings provide an opportunity to probe the fine structure of local kinetic vorticity as a function of azimuthal angle at intermediate and low collision energies by mapping our predictions to the future measurements in experiments., Comment: 8 pages, 3 figures

We study the directed flow of identified particles in Au+Au collisions at $\sqrt{s_\text{NN}}=7.7$ to 62.4 GeV. The Glauber model is extended to include both a tilted deformation of the QGP fireball with respect to the longitudinal direction and a non-zero longitudinal flow velocity gradient in the initial state. By combining this improved initial condition with a (3+1)-dimensional viscous hydrodynamic model calculation, we obtain a satisfactory description of the transverse momentum spectra and the rapidity dependent directed flow coefficient of different hadron species. Our calculation indicates the sensitivity of the hadron directed flow, especially its splitting between protons and antiprotons, to both the initial geometry and the initial longitudinal flow velocity. Therefore, the combination of directed flow of different hadrons can provide a tight constraint on the initial condition of nuclear matter created in heavy-ion collisions. The initial condition extracted from the directed flow is further tested with the global polarization of $\Lambda$ and $\bar{\Lambda}$ within the same theoretical framework, where we obtain a reasonable description of these hyperon polarization observed at different collision energies at RHIC., Comment: 13 pages, 7 figures

We perform a systematic study on the local and global spin polarization of $\Lambda$ and $\overline{\Lambda}$ hyperons in relativistic heavy-ion collisions at beam energy scan energies via the (3+1)-dimensional CLVisc hydrodynamics model with AMPT and SMASH initial conditions. Following the quantum kinetic theory, we decompose the polarization vector as the parts induced by thermal vorticity, shear tensor and the spin Hall effect (SHE). We find that the polarization induced by SHE and the total polarization strongly depends on the initial conditions. At $7.7$GeV, SHE gives a sizeable contribution and even flips the sign of the local polarization along the beam direction for AMPT initial condition, which is not observed for SMASH initial condition. Meanwhile, the local polarization along the out-of-plane direction induced by SHE with AMPT initial condition does not always increase with decreasing collision energies. Next, we find that the polarization along the beam direction is sensitive to the baryon diffusion coefficient, but the local polarization along the out-of-plane direction is not. Our results for the global polarization of $\Lambda$ and $\overline{\Lambda}$ agree well with the STAR data. Interestingly, the global polarization of $\overline{\Lambda}$ is not always larger than that of $\Lambda$ due to various competing effects. Our findings are helpful for understanding the polarization phenomenon and the detailed structure of quark-gluon plasma in relativistic heavy-ion collisions., Comment: 13 pages, 8 figures

Heavy quarks are valuable probes of the electromagnetic field and the initial condition of the quark-gluon plasma (QGP) matter produced in high-energy nuclear collisions. Within an improved Langevin model that is coupled to a (3+1)-dimensional viscous hydrodynamic model, we explore the origin of the directed flow coefficient ($v_1$) of heavy mesons and their decay leptons, and its splitting ($\Delta v_{1}$) between opposite charges. We find that while the rapidity dependence of the heavy flavor $v_1$ is mainly driven by the titled energy density profile of the QGP with respect to the longitudinal direction at the RHIC energy, it is dominated by the electromagnetic field at the LHC energy. The $\Delta v_{1}$ serves as a novel probe of the spacetime evolution profile of the electromagnetic field. Our results of $D$ mesons and their decay electrons are consistent with the available data at RHIC and LHC, and our predictions on the heavy flavor decay muons can be further tested by future measurements., Comment: 12 pages, 10 figures, 1 table. Published: Phys. Rev. C 105, 054907 (2022)

Heavy-ion collisions at varying collision energies provide access to different regions of the QCD phase diagram. In particular collisions at intermediate energies are promising candidates to experimentally identify the postulated first order phase transition and critical end point. While heavy-ion collisions at low and high collision energies are theoretically well described by transport approaches and hydrodynamics+transport hybrid approaches, respectively, intermediate energy collisions remain a challenge. In this work, a modular hybrid approach, the SMASH-vHLLE-hybrid coupling 3+1D viscous hydrodynamics (vHLLE) to hadronic transport (SMASH), is introduced. It is validated and subsequently applied in Au+Au/Pb+Pb collisions between $\sqrt{s_\mathrm{NN}}$ = 4.3 GeV and $\sqrt{s_\mathrm{NN}}$ = 200.0 GeV to study the rapidity and transverse mass distributions of identified particles as well as excitation functions for $\mathrm{dN}/\mathrm{d}y|_{y = 0}$ and $\langle p_\mathrm{T} \rangle$. A good agreement with experimental measurements is obtained, including the baryon stopping dynamics. The transition from a Gaussian rapidity spectrum of protons at lower energies to the double-hump structure at high energies is reproduced. The centrality and energy dependence of charged particle $v_2$ is also described reasonably well. This work serves as a basis for further studies, e.g. systematic investigations of different equations of state or transport coefficients.

We study the origin of the directed flow of charged particles produced in relativistic heavy-ion collisions. Three different initial conditions, Boz$\dot{\textrm{e}}$k-Wyskiel, CCNU and Shen-Alzhrani, of energy density distributions are coupled to the (3+1)-dimensional viscous hydrodynamic model CLVisc, and their effects on the development of the anisotropic medium geometry, pressure gradient and radial flow are systematically compared. By comparing to experimental data at both RHIC and LHC, we find that the directed flow provides a unique constraint on the tilt of the initial medium profile in the plane spanned by the impact parameter and space-time rapidity. Within mid-rapidity, the counter-clockwise tilt is shown to be a crucial source of the positive/negative force by the pressure gradient along the impact parameter ($x$) direction at backward/forward rapidity, which drives a negative slope of the $x$ component of the medium flow velocity with respect to rapidity, and in the end the same feature of the charged particle directed flow., Comment: 12 pages, 10 figures, any comments or suggestions are welcome

We perform the first study on asymmetric longitudinal decorrelations of elliptic, triangular and quadrangular flows in proton-nucleus collisions at the LHC and RHIC energies. To measure the longitudinal flow decorrelations for asymmetric collision systems, we propose a new set of rapidity-asymmetric flow decorrelation functions. Our event-by-event hydrodynamic calculations show that the flow decorrelations in proton-going direction are larger than those in nucleus-going direction. We also find that proton-nucleus collisions at RHIC have larger longitudinal flow decorrelation effects than those at the LHC. Our study opens a new window to probe the longitudinal properties and the origin of flows in relativistic nuclear collisions., Comment: 6 pages, 3 figures

We explore the system size dependence of heavy-quark-QGP interaction by studying the heavy flavor meson suppression and elliptic flow in Pb-Pb, Xe-Xe, Ar-Ar and O-O collisions at the LHC. The space-time evolution of the QGP is simulated using a (3+1)-dimensional viscous hydrodynamic model, while the heavy-quark-QGP interaction is described by an improved Langevin approach that includes both collisional and radiative energy loss inside a thermal medium. Within this framework, we provides a reasonable description of the $D$ meson suppression and flow coefficients in Pb-Pb collisions, as well as predictions for both $D$ and $B$ meson observables in other collision systems yet to be measured. We find a clear hierarchy for the heavy meson suppression with respect to the size of the colliding nuclei, while their elliptic flow coefficient relies on both the system size and the geometric anisotropy of the QGP. Sizable suppression and flow are predicted for both $D$ and $B$ mesons in O-O collisions, which serve as a crucial bridge of jet quenching between large and small collision systems. Scaling behaviors between different collision systems are shown for heavy meson suppression factor and the bulk-eccentricity-rescaled heavy meson elliptic flow as functions of the number of participant nucleons in heavy-ion collisions., Comment: 10 pages, 9 figures

We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor $R_\mathrm{AA}$ and the elliptic flow $v_2$ of $D$ mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the $D$ meson $R_\mathrm{AA}$ and $v_2$, we extract the heavy quark transport parameter $\hat{q}$ and diffusion coefficient $D_\mathrm{s}$ in the temperature range of $1-4~T_\mathrm{c}$, and compare them with the lattice QCD results and other phenomenological studies., Comment: 13 pages, 9 figures

To study the bulk properties of the quark-gluon-plasma (QGP) produced at the beam energy scan (BES) energies at the Relativistic Heavy Ion Collider (RHIC), we extend the (3+1)-dimensional viscous hydrodynamics CLVisc to include net baryon number conservation and Israel-Stewart-like equation for baryon diffusion with the NEOS-BQS equation of state, fluctuating initial conditions from Monte-Carlo Glauber model, and the afterburner SMASH. This integrated framework is shown to provide a good description of identified particle spectra, mean transverse momenta and anisotropic flows for different centralities and over a wide range of collision energies (7.7-62.4 GeV). It is found that the mean momenta of identified particles and anisotropic flows increases mildly with the collision energy due to larger radial flow. We further compute the multiple-particle cumulant ratio $v_2\{4\}/v_2\{2\}$ of elliptic flow across BES energies, and find that the relative fluctuations of elliptic flow are insensitive to the collision energy, consistent with the preliminary STAR data. Our model provides a benchmark for understanding the RHIC-BES data and studying the critical properties and phase structure of hot and dense QCD matter., Comment: 10 pages, 6 figures

Recent experiments have observed large anisotropic collective flows in high multiplicity proton-lead collisions at the Large Hadron Collider (LHC), which indicates the possible formation of mini quark-gluon plasma (QGP) in small collision systems. However, no jet quenching has been confirmed in such small systems so far. To understand this intriguing result, the system size scan experiments have been proposed to bridge the gap between large and small systems. In this work, we perform a systematic study on both heavy and light flavor jet quenching in different collision systems at the LHC energies. Using our state-of-the-art jet quenching model, which combines the next-to-leading-order perturbative QCD framework, a linear Boltzmann transport model and the (3+1)-dimensional viscous hydrodynamics simulation, we provide a good description of nuclear modification factor $R_{\rm AA}$ for charged hadrons and $D$ mesons in central and mid-central Pb+Pb and Xe+Xe collisions measured by CMS collaboration. We further predict the transverse momentum and centrality dependences of $R_{AA}$ for charged hadrons, $D$ and $B$ mesons in Pb+Pb, Xe+Xe, Ar+Ar and O+O collisions at the LHC energies. Our numerical results show a clear system size dependence for both light and heavy flavor hadron $R_{AA}$ across different collision systems. Sizable jet quenching effect is obtained for both heavy and light flavor hadrons in central O+O collisions at the LHC energies. Our study provides a significant bridge for jet quenching from large to small systems, and should be helpful for finding the smallest QGP droplet and the disappearance of QGP in relativistic nuclear collisions., Comment: 12 pages, 8 figures

The charged-particle's final state spectrum is derived from an analytic perturbative solution for the relativistic viscous hydrodynamics. By taking into account the longitudinal acceleration effect in relativistic viscous hydrodynamics, the pseudorapidity spectrum describes well the nucleus-nucleus colliding systems at RHIC and LHC. Based on both the extracted longitudinal acceleration parameters $\lambda^{*}$ and a phenomenological description of the $\lambda^{*}$, the charged-particle's pseudorapidity distributions for $\sqrt{s_{NN}}$ = 5.44 TeV Xe+Xe collisions are computed from the final state expression in a limited space-time rapidity $\eta_{s}$ region., Comment: submitted to Chinese Physics C

We study the longitudinal decorrelations of elliptic, triangular and quadrangular flows in heavy-ion collisions at the LHC and RHIC energies. The event-by-event CLVisc (3+1)-dimensional hydrodynamics model, combined with the fully fluctuating AMPT initial conditions, is utilized to simulate the space-time evolution of the strongly-coupled quark-gluon plasma. Detailed analysis is performed for the longitudinal decorrelations of flow vectors, flow magnitudes and flow orientations. We find strong correlations between final-state longitudinal decorrelations of anisotropic flows and initial-state longitudinal structures and collision geometry: while the decorrelation of elliptic flow shows a non-monotonic centrality dependence due to initial elliptic geometry, typically the longitudinal flow decorrelations are larger in lower energy and less central collisions where the mean lengths of the string structure are shorter in the initial states., Comment: 4 pages, 3 figures, Proceedings of Quark Matter 2018

We perform a systematic study on the decorrelation of anisotropic flows along the pseudorapidity in relativistic heavy-ion collisions at the LHC and RHIC energies. The dynamical evolution of the QGP fireball is simulated via the CLVisc (ideal) (3+1)-dimensional hydrodynamics model, with the fully fluctuating initial condition from A-Multi-Phase-Transport (AMPT) model. Detailed analysis is performed on the longitudinal decorrelations of elliptic, triangular and quadrangular flows in terms of flow vectors, flow magnitudes and flow orientations (event planes). It is found that pure flow magnitudes have smaller longitudinal decorrelation than pure flow orientations, and the decorrelation of flow vectors is a combined effect of both flow magnitudes and orientations. The longitudinal decorrelation of elliptic flow has a strong and non-monotonic centrality dependence due to the initial elliptic collision geometry: smallest decorrelation in mid-central collisions. In contrast, the decorrelations of triangular and quadrangular flows have weak centrality dependence, slightly larger decorrelations in more peripheral collisions. Our numerical results for Pb+Pb collisions at the LHC are in good agreement with the ATLAS data, while our RHIC results predict much larger longitudinal decorrelations as compared to the LHC. We further analyze the longitudinal structures of the AMPT initial conditions and find that the final-state longitudinal decorrelation effects are strongly correlated with the lengths of the initial string structures in the AMPT model. The decorrelation effects are typically larger at lower collision energies and in more peripheral collisions due to shorter lengths of the string structures in the initial states., Comment: 14 pages, 14 figures