Source author record

Jiangyong Jia

Jiangyong Jia appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2026arXiv

Collision energy and system size dependence of longitudinal flow decorrelation in heavy-ion collisions at RHIC energies

In heavy-ion collisions, the initial collision geometry and its fluctuations drive the collective expansion of final-state hadrons in the transverse plane. However, longitudinal fluctuations induce event-plane twist and flow magnitude asymmetries, collectively known as longitudinal flow decorrelation. Using a multi-phase transport (AMPT) model, we systematically investigate the dependence of collision energy and system size of this phenomenon with Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 19.6, 27, 54.4, 200 GeV and isobar collisions (Zr+Zr and Ru+Ru) at $\sqrt{s_{\mathrm{NN}}}$ = 200 GeV. The results reveal two distinct decorrelation components: $r_n(η)$, which includes flow magnitude asymmetry and event-plane twist, and $R_n(η)$ which arises purely from event-plane twist. Both $r_n(η)$ and $R_n(η)$ decrease linearly with $η$ and exhibit a significant dependence on collision energy and the size of the system. Through the slope parameters $F_n$ in the linear parametrization $r_n(η) = 1-2F_nη$, we can quantify the strength of decorrelation. We further observe that both $F_2$ and $F_3$ demonstrate a pronounced power-law scaling behavior with collision energy, following the relation $F_n \propto log \sqrt{s_{NN}}$. These results provide valuable insights into the three-dimensional modeling of the initial stage and the evolution of relativistic heavy-ion collisions.

preprint2022arXiv

Evidence of quadrupole and octupole deformations in $^{96}$Zr+$^{96}$Zr and $^{96}$Ru+$^{96}$Ru collisions at ultra-relativistic energies

In the hydrodynamic model description of heavy ion collisions, the elliptic flow $v_2$ and triangular flow $v_3$ are sensitive to the quadrupole deformation $β_2$ and octupole deformation $β_3$ of the colliding nuclei. The relations between $v_n$ and $β_n$ have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration have just published precision $v_n$ data from isobaric $^{96}$Ru+$^{96}$Ru and $^{96}$Zr+$^{96}$Zr collisions, where they observe large differences in central collisions $v_{2,\mathrm{Ru}}>v_{2,\mathrm{Zr}}$ and $v_{3,\mathrm{Ru}}<v_{3,\mathrm{Zr}}$. Using a transport model simulation, we show that these orderings are a natural consequence of $β_{2,\mathrm{Ru}}\ggβ_{2,\mathrm{Zr}}$ and $β_{3,\mathrm{Ru}}\llβ_{3,\mathrm{Zr}}$. We reproduce the centrality dependence of the $v_2$ ratio qualitatively and $v_3$ ratio quantitatively, and extract values of $β_2$ and $β_3$ that are consistent with those measured at low energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of $^{96}$Zr in heavy ion collisions. Our analysis demonstrates that flow measurements in high-energy heavy ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

preprint2022arXiv

Higher-order transverse momentum fluctuations in heavy-ion collisions

In relativistic heavy-ion collisions, the event-by-event mean transverse momentum fluctuations are sensitive to overlap area and energy density fluctuations in initial state. We present a framework to calculate $p_{\mathrm{T}}$ fluctuations up to $4^{\mathrm{th}}$-order using standard and subevent methods, which is validated using the HIJING model. We observe a power-law dependence for cumulants of all orders as a function of charged particle multiplicity $N_{\mathrm{ch}}$, consistent with a simple independent source picture. The fluctuation in $pp$ collisions is observed to be larger than for $p+$Pb, Pb+Pb and Xe+Xe collisions at the same $N_{\mathrm{ch}}$ due to bias in number of contributing sources. The short-range correlations are greatly suppressed in subevent method in comparison to calculations based on the standard method. This study provides a baseline for transverse momentum fluctuations without the presence of final state effects.

preprint2022arXiv

Impact of event activity variable on the ratio of observables in isobar collisions

