Probability distribution of observables from a Bogoliubov vacuum projected onto good particle number: application to scission configurations of an actinide

This paper proposes and validates a method to compute complete probability distributions for fission observables, such as total kinetic energy, by sampling nucleonic configurations from a Bogoliubov vacuum projected onto good particle number, revealing that significant fluctuations in actinide scission are already captured within the mean-field framework.

The Gist

The Big Picture: Predicting the "Explosion" of an Atom

Imagine a giant, wobbling water balloon (the atomic nucleus) that is about to split in two. This is nuclear fission. When it splits, it creates two smaller balloons (fragments) that fly apart at incredible speeds.

Scientists have been very good at predicting the average behavior of this split. They know, on average, how heavy the pieces will be, how much energy they will release, and how fast they will go.

But here is the problem: Nature isn't just about averages. It's about the fluctuations. Sometimes the pieces fly apart a bit faster, sometimes slower. Sometimes they spin wildly; other times, they barely rotate. These tiny, random variations are crucial for understanding nuclear reactors and the universe, but they are incredibly hard to calculate because they depend on the chaotic dance of billions of tiny particles (protons and neutrons) inside the nucleus.

This paper introduces a new digital simulation method to predict these random fluctuations, specifically for the moment right before the nucleus snaps apart (called "scission").

---

The Old Way vs. The New Way

The Old Way (The "Cloud" View):
Traditionally, physicists treated the nucleus like a smooth, blurry cloud of charge. They calculated the "average density" of the protons and neutrons.

• The Analogy: Imagine looking at a swarm of bees from a mile away. You see a fuzzy, rotating cloud. You can guess the average wind speed of the swarm, but you can't tell if a specific bee is buzzing left or right.
• The Flaw: This smooth cloud view misses the "quantum jitter." It can't explain why the fragments sometimes spin or why their speeds vary so much.

The New Way (The "Individual Bee" View):
The authors developed a method to stop looking at the blurry cloud and start tracking every single bee (every single proton and neutron) individually, even though they are governed by quantum rules.

• The Analogy: Instead of a fuzzy cloud, they use a super-computer to simulate a specific snapshot where every bee has a specific location and a specific spin direction. They do this thousands of times to see all the possible ways the swarm could arrange itself.

How the Method Works: The "Digital Dice"

The nucleus is described by a complex mathematical object called a Bogoliubov vacuum. Think of this as a "recipe" that tells you the probability of finding a particle in a certain spot, but not the exact spot itself.

To get the actual picture, the authors use a technique called Markov Chain Monte Carlo (MCMC).

1. The Setup: Imagine you have a giant, empty room representing the nucleus.
2. The Game: You start by placing all the protons and neutrons randomly in the room.
3. The Rules: You have a "Game Master" (the math) that says, "If you move this neutron here, does it fit the recipe? Yes? Keep it. No? Put it back."
4. The Sampling: You repeat this millions of times, shuffling the particles around. Eventually, the arrangement of particles in the room perfectly matches the "recipe" (the quantum state).
5. The Result: You now have a list of millions of possible "snapshots" of the nucleus. For each snapshot, you can calculate exactly what the fragments would look like if the nucleus split right then.

What They Discovered: The "Neck" is the Key

They applied this method to a heavy atom, Californium-252, right at the moment it is about to split. They looked at two specific ways it splits (called "Standard II" and "Super-Long" modes).

Here are their big findings, explained simply:

1. The Shape Fluctuations are Small

The overall shape of the nucleus (how stretched out it is) doesn't change much from one "snapshot" to the next. It's like a stretched rubber band; it stays pretty consistent.

• Surprise: Even though the shape is stable, the energy released varies a lot. This means the shape isn't the main culprit for the energy variations.

2. The "Neck" is the Wild Card

As the nucleus stretches, a thin "neck" of matter connects the two future fragments.

• The Discovery: The number of particles in this neck is highly unstable. Sometimes there are 2 particles, sometimes 0, sometimes 1.
• The Metaphal: Imagine two people holding hands and pulling apart. The "handshake" (the neck) is the most fragile part. Sometimes they grip tight; sometimes they slip.

3. The Secret Source of Energy Variation

The paper found that the Total Kinetic Energy (TKE)—how fast the fragments fly apart—varies mostly because of the nuclear force in that thin neck.

• The Analogy: Think of the two fragments as magnets. The "neck" is a tiny bridge of magnetic material connecting them.
  • If the bridge is thick (lots of particles in the neck), the magnets hold on tight (strong attraction), and when they finally snap, they fly apart with less speed because they had to fight the pull longer.
  • If the bridge is thin or broken (few particles), the magnets snap apart instantly, flying away with more speed.
