⚛️ phenomenology

Entropy Production and Thermalization in Heavy-Quark Langevin Dynamics

This paper formulates the irreversible entropy balance for heavy quarks in a quark-gluon plasma by separating collisional production from medium-driven effects, deriving closed-form strong- and weak-coupling limits that reveal a maximal dissipation regime at intermediate coupling, and validating these findings numerically for charm and bottom quarks under Bjorken expansion.

Original authors: A. M. Silva

Published 2026-10-19

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: A. M. Silva

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the heart of a heavy-ion collision, where atomic nuclei smash together at nearly the speed of light, a state of matter known as the quark-gluon plasma is born. This is a seething, super-hot soup of the fundamental particles that usually make up protons and neutrons. Within this chaotic environment, some particles are created that are much heavier than the surrounding heat, such as the charm and bottom quarks. Because they are so massive, these heavy particles do not get swept up and thermalized like the lighter particles; instead, they plow through the plasma, acting as independent probes that carry information about the medium they traverse. Physicists often model the movement of these heavy particles as a form of Brownian motion, similar to how pollen grains jitter in a drop of water, but on a subatomic scale. The challenge has been to understand exactly how much information these heavy particles lose to the plasma as they slow down and how the energy of the plasma itself, which is rapidly cooling and expanding, affects this process.

A recent study by A. M. Silva at the University of São Paulo tackles this problem by tracking the flow of disorder, or entropy, associated with these heavy particles. In thermodynamics, entropy is a measure of how much information about a system's initial state has been lost to its surroundings. When a heavy quark moves through the plasma, its momentum changes due to countless tiny collisions. The researcher formulated a precise balance sheet for this process, separating the entropy generated by the collisions themselves from the entropy changes caused by the plasma cooling down around the particle. The study reveals that the total disorder produced is not just a simple sum of friction; it is a complex interplay between the particle trying to catch up to a target that is constantly moving away.

The investigation began with a simplified, mathematical model where the heavy particle's momentum distribution was assumed to be a smooth, bell-shaped curve. In this idealized scenario, the researcher could calculate the exact amount of entropy produced as the particle relaxed toward equilibrium. The analysis showed that the total disorder splits cleanly into two positive parts: one part comes from the particle losing its directed motion, and the other comes from the spread of its momentum adjusting to the temperature of the medium. When the plasma is static and unchanging, the total entropy produced is fixed by how far the particle started from equilibrium. However, the real universe is not static. The plasma created in a collision expands and cools rapidly. The study found that when the plasma cools, the target temperature the particle is trying to reach is constantly shifting. This creates an extra cost in entropy production because the particle is always chasing a moving target.

To test these ideas under realistic conditions, the researcher moved beyond the simple mathematical model to simulate the actual behavior of charm and bottom quarks using a relativistic framework that accounts for the high speeds involved. The simulations followed the particles as they traveled through a cooling environment that mimics the expansion of the early universe. The results confirmed that the total entropy produced depends heavily on how strongly the particle interacts with the plasma. If the interaction is extremely weak, the particle barely notices the plasma and produces very little entropy. If the interaction is extremely strong, the particle adjusts almost instantly, but the cost of constantly chasing the cooling target still generates a specific amount of disorder.

The most striking finding emerged when looking at the middle ground. The study demonstrated that the total entropy production does not simply increase or decrease with the strength of the interaction; instead, it reaches a maximum at an intermediate level of coupling. This means there is a specific "sweet spot" where the heavy quark loses the most information about its initial state. This peak occurs because the particle is slow enough to feel the cooling of the plasma but fast enough to be dragged along by it, creating a perfect storm of irreversible energy loss. The researcher calculated that for charm quarks, this maximum occurs at a specific diffusion parameter, while for the heavier bottom quarks, it happens at a different value. These findings were verified by solving the equations numerically for both types of quarks, showing that the theoretical predictions hold up even when the complex details of relativistic motion are included.

The work also clarified how the expansion of the plasma itself contributes to the process. As the universe expands, the momentum of particles in the direction of the expansion is stretched, a kinematic effect that is distinct from the friction caused by collisions. The study separated this expansion effect from the collisional entropy production, showing that they are different mechanisms that must be accounted for separately. This distinction is crucial for interpreting experimental data from particle colliders, where scientists try to reconstruct the history of the quark-gluon plasma. By understanding exactly how much entropy is produced and where the maximum dissipation occurs, researchers can better determine the properties of the plasma, such as how easily it allows heavy particles to diffuse through it.

Ultimately, this research provides a rigorous diagnostic tool for understanding heavy-quark transport. It establishes that the loss of information by a heavy quark is not a monotonic process but one that peaks at a specific interaction strength. The study confirms that the total entropy produced is the sum of the initial distance from equilibrium and the extra cost incurred by the cooling medium. By calculating these values for both charm and bottom quarks, the work offers a clear, quantitative picture of how these particles thermalize. The results suggest that the maximum dissipation of information is a robust feature of the system, appearing consistently across different coupling strengths and particle types. This insight helps refine the models used to describe the quark-gluon plasma, ensuring that the theoretical tools used to interpret high-energy collision experiments are grounded in a precise understanding of entropy and irreversibility.

Technical Summary: Entropy Production and Thermalization in Heavy-Quark Langevin Dynamics

Problem Statement
Heavy quarks (charm and bottom) produced in relativistic nuclear collisions serve as Brownian probes of the quark-gluon plasma (QGP). While their momentum evolution is typically modeled via Langevin or Fokker–Planck equations, a rigorous formulation of the irreversible entropy balance for these probes—specifically separating production driven by collisions from that driven by the evolving medium—has not been fully detailed in the context of stochastic thermodynamics. Existing studies often focus on reduced density matrices in open quantum systems or total thermodynamic entropy of the bulk medium. This paper addresses the gap by formulating the entropy balance for the classical single-particle momentum density generated by Langevin transport, distinguishing it from the viscous entropy production of the bulk QGP.