The STAR isobar data of $^{96}$Ru+$^{96}$Ru and $^{96}$Zr+$^{96}$Zr collisions at $\sqrt{s_{\mathrm{NN}}}= 200$ GeV show that ratios of observables ($R_{\mathcal{O}}$) such as the multiplicity distribution, $p(N_{\mathrm{ch}})$, and the harmonic flow, $v_n$, deviate from unity, when presented as a function of centrality, $c$. These deviations have been attributed to the differences in the shape and radial profiles between $^{96}$Ru and $^{96}$Zr nuclei. In addition, the ratios $R_{\mathcal{O}}(x)$ depend on the choice of the event activity variable $x$, which could be either $N_{\mathrm{ch}}$ or centrality. We estimate the difference $ΔR$ between these two choices, based on the published $p(N_{\mathrm{ch}})$, as well as those from a multiphase transport (AMPT) model with varied nuclear structure parameters: nuclear radius ($R_0$), surface diffuseness ($a_0$), quadrupole deformation ($β_2$), and octupole deformation ($β_3$). In contrary to $R_{v_n}(c)$, $R_{v_n}(N_{\mathrm{ch}})$ is nearly independent of the analysis approaches, suggesting that nonflow effects are better controlled by $N_{\mathrm{ch}}$ than $c$. The ratios of observables sensitive to the chiral magnetic effect (CME) are also much closer to unity for $x=N_{\mathrm{ch}}$ than $x=c$, indicating that the ratios calculated at the same $N_{\mathrm{ch}}$ provide a better baseline for the non-CME background. According to the AMPT results, the dominant parameter for $ΔR$ is $a_0$, while $R_0$ and $β_n$ are only important in central collisions. The published $p(N_{\mathrm{ch}})$ is also used to estimate $ΔR_{\left\langle p_{\mathrm{T}}\right\rangle}$ for mean transverse momentum, which is non-negligible compared with $R_{\left\langle p_{\mathrm{T}}\right\rangle}-1$.

preprint2022arXiv

Non-flow effects in correlation between harmonic flow and transverse momentum in nuclear collisions

A large anti-correlation signal between elliptic flow $v_2$ and average transverse momentum $[p_{\mathrm{T}}]$ was recently measured in small collision systems, consistent with a final-state hydrodynamic response to the initial geometry. This negative $v_2$-$[p_{\mathrm{T}}]$ correlation was predicted to change to positive correlation for events with very small charged particle multiplicity $N_{\mathrm{ch}}$ due to initial-state momentum anisotropies of the gluon saturation effects. However, the role of non-flow correlations is expected to be important in these systems, which is not yet studied. We estimate the non-flow effects in $pp$, $p$Pb and peripheral PbPb collisions using {\tt Pythia} and {\tt Hijing} models, and compare them with the experimental data. We show that the non-flow effects are largely suppressed using the rapidity-separated subevent cumulant method (details of the cumulant framework are also provided). The magnitude of the residual non-flow is much less than the experimental observation in the higher $N_{\mathrm{ch}}$ region, supporting the final-state response interpretation. In the very low $N_{\mathrm{ch}}$ region, however, the sign and magnitude of the residual non-flow depend on the model details. Therefore, it is unclear at this moment whether the sign change of $v_2$-$[p_{\mathrm{T}}]$ can serve as evidence for initial state momentum anisotropies predicted by the gluon saturation.

preprint2022arXiv

Phenomenological study on correlation between flow harmonics and mean transverse momentum in nuclear collisions

To assess the properties of the quark-gluon plasma formed in nuclear collisions, the Pearson correlation coefficient between flow harmonics and mean transverse momentum, $ρ\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right)$, reflecting the overlapped geometry of colliding atomic nuclei, is measured. $ρ\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right)$ was found to be particularly sensitive to the quadrupole deformation of the nuclei. We study the influence of the nuclear quadrupole deformation on $ρ\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right)$ in $\rm{Au+Au}$ and $\rm{U+U}$ collisions at RHIC energy using $\rm{AMPT}$ transport model, and show that the $ρ\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right)$ is reduced by the prolate deformation $β_2$ and turns to change sign in ultra-central collisions (UCC).

preprint2022arXiv

Probing nuclear quadrupole deformation from correlation of elliptic flow and transverse momentum in heavy ion collisions

In heavy ion collisions, elliptic flow $v_2$ and radial flow, characterized by event-wise average transverse momentum $[p_{\mathrm{T}}]$, are related to the shape and size of the overlap region, which are sensitive to the shape of colliding atomic nuclei. The Pearson correlation coefficient between $v_2$ and $[p_{\mathrm{T}}]$, $ρ_2$, was found to be particularly sensitive to the quadrupole deformation parameter $β$ that is traditionally measured in low energy experiments. Built on earlier insight that the prolate deformation $β>0$ reduces the $ρ_2$ in ultra-central collisions (UCC), we show that the prolate deformation $β<0$ enhances the value of $ρ_2$. As $β>0$ and $β<0$ are the two extremes of triaxiality, the strength and sign of $v_2^2-[p_{\mathrm{T}}]$ correlation can be used to provide valuable information on the triaxiality of the nucleus. Our study provide further arguments for using the hydrodynamic flow as a precision tool to directly image the deformation of the atomic nuclei at extremely short time scale ($<10^{-24}$s).