• The Conclusion: The random fluctuation of just a few particles in that tiny neck is responsible for most of the variation in how fast the fragments fly. The electric repulsion (Coulomb force) is too smooth to cause these wild swings; it's the messy, short-range nuclear force in the neck that does it.

4. The "Spin" Mystery

Fission fragments often spin like tops. But the math used to describe the nucleus (Mean-Field theory) usually predicts a non-spinning, symmetric shape. How do they spin?

• The Explanation: Even if the average shape is symmetric, the individual particles are jittering. This jitter creates tiny, random forces (torques) that twist the fragments as they separate.
• The Metaphor: Imagine two people on ice skates holding a rope. If the rope is perfectly straight, they slide apart. But if the rope is slightly twisted or if they pull at slightly different angles (due to the jitter of their hands), they will start to spin. The paper shows that the "jitter" of the particles in the neck creates these twisting forces.

Why Does This Matter?

1. Better Reactors: Understanding these fluctuations helps engineers design safer and more efficient nuclear reactors.
2. New Physics: It proves that you don't need to invent new physics to explain these variations; the "jitter" of the particles within the standard model is enough.
3. A New Tool: The authors released a free software tool called NucleoScope. It's like a "particle camera" that allows other scientists to take snapshots of nuclei and see the hidden quantum chaos that drives fission.

In a Nutshell

This paper is like upgrading from a weather forecast that only gives the "average temperature" to a simulation that shows you exactly how the wind gusts, raindrops, and clouds move. They found that the "gusts" (random particle movements) in the tiny neck of a splitting atom are the main reason why the explosion varies in speed and spin. It turns out, the chaos of a few particles holds the key to the energy of the whole system.

Technical Summary

1. Problem Statement

Nuclear fission dynamics are increasingly described using Nuclear Energy Density Functional (EDF) frameworks, specifically utilizing Bogoliubov vacua (Hartree-Fock-Bogoliubov states) to model the many-body system. While these mean-field approaches successfully predict average properties (like mass yields and energy balances), they face limitations in describing quantum fluctuations inherent in the fission process.

• The Gap: Current projection techniques allow for the calculation of probability distribution functions (PDFs) for one-body observables linked to spontaneous symmetry breaking (e.g., particle number, angular momentum, mass, and charge).
• The Challenge: Predicting the PDFs for many-body observables (such as total kinetic energy, inter-fragment Coulomb repulsion, or two-body interactions) has remained an unsolved problem. These observables are crucial for understanding the statistical spread of fission outcomes but are difficult to compute directly from a projected Bogoliubov vacuum without encountering the "sign problem" or prohibitive computational costs.

2. Methodology

The authors propose a Monte Carlo sampling method to generate nucleonic configurations consistent with a Bogoliubov vacuum projected onto good particle numbers. This allows for the direct estimation of PDFs for a broad class of observables.

A. Theoretical Framework

• State Representation: The system is described as a tensor product of neutron and proton gases, $|\psi\rangle = |\psi_n\rangle \otimes |\psi_p\rangle$, where each component is a Bogoliubov vacuum projected onto a fixed particle number ($N$ and $Z$).
• Configuration Space: The method samples configurations in the position/intrinsic-spin representation. A configuration is defined by the set of coordinates $\{ \mathbf{r}_1, \dots, \mathbf{r}_N \}$, where $\mathbf{r}_i = (x, y, z, \sigma)$.
• Probability Density: The probability of a specific configuration is proportional to the squared modulus of the wavefunction overlap: $p(\mathbf{r}_1 \dots \mathbf{r}_N) \propto |\langle \mathbf{r}_1 \dots \mathbf{r}_N | \psi \rangle|^2$.
• Determinant Formulation: For projected Bogoliubov states, this overlap squared can be expressed as the determinant of a skew-symmetric matrix $Z$, where elements depend on the canonical single-particle wavefunctions and the pairing amplitudes ($u_k, v_k$).

B. Sampling Algorithm (MCMC)

To handle the high dimensionality of the configuration space (e.g., $>600$ dimensions for an actinide), the authors employ a Markov Chain Monte Carlo (MCMC) approach using the Metropolis algorithm:

1. Initialization: Nucleons are placed in high-density regions of the one-body density with random spins.
2. Proposal Step: A candidate configuration is generated by perturbing one nucleon's position (Gaussian distribution with width $\sigma_{space}$) and/or flipping its spin (probability $P_{flip}$).
3. Acceptance: The move is accepted based on the ratio of the new probability to the old probability.
4. Burn-in and Thinning: To ensure convergence and reduce autocorrelation, a "burn-in" period is discarded, and configurations are recorded only every $N_{jump}$ steps.

