Quantum Simulation of Cranked Zirconium Isotopes: A Fixed-N Approach with a Structured Number-Conserving Ansatz

This paper presents a fixed-particle-number quantum simulation of cranked zirconium isotopes using a structured, number-conserving VQE ansatz that introduces a novel pairing-coherence diagnostic ($\Delta_{\mathrm{coh}}$) to analyze rotational evolution and pairing correlations in a truncated active space.

The Gist

Imagine you are trying to understand how a spinning top made of invisible, quantum Lego bricks behaves. Specifically, you want to know how the shape of this "top" (an atomic nucleus) changes as you spin it faster and faster.

This paper is a report from a team of physicists who used a quantum computer (simulated on a regular computer for now) to solve this puzzle for three specific types of Zirconium atoms (Zirconium-80, 82, and 84).

Here is the breakdown of their work using simple analogies:

1. The Problem: Spinning a Quantum Top

Atomic nuclei are like tiny, dense balls of protons and neutrons. Sometimes, they are perfectly round, but often they are shaped like rugby balls (prolate) or flat pancakes (oblate).

• The Challenge: When you spin a nucleus (like a figure skater pulling in their arms), it tries to align its internal parts with the spin. This creates a tug-of-war:
  • Deformation: The nucleus wants to keep its shape.
  • Pairing: The particles inside like to hold hands in pairs (like dance partners).
  • Spin: The rotation tries to break those pairs and line everyone up.
• The Difficulty: Calculating this tug-of-war is incredibly hard for normal computers because the number of ways the particles can arrange themselves grows explosively (like trying to count every possible way to shuffle a deck of cards, but with trillions of cards).

2. The Solution: A Special Quantum Recipe

The authors didn't just throw a generic quantum algorithm at the problem. They built a customized recipe (called an "ansatz") specifically for this job.

• The "Fixed-N" Rule: In many quantum simulations, scientists allow the number of particles to fluctuate to make the math easier. But in a real nucleus, the number of protons and neutrons is fixed. The authors insisted on keeping the particle count exact.
  • Analogy: Imagine a dance floor where you must always have exactly 10 couples. A standard simulation might let people wander off the floor or bring in strangers to make the math easier. This team said, "No! We must keep exactly 10 couples on the floor at all times."
• The "Structured" Approach: Instead of trying every possible move, they only allowed moves that actually happen in physics:
  • Pair Transfers: Moving a pair of dancers from one spot to another.
  • Coriolis Coupling: Moving a single dancer because the floor is spinning.
  • Why this matters: This is like a chef who only uses fresh, relevant ingredients rather than buying a whole warehouse of random spices. It makes the calculation faster and more accurate for the specific problem.

3. The New Tool: Measuring "Togetherness"

In traditional physics, scientists measure "superconductivity" (how well particles pair up) by looking for a specific signal. But because this team kept the particle count fixed, that old signal disappears (it becomes zero).

• The Innovation: They invented a new way to measure "pairing coherence."
  • Analogy: Imagine you are trying to see if a group of people are holding hands. The old method looked for a "handshake" that changes the number of people in the room (which is impossible here). The new method looks at how much the people lean toward each other and coordinate their movements, even if they don't break the "fixed number" rule. They call this $\Delta_{coh}$.

4. The Results: How the Zirconium Isotopes Behaved

They tested three versions of Zirconium (80, 82, and 84) and found distinct personalities for each:

• Zirconium-80 (The Flat Pancake):

  • No matter how fast they spun it, it stayed flat (oblate). It was stubborn and didn't want to change its shape.
  • Note: This contradicts some real-world experiments that suggest it might be rounder, but the authors admit their model is a simplified "stress test" and not a final answer yet.

• Zirconium-82 (The Shape-Shifter):

  • This one was the most dramatic. It started as a rugby ball (prolate) but as the spin speed increased, it suddenly flattened out and became a pancake. It showed the strongest reaction to the spinning.

• Zirconium-84 (The Stable Rugby Ball):

  • This one stayed a rugby ball the whole time. It was the most "cooperative" with its dance partners (neutrons), maintaining the strongest pairing even while spinning fast.

5. Why This Matters (The "So What?")