preprint2022arXiv

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 $β_2$ 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 $\varepsilon_2$ and the density gradient $d_{\perp}$ in the initial state. Using analytical estimate and a Glauber model, I show that the variances, $\langle\varepsilon_2^2\rangle$ or $\langle(δd_{\perp}/d_{\perp})^2\rangle$, and skewnesses, $\langle\varepsilon_2^2δd_{\perp}/d_{\perp}\rangle$ or $\langle(δd_{\perp}/d_{\perp})^3\rangle$, have a simple analytical form of $a'+b'β_2^2$ and $a'+(b'+c'\cos(3γ))β_2^3$, respectively. From these, I constructed several normalized skewnesses to isolate the $γ$ dependence from that of $β_2$, and show that the correlations between any normalized skewness and any variance can constrain simultaneously the $β_2$ and $γ$. Assuming a linear relation with elliptic flow $v_2$ and mean-transverse momentum $[p_{\mathrm{T}}]$ of final state particles, similar conclusions are also expected for the variances and skewnesses of $v_2$ and $[p_{\mathrm{T}}]$. Our findings motivate a dedicated system scan of high-energy heavy ion collisions to measure triaxiality of atomic nuclei. This is better done by collisions of prolate, $\cos(3γ)=1$, and oblate nuclei, $\cos(3γ)=-1$, with well known $β_2$ values to calibrate the coefficients $b'$ and $c'$, followed by collisions of species of interest especially those with known $β_2$ but unknown $γ$. The results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

preprint2022arXiv

Shape of atomic nuclei in heavy ion collisions

In the hydrodynamic model description of heavy ion collisions, the final-state anisotropic flow $v_n$ are linearly related to the strength of the multi-pole shape of the distribution of nucleons in the transverse plane $\varepsilon_n$, $v_n\propto \varepsilon_n$. The $\varepsilon_n$, for $n=1,2,3,4$, are sensitive to the shape of the colliding ions, characterized by the quadrupole $β_2$, octupole $β_3$ and hexadecapole $β_4$ deformations. This sensitivity is investigated analytically and also in a Monte Carlo Glauber model. One observes a robust linear relation, $\langle\varepsilon_n^2\rangle = a_n'+b_n'β_n^2$, for events in a fixed centrality. The $\langle\varepsilon_1^2\rangle$ has a contribution from $β_3$ and $β_4$, and $\langle\varepsilon_3^2\rangle$ from $β_4$. In the ultra-central collisions, there are little cross contributions between $β_2$ and $\varepsilon_3$ and between $β_3$ and $\varepsilon_2$, but clear cross contributions are present in non-central collisions. Additionally, $\langle\varepsilon_n^2\rangle$ are insensitive to non-axial shape parameters such as the triaxiality. This is good news because the measurements of $v_2$, $v_3$ and $v_4$ can be used to constrain simultaneously the $β_2$, $β_3$, and $β_4$ values. This is best done by comparing two colliding ions with similar mass numbers and therefore nearly identical $a_n'$, to obtain simple equation that relates the $β_n$ of the two species. This opens up the possibility to map the shape of the atomic nuclei at a timescale ($<10^{-24}$s) much shorter than probed by low-energy nuclear structure physics ($<10^{-21}$s), which ultimately may provide information complementary to those obtained in the nuclear structure experiments.

preprint2022arXiv

The impact of nuclear deformation on relativistic heavy-ion collisions: assessing consistency in nuclear physics across energy scales

In the hydrodynamic framework of heavy-ion collisions, elliptic flow, $v_2$, is sensitive to the quadrupole deformation, $β$, of the colliding ions. This enables one to test whether the established knowledge on the low-energy structure of nuclei is consistent with collider data from high-energy experiments. We derive a formula based on generic scaling laws of hydrodynamics to relate the difference in $v_2$ measured between collision systems that are close in size to the value of $β$ of the respective species. We validate our formula in simulations of 238U+238U and 197Au+197Au collisions at top Relativistic Heavy Ion Collider (RHIC) energy, and subsequently apply it to experimental data. Using the deformation of 238U from low-energy experiments, we find that RHIC $v_2$ data implies $0.16 \lesssim |β| \lesssim 0.20$ for 197Au nuclei, i.e., significantly more deformed than reported in the literature, posing an interesting puzzle in nuclear phenomenology.

preprint2021arXiv

Investigation of Experimental Observables in Search of the Chiral Magnetic Effect in Heavy-ion Collisions in the STAR experiment

