Heat’s Ripple Effect: Reshaping Critical Fluctuations in Heavy-Ion Collisions

Author: Denis Avetisyan

New research reveals how temperature gradients fundamentally alter the behavior of the quark-gluon plasma, leading to directional patterns in its evolution.

The study demonstrates how individual, low-energy fluctuation modes contribute to the nonlocal two-point correlation function - measured in MeV and assessed at a fixed spatial point of (5.83, 0) fm with a Matsubara mode of <span class="katex-eq" data-katex-display="false">\lambda = 0</span> and a chemical potential of 240 MeV - revealing patterns determined by combinations of radial and angular momentum quantum numbers. — The study demonstrates how individual, low-energy fluctuation modes contribute to the nonlocal two-point correlation function – measured in MeV and assessed at a fixed spatial point of (5.83, 0) fm with a Matsubara mode of $\lambda = 0$ and a chemical potential of 240 MeV – revealing patterns determined by combinations of radial and angular momentum quantum numbers.

Temperature gradients drive anisotropic correlations in critical fluctuations of the chiral order parameter field, offering a novel pathway to probe the QCD phase transition.

While studies of the QCD phase transition typically assume spatially uniform temperatures, this overlooks the non-equilibrium dynamics of the hot fireball created in heavy-ion collisions. In this work, ‘Critical fluctuation patterns and anisotropic correlations driven by temperature gradients’, we investigate how spatial temperature gradients reshape critical fluctuations of the chiral order parameter, revealing anisotropic correlations within the system. We demonstrate that these gradients induce long-ranged correlations along isotherms, suppressed radially, and result in a superposition of angular momentum modes contributing comparably to observed anisotropic flow. Could azimuthally sensitive observables, therefore, provide a novel pathway for detecting signatures of the QCD phase transition and characterizing the thermal structure of the early Universe?

The Elusive Dance of Quarks: Unveiling the Primordial Transition

The transition between quark-gluon plasma and ordinary hadronic matter, known as the Quantum Chromodynamics (QCD) phase transition, presents a significant hurdle in contemporary theoretical physics. This shift, theorized to have occurred in the early universe microseconds after the Big Bang, involves a dramatic change in the fundamental constituents of matter. Unlike phase transitions observed in everyday life – such as water freezing into ice – the QCD transition occurs at extraordinarily high temperatures and densities, making direct experimental observation incredibly difficult. Current understanding relies heavily on complex simulations and indirect measurements from heavy-ion collisions, attempting to recreate the conditions of the early universe. The challenge lies not only in the extreme conditions but also in the strong interactions governing quarks and gluons, which render traditional perturbative methods ineffective, demanding innovative approaches to model and predict the behavior of matter at these scales.

Investigating the quark-gluon plasma phase transition necessitates exploring matter at temperatures and densities far exceeding those found in everyday experience-conditions akin to those present moments after the Big Bang or within the cores of neutron stars. At these extreme scales, the strong nuclear force, responsible for binding quarks into hadrons, behaves in a fundamentally different manner than predicted by conventional perturbative calculations. These methods, successful in describing weaker interactions, break down because the force becomes so strong that quarks and gluons are no longer distinct particles, but rather participate in a collective, strongly coupled system. Consequently, physicists rely on non-perturbative approaches, such as lattice quantum chromodynamics (QCD) and effective field theories, to model the plasma’s behavior and discern the nature of this phase transition-whether it’s a smooth crossover or a sharp, discontinuous change in state, potentially revealing new states of matter.

The vicinity of the QCD phase transition, known as the critical region, is characterized by dramatic fluctuations in the quark-gluon plasma – akin to a boiling liquid where local order rapidly breaks down and reforms. These intense fluctuations aren’t merely a hindrance to observation; they represent a powerful signal, offering a window into the underlying dynamics of the transition itself. However, extracting meaningful information from this chaotic state requires analytical techniques far beyond standard methods. Researchers employ sophisticated tools, including renormalization group calculations and high-precision lattice simulations, to map the behavior of the plasma and identify universal features – critical exponents and scaling laws – that characterize the transition, regardless of specific details. Effectively, the critical region transforms from a source of complexity into a testing ground for fundamental theories of strong interactions, provided these advanced analytical tools are brought to bear.

