Particle number projected energies at finite temperature

This study incorporates particle number projection into finite-temperature Skyrme density functional calculations to rigorously determine deformation-dependent energies, revealing that while even-odd staggering diminishes near the critical temperature and fission barriers remain similar to unprojected results, the method provides valuable insights into nuclear level densities at both ground states and barriers.

The Gist

Imagine a nuclear physicist trying to understand the behavior of an atomic nucleus not as a solid rock, but as a bustling, chaotic city of protons and neutrons. This city has a strict rule: it must have an exact number of citizens (particles). However, the standard tools physicists use to model these cities often break this rule, allowing the population to fluctuate wildly, like a city where people are constantly appearing and disappearing out of thin air.

This paper, titled "Particle number projected energies at finite temperature," introduces a new, more rigorous way to fix this counting error, especially when the nuclear city is hot and vibrating.

Here is a breakdown of the paper's concepts using everyday analogies:

1. The Problem: The "Leaky Bucket" of Physics

In the world of nuclear physics, scientists use a powerful tool called Density Functional Theory (DFT) to predict how heavy atoms behave. It's like a weather forecast for the nucleus.

• The Issue: To make the math work, the standard model treats the nucleus as if it's in a "grand canonical ensemble." Think of this like a bucket with a hole in the bottom. You can control the average amount of water (particles) inside, but the actual amount fluctuates.
• Why it matters: In reality, a specific nucleus (like Uranium) has a fixed number of protons and neutrons. It doesn't have a "fuzzy" population. When the nucleus gets hot (high temperature), these fluctuations get worse, and the model starts to lose important details, like the "odd-even staggering" effect (a subtle rhythm in how stable nuclei are based on whether they have an even or odd number of particles).

2. The Solution: The "Strict Bouncer" (Particle Number Projection)

The authors developed a method called Particle Number Projection (PNP).

• The Analogy: Imagine the standard model is a party where people wander in and out freely. The PNP method is like hiring a strict bouncer at the door. No matter how chaotic the party gets inside, the bouncer ensures that exactly 100 people are in the room at all times.
• The Innovation: While physicists have used this "bouncer" for cold, quiet nuclei (zero temperature), doing it for hot, vibrating nuclei is incredibly difficult. It's like trying to count exactly 100 people in a room while they are all dancing, jumping, and running around. This paper successfully figured out the math to make that count accurate even when the nucleus is "hot."

3. What Happens When the Nucleus Heats Up?

The researchers applied this method to heavy nuclei like Uranium-238 and the super-heavy Flerovium-292.

• The "Staggering" Disappears: At low temperatures, the nucleus acts like a superfluid (a frictionless fluid) where the even/odd particle count matters a lot. It's like a dance floor where couples (pairs) move in perfect sync. As the temperature rises, the heat breaks these couples apart. The "odd-even" rhythm fades away, and the particle distribution becomes a smooth, bell-shaped curve (Gaussian).
• The "Fission Barrier" (The Hill to Climb): One of the main goals was to see how hard it is for a nucleus to split apart (fission). This is measured by the "fission barrier"—a hill the nucleus must climb to break.
  • The Surprise: Even though the "exact" energy calculations changed significantly (the "bouncer" corrected the energy levels), the height of the hill (the fission barrier) didn't change much compared to the old, leaky-bucket model. This suggests that for predicting when a nucleus will split at high temperatures, the simpler, older models are actually "good enough," even if they aren't perfectly precise.

4. Counting the "States" (Level Density)

The paper also looked at Level Density, which is essentially counting how many different ways the nucleus can vibrate or arrange itself at a given energy.

• The Analogy: Think of a piano. At low energy, there are only a few notes you can play. As you add energy, the number of possible chords and melodies explodes.
• The Result: The authors compared their new "strict bouncer" method against an older approximation (the Discrete Gaussian method). They found that at low energies, the old method was wrong (it missed the subtle odd-even rhythms), but at high energies, both methods agreed. This gives scientists a better "map" of the nuclear landscape, which is crucial for predicting how long super-heavy elements survive before they decay.