The chiral magnetic effect (CME) is a novel transport phenomenon, arising from the interplay between quantum anomalies and strong magnetic fields in chiral systems. In high-energy nuclear collisions, the CME may survive the expansion of the quark-gluon plasma fireball and be detected in experiments. Over the past decade, the experimental searches for the CME have aroused extensive interest at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). The main goal of this article is to investigate three pertinent experimental approaches: the $γ$ correlator, the $R$ correlator and the signed balance functions. We will exploit both simple Monte Carlo simulations and a realistic event generator (EBE-AVFD) to verify the equivalence in the kernel-component observables among these methods and to ascertain their sensitivities to the CME signal for the isobaric collisions at RHIC.

preprint2020arXiv

Centrality fluctuations and decorrelations in heavy-ion collisions

The centrality or the number of initial-state sources $V$ of the system produced in heavy ion collision is a concept that is not uniquely defined and subject to significant theoretical and experimental uncertainties. We argue that a more robust connection between the initial-state sources with final-state multiplicity could be established from the event-by-event multiplicity correlation between two subevents separated in pseudorapidity, $N_a$ vs $N_b$. This correlation is sensitive to two main types of centrality fluctuations (CF): 1) particle production for each source $p(n)$ which smears the relation between $V$ and $N_a$ used for experimental centrality, and 2) decorrelations between the sources in the two subevents $V_b$ and $V_a$. The CF is analyzed in terms of cumulants of $V_b$ and $N_b$ as a function of $N_a$, i.e. experimental centrality is defined with $N_a$. We found that the mean values $\langle V_b\rangle_{N_a}$ and $\langle N_b\rangle_{N_a}$ increase linearly with $N_a$ in mid-central collisions, but flatten out in ultra-central collisions. Such non-linear behavior is sensitive to the centrality resolution of $N_a$. In the presence of centrality decorrelations, the scaled variances $\langle(δV_b)^2\rangle/\langle V_b\rangle$ and $\langle(δN_b)^2\rangle/\langle N_b\rangle$ are found to decrease linearly with $N_a$ in mid-central collisions, while the $p(n)$ leads to another sharp decrease in the ultra-central region. The higher-order cumulants of $V_b$ and $N_b$ show interesting but rather complex behaviors which deserve further studies. Our results suggest that one can use the cumulants of the two-dimensional multiplicity correlation, especially the mean and variance, to constrain the particle production mechanism as well as the longitudinal fluctuations of the initial-state sources.

preprint2020arXiv

Disentangle contributions to small-system collectivity via scans of light nucleus-nucleus collisions

The observation of multi-particle azimuthal correlations in high-energy small-system collisions has led to intense debate on its origin and the possible coexistence from two competing theoretical scenarios: one based on initial-state intrinsic momentum anisotropy (ISM), and the other based on final-state collective response to the collision geometry (FSM). To complement the previous scan of asymmetric collision systems ($p$+Au, $d$+Au and He+Au), we propose a scan of small symmetric collision systems at RHIC, such as C+C, O+O, Al+Al and Ar+Ar $\sqrt{s_{\mathrm{NN}}}=0.2$ TeV, to provide further insights in disentangling contributions from these two scenarios. These symmetric small systems have the advantage of providing a better controlled initial geometry dominated by the average shape of the overlap region, as opposed to fluctuation-driven geometries in asymmetric systems. A transport model is employed to investigate the expected geometry response in the FSM scenario. Different trends of elliptic flow with increasing charge particle multiplicity are observed between symmetric and asymmetric systems, while triangular flow appears to show a similar behavior. Furthermore, a comparison of O+O collisions at $\sqrt{s_{\mathrm{NN}}}=0.2$ TeV and at $\sqrt{s_{\mathrm{NN}}}=2.76-7$ TeV, as proposed at the LHC, provides a unique opportunity to disentangle the collision geometry effects at nucleon level from those arising from subnucleon fluctuations.

preprint2020arXiv

Longitudinal eccentricity decorrelations in heavy ion collisions

In heavy-ion collisions, the harmonic flow $V_n$ of final-state particles are driven by the eccentricity vector ${\mathcal{E}}_n$ that describe the shape of the initial fireball projected in the transverse plane. It is realized recently that the structure and shape of the fireball, and consequently the ${\mathcal{E}}_n$, fluctuate in pseudorapidity $η$ in a single event, ${\mathcal{E}}_n(η)$. This leads to eccentricity decorrelation between different $η$, driving the longitudinal flow decorrelations observed in the experiments. Using a Glauber model with a paramerterized longitudinal structure, we have estimated the eccentricity decorrelations and related them to the measured flow decorrelation coefficients for elliptic flow $n=2$ and triangular flow $n=3$. We investigated the dependence of eccentricity decorrelations on the choice of collision system in terms of the size, asymmetry and deformation of the nuclei. We found that these nuclear geometry effects lead to significant and characteristic patterns on the eccentricity decorrelations, which describe the measured ratios of the flow decorrelations between Xe+Xe and Pb+Pb collisions. These patterns can be searched for using existing experimental data at RHIC and the LHC, and if confirmed, they will provide a mean to improve our understanding of the initial state of the heavy-ion collisions.