The correlation <span class="katex-eq" data-katex-display="false">\langle\tilde{\sigma}({\bm{r}}, \tau)\tilde{\sigma}({\bm{r}^{\prime}}, \tau)\rangle d_{z}</span> in MeV, evaluated at <span class="katex-eq" data-katex-display="false">{\bm{r}^{\prime}}=(2.5,0)</span>, <span class="katex-eq" data-katex-display="false">(6,0)</span>, and <span class="katex-eq" data-katex-display="false">(7.5,0)</span> fm, reveals dependencies on the energy cutoff (250 and 450 MeV) and angular momentum (up to <span class="katex-eq" data-katex-display="false">l=6</span>) for radial modes <span class="katex-eq" data-katex-display="false">n=1</span> and <span class="katex-eq" data-katex-display="false">n=2</span> at a chemical potential of 240 MeV and Matsubara mode <span class="katex-eq" data-katex-display="false">\lambda=0</span>. — The correlation $\langle\tilde{\sigma}({\bm{r}}, \tau)\tilde{\sigma}({\bm{r}^{\prime}}, \tau)\rangle d_{z}$ in MeV, evaluated at ${\bm{r}^{\prime}}=(2.5,0)$ , $(6,0)$ , and $(7.5,0)$ fm, reveals dependencies on the energy cutoff (250 and 450 MeV) and angular momentum (up to $l=6$ ) for radial modes $n=1$ and $n=2$ at a chemical potential of 240 MeV and Matsubara mode $\lambda=0$ .

Navigating the Complexity: Theoretical Frameworks for Strong Interactions

Lattice Quantum Chromodynamics (QCD) represents a non-perturbative, first-principles approach to solving QCD equations. This method discretizes spacetime into a four-dimensional lattice, transforming differential equations into algebraic equations that can be solved numerically. Calculations involve evaluating functional integrals over gauge fields and quark/antiquark fields on this lattice, approximating the path integral formulation of QCD. While offering a theoretically sound method for addressing phenomena inaccessible to perturbative calculations – such as hadron masses, decay constants, and the quark-gluon plasma equation of state – Lattice QCD demands substantial computational resources. The required computing power scales rapidly with decreasing lattice spacing (to approach the continuum limit) and increasing volume, necessitating the use of high-performance computing facilities and advanced algorithms to manage the associated computational cost.

Effective models in strong interaction theory utilize simplification to achieve analytical solutions where direct computation of Quantum Chromodynamics (QCD) is intractable. These models, exemplified by the Ising Model, reduce the complexity of $SU(N_c)$ gauge theory by focusing on relevant degrees of freedom, often associated with the chiral order parameter, which describes spontaneous chiral symmetry breaking. This parameter quantifies the condensation of quark-antiquark pairs and dictates the properties of hadrons. By concentrating on these collective behaviors, rather than individual quarks and gluons, effective models allow for the investigation of phenomena like confinement and hadronization through techniques like mean-field theory and renormalization group analysis, providing qualitative and, in some cases, semi-quantitative insights into the behavior of strongly interacting matter.

Functional methods in quantum chromodynamics (QCD) utilize integral representations, specifically functional integrals, to study the behavior of strong interactions when perturbative approaches fail. These methods bypass the limitations of asymptotic expansions by summing over all possible field configurations, weighted by the QCD action, and effectively include non-perturbative effects. Techniques within this framework include the Dyson-Schwinger equations and the Functional Renormalization Group (FRG), which allow for the calculation of Green’s functions and the running of coupling constants beyond the perturbative regime. By exploring the functional dependence of quantities like the quark propagator and gluon propagator, functional methods aim to reveal the dynamical chiral symmetry breaking and confinement mechanisms inherent in QCD, providing insights into hadron structure and properties that are inaccessible through traditional perturbation theory. The resulting equations are typically solved numerically or through variational approximations due to their inherent complexity.

Solutions to equation (15) demonstrate that the spatial temperature distribution and base field <span class="katex-eq" data-katex-display="false">\sigma_c</span> vary with radius ρ depending on the chemical potential, contrasting with the constant-volume condition <span class="katex-eq" data-katex-display="false">\partial_{\sigma}\mathcal{V}=0</span>. — Solutions to equation (15) demonstrate that the spatial temperature distribution and base field $\sigma_c$ vary with radius ρ depending on the chemical potential, contrasting with the constant-volume condition $\partial_{\sigma}\mathcal{V}=0$ .

Decoding the Signals: Probing Fluctuations and Correlations

The fluctuation spectrum quantifies the amplitude of disturbances across various spatial frequencies, or modes, within a system. Analysis of this spectrum reveals how the system responds to external or internal perturbations; a broader spectrum indicates greater sensitivity and a wider range of responsive modes, while a narrower spectrum suggests a more constrained response. Specifically, the energy distribution across these modes – typically represented as a power spectral density – provides a detailed map of the system’s dynamic properties. Examining the spectral shape allows for the identification of dominant modes and the characterization of the system’s preferred patterns of instability or relaxation. $S(q, \omega)$ typically represents the fluctuation spectrum, where $q$ denotes the wavevector and ω the frequency of the fluctuations.

