Effects of Subnucleonic Fluctuations on the Longitudinal Dynamics of Heavy-Ion Collisions

Basic Information

• Paper ID: 2510.14118
• Title: Effects of Subnucleonic Fluctuations on the Longitudinal Dynamics of Heavy-Ion Collisions
• Authors: Oscar Garcia-Montero, Sören Schlichting, Jie Zhu
• Classification: nucl-th (Nuclear Theory), hep-ph (High Energy Physics - Phenomenology)
• Publication Date: October 15, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.14118v1

Abstract

This study provides an in-depth investigation of the effects of subnucleonic fluctuations on the longitudinal dynamics in the initial state of heavy-ion collisions. The authors employ a hybrid model combining a three-dimensionally resolved McDipper initial state model with the (3+1)-dimensional viscous hydrodynamic framework CLVisc to study the emergence of long-range rapidity structures in Pb+Pb collisions. By incorporating (sub)nucleonic-scale fluctuations from different sources, the authors analyze the effects of these fluctuations on the longitudinal structure of observables including charged hadron multiplicity, baryon stopping, directed flow, and flow decorrelation. The study finds significant agreement with experimental data regarding directed flow and rapidity-resolved charged particle multiplicity.

Research Background and Motivation

Core Issues

1. Inadequacy in Longitudinal Dynamics Description: Traditional theoretical models for heavy-ion collisions predominantly employ the approximate boost invariance assumption, which is appropriate in the central rapidity region at high collision energies but proves insufficient in describing the longitudinal structures revealed by experiments, particularly for small systems (O+O, Ne+Ne collisions) and asymmetric systems (p+A, d+A).
2. Importance of Subnucleonic Fluctuations: While nucleonic-scale fluctuations are crucial for explaining collective phenomena and centrality dependence in heavy-ion collisions, subnucleonic fluctuations are necessary for describing fluctuation-driven observables or smaller collision systems.

Research Significance

• Experimental Observation Requirements: Observables such as baryon stopping, rapidity-dependent anisotropic flow, and flow decorrelation across pseudorapidity indicate the need for more realistic three-dimensional initial conditions
• Theoretical Refinement: Most studies focus primarily on the effects of subnucleonic fluctuations in the transverse plane, with relatively limited attention to their effects on longitudinal decorrelation observables

Limitations of Existing Approaches

• Boost Invariance Assumption: Fails in non-central rapidity regions and small system collisions
• Lack of 3D Description: Cannot accurately describe longitudinal structures and long-range correlations
• Incomplete Fluctuation Sources: Lacks systematic investigation of subnucleonic-level fluctuations

Core Contributions

1. Established a Complete 3D Hybrid Model: Combined McDipper v1.2 initial state model with CLVisc (3+1)D viscous hydrodynamic framework
2. Systematically Studied Subnucleonic Fluctuation Effects: Compared five different fluctuation configurations and their effects on longitudinal observables
3. Achieved Good Agreement with Experimental Data: Demonstrated significant consistency with directed flow and rapidity-resolved charged particle multiplicity
4. Revealed the Importance of Subnucleonic Structure: Demonstrated the relevance of nucleonic substructure in shaping the longitudinal dynamics of heavy-ion collisions

Methodology Details

Task Definition

Study the longitudinal dynamics of Pb+Pb collisions at √sNN = 2.76 TeV by incorporating different levels of subnucleonic fluctuations to predict and explain experimentally observed rapidity-dependent observables.

Model Architecture

McDipper Initial State Model

• Theoretical Foundation: Based on the Color Glass Condensate (CGC) kT-factorization framework
• Gluon Production: Described through single gluon production formula
• Quark Energy and Charge Deposition: Derived from multiple scattering of collinear quarks and stopping
• Subnucleonic Fluctuation Implementation:
  • Decompose nucleons into Gaussian hotspots with event-by-event position fluctuations
  • Implement thickness fluctuations by sampling normalized nucleon (hotspot) intensities from log-normal distributions

Five Fluctuation Configurations

1. Smooth Nucleon: No subnucleonic fluctuations
2. Three Hotspots: Includes position fluctuations
3. Nucleon Fluctuations (σ=0.637): Nucleonic-scale thickness fluctuations
4. Hotspots + Fluctuations (σ=0.637): Three hotspots + moderate thickness fluctuations
5. Hotspots + Strong Fluctuations (σ=1.2): Three hotspots + strong thickness fluctuations

CLVisc Hydrodynamic Evolution

• Dimensionality: (3+1)D viscous hydrodynamics
• Conservation Equations: Energy-momentum and baryon number conservation, including second-order Israel-Stewart-type terms
• Transport Coefficients: Selected based on Bayesian analysis of LHC data
• Equation of State: NEOS-B, including baryon chemical potential effects
• Particlization: Using Cooper-Frye prescription at energy density εfrz = 0.266 GeV/fm³

Technical Innovations

1. Complete 3D Description: Breaks through traditional boost invariance assumption, achieving genuine three-dimensional evolution
2. Multi-level Fluctuation Modeling: Systematically incorporates different fluctuation sources from nucleonic to subnucleonic scales
3. End-to-End Framework: Complete description from initial to final state without requiring additional parameter tuning

Experimental Setup

Collision System

• Collision Type: Pb+Pb
• Collision Energy: √sNN = 2.76 TeV
• Centrality: 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%

Evaluation Metrics

1. Charged Hadron Multiplicity: dN/dη distribution
2. Net Proton Multiplicity: dN/dy distribution
3. Directed Flow: Odd and even components of v₁
4. Flow Decorrelation: Forward-backward ratio rn(ηa, ηb)

Comparison Data

