Improved quasiparticle nuclear Hamiltonians for quantum computing

This paper presents a systematic improvement to quasiparticle nuclear Hamiltonians for quantum computing by applying Brillouin-Wigner perturbation theory and a mean-field Hartree-Fock approximation to accurately describe open-shell nuclei in the $sd$ shell with minimal error, thereby enabling efficient quantum simulations on near-term devices.

The Gist

Imagine you are trying to predict the weather for a specific city. To do this perfectly, you would need to track every single air molecule, every drop of water, and how they all interact with each other. In the world of physics, this is what scientists try to do when they study the atomic nucleus (the core of an atom). They want to know exactly how protons and neutrons dance together.

The problem? There are so many particles, and they interact in such complex ways, that even the world's most powerful supercomputers get overwhelmed. It's like trying to solve a puzzle with a billion pieces where every piece changes shape every second.

This is where Quantum Computing comes in. Instead of using a regular computer, scientists use quantum computers, which are built to speak the same "language" as these tiny particles. However, there's a catch: translating the physics of the nucleus into a language a quantum computer can understand is like trying to translate a novel into a language that only has 10 words. You lose a lot of detail, or the translation becomes too long to fit on the computer.

The "Pairing" Shortcut

In this paper, the authors (Emanuele Costa and Javier Menéndez) are working on a specific translation method called Quasiparticle Pairing.

Think of the nucleus as a crowded dance floor.

• The Old Way: You try to track every single dancer individually. This is accurate but requires a massive amount of space (qubits) on the quantum computer.
• The Quasiparticle Way: You notice that most dancers are holding hands in pairs. So, instead of tracking 100 individuals, you just track 50 couples. This cuts the "space" needed in half and makes the dance floor much easier to manage.

This method works beautifully for "semimagic" nuclei (nuclei where the dance floor is perfectly organized). But for "open-shell" nuclei (where the dance floor is messy, with unpaired dancers running around), this shortcut fails. It's like assuming everyone is in a couple when, in reality, some people are dancing solo or in groups of three. The prediction becomes inaccurate.

The "Brillouin-Wigner" Fix

The authors wanted to keep the efficiency of the "couple" shortcut but fix the accuracy for the messy dance floors. They used a mathematical tool called Brillouin-Wigner (BW) Perturbation Theory.

The Analogy:
Imagine you are trying to predict the outcome of a game of chess.

1. The Shortcut: You only look at the pieces currently on the board.
2. The Problem: You miss the fact that a piece could move from the back row to attack you in three turns.
3. The Fix (BW Theory): You don't just look at the current board; you simulate all the possible future moves (virtual excitations) that aren't currently happening but could influence the game. You then fold the effect of these "ghost moves" back into your current board calculation.

In the paper, they use this to calculate how the "unpaired" dancers (protons and neutrons) affect the "paired" couples. This allows them to get a much more accurate picture without needing to track every single particle individually.

The "Hartree-Fock" Simplification

Here is the tricky part: Calculating all those "ghost moves" is still too hard for today's quantum computers (which are still in their "childhood" phase).

To solve this, the authors introduced a clever approximation called the Hartree-Fock (HF) Ansatz.

The Analogy:
Imagine you want to know the average temperature of a whole city.

• The Exact Way: Measure the temperature of every single house, street, and park. (Too hard).
• The HF Way: You pick one "average" house that represents the whole city. You assume the rest of the city behaves roughly like this one house. It's not perfect, but it's a very good guess that is easy to calculate.

By using this "average house" method to simplify the complex math, they created a new, simplified Hamiltonian (a set of rules for the quantum computer).

The Results

The authors tested this new method on a group of nuclei (the "sd shell," which includes elements like Neon, Magnesium, and Sulfur).

• The Old Shortcut: Made errors of up to 15% for messy nuclei.
• The New Method: Reduced the error to less than 2%.

This is a huge deal. It means they can now use near-future quantum computers to study complex, messy nuclei with high accuracy. They managed to keep the "shortcut" (saving space) while fixing the "accuracy" problem.

Why This Matters

This paper is like building a better bridge between the messy reality of atomic nuclei and the limited tools of today's quantum computers.

• For Scientists: It opens the door to simulating heavier, more complex atoms that were previously impossible to study with quantum methods.
• For the Future: It paves the way for understanding nuclear reactions, which could help in everything from developing better nuclear energy to understanding how stars explode.

In short, the authors took a rough, fast approximation, added a layer of smart math to fix its mistakes, and simplified it just enough so that today's experimental quantum computers can actually run it. It's a major step forward in making quantum physics practical.

Technical Summary

1. Problem Statement

Quantum computing offers a promising avenue for simulating nuclear structure, potentially overcoming the exponential scaling of classical diagonalization methods (Nuclear Shell Model, NSM). However, current approaches face significant hurdles:

• Encoding Overhead: Standard fermion-to-qubit mappings (e.g., Jordan-Wigner, Bravyi-Kitaev) require non-local operator strings, increasing circuit depth and gate counts.
• Quasiparticle Limitations: A more efficient encoding based on collective nucleon pairing modes (quasiparticles) reduces qubit requirements by half and yields local operators. However, this approximation is accurate only for semimagic nuclei (where pairing dominates). It fails for open-shell nuclei (where proton and neutron numbers are similar), as it neglects crucial proton-neutron correlations beyond the pairing channel, leading to energy errors exceeding 10%.

The Core Challenge: How to systematically improve the accuracy of the quasiparticle framework for open-shell nuclei without sacrificing its computational advantages (locality and reduced qubit count) for near-term quantum devices.

2. Methodology

