Clustered Nature of Hot and Dense Nuclear Matter: A quantum statistical approach

Technical Summary: Clustered Nature of Hot and Dense Nuclear Matter: A quantum statistical approach

Problem Statement
The equilibrium abundances of light nuclear clusters (deuteron $^2$H, triton $^3$H, helion $^3$He, and $\alpha$-particle $^4$He) in hot, dense nuclear matter near saturation density ($n_B \approx 0.16 \text{ fm}^{-3}$) are critical for nuclear physics and astrophysical applications, such as binary neutron star mergers. However, theoretical approaches yield diverging predictions. Specifically, there is a significant discrepancy between the Quantum Statistical (QS) approach and the recently proposed Phase-Space Excluded-Volume (PSEV) approach regarding the survival and abundance of $\alpha$-particles at densities approaching and exceeding saturation. While QS calculations predict a sharp decrease in cluster abundances as density increases, PSEV calculations suggest unexpectedly high $\alpha$-clustering ($X_\alpha \approx 0.18 - 0.26$) at densities of 1–2 times saturation density for temperatures between 20 and 40 MeV.

Methodology
The authors conduct a comparative analysis of the QS approach and the PSEV approach to identify the physical origins of these discrepancies. The study focuses on two primary ingredients:

1. The Mott Momentum ($P^{\text{Mott}}_\nu$): The critical momentum above which a bound state can exist in the medium. Below this momentum, the bound state dissolves into the continuum due to Pauli blocking.
2. Momentum Distribution Functions: The statistical distribution of clusters as a function of momentum and density.

The QS approach utilizes a Green's function formalism leading to an in-medium Schrödinger equation. It accounts for single-particle self-energy shifts, Pauli blocking, and continuum correlations via a generalized Beth-Uhlenbeck formula. The PSEV approach approximates Pauli blocking by introducing an overlap integral between the bound state wave function and the phase-space occupation numbers of the medium, governed by an empirical cutoff parameter ($F^{\text{cut}}_A$).

The authors derive exact relations for the deuteron Mott momentum using perturbation theory and compare them with the PSEV approximation. They also analyze the momentum distribution functions, contrasting calculations that use free-space binding energies against those that incorporate medium-modified quasiparticle energies (Pauli blocking shifts).

Key Contributions and Results

• Analysis of the Mott Momentum:

  • The PSEV approach determines the Mott momentum using an overlap integral condition (Eq. 4) with an empirical cutoff parameter $F^{\text{cut}}_A$. The authors demonstrate that the choice of this parameter significantly alters the density range where bound states are predicted to exist.
  • For deuterons, a cutoff of $F^{\text{cut}}_2 \approx 0.1$ aligns the PSEV Mott density with the QS result. However, for $\alpha$-particles, the PSEV parameter $F^{\text{cut}}_4 = 0.35$ (used in recent literature) allows bound states to exist at much higher densities than the QS approach predicts.
  • The authors show that using more realistic wave functions (e.g., Jastrow or Hulthen forms) rather than simple Gaussian approximations reduces the predicted existence region of bound states, bringing PSEV results closer to QS predictions, though the empirical cutoff remains a source of uncertainty.

• Impact of Pauli Blocking on Binding Energies:

  • A critical finding is that the PSEV approach often utilizes free-space binding energies for the cluster energy $\epsilon_\nu(P)$, neglecting the medium-induced shift (Pauli blocking shift) in the binding energy itself.
  • When the authors incorporate the Pauli blocking shift into the distribution function (Eq. 8), the calculated mass fractions drop drastically. For example, at $T=20$ MeV and $n_B = 0.03 \text{ fm}^{-3}$, the $\alpha$-particle mass fraction $X_\alpha$ decreases from 0.2212 to 0.01431 when the shift is included.
  • This omission in the PSEV approach leads to a significant overestimation of light cluster abundances, particularly $\alpha$-particles, at high densities.

• Continuum Correlations:

  • The QS approach treats the dissolution of bound states as a transition to resonance states in the continuum, which still contribute to the density. The PSEV approach, in its current formulation, does not provide a general method for treating these continuum correlations, often assuming a sharp cutoff (Heaviside step function) in the distribution function. The authors argue that the distribution function should approach zero smoothly as momentum falls below the Mott limit.

• Comparison with Generalized RMF:

  • The authors compare their results with the Generalized Relativistic Mean-Field (gRMF) approach. The gRMF results, which include medium shifts for clusters, show reasonable agreement with the QS approach, further validating the QS prediction of low $\alpha$-abundance near saturation density and highlighting the discrepancy with the unshifted PSEV results.

Significance and Conclusions
The paper concludes that the primary reasons for the diverging results between the QS and PSEV approaches are:

1. The determination of the Mott momentum, specifically the reliance on empirical cutoff parameters in PSEV versus the rigorous solution of the in-medium Schrödinger equation in QS.
2. The neglect of Pauli blocking shifts in the binding energies within the PSEV distribution functions, which artificially sustains bound states at higher densities.

The authors assert that while the PSEV approach offers a computationally simpler method to account for momentum-dependent Pauli blocking, it requires further improvements to reproduce the systematic results of the quantum-statistical approach. Specifically, it must account for medium-modified binding energies and continuum correlations.

Regarding the high $\alpha$-clustering observed in some heavy-ion collision (HIC) analyses, the authors caution that confirming abundant $\alpha$-clustering above saturation density requires more fundamental investigations. They suggest that experimental "excesses" of $\alpha$-particles could be attributed to feed-down from the decay of excited $A=5-8$ fragments or short-range correlations, rather than equilibrium $\alpha$-clustering. The paper emphasizes that local thermodynamic equilibrium properties are a prerequisite for understanding non-equilibrium HIC evolution, but current PSEV formulations may overestimate cluster survival in dense matter. Future work is needed to extend these equilibrium analyses to non-equilibrium conditions relevant to HIC experiments and astrophysical scenarios.