• ALICE Experimental Data: Charged hadron multiplicity, directed flow
• CMS Experimental Data: Flow decorrelation

Implementation Details

• Initialization Time: τ = 0.6 fm
• Reference Plane: Spectator plane
• Reference Rapidity: 4.4 < ηb < 5 (flow decorrelation analysis)

Experimental Results

Main Findings

Charged Hadron Multiplicity

• Overall Agreement: McDipper hybrid model shows significant consistency with experimental data across most rapidity ranges
• Large Rapidity Improvement: In the |η| > 3 region, agreement with data improves after including subnucleonic fluctuations
• Physical Interpretation: For smaller (lower density) systems, subnucleonic scales represent relevant degrees of freedom

Net Proton Multiplicity

• Fluctuation Effects: Including subnucleonic fluctuations reduces total baryon charge deposition
• Stopping Power: Does not affect stopping power (peak position), primarily dependent on collision energy
• Mechanism Understanding: Results from increased initial state granularity

Directed Flow v₁

• Odd Component: Reflects tilted collision geometry, well reproduced in all cases
• Even Component: Originates from event-by-event asymmetries
  • Configuration (i) is underestimated
  • Significantly enhanced after including fluctuations, improving agreement with experiment

Flow Decorrelation

• Baseline Issue: Configuration (i) underestimates rn observed in CMS data
• Fluctuation Improvement: Including subnucleonic fluctuations increases decorrelation, particularly in central collisions
• Residual Insufficiency: Still falls below experimental results in peripheral collisions, attributed to missing fluctuation sources in the quark sector

Ablation Studies

By comparing five different fluctuation configurations, the authors systematically analyze the contributions of various fluctuation sources:

• Position fluctuations (hotspots) primarily affect small rapidity regions
• Thickness fluctuations show significant effects across all rapidity regions
• Increased fluctuation strength (σ values) further enhances effects

Related Work

Main Research Directions

1. CGC Framework: Application of Color Glass Condensate theory in describing heavy-ion collision initial states
2. Hydrodynamic Simulations: Application of (3+1)D viscous hydrodynamics in QGP evolution
3. Fluctuation Effects: Effects of nucleonic and subnucleonic fluctuations on collective flow

Advantages of This Work

• Complete 3D Description: Significant improvement over boost-invariant models
• Systematic Investigation: Systematic comparison of different fluctuation sources
• Experimental Verification: Detailed comparison with multiple experimental observables

Conclusions and Discussion

Main Conclusions

1. Importance of Subnucleonic Fluctuations: Subnucleonic fluctuations have clear effects on the longitudinal dynamics of heavy-ion collisions
2. Experimental Agreement: The framework provides good descriptions of observables without requiring adjustments to transport coefficients
3. Relevance of Nucleonic Substructure: Nucleonic substructure plays an important role in shaping the longitudinal dynamics of heavy-ion collisions

Limitations

1. Absence of Hadronic Afterburner: May affect quantitative accuracy
2. Quark Sector Fluctuations: Missing quark fluctuation sources most relevant in fragmentation regions
3. Pre-equilibrium Stage: Lacks realistic pre-equilibrium stage description

Future Directions

1. Valence Quark Fluctuations: Incorporate valence quark fluctuations
2. Realistic Pre-equilibrium Stage: Implement more realistic pre-equilibrium dynamics
3. Hadronic Afterburner: Include hadronic afterburner for improved quantitative accuracy

In-Depth Evaluation

Strengths

1. Methodological Innovation:
   • First systematic study of subnucleonic fluctuation effects on longitudinal dynamics
   • Established complete 3D hybrid model framework
2. Experimental Sufficiency:
   • Detailed comparison with multiple experimental observables
   • Covers different centralities and rapidity ranges
3. Result Convincingness:
   • Achieves good agreement with experiments on key observables
   • Provides clear physical picture and mechanism explanations
4. Writing Clarity:
   • Clear structure with rigorous logic
   • Accurate technical detail descriptions

Weaknesses

1. Model Completeness:
   • Absence of hadronic afterburner may affect final-state observables
   • Simplified treatment of pre-equilibrium stage
2. Parameter Tuning:
   • Some parameters (e.g., gluon energy normalization) require matching with experimental data
   • Fluctuation parameter selection lacks first-principles guidance
3. Systematic Uncertainties:
   • Insufficient discussion of theoretical uncertainties
   • Lacks systematic analysis of model parameter sensitivity

Impact

1. Field Contribution:
   • Provides new theoretical framework for heavy-ion collision longitudinal dynamics research
   • Deepens understanding of subnucleonic fluctuation effects
2. Practical Value:
   • Applicable to other collision systems and energies
   • Provides theoretical support for experimental data analysis
3. Reproducibility:
   • Clear description of model parameters and settings
   • Used code frameworks (McDipper, CLVisc) are publicly available

Applicable Scenarios

1. Heavy-Ion Collisions: Particularly suitable for longitudinal dynamics studies in large system collisions
2. Small System Collisions: For small system collisions requiring consideration of subnucleonic fluctuations
3. Energy Scan: Extensible to studies at different collision energies
4. New Observable Predictions: Can be used to predict new longitudinal correlation observables

References

The paper cites 16 important references covering key areas including CGC theory, hydrodynamic simulations, and experimental data, providing a solid theoretical and experimental foundation for the research.

Overall Assessment: This is a high-quality theoretical physics research paper that achieves important progress in the study of longitudinal dynamics in heavy-ion collisions. The authors established a complete 3D theoretical framework, systematically investigated the effects of subnucleonic fluctuations, and conducted detailed comparisons with experimental data. Despite some technical limitations, it makes significant contributions to the development of this field.