Methodology
The author employs a multi-stage analytical and numerical approach:

Linear Langevin Framework: The study begins with a homogeneous, isotropic $d$ -dimensional Ornstein–Uhlenbeck process. The momentum distribution is assumed Gaussian, allowing for exact solutions for the mean momentum and covariance. The entropy balance is derived by separating the change in distribution entropy into irreversible production ( $\Pi$ ) and entropy flow ( $\Phi$ ).
Static vs. Driven Baths:
- Static Bath: The paper establishes an exact H-theorem for a static thermal bath, showing that the integrated entropy production equals the initial relative entropy (Kullback–Leibler divergence) between the initial state and the Maxwell distribution.
- Driven Bath (Cooling): The framework is extended to a cooling medium, specifically utilizing the Bjorken temperature protocol ( $T(t) \propto t^{-1/3}$ ). Here, the thermal target moves, introducing a driving term to the relative entropy balance.
Asymptotic Analysis: Closed-form expressions for total entropy production are derived in the strong-coupling (fast relaxation) and weak-coupling (slow relaxation) limits. The strong-coupling limit utilizes a quasistatic expansion, while the weak-coupling limit assumes a frozen distribution.
Relativistic Extension: To address physical relevance, the non-relativistic Gaussian model is extended to a relativistic Fokker–Planck equation with a Jüttner stationary state. This involves a detailed-balance generator with momentum-dependent diffusion coefficients.
Numerical Validation: The relativistic equations are solved numerically for charm and bottom quarks using a conservative radial grid and Scharfetter–Gummel fluxes. Linear-response coefficients are computed from the static generator to verify the analytical limits.

Key Contributions

Decomposition of Entropy Production: For linear dynamics, the irreversible entropy production is shown to split exactly into two positive contributions: one arising from the relaxation of the directed momentum (drift) and another from the relaxation of the momentum width (variance) toward its thermal value.
The Driven H-Theorem: The paper demonstrates that in a cooling bath, the integrated entropy production is not solely determined by the initial relative entropy. An explicit driving term arises from the motion of the thermal reference state.
Coupling Dependence and Maximal Dissipation:
- Strong Coupling: The total production exceeds the initial relative entropy by a term linear in the relaxation time. This term comprises a "thermodynamic length" contribution (quasistatic tracking) and a transient cost associated with relaxing against a moving target.
- Weak Coupling: The production vanishes inversely with the relaxation time.
- Intermediate Coupling: The combination of these limits proves that total dissipation reaches a maximum at intermediate coupling strengths, provided the strong-coupling slope is positive.
Longitudinal Expansion: The work distinguishes between collisional entropy production and kinematic drives. Longitudinal Bjorken expansion is identified as a separate kinematic drive (phase-space contraction) that modifies the entropy balance but does not contribute to the irreversible collisional production term.
Relativistic Benchmark: A relativistic detailed-balance generator is solved numerically, confirming that the linear-response coefficients derived from the static generator accurately reproduce the strong- and weak-coupling limits for both charm and bottom quarks.

Results

Analytical Limits: For the Bjorken temperature protocol, the total entropy production $\Sigma(t_f)$ $Σ (t_{f})$ is derived as:
- $\Sigma \approx H(0) + \epsilon \cdot (\text{transient} + \text{tracking})$ for strong coupling ( $\epsilon \to 0$ ), where $\epsilon$ is the relaxation-time ratio.
- $\Sigma \approx \text{const}/\epsilon$ for weak coupling ( $\epsilon \to \infty$ ).
Numerical Verification: Simulations for charm ( $M_c = 1.5$ $M_{c} = 1.5$ GeV) and bottom ( $M_b = 4.8$ $M_{b} = 4.8$ GeV) quarks under a temperature drop from 0.40 GeV to 0.155 GeV confirm the theoretical predictions.
- The numerical entropy balance is satisfied to high accuracy, with maximum absolute residuals of $5.9 \times 10^{-6}$ for charm and $1.6 \times 10^{-5}$ for bottom.
- The total production exceeds the initial relative entropy at strong coupling.
- A distinct maximum in total production is observed at intermediate diffusion strengths ( $c_D \approx 7.5$ for charm, $c_D \approx 2.9$ for bottom).
Transient Dominance: In the relativistic benchmark, the transient cost of initial relaxation against the cooling bath is found to be a larger source of excess dissipation than the quasistatic tracking term.

Significance and Claims
The paper claims to provide a "clean diagnostic" for transport calculations in heavy-flavor phenomenology. By isolating the momentum-space irreversibility, the derived entropy balance serves as a rigorous test for numerical Langevin codes, independent of spatial propagation or hadronization models.

The author emphasizes that this quantity measures the irreversible loss of information regarding the initial momentum state under a specific medium protocol. It is distinct from the total thermodynamic entropy of the QGP. The primary significance lies in establishing that the dissipation of heavy quarks is not monotonic with coupling strength; rather, it is maximized at intermediate coupling. This structure holds for both non-relativistic Gaussian approximations and relativistic detailed-balance generators, suggesting that the interplay between relaxation rates and the cooling of the medium is a fundamental feature of heavy-quark thermalization. The work does not propose new experimental setups but offers a theoretical framework to interpret transport coefficients and entropy budgets in existing hydrodynamic simulations.

More like this