Prospects for Quantum Computation in Propellant Design: Assessing the Stability of Cyclic Ozone in Nanoscale Confinement

Technical Summary: Prospects for Quantum Computation in Propellant Design

Problem Statement
The paper addresses the challenge of designing advanced rocket propellants, specifically focusing on the potential of cyclic ozone ($O_3(D_{3h})$) as a high-energy-density additive. While cyclic ozone offers superior energy density compared to its bent isomer ($O_3(C_{2v})$), it is inherently unstable and has not been isolated in bulk. The authors propose that encapsulating cyclic ozone within fullerene matrices (endofullerenes) could stabilize the isomer through nanoscale confinement. However, synthesizing and characterizing these systems requires predicting reaction pathways, isomerization rates, and thermodynamic stability with high accuracy.

The core computational bottleneck is the electronic structure of these systems. Ozone exhibits strong non-adiabatic effects and pronounced multireference character, rendering standard classical quantum chemical methods (like Density Functional Theory or single-reference coupled cluster) insufficient for high-precision predictions. Furthermore, the system size (e.g., $(12 O_3)@C_{180}$) exceeds the tractable limits of classical Full Configuration Interaction (FCI), which scales exponentially with system size. The authors posit that fault-tolerant quantum computation (FTQC) is required to achieve the necessary accuracy (FCI-quality) for these large, strongly correlated systems.

Methodology
The authors present an end-to-end analysis of a computational workflow designed to assess the stability of cyclic ozone in fullerene cages. The workflow integrates classical and quantum components:

1. Classical Pre-processing:

   • Geometry Generation: Initial geometries are generated using sphere packing algorithms and relaxed via semi-empirical methods (XTB).
   • Molecular Dynamics: Finite-temperature configurations are sampled using ab initio molecular dynamics (MD) with tight-binding DFT and a Langevin thermostat.
   • Reaction Pathway: The isomerization pathway between bent and cyclic ozone is mapped using the Climbing Image Nudged Elastic Band (CI-NEB) method with DFT.
   • Variational Transition State Theory (VTST): The workflow uses VTST to estimate reaction rate constants. This requires high-accuracy energy corrections ($\Delta V$) at specific points along the reaction coordinate (reactants, products, transition state, and an auxiliary point) to correct lower-level DFT energies.

2. Quantum Algorithm Framework:

   • Algorithm: The study utilizes Qubitized Quantum Phase Estimation (QPE) to determine ground-state energies with FCI-level accuracy.
   • Hamiltonian Encodings: Two distinct encodings are evaluated to assess resource overheads:
     • Local Basis (Double-Factorized): Uses an active space approach with a double-factorized (DF) Hamiltonian encoding to reduce the number of terms in the electronic structure Hamiltonian.
     • Semilocal Basis (Dual Plane-Wave): Uses a Dual Plane-Wave (DPW) basis, treating the isolated molecule within a periodic supercell with vacuum padding.
   • Electrostatic Embedding: To reduce computational cost, the authors explore an embedding strategy where the fullerene cage is treated as a fixed electrostatic environment (point charges or frozen potential), allowing the quantum calculation to focus solely on the oxygen valence degrees of freedom.

3. Resource Estimation:

   • The authors employ three independent toolkits to estimate logical and physical resource requirements: Azure Quantum Resource Estimator (AzureQRE), PennyLane's resource estimation framework, and MIT Lincoln Laboratory's pyLIQTR.
   • Estimates are calculated for specific system sizes, including $(12 O_3)@C_{180}$ with varying active space sizes (from 2 to 192 orbitals) and $O_3@C_{60}$.
   • The analysis assumes a target energy accuracy of 1 mHa (below the standard 1.6 mHa chemical accuracy) and models fault-tolerant hardware using a surface code implementation with optimistic superconducting qubit parameters.

Key Results
The resource estimation results delineate the scale of hardware required for practical application:

• Logical Resources: For a large system like $(12 O_3)@C_{180}$ with a 192-orbital active space, the study estimates a requirement of approximately 4,330 logical qubits (AzureQRE) to 5,642 logical qubits (pyLIQTR with DPW basis).
• Gate Counts: The temporal overhead is substantial. The largest active space scenarios require between $2.6 \times 10^{13}$ and $1.2 \times 10^{14}$ T-gates.
• Physical Resources: Translating logical requirements to physical hardware (assuming a distance-$d=19$ surface code), the study estimates a need for approximately $6 \times 10^6$ physical qubits.
• Runtime: Under current optimistic assumptions, a single-shot calculation for a 64-orbital active space is projected to take roughly 13.43 days, while the largest active space (192 orbitals) would require approximately three years.
• Impact of Embedding: The electrostatic embedding strategy significantly reduces the number of correlated orbitals, thereby lowering the logical qubit footprint and gate complexity, though the absolute resource requirements remain high.
• Algorithmic Improvements: The authors note that state-of-the-art encodings (e.g., DFTHC-BLISS-SA) could potentially reduce T-gate counts by factors ranging from 7 to 195, though the systems remain intractable for classical methods.

Significance and Claims
The paper does not claim that cyclic ozone propellants are immediately realizable or that current quantum hardware can solve these problems. Instead, its significance lies in providing a realistic, end-to-end resource assessment for a practical, industrially relevant molecular engineering problem.

• Defining the Scale: The work establishes a "plausible scale" for fault-tolerant quantum computation in molecular design, moving beyond asymptotic complexity to report concrete logical and physical resource estimates.
• Benchmarking: By cross-checking estimates across three independent toolkits (AzureQRE, PennyLane, pyLIQTR), the authors provide a robust bracketing of realistic overheads, highlighting the gap between current capabilities and the requirements for FCI-quality calculations on large endofullerenes.
• Strategic Workflow: The paper advocates for a hybrid approach where quantum computation is used selectively (e.g., for specific points in a VTST workflow) to correct classical approximations, rather than attempting to simulate entire reaction surfaces solely with quantum resources.
• Feasibility Outlook: While the resource requirements are currently beyond reach, the authors argue that algorithmic improvements (such as tensor hypercontraction and compressed double-factorization) could reduce overhead by orders of magnitude, potentially bringing these calculations within the realm of "entry-level utility" for quantum computers in the future.

In summary, the paper serves as a feasibility study demonstrating that while quantum computation holds the theoretical potential to solve the electronic structure problems preventing the use of cyclic ozone propellants, realizing this potential requires significant advancements in both hardware scale and algorithmic efficiency.