Sensitivity of Heavy-Quark Dipolar Flow to its Initial Spatial Distributions in Cu+Au Collisions
Pith reviewed 2026-07-01 16:28 UTC · model grok-4.3
The pith
Heavy-quark directed flow in Cu+Au collisions is an order of magnitude larger than light-hadron flow and sensitive to initial positions plus drag coefficient.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In Cu+Au collisions at top RHIC energy the intrinsic asymmetry of the colliding nuclei produces a spatially lopsided initial energy-density profile that generates a dipolar flow structure at midrapidity. Charm quarks propagating through this medium via Langevin dynamics acquire a finite directed flow v1 whose pT-integrated value is approximately an order of magnitude larger than that of charged hadrons; the pT-differential v1 exhibits strong sensitivity to the assumed initial spatial distribution of the heavy quarks and to the temperature dependence of the drag coefficient.
What carries the argument
Langevin dynamics for charm quarks embedded inside a realistic hydrodynamic background, with the temperature-dependent drag coefficient acting as the principal medium-interaction term that converts the initial dipolar asymmetry into final-state directed flow v1.
If this is right
- Precision measurements of heavy-flavor directed flow can constrain the temperature-dependent drag coefficient and thereby improve Langevin descriptions of heavy-quark transport.
- The strong dependence on initial heavy-quark positions demonstrates that pre-equilibrium dynamics must be included to predict final-state anisotropies accurately.
- The same framework can be used to compare predictions across different asymmetric collision systems once the initial spatial distributions are specified.
- Directed flow of heavy quarks supplies an independent handle on medium interactions that is complementary to elliptic flow or nuclear modification factor observables.
Where Pith is reading between the lines
- If the reported sensitivity persists, heavy-quark v1 could serve as a cross-check on initial-state models that are currently tuned mainly on light-hadron data.
- A natural next step would be to repeat the calculation for bottom quarks to test whether the order-of-magnitude enhancement and the sensitivity pattern survive the increase in mass.
- Comparison of the predicted v1 with existing data from symmetric collisions could quantify how much of the observed directed flow is driven by geometry versus fluctuations.
Load-bearing premise
The hydrodynamic background is assumed to be realistic and the Langevin dynamics with a temperature-dependent drag coefficient is assumed to capture the dominant medium interactions without significant missing contributions from coalescence, hadronization, or pre-equilibrium evolution beyond the initial heavy-quark spatial distribution.
What would settle it
An experimental measurement in which the pT-integrated heavy-quark v1 is not approximately an order of magnitude larger than the v1 of charged hadrons, or in which the pT-differential v1 shows no appreciable change when different initial spatial distributions for the heavy quarks are assumed.
Figures
read the original abstract
We investigate charm-quark dynamics in asymmetric Cu+Au collisions at top RHIC energy using a Langevin approach embedded in a realistic hydrodynamic background. The intrinsic asymmetry of the colliding nuclei leads to a spatially lopsided initial energy-density profile, which generates a dipolar flow structure in the transverse plane even at midrapidity. As charm quarks propagate through this medium, they acquire a finite directed flow, $v_1$. We find that the $p_T$-integrated heavy-quark $v_1$ is approximately an order of magnitude larger than that of charged hadrons. In addition, the $p_T$-differential $v_1$ exhibits strong sensitivity to the initial spatial distribution of heavy quarks, emphasizing the importance of pre-equilibrium dynamics in determining final-state anisotropies. Beyond geometric effects, $v_1$ also provides direct sensitivity to medium interactions through the temperature-dependent drag coefficient. Its pronounced dependence on this transport input indicates that precision measurements of heavy-flavor directed flow could place meaningful constraints on heavy-quark transport coefficients, thereby improving Langevin-based descriptions and predictive power for heavy-flavor observables in heavy-ion collisions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates charm-quark directed flow (v1) in asymmetric Cu+Au collisions at top RHIC energy using a Langevin approach embedded in a realistic hydrodynamic background. It reports that the p_T-integrated heavy-quark v1 is approximately an order of magnitude larger than that of charged hadrons and that the p_T-differential v1 exhibits strong sensitivity to the initial spatial distribution of heavy quarks as well as to the temperature-dependent drag coefficient, with implications for constraining heavy-quark transport coefficients.
Significance. If the quantitative results hold after addressing model assumptions, the work would be significant for highlighting heavy-quark v1 in asymmetric collisions as a probe of initial geometry, pre-equilibrium dynamics, and medium interactions. The focus on sensitivity to initial heavy-quark placement and the drag coefficient could strengthen constraints on transport inputs in Langevin-based models, provided the hydrodynamic background and evolution are shown to be robust.
major comments (2)
- [Abstract] Abstract: The central claim that p_T-integrated heavy-quark v1 is an order of magnitude larger than charged-hadron v1, along with the stated sensitivities of p_T-differential v1, rests on the assumption that Langevin dynamics with the given temperature-dependent drag fully determines the final v1. The text does not quantify how coalescence at hadronization or additional pre-equilibrium effects beyond initial placement would modify these results; if such channels contribute at the 20-30% level, both the magnitude comparison and the sensitivity claims would be diluted.
- [Model and results sections] Model and results sections: The drag coefficient is described as temperature-dependent and is a standard fitted input. Without explicit demonstration that the reported v1 values do not reduce to the fit itself (e.g., via variation studies or comparison to a constant-drag baseline), the claimed direct sensitivity to this transport input cannot be assessed as independent of the model tuning.
minor comments (1)
- [Abstract] The abstract refers to a 'realistic hydrodynamic background' without specifying the particular hydrodynamic model, initial conditions, or viscosity parameters used; adding these details would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each major point below, providing clarifications on the scope of our Langevin-based study of charm-quark directed flow in Cu+Au collisions and indicating revisions where they strengthen the presentation without altering the core findings.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central claim that p_T-integrated heavy-quark v1 is an order of magnitude larger than charged-hadron v1, along with the stated sensitivities of p_T-differential v1, rests on the assumption that Langevin dynamics with the given temperature-dependent drag fully determines the final v1. The text does not quantify how coalescence at hadronization or additional pre-equilibrium effects beyond initial placement would modify these results; if such channels contribute at the 20-30% level, both the magnitude comparison and the sensitivity claims would be diluted.
Authors: Our study computes the directed flow acquired by charm quarks during their propagation through the hydrodynamic medium using the Langevin equation; the reported v1 therefore refers to the heavy quarks themselves prior to hadronization. The order-of-magnitude enhancement relative to charged-hadron v1 arises directly from the asymmetric initial geometry and is a feature of the quark-level dynamics. We vary the initial spatial distribution of heavy quarks precisely to probe sensitivity to pre-equilibrium placement, which is the dominant pre-equilibrium effect under consideration. Coalescence is outside the present scope, as the manuscript focuses on the quark transport stage; a quantitative assessment of its impact would require a separate hadronization model. We have added a clarifying sentence in the conclusions to state this limitation explicitly. revision: partial
-
Referee: [Model and results sections] Model and results sections: The drag coefficient is described as temperature-dependent and is a standard fitted input. Without explicit demonstration that the reported v1 values do not reduce to the fit itself (e.g., via variation studies or comparison to a constant-drag baseline), the claimed direct sensitivity to this transport input cannot be assessed as independent of the model tuning.
Authors: The temperature-dependent drag is indeed a standard parametrization, but the sensitivity we report stems from its interplay with the space-time temperature profile of the evolving medium in asymmetric collisions. To make this explicit, we have added a new figure and accompanying text comparing results obtained with the temperature-dependent drag against a constant-drag baseline (using the same average value). The comparison shows that the p_T-differential v1 shape and magnitude differ in a manner not reproducible by a simple rescaling of the constant case, confirming that the temperature dependence introduces genuine additional sensitivity beyond the overall normalization of the fit. revision: yes
Circularity Check
No significant circularity; forward model predictions are independent of inputs
full rationale
The paper computes heavy-quark v1 via Langevin evolution in a fixed hydrodynamic background, reporting its magnitude and sensitivity to initial heavy-quark spatial distributions plus the temperature-dependent drag coefficient. These outputs are generated by propagating the model equations forward from the stated inputs; no quoted equations, self-citations, or reductions show v1 being redefined as the drag fit itself or forced by construction. The derivation chain remains self-contained as a standard sensitivity study without load-bearing self-referential steps.
Axiom & Free-Parameter Ledger
free parameters (1)
- temperature-dependent drag coefficient
axioms (1)
- domain assumption Hydrodynamic background accurately represents the medium evolution in Cu+Au collisions
Reference graph
Works this paper leans on
-
[1]
H. Song, S. A. Bass, U. Heinz, T. Hirano, and C. Shen, Phys. Rev. Lett.106, 192301 (2011), [Erratum: Phys.Rev.Lett. 109, 139904 (2012)], arXiv:1011.2783 6 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[2]
Heavy Ion Collisions: The Big Picture, and the Big Questions
W. Busza, K. Rajagopal, and W. van der Schee, Ann. Rev. Nucl. Part. Sci.68, 339 (2018), arXiv:1802.04801 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[3]
P. F. Kolb, J. Sollfrank, and U. W. Heinz, Phys. Rev. C 62, 054909 (2000), arXiv:hep-ph/0006129
work page internal anchor Pith review Pith/arXiv arXiv 2000
-
[4]
Collective flow and viscosity in relativistic heavy-ion collisions
U. Heinz and R. Snellings, Ann. Rev. Nucl. Part. Sci.63, 123 (2013), arXiv:1301.2826 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[5]
R. Rapp, D. Blaschke, and P. Crochet, Prog. Part. Nucl. Phys.65, 209 (2010), arXiv:0807.2470 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[6]
M. He, R. J. Fries, and R. Rapp, Phys. Lett. B735, 445 (2014), arXiv:1401.3817 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[7]
A. Andronic, P. Braun-Munzinger, M. K. K¨ ohler, A. Mazeliauskas, K. Redlich, J. Stachel, and V. Vislavi- cius, JHEP07, 035, arXiv:2104.12754 [hep-ph]
-
[8]
Extraction of Heavy-Flavor Transport Coefficients in QCD Matter
A. Beraudoet al., Nucl. Phys. A979, 21 (2018), arXiv:1803.03824 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[9]
S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C88, 044907 (2013), arXiv:1308.0617 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[10]
Towards the determination of heavy-quark transport coefficients in quark-gluon plasma
S. Caoet al., Phys. Rev. C99, 054907 (2019), arXiv:1809.07894 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[11]
F. Scardina, S. K. Das, V. Minissale, S. Plumari, and V. Greco, Phys. Rev. C96, 044905 (2017), arXiv:1707.05452 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[12]
S. K. Das, J.-e. Alam, and P. Mohanty, Phys. Rev. C80, 054916 (2009), arXiv:0908.4194 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[13]
G. D. Moore and D. Teaney, Phys. Rev. C71, 064904 (2005), arXiv:hep-ph/0412346
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[14]
S. L. Olsen, T. Skwarnicki, and D. Zieminska, Rev. Mod. Phys.90, 015003 (2018), arXiv:1708.04012 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [15]
-
[16]
J. Adamet al.(STAR), Phys. Rev. C99, 034908 (2019), arXiv:1812.10224 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[17]
S. K. Das, F. Scardina, S. Plumari, and V. Greco, Phys. Lett. B747, 260 (2015), arXiv:1502.03757 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[18]
W. M. Alberico, A. Beraudo, A. De Pace, A. Molinari, M. Monteno, M. Nardi, F. Prino, and M. Sitta, Eur. Phys. J. C73, 2481 (2013), arXiv:1305.7421 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[19]
S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C92, 024907 (2015), arXiv:1505.01413 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[20]
M. He, R. J. Fries, and R. Rapp, Phys. Rev. C86, 014903 (2012), arXiv:1106.6006 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[21]
Directed and elliptic flow of charged pions and protons in Pb+Pb collisions at 40 and 158A GeV
C. Altet al.(NA49), Phys. Rev. C68, 034903 (2003), arXiv:nucl-ex/0303001
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[22]
C. Shen and S. Alzhrani, Phys. Rev. C102, 014909 (2020), arXiv:2003.05852 [nucl-th]
-
[23]
Magnetohydrodynamics, charged currents and directed flow in heavy ion collisions
U. Gursoy, D. Kharzeev, and K. Rajagopal, Phys. Rev. C89, 054905 (2014), arXiv:1401.3805 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[24]
Beam-Energy Dependence of Directed Flow of Protons, Antiprotons and Pions in Au+Au Collisions
L. Adamczyket al.(STAR), Phys. Rev. Lett.112, 162301 (2014), arXiv:1401.3043 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [25]
-
[26]
S. K. Das, S. Plumari, S. Chatterjee, J. Alam, F. Scar- dina, and V. Greco, Phys. Lett. B768, 260 (2017), arXiv:1608.02231 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [27]
-
[28]
S. Chatterjee and P. Bozek, Phys. Lett. B798, 134955 (2019), arXiv:1804.04893 [nucl-th]
- [29]
-
[30]
S. Chatterjee and P. Bo˙ zek, Phys. Rev. Lett.120, 192301 (2018), arXiv:1712.01189 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[31]
S. K. Das, M. Ruggieri, F. Scardina, S. Plumari, and V. Greco, J. Phys. G44, 095102 (2017), arXiv:1701.05123 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [32]
-
[33]
S. K. Das, M. Ruggieri, S. Mazumder, V. Greco, and J.- e. Alam, J. Phys. G42, 095108 (2015), arXiv:1501.07521 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2015
- [34]
-
[35]
Elliptic Flow from Non-equilibrium Initial Condition with a Saturation Scale
M. Ruggieri, F. Scardina, S. Plumari, and V. Greco, Phys. Lett. B727, 177 (2013), arXiv:1303.3178 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[36]
Directed flow in ultrarelativistic heavy-ion collisions
P. Bozek and I. Wyskiel, Phys. Rev. C81, 054902 (2010), arXiv:1002.4999 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[37]
Directed flow at midrapidity in heavy-ion collisions
M. Luzum and J.-Y. Ollitrault, Phys. Rev. Lett.106, 102301 (2011), arXiv:1011.6361 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[38]
Charge-dependent directed flow in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV
L. Adamczyket al.(STAR), Phys. Rev. Lett.118, 012301 (2017), arXiv:1608.04100 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [39]
-
[40]
J. Jia, J. Phys. G41, 124003 (2014), arXiv:1407.6057 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[41]
J. Y. Zhang, Y. J. Song, L. Ma, W. H. Ma, and H. Z. Huang, Phys. Rev. C112, 044901 (2025)
2025
-
[42]
T. Parida and S. Chatterjee, Phys. Rev. C106, 044907 (2022), arXiv:2204.02345 [nucl-th]
-
[43]
Charge-dependent directed flow in asymmetric nuclear collisions
V. Voronyuk, V. D. Toneev, S. A. Voloshin, and W. Cass- ing, Phys. Rev. C90, 064903 (2014), arXiv:1410.1402 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[44]
A. Adareet al.(PHENIX), Phys. Rev. Lett.121, 222301 (2018), arXiv:1807.11928 [nucl-ex]
-
[45]
K. Nakamura, T. Miyoshi, C. Nonaka, and H. R. Takahashi, Phys. Rev. C107, 014901 (2023), arXiv:2209.00323 [nucl-th]
-
[46]
Event-by-event viscous hydrodynamics for Cu-Au collisions at 200GeV
P. Bo˙ zek, Phys. Lett. B717, 287 (2012), arXiv:1208.1887 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[47]
A. Adareet al.(PHENIX), Phys. Rev. C94, 054910 (2016), arXiv:1509.07784 [nucl-ex]
-
[48]
Triangularity and Dipole Asymmetry in Heavy Ion Collisions
D. Teaney and L. Yan, Phys. Rev. C83, 064904 (2011), arXiv:1010.1876 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[49]
F. G. Gardim, F. Grassi, Y. Hama, M. Luzum, and J.-Y. Ollitrault, Phys. Rev. C83, 064901 (2011), arXiv:1103.4605 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[50]
K. Shafi and S. Chatterjee, Eur. Phys. J. C86, 93 (2026), arXiv:2505.11713 [nucl-th]
- [51]
- [52]
-
[53]
Heavy-Quark Probes of the Quark-Gluon Plasma at RHIC
H. van Hees, V. Greco, and R. Rapp, Phys. Rev. C73, 034913 (2006), arXiv:nucl-th/0508055
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[54]
Thermalization of charm quarks in infinite and finite QGP matter
S. Cao and S. A. Bass, Phys. Rev. C84, 064902 (2011), arXiv:1108.5101 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[55]
M. He, H. van Hees, P. B. Gossiaux, R. J. Fries, and R. Rapp, Phys. Rev. E88, 032138 (2013), arXiv:1305.1425 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[56]
QCD Predictions for Charm and Bottom Production at RHIC
M. Cacciari, P. Nason, and R. Vogt, Phys. Rev. Lett.95, 7 122001 (2005), arXiv:hep-ph/0502203
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[57]
Theoretical predictions for charm and bottom production at the LHC
M. Cacciari, S. Frixione, N. Houdeau, M. L. Mangano, P. Nason, and G. Ridolfi, JHEP10, 137, arXiv:1205.6344 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv
-
[58]
3+1D hydrodynamic simulation of relativistic heavy-ion collisions
B. Schenke, S. Jeon, and C. Gale, Phys. Rev. C82, 014903 (2010), arXiv:1004.1408 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[59]
Elliptic and triangular flow in event-by-event (3+1)D viscous hydrodynamics
B. Schenke, S. Jeon, and C. Gale, Phys. Rev. Lett.106, 042301 (2011), arXiv:1009.3244 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[60]
The production of photons in relativistic heavy-ion collisions
J.-F. Paquet, C. Shen, G. S. Denicol, M. Luzum, B. Schenke, S. Jeon, and C. Gale, Phys. Rev. C93, 044906 (2016), arXiv:1509.06738 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[61]
Directed Flow of Identified Particles in Au + Au Collisions at $\sqrtsNN = 200$ GeV at RHIC
L. Adamczyket al.(STAR), Phys. Rev. Lett.108, 202301 (2012), arXiv:1112.3930 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2012
- [62]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.