preprint2020arXiv

Robustness of principal component analysis on harmonic flow in heavy ion collisions

The principal component analysis (PCA), a mathematical tool commonly used in statistics, has recently been employed to interpret the $p_T$-dependent fluctuations of harmonic flow $v_n$ in terms of leading and subleading flow modes in heavy ion collisions. Using simulated data from AMPT and HIJING models, we show that the PCA modes are not fixed, but depend on the choice of the particle weight and the $p_T$ range. Furthermore, the shape of the leading mode is affected by the presence of non-flow correlations, and fake subleading mode may arise from the mixing of non-flow correlations with leading flow mode with a magnitude that could be larger than the genuine subleading flow mode. Therefore, the meaning of PCA modes and their relations to physical leading and subleading flow modes associated initial state eccentricities need to be further clarified/validated in realistic model simulations.

preprint2019arXiv

Medium-Induced Transverse Momentum Broadening via Forward Dijet Correlations

Dijet azimuthal angle correlation is arguably one of the most direct probes of the medium-induced broadening effects. The evidence for such broadening, however, is not yet clearly observed within the precision of current mid-rapidity measurements at RHIC and the LHC. We show that the dijet correlation in forward rapidity from the future LHC RUN3, aided by forward detector upgrades, can reveal this broadening thanks to the steeper jet spectra, suppressed vacuum radiations and lower underlying event background, with a sensitivity comparable to that of the future high-luminosity Au+Au run at RHIC. Dijet correlation measurements at the two facilities together can provide powerful constraints on the temperature dependence of medium transport properties.

preprint2016arXiv

Forward-backward multiplicity fluctuation and longitudinal harmonics in high-energy nuclear collisions

An analysis method is proposed to study the forward-backward (FB) multiplicity fluctuation in high-energy nuclear collisions, built on the earlier work of Bzdak and Teaney. The method allows the decomposition of the centrality dependence of average multiplicity from the dynamical event-by-event (EbyE) fluctuation of multiplicity in pseudorapidity. Application of the method to AMPT and HIJING models shows that the long-range component of the FB correlation is captured by a few longitudinal harmonics, with the first component driven by the asymmetry in the number of participating nucleons in the two colliding nuclei. The higher-order longitudinal harmonics are found to be strongly damped in AMPT compare to HIJING, due to weaker short-range correlations as well as the final-state effects present in the AMPT model. Two-particle pseudorapidity correlation reveals interesting charge-dependent short-range structures that are absent in HIJING model. The proposed method opens an avenue to elucidate the particle production mechanism and early time dynamics in heavy-ion collisions. Future analysis directions and prospects of using the pseudorapidity correlation function to understand the centrality bias in $p$+$p$, $p$+A and A+A collisions are discussed.

preprint2015arXiv

Limitation of multi-particle correlations for studying the event-by-event distribution of harmonic flow in heavy-ion collisions

The sensitivity of flow harmonics from cumulants on the event-by-event flow distribution $p(v_n)$ is investigated using a simple central moment expansion approach. For narrow distribution whose width is much smaller than the mean $σ_n\ll < v_n >$, the difference between the first three higher-order cumulant estimates $v_n\{4\}$, $v_n\{6\}$ and $v_n\{8\}$ are not very sensitive to the shape of $p(v_n)$. For broad distribution $σ_n\gtrsim < v_n >$, the higher-order cumulant estimates differ from each other but may change sign and become ill-defined. This sign change arises from the choice of $p(v_n)$, without the need to invoke non-flow effects. Direct extraction of $p(v_n)$ via a data-driven unfolding method is a more preferred approach for flow distribution measurement.

preprint2015arXiv

The Hot QCD White Paper: Exploring the Phases of QCD at RHIC and the LHC

The past decade has seen huge advances in experimental measurements made in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and more recently at the Large Hadron Collider (LHC). These new data, in combination with theoretical advances from calculations made in a variety of frameworks, have led to a broad and deep knowledge of the properties of thermal QCD matter. Increasingly quantitative descriptions of the quark-gluon plasma (QGP) created in these collisions have established that the QGP is a strongly coupled liquid with the lowest value of specific viscosity ever measured. However, much remains to be learned about the precise nature of the initial state from which this liquid forms, how its properties vary across its phase diagram and how, at a microscopic level, the collective properties of this liquid emerge from the interactions among the individual quarks and gluons that must be visible if the liquid is probed with sufficiently high resolution. This white paper, prepared by the Hot QCD Writing Group as part of the U.S. Long Range Plan for Nuclear Physics, reviews the recent progress in the field of hot QCD and outlines the scientific opportunities in the next decade for resolving the outstanding issues in the field.