Eigenmodes represent the natural spatial patterns of fluctuations within the system and are mathematically described by their associated angular momentum. These modes are not simply static shapes; their frequencies and spatial extents are directly influenced by the dynamics of the phase transition. Specifically, changes in the transition’s velocity or the presence of inhomogeneities alter the eigenmode spectrum, leading to shifts in peak locations and broadening of spectral lines. Analyzing these changes provides a means to probe the underlying mechanisms driving the transition, as the sensitivity of eigenmodes allows for the detection of subtle variations in the system’s order parameter and correlation length. The angular momentum component further defines the rotational symmetry of each mode, impacting how fluctuations propagate and interact within the system.

The presence of a temperature gradient introduces a finite energy gap in the observed fluctuation spectrum, a deviation from the expectation of a zero-energy gap in systems maintained at a uniform temperature. This energy gap represents a minimum energy required to excite fluctuations within the system and is directly attributable to the imposed thermal disequilibrium. Measurements indicate that the magnitude of this gap is proportional to the temperature gradient; larger gradients result in wider gaps, suggesting a direct relationship between the driving force and the energetic cost of fluctuation. This finding demonstrates that spatial inhomogeneities, specifically temperature gradients, fundamentally alter the energetic landscape of fluctuations, impacting their dynamics and observable characteristics.

Spatial correlations within the system demonstrate how fluctuations at one location are linked to those at others, a relationship directly affected by the imposed temperature gradient. These correlations are not isotropic; instead, they exhibit anisotropy, meaning their strength and extent vary depending on the direction relative to the temperature gradient. Quantification of these anisotropic correlations is achieved through the calculation of correlation functions in multiple spatial dimensions, revealing elongated or directional patterns of interconnected fluctuations. The degree of anisotropy provides a measurable characteristic of how the temperature gradient modifies the system’s collective behavior, moving it away from the isotropic correlations observed in systems at uniform temperature.

Analysis demonstrates that fluctuations in the system are not uniformly distributed, but are instead concentrated in the immediate vicinity of the phase transition interface. This localization is a direct consequence of the imposed temperature gradient, which disrupts the expected fluctuation behavior observed in homogeneous systems. Specifically, the non-uniform thermal profile creates a region of heightened sensitivity at the interface, leading to an increase in the amplitude of fluctuations in that area and a corresponding decrease further away. This represents a discernible shift in the typical fluctuation patterns, moving from a broad distribution to one strongly peaked at the phase transition boundary.

Three- and four-point correlation functions reveal dependencies on density ρ and angular momentum mode at μ values of 120, 240, and 360 MeV.

Modeling the Transient: Non-Equilibrium Dynamics in Detail

Langevin and Fokker-Planck equations are employed to describe the probabilistic behavior of systems not in thermodynamic equilibrium. The Langevin equation, a stochastic differential equation, models the time evolution of a system’s state by incorporating a deterministic force alongside a random, fluctuating force and a dissipative force proportional to the system’s velocity. The Fokker-Planck equation, derived from the Kramers-Moyal expansion of the master equation, provides a complementary approach by directly describing the time evolution of the probability density function $P(x,t)$ for the system’s state $x$ at time $t$ . Both equations are particularly useful for analyzing systems subject to thermal fluctuations and provide a framework for understanding phenomena like Brownian motion, diffusion, and relaxation processes in complex systems where a full microscopic description is impractical.

Langevin and Fokker-Planck equations incorporate stochastic forces and diffusion terms to model the probabilistic evolution of dynamical systems away from equilibrium. Stochastic forces, often represented as $\xi(t)$ , account for random fluctuations originating from numerous microscopic degrees of freedom not explicitly included in the macroscopic description. Diffusion, characterized by a diffusion coefficient $D$ , describes the spreading of probability due to these random fluctuations, effectively representing the system’s tendency toward increased entropy. These equations do not provide deterministic trajectories, but rather describe the time evolution of probability distributions, allowing for the calculation of ensemble averages and statistical properties of the system’s dynamics despite the inherent randomness.

Finite-size effects arise in modeling non-equilibrium dynamics due to the inherent limitations imposed by simulating or analyzing systems with a restricted number of degrees of freedom or spatial extent. These effects manifest as deviations from the behavior predicted by the corresponding infinite-size or continuum limit. Specifically, boundary conditions and correlations between system elements become more pronounced, influencing quantities such as diffusion coefficients, relaxation times, and correlation functions. Consequently, observed dynamics may exhibit artificial constraints or altered scaling behavior not representative of the true, macroscopic limit, necessitating careful consideration and potentially requiring scaling analyses or extrapolation techniques to accurately characterize the system’s behavior.

