Finite-size effects and scaling properties of chiral and baryon-number fluctuations

Technical Summary: Finite-size effects and scaling properties of chiral and baryon-number fluctuations

Problem Statement
Locating the critical endpoint (CEP) in the phase diagram of strongly interacting matter is a primary objective in heavy-ion physics. While theoretical effective models in the thermodynamic limit can describe critical behaviors at high baryon density, experimental systems in heavy-ion collisions possess finite linear extents (typically a few to tens of femtometers). In such finite systems, true criticality is precluded by finite-size rounding, critical slowing down, and finite time effects. Consequently, theoretical frameworks must incorporate finite-size effects to interpret fluctuation observables. Existing approaches face limitations: Lattice QCD is restricted to vanishing baryon chemical potential ($\mu_B$) due to the sign problem, and Chiral Perturbation Theory (ChPT) is valid only at low temperatures. Furthermore, previous mean-field studies incorporating momentum-space constraints often fail to remove critical divergences or exhibit trivial finite-size scaling. This work addresses the need for a reliable modeling of finite-volume systems valid at high density to capture signatures of a potential critical point.

Methodology
The authors introduce an effective chiral model based on a quark-meson framework ($N_f=1$) to study finite-size effects. The core methodology involves a non-perturbative integration over the constant (zero) mode of the scalar field, $\phi_0$, while treating the system in a finite cubic volume $V=L^3$. The approach is structured in three progressive layers of approximation:

1. Zero-Mode Integration: The partition function is reduced to an integral over the spatially constant mode $\bar{\phi}$, governed by a volume-dependent statistical weight $e^{-\beta V U_{eff}(\bar{\phi})}$. This isolates the effect of statistical weighting in finite volumes, effectively reducing the theory to a zero-dimensional field theory.
2. Momentum-Space Discretization: The model incorporates constraints on the fermionic momentum modes due to the finite volume (using periodic boundary conditions). This introduces an intrinsic size dependence into the effective potential, specifically modifying the fermionic vacuum and matter contributions.
3. Gradient Effects: To account for spatial correlations beyond the zero-mode approximation, the authors employ an ansatz for a static, spatially varying classical field $\phi(\vec{x}) = \bar{\phi} f(\vec{x})$. This incorporates gradient terms into the effective action without solving the full equation of motion, allowing for an assessment of kinetic contributions.

The study utilizes the Ising scaling ansatz to analyze the singular part of the free energy density, mapping the model parameters ($T, \mu_B, h$) to Ising scaling fields ($\tau, \eta$). The analysis focuses on the chiral condensate, chiral susceptibility, and net-baryon number fluctuations (specifically the Binder cumulant and kurtosis) along the phase boundary and the freeze-out line.

Key Contributions and Results

• Finite-Size Rounding and Scaling: The integration over the zero mode naturally leads to finite-size rounding of phase transitions. For a finite but large volume, non-critical quantities remain largely unaffected, while fluctuations near the critical point follow expected finite-size scaling laws. However, for sufficiently small systems ($L \lesssim 10-20$ fm), genuine finite-size effects emerge, altering the thermodynamics and phase diagram significantly.
• Chiral Condensate and Susceptibility:
  • Using a volume-independent effective potential ($U_{eff}^\infty$), the model predicts a decrease in the chiral condensate and a softening of the transition as $L$ decreases, leading to a reduction in the pseudocritical temperature.
  • Including momentum-space discretization reverses this trend at low temperatures, causing the condensate to increase and the transition to steepen, consistent with some lattice results but differing in the behavior of the pseudocritical temperature.
  • The inclusion of gradient terms further reduces the condensate and softens the transition, significantly lowering the pseudocritical temperature for $L \lesssim 10-20$ fm.
• Scaling Behavior:
  • Critical Point: The singular parts of the chiral condensate and susceptibility exhibit standard finite-size scaling for $L \gtrsim 20$ fm.
  • First-Order Transition: In the vicinity of the first-order transition, the susceptibility scales as $\chi \propto L^d$ due to phase coexistence. This coexistence prevents the identification of an "apparent critical endpoint" merely by maximizing susceptibility along the phase boundary, a limitation noted in functional methods relying solely on momentum constraints.
• Fluctuations and Observables:
  • Binder Cumulant: The Binder cumulant ($\kappa_B$) is shown to be size-independent ($L^0$) along the phase boundary, crossing at the critical point. This provides a robust method to locate the CEP even when the thermodynamic limit is not reached.
  • Freeze-out Line: Along the freeze-out line (outside the scaling region), the kurtosis ratio ($R_{42} = \chi_4/\chi_2$) exhibits a dip-peak structure. The authors find that this non-monotonic behavior arises from the separation between the freeze-out line and the phase boundary, not necessarily from the proximity of a critical point.
  • Finite-Size Impact on Observables: The dip-peak structure in $R_{42}$ is sensitive to finite-size effects, becoming suppressed and shifting to lower collision energies ($\sqrt{s}$) for $L \lesssim 40$ fm. Momentum discretization and gradient terms alter the magnitude and position of these features but do not eliminate the finite-volume modifications.

Significance and Claims
The paper claims to provide a qualitative theoretical picture of how finite spatial extent affects phase transitions and fluctuation observables. Its primary significance lies in distinguishing between "scaling" effects (valid for large volumes) and "genuine" finite-size effects (dominant at small volumes, $L \lesssim 10-20$ fm).

The authors emphasize that:

1. Coexistence Dominance: In finite volumes, the coexistence of phases at a first-order transition leads to a monotonic increase in susceptibility along the phase boundary, obscuring the CEP. Therefore, standard extremum searches for the CEP are insufficient; ratios like the Binder cumulant are required.
2. Interpretation of Non-Monotonicity: The dip-peak structure in baryon number kurtosis along the freeze-out line is not an unambiguous signal of a nearby critical point in finite systems; it can be generated solely by the geometric separation of the freeze-out and phase boundaries.
3. Model Dependence: The specific behavior of finite-size effects (e.g., the shift in pseudocritical temperature) depends heavily on the treatment of momentum constraints and gradient terms. The gradient term, often neglected in mean-field studies, is shown to be substantial for small systems.

The work concludes that a comprehensive theoretical description including finite-size effects is mandatory for comparing experimental fluctuation data with theoretical predictions to locate the CEP. The authors note that their framework is schematic and qualitative, intended to isolate mechanisms rather than provide quantitative predictions for specific experimental energies, and that future work should address canonical ensemble effects and self-consistent solutions for spatially varying fields.