preprint2014arXiv

A method for studying the rapidity fluctuation and decorrelation of harmonic flow in heavy-ion collisions

An "event-shape twist" technique is proposed to study the longitudinal dynamics of harmonic flow, in particular the effects of rapidity fluctuation and event-plane decorrelation. This technique can distinguish between two types of rapidity decorrelation effects: a systematic rotation versus a random fluctuation of flow angles along the rapidity direction. The technique is demonstrated and the magnitude of the two decorrelation effects is predicted using the AMPT model via a single particle analysis and two-particle correlation analysis. An observed decorrelation can be attributed to a systematic rotation of event-plane angle along the pseudorapidity, consistent with a collective response to an initial state twist of the fireball proposed by Bozek {\it et.al.}. This rotation is also observed for several higher-order harmonics with the same sign and similar magnitudes.

preprint2014arXiv

Elucidating the event-by-event flow fluctuations in heavy-ion collisions via the event shape selection technique

The presence of large event-by-event flow fluctuations in heavy ion collisions at RHIC and the LHC provides an opportunity to study a broad class of flow observables. This paper explores the correlations among harmonic flow coefficients $v_n$ and their phases $Φ_n$, and the rapidity fluctuation of $v_n$. The study is carried out usin Pb+Pb events generated by the AMPT model with fixed impact parameter. The overall ellipticity/triangularity of events is varied by selecting on the eccentricities $ε_n$ or the magnitudes of the flow vector $q_n$ for n=2 and 3, respectively. The responses of the $v_n$, the event-plane correlations, and the rapidity fluctuations, to the change in $ε_n$ and $q_n$ are then systematized. Strong positive correlations are observed among all even harmonics $v_2, v_4$, and $v_6$ (all increase with $q_2$), between $v_2$ and $v_5$ (both increase with $q_2$) and between $v_3$ and $v_5$ (both increase with $q_3$), consistent with the effects of nonlinear collective response. In contrast, an anti-correlation is observed between $v_2$ and $v_3$ similar to that seen between $ε_2$ and $ε_3$. These correlation patterns are found to be independent of whether selecting on $q_n$ or $ε_n$, validating the ability of $q_n$ in selecting the initial geometry. A forward/backward asymmetry of $v_n(η)$ is observed for events selected on $q_n$ but not on $ε_n$, reflecting dynamical fluctuations exposed by the $q_n$ selection. Many event-plane correlators show good agreement between $q_n$ and $ε_n$ selections, suggesting that their variations with $q_n$ are controlled by the change of $ε_n$ in the initial geometry. Hence these correlators may serve as promising observables for disentangling the fluctuations generated in various stages of the evolution of the matter created in heavy ion collisions.

preprint2014arXiv

Event-shape fluctuations and flow correlations in ultra-relativistic heavy-ion collisions

I review recent measurements of a large set of flow observables associated with event-shape fluctuations and collective expansion in heavy ion collisions. First, these flow observables are classified and experiment methods are introduced. The experimental results for each type of observables are then presented and compared to theoretical calculations. A coherent picture of initial condition and collective flow based on linear and non-linear hydrodynamic responses is derived, which qualitatively describe most experimental results. I discuss new types of fluctuation measurements that can further our understanding of the event-shape fluctuations and collective expansion dynamics.

preprint2014arXiv

Forward-backward eccentricity and participant-plane angle fluctuations and their influences on longitudinal dynamics of collective flow

We argue that the transverse shape of the fireball created in heavy ion collision is controlled by event-by-event fluctuations of the eccentricity vectors for the forward-going and backward-going wounded nucleons: $\vecε_n^{\mathrm{F}}\equiv ε_n^{\mathrm{F}} e^{i nΦ_n^{\mathrm{*F}}}$ and $\vecε_n^{\mathrm{B}}\equiv ε_n^{\mathrm{B}} e^{i nΦ_n^{\mathrm{*B}}}$. Due to the asymmetric energy deposition of each wounded nucleon along its direction of motion, the eccentricity vector of the produced fireball is expected to interpolate between $\vecε_n^{\mathrm{F}}$ and $\vecε_n^{\mathrm{B}}$ along the pseudorapidity, and hence exhibits sizable forward-backward(FB) asymmetry ($ε_n^{\rm B}\neqε_n^{\rm F}$) and/or FB-twist ($Φ_n^{\mathrm{*F}}\neqΦ_n^{\mathrm{*B}}$). A transport model calculation shows that these initial state longitudinal fluctuations for $n=2$ and 3 survive the collective expansion, and result in similar FB asymmetry and/or a twist in the final state event-plane angles. These novel EbyE longitudinal flow fluctuations should be accessible at RHIC and the LHC using the event-shape selection technique proposed in earlier papers. If these effects are observed experimentally, it could improve our understanding of the initial state fluctuations, particle production and collective expansion dynamics of the heavy ion collision.

