🔬 physics

Radical-Fragment Many-Body Expansion for Linear Alkane Quantum Chemistry

This paper introduces a radical-fragment many-body expansion (MBE2) that utilizes homolytic bond cleavage and restricted open-shell Hartree-Fock to reconstruct linear alkane energies from just four unique fragment calculations, achieving significant qubit reductions and demonstrating the viability of fragmentation-based quantum chemistry for scaling quantum solvers to large molecular systems.

Original authors: Daniel Sierra-Sosa, Jorge Saavedra, Santiago Solares, Gregorio Toscano-Pulido

Published 2026-07-01

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Daniel Sierra-Sosa, Jorge Saavedra, Santiago Solares, Gregorio Toscano-Pulido

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to solve a massive, intricate jigsaw puzzle representing a long chain of carbon atoms (a linear alkane, like a piece of plastic or wax). In the world of quantum chemistry, solving this puzzle means calculating exactly how the electrons behave.

The problem is that as the chain gets longer, the puzzle becomes impossibly huge. For a very long chain (hexacosane), the number of pieces required to describe it perfectly would need a supercomputer that doesn't exist yet, let alone a quantum computer. It's like trying to lift a skyscraper with a single pair of tweezers.

The Big Idea: Break it Down
The authors of this paper propose a clever strategy: instead of trying to lift the whole skyscraper at once, break it into small, manageable bricks. They call this the "Radical-Fragment Many-Body Expansion."

Here is how their specific method works, using simple analogies:

1. The "Radical" Cut (Homolytic Cleavage)
Most methods try to cut the molecule and then "cap" the broken ends with hydrogen atoms to make them look stable again (like putting a plastic cap on a broken pipe).

This paper's approach: They cut the molecule in half right down the middle of the bond, giving one electron to each side. This creates "radical" pieces—fragments with an unpaired electron, like a magnet with one pole sticking out.
Why it matters: They treat these open, un-capped pieces exactly as they are. They don't add fake atoms to hide the cut. They calculate the energy of these "wild" pieces in isolation.

2. The "LEGO" Assembly (Many-Body Expansion)
Once they have the energy of the small pieces, they don't just add them up. They realize that when two pieces snap together, there is a tiny bit of extra energy from the connection itself.

The Formula: Total Energy = (Energy of all individual pieces) + (Energy of the connections between neighbors).
The Magic Trick: Because these chains are so repetitive (like a row of identical LEGO bricks), the authors found they only need to calculate four unique types of pieces to solve the puzzle for any length of chain, from a short one to a very long one.
- Two types of single pieces (the ends and the middle).
- Two types of connected pairs (how the ends connect to the middle, and how middles connect to middles).