• It's a Proof of Concept: The authors aren't claiming their quantum computer is better than a supercomputer yet (they admit their simulation was run on a classical computer).
• The Real Win: They proved that you can build a quantum algorithm that respects the strict rules of nature (like keeping the particle count fixed) without breaking the math.
• Future Potential: As real quantum computers get bigger, this specific "recipe" will allow scientists to study heavy, complex nuclei that are currently impossible to simulate accurately.

Summary

Think of this paper as a blueprint for a new type of quantum microscope. The authors built a specialized lens that can look at spinning atomic nuclei without losing track of the particles. They tested it on three Zirconium atoms, found interesting differences in how they spin and pair up, and showed that this method is ready to be used on real quantum hardware in the future to solve even bigger nuclear mysteries.

Technical Summary

1. Problem Statement

The paper addresses the challenge of simulating the rotational response of open-shell atomic nuclei, specifically the neutron-deficient zirconium isotopes ($^{80,82,84}\text{Zr}$). These nuclei exhibit complex behavior governed by the competition between:

• Deformation: The shape of the nucleus (prolate vs. oblate).
• Pairing: Correlations between nucleons forming pairs.
• Single-particle alignment: The alignment of individual angular momenta with the rotational axis under "cranking" (rotation).

Traditional mean-field methods (like BCS or HFB) often break particle-number symmetry to handle pairing, which can be problematic for exact spectroscopic predictions. Conversely, full Configuration Interaction (CI) methods scale combinatorially and become intractable for larger active spaces. The authors aim to demonstrate a Variational Quantum Eigensolver (VQE) approach that:

1. Preserves exact particle number ($N$).
2. Handles the cranking term (Coriolis coupling) and genuine two-body pairing interactions simultaneously.
3. Operates within a controlled, truncated active space to demonstrate feasibility on quantum hardware (simulated here).

2. Methodology

A. Hamiltonian and Model

• Model: The study uses a Nilsson + Pairing Hamiltonian on a fixed deformation grid.
• Routhian: The many-body Routhian is defined as:
  $$\hat{H}' = \sum (\epsilon_k - \lambda_F) \hat{n}_k - G \sum P^\dagger_k P_l - \omega \hat{J}_x + \lambda_p(\hat{N} - N_{act})^2$$
  Where $\omega \hat{J}_x$ is the cranking term, $G$ is the pairing strength, and $\hat{J}_x$ matrix elements are reconstructed from deformation-dependent Nilsson eigenvectors.
• Active Space: A window of $M=8$ doubly degenerate Nilsson orbitals around the Fermi surface is selected. This maps to 16 qubits via the Jordan-Wigner transformation.

B. Structured Number-Conserving Ansatz

Instead of using generic hardware-efficient circuits or full Unitary Coupled Cluster (UCC), the authors propose a structured ansatz tailored to the physics of the problem:

• Reference State: A Hartree-Fock-like determinant with the lowest active pair orbitals occupied.
• Excitations:
  • Double Excitations: Restricted to pair-transfer operations ($P^\dagger_k P_l$). This directly implements the pairing interaction.
  • Single Excitations: Restricted strictly to the non-zero Coriolis coupling graph ($j_x$) of the active basis. This ensures the ansatz only explores excitations relevant to the cranking term.
• Parameter Count: For $M=8$, this results in 42 variational parameters (28 double + 14 single), significantly fewer than a generic UCCSD ansatz ($O(M^4)$).
• Key Feature: The ansatz exactly conserves particle number, avoiding the symmetry-breaking issues of BCS.

C. Pairing Coherence Diagnostic

Since the ansatz conserves particle number, the standard BCS pairing gap ($\Delta_\kappa \propto |\langle P_k \rangle|$) vanishes identically because $\langle P_k \rangle = 0$ in a fixed-$N$ sector.

• Solution: The authors introduce a new scalar diagnostic, $\Delta_{coh}$ (Pairing Coherence):
  $$\Delta_{coh} = G \sqrt{\sum_{k \neq l} |\langle P^\dagger_k P_l \rangle|}$$
• This quantity measures off-diagonal pair coherence (pair density matrix elements) and serves as a proxy for pairing correlations in fixed-$N$ systems without requiring symmetry breaking.

D. Computational Setup

• Solver: VQE optimized using the L-BFGS-B algorithm.
• Simulation: Statevector simulation (classical simulation of the quantum circuit) on a 16-qubit register.
• Baseline: A classical cranked-BCS calculation on the same mesh is used as a qualitative mean-field reference.