preprint2013arXiv

A method for studying initial geometry fluctuations via event plane correlations in heavy ion collisions

A method is proposed to measure the relative azimuthal angle distributions involving two or more event planes of different order in heavy ion collisions using a Fourier analysis technique. The analysis procedure is demonstrated for correlations involving two and three event planes (Phi_n, Phi_m and Phi_h). The Fourier coefficients of these distributions are found to coincide with previously proposed correlators, such as cos(6Phi_2-6Phi_3) and cos(Phi_1+2Phi_2-3Phi_3) etc, hence the method provides a natural framework for studying these correlators at the same time. Using a Monte Carlo Glauber model to simulate Au+Au collisions with fluctuating initial geometry, we are able to identify several new two- or three-plane correlators that have sizable magnitudes and should be measured experimentally.

preprint2013arXiv

A study of the anisotropy associated with dipole asymmetry in heavy ion collisions

The anisotropy associated with the initial dipole asymmetry in heavy ion collisions is studied via the two-particle relative azimuthal azimuthal angle (Δϕ=ϕ^a-ϕ^b) correlations, within a multi-phase transport model. For a broad selection of centrality, transverse momenta (pT^{a,b}) and pseudorapidity (η^{a,b}), a fitting method is used to decompose the anisotropy into a rapidity-even component, characterized by the Fourier coefficient v1, and a global momentum conservation component. The extracted v1 values are negative for pT<=0.7-0.9 GeV, reach a maximum at 2-3 GeV, and decreases at higher pT. The v1 values vary weakly with ηand centrality, but increases with collision energy and parton cross-section. The extracted global momentum conservation component is found to depend on Δη= η^a-η^b for Δη<3.

preprint2013arXiv

Disentangling flow and nonflow correlations via Bayesian unfolding of the event-by-event distributions of harmonic coefficients in ultrarelativistic heavy-ion collisions

The performance of the Bayesian unfolding method in extracting the event-by-event (EbyE) distributions of harmonic flow coefficients v_n is investigated using a toy model simulation, as well as simulations based on the HIJING and AMPT models. The unfolding method is shown to recover the input v_2-v_4 distributions for multiplicities similar to those observed in Pb+Pb collisions at the LHC. The effects of the nonflow are evaluated using HIJING simulation with and without a flow afterburner. The probability distribution of v_n resulting only from nonflow in HIJING is nearly a Gaussian and can be largely suppressed in the data-driven unfolding method used by the ATLAS Collaboration. The residual nonflow effects have no appreciable impact on the v_3 distributions, but are found to affect the tails of the v_2 and v_4 distributions; these effects manifest as a small simultaneous change in the mean and standard deviation of the $v_n$ distributions. For the AMPT model, which contains both flow fluctuations and nonflow effects, the reduced shape of the extracted v_n distributions is found to be independent of pT in the low pT region, similar to what is observed in the ATLAS data. The prospect of using the EbyE distribution of the harmonic spectrum aided by the unfolding technique as a general tool to study azimuthal correlations in high energy collisions is also discussed.

preprint2012arXiv

Azimuthal anisotropy in a jet absorption model with fluctuating initial geometry in heavy ion collisions

The azimuthal anisotropy due to path-length dependent jet energy loss is studied in a simple jet absorption model that include event by event fluctuating Glauber geometry. Significant anisotropy coefficients v_n are observed for n=1,2 and 3, but they are very small for n>3. These coefficients are expected to result in a "ridge" for correlations between two independently produced jets. The correlations between the orientation of the n^{th}-order anisotropy induced by jet absorption (Φ_n^{QP}) and the n^{th}-order participant plane (Φ_n^{PP}) responsible for harmonic flow are studied. Tight correlations are observed for n=2 in mid-central collisions, but they weaken significantly for n\neq2. The correlations are positive for n>=3, but become negative in central collisions for n>3. The dispersion between Φ_n^{QP} and Φ_n^{PP} is expect to break the factorization of the Fourier coefficients from two-particle correlation v_{n,n} into the single particle v_n, and has important implications for the high-pT ridge phenomena.

preprint2012arXiv

