pith. sign in

arxiv: 2605.26766 · v1 · pith:D4UWIXQJnew · submitted 2026-05-26 · ⚛️ nucl-th · hep-ex· hep-ph· nucl-ex

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

classification ⚛️ nucl-th hep-exhep-phnucl-ex
keywords heavy-quark directed flowCu+Au collisionsLangevin dynamicsdipolar flowhydrodynamic backgroundtemperature-dependent draginitial spatial distributiontransport coefficients
0
0 comments X

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.

The paper examines charm-quark motion through the lopsided energy-density profile created when a copper nucleus collides with a gold nucleus. This geometric asymmetry produces a dipolar flow pattern that imparts a net directed flow v1 to the heavy quarks as they traverse the expanding medium. Calculations using Langevin dynamics show that the resulting v1 is roughly ten times larger than the v1 measured for ordinary charged hadrons, and that the shape of the v1 spectrum changes markedly when the starting locations of the charm quarks or the temperature dependence of the drag force are altered. A reader should care because the result identifies a concrete observable whose magnitude and shape can be measured to test both the earliest spatial arrangement of heavy quarks and the strength of their coupling to the medium.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.26766 by Ankit Kumar Panda, Tribhuban Parida.

Figure 1
Figure 1. Figure 1: FIG. 1. Initial energy density profile for a Cu+Au collision [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Initial heavy-quark distributions along the impact [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

The central claims rest on a hydrodynamic background whose realism is asserted without independent verification in the abstract and on a temperature-dependent drag coefficient whose functional form is a standard but adjustable input in Langevin models.

free parameters (1)
  • temperature-dependent drag coefficient
    Central transport input whose temperature dependence is used to generate the reported v1 sensitivities; functional form and normalization are model inputs.
axioms (1)
  • domain assumption Hydrodynamic background accurately represents the medium evolution in Cu+Au collisions
    The Langevin propagation is performed inside this background.

pith-pipeline@v0.9.1-grok · 5743 in / 1204 out tokens · 31511 ms · 2026-07-01T16:28:04.884463+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

62 extracted references · 61 canonical work pages · 43 internal anchors

  1. [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]

  2. [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]

  3. [3]

    P. F. Kolb, J. Sollfrank, and U. W. Heinz, Phys. Rev. C 62, 054909 (2000), arXiv:hep-ph/0006129

  4. [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]

  5. [5]

    R. Rapp, D. Blaschke, and P. Crochet, Prog. Part. Nucl. Phys.65, 209 (2010), arXiv:0807.2470 [hep-ph]

  6. [6]

    M. He, R. J. Fries, and R. Rapp, Phys. Lett. B735, 445 (2014), arXiv:1401.3817 [nucl-th]

  7. [7]

    Andronic, P

    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. [8]

    Extraction of Heavy-Flavor Transport Coefficients in QCD Matter

    A. Beraudoet al., Nucl. Phys. A979, 21 (2018), arXiv:1803.03824 [nucl-th]

  9. [9]

    Heavy quark dynamics and hadronization in ultra-relativistic heavy-ion collisions: collisional versus radiative energy loss

    S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C88, 044907 (2013), arXiv:1308.0617 [nucl-th]

  10. [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]

  11. [11]

    Estimating the Charm Quark Diffusion Coefficient and thermalization time from D meson spectra at RHIC and LHC

    F. Scardina, S. K. Das, V. Minissale, S. Plumari, and V. Greco, Phys. Rev. C96, 044905 (2017), arXiv:1707.05452 [nucl-th]

  12. [12]

    S. K. Das, J.-e. Alam, and P. Mohanty, Phys. Rev. C80, 054916 (2009), arXiv:0908.4194 [nucl-th]

  13. [13]

    G. D. Moore and D. Teaney, Phys. Rev. C71, 064904 (2005), arXiv:hep-ph/0412346

  14. [14]

    S. L. Olsen, T. Skwarnicki, and D. Zieminska, Rev. Mod. Phys.90, 015003 (2018), arXiv:1708.04012 [hep-ph]

  15. [15]

    J. W. Harris and B. M¨ uller, Eur. Phys. J. C84, 247 (2024), arXiv:2308.05743 [hep-ph]

  16. [16]
  17. [17]

    S. K. Das, F. Scardina, S. Plumari, and V. Greco, Phys. Lett. B747, 260 (2015), arXiv:1502.03757 [nucl-th]

  18. [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]

  19. [19]

    Energy loss, hadronization and hadronic interactions of heavy flavors in relativistic heavy-ion collisions

    S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C92, 024907 (2015), arXiv:1505.01413 [nucl-th]

  20. [20]

    M. He, R. J. Fries, and R. Rapp, Phys. Rev. C86, 014903 (2012), arXiv:1106.6006 [nucl-th]

  21. [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

  22. [22]

    Shen and S

    C. Shen and S. Alzhrani, Phys. Rev. C102, 014909 (2020), arXiv:2003.05852 [nucl-th]

  23. [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]

  24. [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]

  25. [25]

    A. K. Panda, R. Gangadharan, and V. Roy, J. Phys. G 50, 075102 (2023), arXiv:2301.00632 [nucl-th]

  26. [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]

  27. [27]

    Y. Sun, S. Plumari, and S. K. Das, Phys. Lett. B843, 138043 (2023), arXiv:2304.12792 [nucl-th]

  28. [28]

    Chatterjee and P

    S. Chatterjee and P. Bozek, Phys. Lett. B798, 134955 (2019), arXiv:1804.04893 [nucl-th]

  29. [29]

    A. K. Panda, Pooja, M. L. Sambataro, S. Plumari, and S. K. Das, Directed Flow of D and B Mesons in an Elec- trically and Chirally Conductive QGP at LHC Energies (2026), arXiv:2603.09636 [hep-ph]

  30. [30]

    Large directed flow of open charm mesons probes the three dimensional distribution of matter in heavy ion collisions

    S. Chatterjee and P. Bo˙ zek, Phys. Rev. Lett.120, 192301 (2018), arXiv:1712.01189 [nucl-th]

  31. [31]

    S. K. Das, M. Ruggieri, F. Scardina, S. Plumari, and V. Greco, J. Phys. G44, 095102 (2017), arXiv:1701.05123 [nucl-th]

  32. [32]

    Singh, M

    M. Singh, M. Kurian, B. Schenke, S. Jeon, and C. Gale, Phys. Rev. C113, 024904 (2026), arXiv:2509.18647 [nucl-th]

  33. [33]

    S. K. Das, M. Ruggieri, S. Mazumder, V. Greco, and J.- e. Alam, J. Phys. G42, 095108 (2015), arXiv:1501.07521 [nucl-th]

  34. [34]

    Y. Sun, G. Coci, S. K. Das, S. Plumari, M. Rug- gieri, and V. Greco, Phys. Lett. B798, 134933 (2019), arXiv:1902.06254 [nucl-th]

  35. [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]

  36. [36]

    Directed flow in ultrarelativistic heavy-ion collisions

    P. Bozek and I. Wyskiel, Phys. Rev. C81, 054902 (2010), arXiv:1002.4999 [nucl-th]

  37. [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]

  38. [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]

  39. [39]

    M. I. Abdulhamidet al.(STAR), Phys. Rev. X14, 011028 (2024), arXiv:2304.03430 [nucl-ex]

  40. [40]

    J. Jia, J. Phys. G41, 124003 (2014), arXiv:1407.6057 [nucl-ex]

  41. [41]

    J. Y. Zhang, Y. J. Song, L. Ma, W. H. Ma, and H. Z. Huang, Phys. Rev. C112, 044901 (2025)

  42. [42]

    Parida and S

    T. Parida and S. Chatterjee, Phys. Rev. C106, 044907 (2022), arXiv:2204.02345 [nucl-th]

  43. [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]

  44. [44]

    Adareet al.(PHENIX), Phys

    A. Adareet al.(PHENIX), Phys. Rev. Lett.121, 222301 (2018), arXiv:1807.11928 [nucl-ex]

  45. [45]

    Nakamura, T

    K. Nakamura, T. Miyoshi, C. Nonaka, and H. R. Takahashi, Phys. Rev. C107, 014901 (2023), arXiv:2209.00323 [nucl-th]

  46. [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]

  47. [47]

    Adareet al.(PHENIX), Phys

    A. Adareet al.(PHENIX), Phys. Rev. C94, 054910 (2016), arXiv:1509.07784 [nucl-ex]

  48. [48]

    Triangularity and Dipole Asymmetry in Heavy Ion Collisions

    D. Teaney and L. Yan, Phys. Rev. C83, 064904 (2011), arXiv:1010.1876 [nucl-th]

  49. [49]

    F. G. Gardim, F. Grassi, Y. Hama, M. Luzum, and J.-Y. Ollitrault, Phys. Rev. C83, 064901 (2011), arXiv:1103.4605 [nucl-th]

  50. [50]

    Shafi and S

    K. Shafi and S. Chatterjee, Eur. Phys. J. C86, 93 (2026), arXiv:2505.11713 [nucl-th]

  51. [51]

    J. S. Moreland, J. E. Bernhard, and S. A. Bass, Estimat- ing initial state and quark-gluon plasma medium prop- erties using a hybrid model with nucleon substructure calibrated top-Pb and Pb-Pb collisions at √sNN = 5.02 TeV (2018), arXiv:1808.02106 [nucl-th]

  52. [52]

    M. He, H. van Hees, and R. Rapp, Prog. Part. Nucl. Phys. 130, 104020 (2023), arXiv:2204.09299 [hep-ph]

  53. [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

  54. [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]

  55. [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]

  56. [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

  57. [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]

  58. [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]

  59. [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]

  60. [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]

  61. [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]

  62. [62]

    Jiang, S

    Z.-F. Jiang, S. Cao, W.-J. Xing, X.-Y. Wu, C. B. Yang, and B.-W. Zhang, Phys. Rev. C105, 054907 (2022), arXiv:2202.13555 [nucl-th]