3. Key Results

The study analyzes $^{80}\text{Zr}$, $^{82}\text{Zr}$, and $^{84}\text{Zr}$ across a cranking frequency range ($0 \le \hbar\omega \le 1.0$ MeV).

• Shape Evolution:

  • $^{80}\text{Zr}$: Exhibits a stable oblate minimum ($\delta^* \approx -0.25$) throughout the cranking range. (Note: This contrasts with some experimental interpretations favoring prolate shapes, highlighting the model's sensitivity).
  • $^{82}\text{Zr}$: Shows the strongest rotational evolution. It starts prolate but softens to a spherical shape ($\delta^* = 0$) at high frequencies, coinciding with the largest angular momentum alignment.
  • $^{84}\text{Zr}$: Retains a robust prolate minimum ($\delta^* \approx +0.25$) throughout the scan.

• Rotational Alignment ($J_x$):

  • $^{82}\text{Zr}$ achieves the highest total alignment ($J_x \approx 19.38$ at $\hbar\omega=1.0$ MeV), driven by the added neutron pair.
  • $^{80}\text{Zr}$ shows symmetric proton-neutron alignment ($J_x \approx 12.02$).
  • $^{84}\text{Zr}$ shows moderate alignment ($J_x \approx 12.61$) with neutron dominance.

• Pairing Coherence ($\Delta_{coh}$):

  • The $\Delta_{coh}$ metric remains finite for all isotopes, successfully capturing the suppression of pairing with increasing rotation.
  • $^{84}\text{Zr}$ exhibits the strongest neutron pairing coherence, consistent with the stability of its prolate shape.
  • The fixed-$N$ quantum results show a "cleaner" pairing signal compared to the classical BCS baseline, which collapses more rapidly at high frequencies.

• Sensitivity Analysis:

  • Comparisons between $M=6$ and $M=8$ active spaces show stable qualitative trends but visible quantitative shifts. This confirms the results are not artifacts of the smallest truncation but acknowledges that full convergence is not yet claimed.

4. Key Contributions

1. Structured Ansatz Design: The development of a compact, physically motivated ansatz that combines pair-transfer doubles and Coriolis-driven singles. This reduces the parameter space while preserving exact particle number.
2. Fixed-$N$ Pairing Diagnostic: The introduction of $\Delta_{coh}$ as a viable observable for pairing correlations in number-conserving quantum simulations, solving the problem of the vanishing BCS gap in fixed-$N$ manifolds.
3. Methodological Framework: A demonstration that correlated quantum cranking can be formulated directly in a qubit register, retaining genuine two-body interactions and Coriolis couplings without mean-field decoupling.
4. Isotope Trends: Identification of distinct rotational and pairing behaviors across the Zr isotopic chain, specifically the unique stability of the oblate $^{80}\text{Zr}$ and the strong alignment of $^{82}\text{Zr}$ within the truncated model space.

5. Significance and Limitations

• Significance:

  • The work bridges the gap between nuclear structure physics and quantum computing by addressing a specific, difficult many-body problem (cranking + pairing) with a tailored quantum algorithm.
  • It highlights the scalability motivation: while the current 16-qubit simulation is not faster than classical mean-field methods, the exact fixed-$N$ Hilbert space dimension grows combinatorially ($\binom{2M}{M}$), reaching ~12,870 for $M=8$. This demonstrates the potential necessity of quantum resources for larger active spaces where classical exact diagonalization fails.
  • It provides a practical framework for analyzing pairing without symmetry breaking.

• Limitations:

  • Not a Spectroscopic Prediction: The results are presented as trends within a truncated active space and a specific Nilsson model, not as converged experimental predictions. The oblate minimum in $^{80}\text{Zr}$ contradicts some experimental data, emphasizing the model's sensitivity to the active space choice.
  • No Quantum Advantage: The calculations are performed via statevector simulation on classical hardware. The paper does not claim quantum advantage but rather validates the formulation for future hardware.
  • Convergence: Full convergence with respect to the active space size ($M$) is not achieved; $M=8$ is a representative truncation.

In conclusion, this paper serves as a methodological proof-of-concept for simulating rotating, pairing nuclei on quantum computers using a symmetry-preserving, structured ansatz, offering a new diagnostic tool ($\Delta_{coh}$) for pairing correlations in fixed-$N$ systems.