The authors propose a systematic improvement using Brillouin-Wigner (BW) perturbation theory to derive an effective Hamiltonian within the quasiparticle subspace.

A. Theoretical Framework

• Partitioning: The full Hilbert space is partitioned into a Model Space ($Q$) spanned by quasiparticle pair states and a Complementary Space ($R$) containing virtual excitations (e.g., unpaired nucleons).
• Effective Hamiltonian: Using BW theory, the effect of the excluded space $R$ is folded into an energy-dependent effective Hamiltonian acting solely within $Q$:
  $$H_{\text{eff}}(E) = H_Q + H_{QR} (E I_R - H_{RR})^{-1} H_{RQ}$$
  Unlike Rayleigh-Schrödinger theory, the BW resolvent is evaluated at the exact eigenvalue $E$, requiring a self-consistent iterative solution.
• Iterative Solution: The resolvent is expanded as a geometric series. The method iterates between calculating the effective correction $\Delta H(E)$ and diagonalizing the resulting Hamiltonian until the energy converges.

B. Truncation for Quantum Hardware (HF-BW)

Direct implementation of the full BW method on quantum devices is infeasible because the perturbative expansion generates high-body interactions (up to $(N_n+Z_p)/2$-body terms), which drastically increase circuit depth. To make this suitable for Near-Term Quantum Devices (NISQ), the authors introduce two key approximations:

1. Two-Body Truncation: The method is restricted to the subspace of one quasiparticle-proton and one quasiparticle-neutron ($N_n=Z_p=2$). This limits the effective interaction to two-body terms.
2. Hartree-Fock (HF) Ansatz: Instead of diagonalizing the full $H_{RR}$ matrix (which is computationally expensive), the authors approximate the resolvent using a mean-field Hartree-Fock reference state. They retain only the ground-state pole of the spectral representation.
   • This replaces the complex correlated propagator with a simple mean-field energy denominator.
   • The resulting effective Hamiltonian ($H_{\text{eff, HF}}$) contains only local two-body quasiparticle interactions, preserving the qubit efficiency of the original mapping.

C. Pauli Blocking

To address the fact that the truncated method ignores the Pauli exclusion principle for the remaining nucleons, the authors introduce a statistical Pauli blocking factor. This factor estimates the probability that single-particle states involved in a transition are unoccupied, preventing unphysical transitions into already occupied states.

3. Key Contributions

1. Systematic Improvement via BW Theory: The first application of Brillouin-Wigner perturbation theory to systematically correct quasiparticle nuclear Hamiltonians for open-shell nuclei.
2. HF-Truncated Effective Hamiltonian: Development of a practical, mean-field-based truncation scheme that reduces the effective Hamiltonian to strictly two-body, local qubit operators. This bridges the gap between high-accuracy many-body theory and hardware constraints.
3. Pauli Blocking Correction: Introduction of a statistical blocking mechanism to mitigate errors arising from the truncation of the many-body space.
4. Benchmarking: Comprehensive validation against exact shell-model results for the $sd$-shell, demonstrating that the method recovers open-shell physics with high fidelity.

4. Results

The authors benchmarked their methods on even-even $sd$-shell nuclei (from $^{20}\text{Ne}$ to $^{36}\text{Ar}$) using the USDB Hamiltonian.

• Full BW Method (Classical):
  • Achieves a relative energy error ($\Delta r_e$) below 0.2% ($\sim 10^{-4}$) for most nuclei compared to the exact shell model.
  • Convergence is generally fast, though heavier or more deformed nuclei (e.g., $^{24}\text{Mg}$) require more iterations.
• Truncated HF-BW Method (Quantum-Suitable):
  • Accuracy: The ground-state energy errors are typically within 2% of the exact result. The average error is $\sim 1.6\%$.
  • Exceptions: Nuclei with strong non-pairing correlations (e.g., $^{22}\text{Ne}$, $^{24}\text{Mg}$) show slightly higher errors ($\sim 2.6\%$ and $5\%$ respectively), indicating a need for higher-order truncations for these specific cases.
  • Fidelity: The ground states of the truncated Hamiltonians have fidelities ($F$) exceeding 0.88 (average) with respect to the full BW reference state.
  • Operator Distance: The normalized operator distance between the truncated and full effective Hamiltonians is small ($< 3.2\%$), confirming the validity of the truncation.
• Importance of Pauli Blocking:
  • Methods without Pauli blocking showed large, systematic variational violations (errors up to -8.6%) and poor physical consistency.
  • Including Pauli blocking restored accuracy to within the 2% threshold, proving it is essential for physical meaningfulness.

5. Significance and Outlook

• Overcoming the "Open-Shell" Barrier: This work extends the applicability of efficient quasiparticle encodings from semimagic nuclei to the more complex open-shell systems, which are crucial for understanding nuclear deformation and structure.
• NISQ Readiness: By reducing the effective Hamiltonian to local two-body terms, the method makes high-accuracy nuclear simulations feasible on current and near-future quantum hardware (e.g., via VQE or QPE algorithms).
• Scalability: The approach offers a systematic path to improve accuracy (by increasing the truncation order or refining the HF ansatz) while maintaining polynomial classical preprocessing costs and linear qubit scaling.
• Future Directions: The authors suggest testing these effective Hamiltonians with advanced quantum algorithms (e.g., Unitary Coupled Cluster, ADAPT-VQE) and exploring further truncations (e.g., gadget Hamiltonians) to optimize resource usage on quantum devices.

In summary, the paper provides a robust theoretical and practical framework to enhance quasiparticle nuclear models, making them accurate enough for open-shell nuclei while remaining compatible with the constraints of near-term quantum computing.