Summary: Why Does This Matter?

This paper is a technical triumph in mathematical housekeeping.

1. It fixed the counting: It created a rigorous way to ensure nuclei have the exact right number of particles, even when they are hot and chaotic.
2. It validated the shortcuts: It showed that while the exact energy calculations are different, the big-picture predictions (like fission barriers) from older, simpler models are surprisingly robust at high temperatures.
3. It helps the future: By providing more accurate "level density" numbers, this work helps scientists better predict the stability of super-heavy elements (the ones at the very edge of the periodic table), which is essential for understanding how to create new elements in the lab.

In short, the authors built a better ruler to measure the heat and stability of the atomic world, confirming that while the details are complex, the big picture remains surprisingly stable.

Technical Summary

1. Problem Statement

Nuclear Density Functional Theory (DFT) is a powerful tool for describing heavy and superheavy nuclei, particularly regarding fission and fusion dynamics. However, standard DFT calculations (such as Hartree-Fock-Bogoliubov or HF+BCS) suffer from symmetry breaking at the mean-field level. Specifically:

• Gauge Symmetry Breaking: The BCS theory and Bogoliubov transformation break particle number conservation, leading to fluctuations in the particle number.
• Finite Temperature Complications: At finite temperatures, thermal excitations are typically described using the grand canonical ensemble, which inherently introduces particle number fluctuations. As temperature increases, pairing correlations and quantum effects are "washed out," and the lack of particle number conservation becomes a significant issue for finite nuclei.
• Limitations of Existing Methods: While particle number projection (PNP) is well-established at zero temperature, its implementation at finite temperature is complex. Previous finite-temperature PNP approaches often rely on approximations (e.g., saddle-point or discrete Gaussian approximations) or Monte Carlo methods, lacking exact calculations of projected energies within self-consistent DFT frameworks. Furthermore, calculating the entropy exactly via $Tr(\hat{D}_p \ln \hat{D}_p)$ is computationally prohibitive for heavy nuclei.

The core problem addressed is the lack of a rigorous, self-consistent framework to calculate exact particle-number-projected energies and free energies for compound nuclei at finite temperatures within Skyrme DFT.

2. Methodology

The authors developed a formalism to incorporate exact particle number projection into self-consistent Skyrme-Hartree-Fock+BCS calculations at finite temperatures.

• Theoretical Framework:

  • Projected Density Matrix: The system is described by a projected density matrix $\hat{D}_p = \frac{1}{Z_p} \hat{P} \hat{D} \hat{P}$, where $\hat{P}$ is the projection operator defined via an integral over the gauge angle $\theta$.
  • Exact Energy Calculation: The projected energy $E$ is calculated using the trace formula:
    $$E = \frac{1}{Z_p} \frac{1}{2\pi} \int d\theta e^{-i\theta N} \text{Tr}(e^{i\theta \hat{N}} \hat{D} \hat{H})$$
  • Generalized Wick's Theorem: The authors derived expressions for the single-particle density matrices ($\tilde{\rho}, \tilde{\lambda}, \tilde{\kappa}$) in terms of the gauge angle $\theta$, quasi-particle energies ($E_p$), and BCS occupation numbers ($u_p, v_p$). This allows the construction of the energy functional $E(\theta)$ without relying on the saddle-point approximation for the energy itself.
  • Entropy Approximation: Since exact entropy calculation is too costly, they adopted the method from Fanto et al., calculating entropy via the partition function: $S_p = \ln(Z_p) - \beta \frac{\partial \ln(Z_p)}{\partial \beta}$. This is noted to be accurate above the critical temperature ($T_c$) where the pairing gap vanishes.
  • Level Density: The level density $\rho_N(E^*)$ is calculated using the saddle-point approximation based on the projected partition function.

• Computational Implementation:

  • Code: Calculations were performed using the SkyrmeAx code in axial-symmetric coordinate spaces.
  • Interactions: The SkM* force was used for the particle-hole channel (suitable for large deformations), and a density-dependent $\delta$ interaction was used for the pairing channel.
  • Variational Approach: The study implements "Projection After Variation" (PAV). While "Variation After Projection" (VAP) is theoretically more accurate, it is practically impossible at finite temperature in self-consistent DFT due to the inability to apply Wick's theorem to the entropy matrix elements.