C. Observable Estimation

The method is restricted to observables that are diagonal in the position/spin representation (e.g., particle number, multipole moments, local two-body interactions).

• For a diagonal observable $\hat{O}$, the expectation value is the average of the kernel $O(\mathbf{r}_1 \dots \mathbf{r}_N)$ over the sampled configurations.
• Crucially, the full PDF of the observable is obtained directly by histogramming the values of the kernel calculated for each sampled configuration.

3. Key Contributions

1. Novel Sampling Technique: Development of a method to sample nucleon positions and spins from a projected Bogoliubov vacuum, enabling the calculation of PDFs for many-body observables previously inaccessible in EDF frameworks.
2. Open-Source Implementation: The authors released the code NucleoScope (C++) to facilitate these calculations.
3. Validation: Rigorous benchmarking against deterministic calculations (in a harmonic oscillator basis) for light ($^{20}\text{Ne}$) and heavy ($^{252}\text{Cf}$) nuclei, demonstrating convergence and accuracy within statistical uncertainties.
4. Application to Scission: Application of the method to two specific scission configurations of $^{252}\text{Cf}$ (Standard II and Super-Long modes) to analyze fluctuations in geometry, interaction energies, and torques.

4. Key Results

The study was applied to two scission configurations of $^{252}\text{Cf}$: the Standard II (SII) mode (asymmetric) and the Super-Long (SL) mode (symmetric).

A. Geometric Fluctuations

• Compound System: The elongation ($\beta_{20}$) and octupole ($\beta_{30}$) moments of the compound nucleus show nearly Gaussian distributions with small relative fluctuations ($\sim 2\%$ for $\beta_{20}$).
• Neck Region: The number of particles in the neck ($Q_{neck}$) exhibits large fluctuations (relative fluctuation $\approx 100\%$) and a non-Gaussian distribution, often peaking at zero particles.
• Fragments: The shapes of the pre-fragments fluctuate by $0.05$–$0.1$ units. These shape fluctuations imply an excitation energy fluctuation of only $\sim 1\text{--}2$ MeV, which is insufficient to explain the full experimental fluctuation in prompt neutron multiplicity.

B. Inter-Fragment Interaction Energies

• Total Kinetic Energy (TKE) Fluctuation: The study found that a significant fraction of the experimental fluctuation in the Total Kinetic Energy (TKE) of fission fragments originates from the fluctuation of the inter-fragment interaction potential at the scission point.
• Nuclear vs. Coulomb: Surprisingly, the residual nuclear interaction (short-range) is the dominant source of this fluctuation, not the Coulomb repulsion.
  • The nuclear interaction fluctuation is strongly correlated ($\rho \approx 0.75$) with the number of particles in the neck.
  • The Coulomb interaction is less sensitive to neck particles ($\rho \approx 0.1$).
• Mechanism: The fluctuation arises because the positions of a few nucleons in the neck region vary significantly, altering the short-range nuclear force between the separating fragments.

C. Torques and Angular Momentum

• Residual Torques: The study calculated the PDF of torques acting on fragments due to inter-fragment interactions.
• Origin: Similar to energy, the nuclear interaction is the primary driver of torque fluctuations, while Coulomb torques follow a normal distribution.
• Implication: The residual torque is mostly oriented perpendicular to the fission axis. While the authors estimate this could contribute $\sim 8\hbar$ to the angular momentum, they note this is likely an overestimate due to classical assumptions, suggesting the mean-field picture already captures a significant portion of the spin generation.

5. Significance and Conclusion

• Mean-Field Sufficiency: The results suggest that a significant portion of the fluctuations in measured fission observables (specifically TKE and angular momentum) is already present within the mean-field picture when quantum fluctuations of nucleon positions are properly accounted for via projection.
• Neck Dynamics: The study highlights the critical role of the nuclear neck and the few nucleons residing there. Their stochastic positioning drives the variability in nuclear forces, which in turn dictates the final kinetic energy and spin of the fragments.
• Future Directions: The method provides a pedagogical tool to visualize the "nucleon cloud" and offers a pathway to study non-diagonal observables (like kinetic energy operators directly) by extending the approach to non-diagonal Monte Carlo techniques.

In summary, this work bridges the gap between static mean-field descriptions and the statistical nature of fission outcomes, demonstrating that quantum fluctuations in nucleon positions at scission are a primary driver of the observed variability in fission fragment properties.