Study on initial geometry fluctuations via participant plane correlations in heavy ion collisions: part II

Further investigation of the participant plane correlations within a Glauber model framework is presented, focusing on correlations between three or four participant planes of different order. A strong correlation is observed for $\cos(2Φ_{2}^*+3Φ_{3}^*-5Φ_{5}^*)$ which is a reflection of the elliptic shape of the overlap region. The correlation between the corresponding experimental reaction plane angles can be easily measured. Strong correlations of similar geometric origin are also observed for $\cos(2Φ_{2}^*+4Φ_{4}^*-6Φ_{6}^*)$, $\cos(2Φ_2^*-3Φ_3^*-4Φ_4^*+5Φ_5^*)$, $\cos(6Φ_2^*+3Φ_3^*-4Φ_4^*-5Φ_5^*)$, $\cos(Φ_1^*-2Φ_2^*-3Φ_3^*+4Φ_4^*)$, $\cos(Φ_1^*+6Φ_2^*-3Φ_3^*-4Φ_4^*)$, and $\cos(Φ_1^*+2Φ_2^*+3Φ_3^*-6Φ_6^*)$, which are also measurable. Experimental measurements of the corresponding reaction plane correlators in heavy ion collisions at RHIC and the LHC may improve our understanding of the physics underlying the measured higher order flow harmonics.

preprint2011arXiv

A study of the correlations between jet quenching observables at RHIC

Focusing on four types of correlation plots, $R_{\rm AA}$ vs. $v_2$, $R_{\rm AA}$ vs. $I_{\rm AA}$, $I_{\rm AA}$ vs. $v_2^{I_{\rm AA}}$ and $v_2$ vs.\ $v_2^{I_{\rm AA}}$, we demonstrate how the centrality dependence of \emph{correlations} between multiple jet quenching observables provide valuable insight into the energy loss mechanism in a quark-gluon plasma. In particular we find that a qualitative energy loss model gives a good description of $R_{\rm AA}$ vs.\ $v_2$ only when we take $ΔE\sim l^3$ and a medium geometry generated by a model of the Color Glass Condensate. This same $ΔE\sim l^3$ model also qualitatively describes the trigger $p_T$ dependence of $R_{\rm AA}$ vs.\ $I_{\rm AA}$ data and makes novel predictions for the centrality dependence for this $R_{\rm AA}$ vs.\ $I_{\rm AA}$ correlation. Current data suggests, albeit with extremely large uncertainty, that $v_2^{I_{\rm AA}}\gg v_2$, a correlation that is difficult to reproduce in current energy loss models.

preprint2010arXiv

Dissecting the role of initial collision geometry for jet quenching observables in relativistic heavy ion collisions

The observation of large azimuthal anisotropy or $v_2$ for hadrons above $p_T>5$ GeV/$c$ in Au+Au collisions at $\sqrt{s_{\rm nn}}=200$ GeV has been a longstanding challenge for jet quenching models based on perturbative QCD (pQCD). Using a simple jet absorption model, we seek to clarify the situation by exploring in detail how the calculated $v_2$ varies with choices of the collision geometry as well as choices of the path length dependence and thermalization time $τ_0$ in the energy loss formula. Besides the change of eccentricity due to distortion from gluon saturation or event-by-event fluctuation, we find that the $v_2$ is also sensitive to the centrality dependence of multiplicity and the relative size between the matter profile and the jet profile. We find that the $v_2$ calculated for the naive quadratic path length dependence of energy loss, even including eccentricity fluctuation and the gluon saturation, is not enough to describe the experimental data at high $p_T$ ($\sim$ 6 GeV/$c$) in Au+Au collisions. However, it can match the full centrality dependence of $v_2$ data if higher power path length dependence of energy loss is allowed. We also find that the calculated $v_2$ is sensitive to the assumption of the early time dynamics but generally increases with $τ_0$, opposite to what one expects for elliptic flow. This study attests to the importance of confining the initial geometry, possibly by combining jet quenching $v_2$ with elliptic flow and other jet quenching observables, for proper interpretation of the experimental data.

preprint2006arXiv

Jet properties from $π^{\pm}$ - $h^{\pm}$ correlation in $p+p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

We review recent results on the charged pion - charged hadron correlation in $p+p$ and $d$ + Au collisions as measured by the PHENIX Collaboration. Properties of di-jet system, such as the jet shape, associated hadron yield per trigger pion, and the underlying event are extracted statistically from the $π^{\pm}-h^{\pm}$ correlation function in $Δϕ$ and $Δη$. For jet triggered with high $p_T$ pions ($p_T>5$ GeV/c), no apparent differences in the jet properties are seen between $p+p$ and $d$ + Au.