3. Key Contributions

1. Derivation of Exact PNP Formalism at Finite Temperature: The paper provides the first rigorous derivation and implementation of exact particle number projected energies for compound nuclei within a self-consistent Skyrme DFT framework at finite temperatures.
2. Analysis of Partition Function Ratios: The study analyzes the ratio of the projected partition function to the grand canonical partition function ($Z_p/\Xi$), demonstrating the transition from odd-even staggering at low temperatures to a Gaussian distribution at high temperatures.
3. Comparison of Approximations: It systematically compares exact PNP energies with approximate canonical energies (derived from partial derivatives of the partition function) and standard unprojected FT-BCS results.
4. Level Density Parameters: The authors extracted level density parameters ($a_0$ at the ground state and $a_f$ at the fission barrier) for $^{238}\text{U}$ and $^{292}\text{Fl}$, comparing the PNP method against the Discrete Gaussian (DG) approximation.

4. Results

• Partition Function & Staggering:

  • At low temperatures ($T < 0.4$ MeV), the system is dominated by even particle number components, showing clear odd-even staggering.
  • As temperature increases, the contribution of odd components rises, and the staggering diminishes.
  • Above $T \approx 0.5$ MeV, the distribution of particle number components becomes Gaussian-like, driven by thermal effects.

• Fission Barriers:

  • Energy Trends: Exact PNP energies are systematically lower (by $\sim 2.0$ MeV) than unprojected FT-BCS energies due to the inclusion of additional correlations.
  • Barrier Heights: Despite the shift in absolute energy, the fission barrier heights at high temperatures ($T=1.0$ MeV) are very similar between projected and unprojected calculations (5.47 MeV vs. 5.20 MeV). This suggests that unprojected FT-BCS is sufficient for estimating barrier heights at high temperatures.
  • Zero Temperature: At $T=0$, the difference is more pronounced (7.38 MeV with PNP vs. 7.69 MeV without), highlighting the importance of PNP when pairing correlations are strong.
  • Approximate Energies: Approximate canonical energies (Eq. 13) deviate significantly from exact PNP energies, particularly at the fission barrier, validating the necessity of exact calculations.

• Level Densities:

  • The calculated level densities for $^{238}\text{U}$ using the PNP method agree well with experimental data when collective enhancement factors are included.
  • The Discrete Gaussian (DG) method shows deviations at low excitation energies ($E_x < 1$ MeV) because it fails to account for the odd-even staggering induced by pairing interactions at low temperatures.
  • Parameters: The extracted level density parameter ratio $a_f/a_0$ (barrier to ground state) is found to be 1.0718 for $^{238}\text{U}$ and 1.0625 for $^{292}\text{Fl}$ at high excitation energies. These values are close to the commonly used empirical value of 1.1.
  • The ground state parameter $a_0$ is found to be approximately $A/9.97$ and $A/10.38$ for the respective nuclei, providing microscopic guidance for statistical models.

5. Significance

• Theoretical Advancement: This work bridges a critical gap in nuclear theory by enabling exact particle number projection in finite-temperature DFT, moving beyond saddle-point approximations for energy calculations.
• Superheavy Nuclei Synthesis: The results provide crucial inputs (specifically the level density parameters $a_0$ and $a_f/a_0$) for statistical models used to calculate the survival probabilities of superheavy compound nuclei. Accurate level densities are essential for predicting the cross-sections of synthesizing new elements.
• Validation of Approximations: The study clarifies the regimes where unprojected FT-BCS is sufficient (high-temperature fission barriers) versus where exact projection is necessary (low-temperature properties and precise barrier energetics).
• Methodological Utility: The derived formalism for projected densities and energies can be extended to calculate other observables at finite temperature, offering a robust tool for future investigations into nuclear dynamics and quantum computing applications involving pairing.