Analysis of non-equilibrium dynamics demonstrates that various angular momentum modes contribute with comparable magnitudes to the overall correlation strength. This finding diverges from the behavior observed in homogeneous systems, where a dominant zero-mode typically dictates correlation behavior. Specifically, the research indicates that higher-order angular momentum modes are not negligible, and their collective influence is essential for accurately modeling the system’s response. This suggests that models relying solely on the zero-mode approximation may underestimate the full extent of correlations in these non-equilibrium scenarios, and a multi-mode approach is necessary for comprehensive analysis.

From Theory to Observation: Bridging the Gap with Experiment

To experimentally investigate the predicted phase transition from ordinary matter to the quark-gluon plasma (QGP), scientists utilize heavy-ion collisions at relativistic energies. These collisions, typically involving gold or lead nuclei, generate extraordinarily high temperatures and energy densities – exceeding those found in the cores of supernovae or even moments after the Big Bang. Such extreme conditions momentarily liberate quarks and gluons, effectively recreating the state of matter theorized to have existed in the early universe. The resulting QGP is not merely hot matter; it behaves as an almost-perfect fluid, displaying properties fundamentally different from those of ordinary hadronic matter. By carefully analyzing the debris produced in these collisions, physicists can indirectly probe the characteristics of this fleeting, exotic state and test the predictions of Quantum Chromodynamics (QCD), the theory governing the strong nuclear force.

The incredibly energetic heavy-ion collisions don’t just produce heat; they forge a rapidly evolving temperature gradient across the newly created matter. This gradient isn’t uniform; instead, it fuels fluctuations in density and pressure, leading to correlated movements of particles. These correlations aren’t random; they exhibit a directional preference, resulting in what physicists term anisotropic flow. Essentially, the particles don’t expand equally in all directions, but rather preferentially along the direction of the initial temperature gradient. The strength and pattern of this anisotropic flow – often described by coefficients like $v_2$ – provides a crucial window into the early stages of the collision and the properties of the quark-gluon plasma formed, allowing researchers to map the behavior of matter under extreme conditions.

Anisotropic flow, a collective behavior observed in heavy-ion collisions, serves as a powerful diagnostic tool for understanding the quark-gluon plasma (QGP). This phenomenon, where particles preferentially stream along certain directions, arises from the initial temperature gradient and subsequent pressure gradients established in the extremely hot and dense medium. By meticulously analyzing the patterns of this flow – specifically, the coefficients characterizing its different harmonic components $v_n$ – experimentalists can reconstruct the transport properties of the QGP, such as its viscosity and equation of state. These measurements provide stringent tests of theoretical predictions derived from lattice QCD calculations and hydrodynamic models, allowing scientists to map out the phase diagram of strongly interacting matter and refine their understanding of the fundamental forces governing the universe at its most extreme conditions.

The study meticulously details how temperature gradients introduce anisotropy into the critical fluctuations surrounding the QCD phase transition. This reshaping of the order parameter field, and the resultant anisotropic correlations, speaks to a deeper principle: that form follows function, and that even seemingly random fluctuations reveal underlying order when examined with care. As Georg Wilhelm Friedrich Hegel observed, “The truth is the whole.” This research doesn’t simply document a phenomenon; it seeks the complete picture of non-equilibrium dynamics, demonstrating how external forces sculpt the very fabric of the critical transition and ultimately, influence observable signatures in heavy-ion collisions. The elegance lies in recognizing this holistic connection.

Beyond the Horizon

The presented framework, while offering a compelling mechanism for temperature-gradient-induced anisotropy in critical fluctuations, merely sketches the contours of a far more complex landscape. The assumption of a simple hydrodynamic background, a necessary simplification, feels increasingly provisional. One suspects the true non-equilibrium dynamics governing these heavy-ion collisions are textured with subtle, long-range correlations that this approach, in its current form, cannot fully capture. A consistent treatment incorporating genuine memory effects – a recognition that the system’s past influences its present behavior – feels critical, yet remains elusive.

Further refinement demands a rigorous investigation into the interplay between the chiral order parameter field and other relevant degrees of freedom. The model presently focuses on the effect of the temperature gradient on the critical fluctuations, but it would be prudent to explore how these fluctuations, in turn, modify the thermal transport itself. Such feedback loops, while mathematically challenging, may well be essential for a predictive theory. One hopes future work will embrace these complexities, acknowledging that elegance often resides not in simplicity, but in a faithful rendering of nature’s inherent intricacy.

Ultimately, the true test lies in confronting experimental data. Anisotropic flow measurements, as the authors suggest, offer a promising avenue for exploration. However, disentangling the signatures of these critical fluctuations from other sources of anisotropy will require a level of precision and theoretical control that currently remains a significant hurdle. The pursuit, nonetheless, is worthwhile, for it promises a deeper understanding of matter under extreme conditions – a glimpse into the very fabric of the early universe.

Original article: https://arxiv.org/pdf/2603.04923.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/