3. The Quantum Advantage
To solve the energy of these four unique pieces, they used quantum computers (specifically IBM's hardware).

The Old Way: To solve the whole long chain, you would need a quantum computer with 368 qubits (quantum bits). Current machines can't handle that.
The New Way: Because they only need to solve for the small pieces, the biggest calculation they ever need is just 30 qubits.
The Result: They reduced the required power by 12.3 times. It's like going from needing a fleet of trucks to deliver a package, to just needing a bicycle.

4. Testing the Theory
The team tested this on 11 different chains of varying lengths. They compared their "broken-down" method against:

Classical Computers: Using standard math (RHF and CCSD) to see if the "broken" method was accurate.
Quantum Solvers: Using three different quantum algorithms (VQE, ADAPT-VQE, and SQD) to solve the small pieces.

What They Found

Accuracy: The "broken-down" method was very accurate. The small errors introduced by breaking the chain were consistent and predictable. The quantum computers successfully solved the small pieces and reconstructed the total energy very close to the classical reference.
Scalability: No matter how long the chain got, the quantum computer never needed more than 30 qubits. The work didn't get harder; it just got longer, but since the pieces were identical, they only had to do the math four times.
Hardware Success: They successfully ran these calculations on real IBM quantum hardware, proving that this "break it down" strategy works on today's noisy, imperfect machines.

In Summary
This paper shows that you don't need a giant, perfect quantum computer to study huge molecules. Instead, you can cut the molecule into small, radical pieces, solve those tiny pieces on a small quantum computer, and then use a smart formula to snap the answers back together. It turns an impossible task into a manageable one, making quantum chemistry for large molecules possible on current technology.

Technical Summary: Radical-Fragment Many-Body Expansion for Linear Alkanes

Problem Statement

Linear alkanes ( $C_nH_{2n+2}$ ) serve as fundamental testbeds for electronic structure methods due to their repetitive $sp^3$ -hybridized carbon backbones. However, their computational cost scales combinatorially with system size. Under the Jordan–Wigner transformation with the STO-3G minimal basis, a full-molecule treatment requires $q = 14n + 4$ qubits. For hexacosane ( $C_{26}H_{54}$ ), this necessitates 368 qubits, a scale far exceeding current quantum hardware capabilities. Standard fragmentation methods, such as the Fragment Molecular Orbital (FMO) approach, typically rely on heterolytic bond cleavage with hydrogen capping atoms and electrostatic embedding, which introduces complexity in self-consistent charge iterations and closed-shell constraints.

Methodology

The authors propose a Radical-Fragment Many-Body Expansion at the two-body level (MBE2) tailored for linear alkanes. The core methodology diverges from standard FMO in three key aspects:

Homolytic Cleavage and Radical Fragments: Instead of heterolytic cleavage with capping atoms, the method employs homolytic C–C bond cleavage. This produces open-shell, electrically neutral radical fragments: the methyl radical ( $CH_3^\bullet$ ) and the methylene diradical ( $^\bullet CH_2^\bullet$ ). These fragments are treated as isolated systems using Restricted Open-Shell Hartree–Fock (ROHF), eliminating the need for electrostatic embedding or capping atoms.
Two-Body Assembly Formula: The total energy is reconstructed using a many-body expansion truncated at the two-body level (MBE2). The formula sums monomer energies and pairwise bonded-dimer corrections:
$E_{MBE2} = \sum_I E_I + \sum_{bonded I<J} \Delta E_{IJ}$
where $\Delta E_{IJ} = E_{IJ} - E_I - E_J$ .
Symmetry Exploitation: Due to the translational symmetry of linear alkanes, the decomposition requires only four unique fragment calculations regardless of chain length:
- Two monomer types: $CH_3^\bullet$ and $^\bullet CH_2^\bullet$ .
- Two bonded-dimer types: $CH_3-CH_2$ and $CH_2-CH_2$ .
  This reduces the computational scaling from $O(n)$ to $O(1)$ regarding the number of unique quantum calculations and fixes the maximum qubit requirement.

Quantum Solvers and Benchmarks:
The framework was benchmarked against five energy solvers across 11 linear alkanes (from butane to hexacosane):

Classical References: Restricted Hartree-Fock (RHF) and Coupled Cluster Singles and Doubles (CCSD).
Quantum Solvers:
- VQE: Variational Quantum Eigensolver using a hardware-efficient ansatz (EfficientSU2) on a statevector simulator.
- ADAPT-VQE: Adaptive Derivative-Assembled Pseudo-Trotter VQE using a Generalized Singles and Doubles (GSD) operator pool.
- SQD: Sample-based Quantum Diagonalization executed on IBM's ibm_pittsburgh processor, utilizing a Local Unitary Cluster Jastrow (LUCJ) ansatz initialized from CCSD amplitudes.

For quantum calculations, fragments were reduced to an 8-qubit active space (CAS(n,4)) via CASSCF and Cholesky decomposition, except for SQD, which operated on the full fragment Hamiltonians (up to 30 qubits) on real hardware.

Key Results

1. Qubit Reduction and Scaling

The MBE2 approach achieves a dramatic reduction in quantum resource requirements:

Qubit Count: For hexacosane ( $C_{26}H_{54}$ ), the qubit requirement drops from 368 (full molecule) to a maximum of 30 qubits (for the $CH_3-CH_2$ dimer), representing a 12.3× reduction.
Calculation Count: By exploiting symmetry, the number of unique calculations is reduced from 51 (brute-force fragmentation) to 4, a 12.8× reduction.
Fixed Cost: The maximum fragment size (30 qubits) and the number of unique calculations remain constant regardless of chain length, making the approach scalable to arbitrarily long chains without additional quantum resources.

2. Accuracy and Error Analysis

Fragmentation Error Dominance: The primary source of error in the MBE2 framework is the "bonded-pair-only" approximation (neglecting non-bonded pairs and three-body terms), not the quantum solver itself.
- MBE2-RHF and MBE2-CCSD show errors growing linearly with chain length, converging to a per-carbon penalty of ~14–15 kcal/mol.
- The quantum solvers closely track their respective classical MBE2 references.
Solver Performance:
- VQE: Tracks MBE2-RHF with high fidelity. The per-CH $_2$ deviation is only +0.13 kcal/mol relative to the ROHF baseline.
- SQD: Executed on real hardware, SQD tracks MBE2-CCSD with a hardware penalty of only +1.4 to +5.2 kcal/mol across the series, demonstrating viability on noisy intermediate-scale quantum (NISQ) devices.
- ADAPT-VQE: While converging closer to the exact CASSCF ground state for individual fragments, it exhibited an imbalance between monomer and dimer convergence. This led to an accumulation of error in the total energy for long chains (up to +457.5 kcal/mol for hexacosane), highlighting the sensitivity of the MBE2 assembly step to uniform convergence across fragment types.

Significance and Claims

The paper claims that the radical-fragment MBE2 approach renders the quantum chemistry of large molecular systems tractable on near-term quantum hardware. Its significance lies in:

Resource Efficiency: It decouples the quantum hardware requirements from the system size, reducing the qubit count for large alkanes by over an order of magnitude.
Chemical Naturalness: The use of homolytic cleavage to generate open-shell radicals is physically natural for non-polar C–C bonds and avoids the artificial charges or capping atoms required by heterolytic methods.
Solver Agnosticism: The framework is compatible with various quantum algorithms (VQE, ADAPT-VQE, SQD) and classical references, provided the underlying solvers can handle the reduced active spaces or full fragment Hamiltonians.
Scalability: The constant-size fragment advantage suggests the method can be extended to arbitrarily long alkane chains and potentially other polymer systems without increasing the maximum qubit requirement.

The authors conclude that while the current accuracy is bounded by the two-body approximation, the framework provides a viable pathway for scaling quantum solvers. Future work aims to incorporate non-bonded pair corrections and MBE3 terms, as well as applying the method to branched hydrocarbons and heteroatom-containing systems.