Quantum Algorithm for Open-System Battery Cathodes by Modeling Multiple Strongly Coupled Holstein Polarons with Chain-Mapped Caldeira-Leggett Dynamics

Abstract: Cathode lithiation occupies a chemical regime of tightly localized orbitals, narrow bandwidths, and strong electron-lattice coupling. The defining electrochemical observables (open-circuit voltage and differential capacity) are open-system, reservoir-equilibration quantities that closed-Hamiltonian quantum simulation cannot produce, set by exchange with electron, Li$<sup>+$, and phonon baths. We present a fault-tolerant quantum algorithm that recovers them through a unitary chain-mapped Caldeira-Leggett embedding, rendering the baths Trotterizable. The resulting fourth-order Trotter step has a T-gate count polynomial in system size, validating its open-system dynamics against hierarchical equations of motion (HEOM) at strong coupling and the Lindblad limit at weak coupling. For single-carrier olivine LiFePO$_4$, a single voltage anchor on an otherwise DFT-fixed Hamiltonian places the differential-capacity peak within the $\pm5$ mV reproducibility of the experimental plateau. For multi-carrier spinel LiMn$_2$O$_4$, whose $1{:}1$ Mn$<sup>{3+}$/Mn$<sup>{4+}$ filling makes the inter-site Coulomb repulsion dynamically active, the same kernel yields a two-plateau voltage curve with a $125$ mV split, within $17\%$ of the observed $150$ mV. We deliver an end-to-end fault-tolerant resource estimate for such a multi-carrier, three-reservoir observable: $368$ logical qubits and $\sim3\times10<sup>5$ T-gates per step, or $\sim1.7\times10<sup>{12}$ T-gates for a full voltage curve (parallelizable over $\sim10<sup>3$ trajectories), leaving the production-scale dynamical run as a milestone for future hardware. The same kernel reproduces macroscopic quantum coherence, two-band superconductivity, and the Mikheyev-Smirnov-Wolfenstein resonance without modification, placing dynamical battery chemistry and similar Hamiltonians within scope for fault-tolerant quantum simulation.

• The paper presents a fault-tolerant quantum algorithm that models multi-polaron dynamics in battery cathodes using a chain-mapped Caldeira–Leggett embedding.
• It employs a palindromic fourth-order Suzuki decomposition and reversible arithmetic to efficiently simulate open-system dynamics with polynomial gate complexity.
• Resource analysis shows that simulating materials like LiFePO₄ and LiMn₂O₄ yields precise voltage predictions, making the approach feasible for future quantum hardware.

Quantum Simulation of Open-System Battery Cathode Dynamics via Multi-Polaron Chain-Mapped Caldeira–Leggett Embedding

Overview and Motivation

This work presents a fault-tolerant quantum algorithm for simulating electrochemical observables of battery cathodes, specifically targeting the regime of strong electron-lattice coupling, multiple charge carriers (Holstein polarons), and finite temperature—features central to commercial chemistries such as olivine LiFePO$_4$ and spinel LiMn$_2$O$_4$. The methodology addresses a critical open-system modeling gap: observables relevant to battery design (open-circuit voltage, differential capacity, impedance, etc.) fundamentally require equilibration with external electron, Li$^+$, and phonon reservoirs. Existing quantum simulation frameworks typically treat only single-carrier, weak-coupling, or closed-system limits, which are insufficient for the strongly correlated, open-system behaviors in real battery materials.

The central innovation is a unitary chain-mapped Caldeira–Leggett (CL) embedding of these reservoirs, implemented via the Tamascelli finite-temperature extension, enabling Trotterizable simulation of the full open-system dynamics on quantum circuits with polynomial gate complexity in system size. This enables, for the first time, rigorous resource estimates for dynamical, multi-carrier, three-reservoir observables at scale, including end-to-end counts for logical qubits and T-gates necessary for execution on fault-tolerant quantum hardware.

Hamiltonian Framework and Open-System Mapping

The cathode Hamiltonian incorporates carrier self-trapping, strong electron-phonon coupling (Holstein), inter-site Coulomb interactions (multi-polaron sector), and all relevant reservoirs, giving:

$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$

where all system parameters are taken from DFT or experimental data; $V_{\mathrm{wb}}$ is the Li washboard potential, $V_{\mathrm{LiP}}$ the screened polaron–Li interaction, and the phonon bath is Ohmic with a hard cutoff, as appropriate for cathode materials.

Coupling to reservoirs is executed by the chain-mapped CL transformation:

$$H_{\text{tot}} = H_S + \sum_k \tfrac{1}{2} (p_k^2 + \omega_k^2 q_k^2) + q \sum_k c_k q_k + q^2 \sum_k \frac{c_k^2}{2\omega_k^2}$$

Chain mapping (via the Stieltjes recurrence) renders the baths as 1D chains of modes, Trotterizable and amenable to digital quantum simulation, with finite-temperature effects captured by Tamascelli doubling (auxiliary chain for negative frequencies). The explicit inclusion of the CL counter-term ensures bath-induced static renormalization is canceled and the polaron binding energy remains invariant with bath discretization (chain length $K$).

Quantum Algorithm Design and Resource Analysis

The simulation algorithm, REQWIEM, is a first-quantized, product-formula approach—entirely within Clifford$+$T—using a palindromic (time-reversible) fourth-order Suzuki decomposition for real-time propagation. Each Trotter step involves:

• Grid-based representation of all continuous degrees of freedom (DVR grid for Li position, bosonic chain modes, system phonons)
• Jordan–Wigner encoding for fermionic carriers
• Locality-collapsed, polynomial on-diagonal multiplexing for potential energy (no $_2$0 channel scaling)
• Controlled rotation synthesis (Walsh–Hadamard/Gray transform via Möttönen UCR circuits)
• Non-polynomial diagonal terms (e.g., soft Coulomb) computed by reversible arithmetic with $_2$184-qubit ancilla workspace
• Chain-reservoir layers constructed from precomputed chain coefficients (analytic moment-matching and spectral density control)

The step complexity is polynomial in $_2$2 (number of active sites), and chain register sizes scale logarithmically with the relevant energy-to-frequency ratios ($_2$3), with all truncation errors (Fock tail, register aliasing, chain recurrence) tightly controlled and tabulated.

Model Validation and Benchmarking

Validation is performed at several levels:

• Internal Consistency & Operator Equivalence: Gate-level statevector norm equivalence between circuit and classical propagation ($_2$4) for all primitives up to $_2$5 qubits.
• Regime Benchmarking: For single-site/small-$_2$6 instances, chain-mapped open-system dynamics reproduce numerically exact HEOM reference results to $_2$72% RMS population error at strong coupling, where Lindblad descriptions fail. Multi-site, two-bath excitonic dimer benchmarks reach $_2$83.5% RMS agreement (see Figure 1).
• Observable Accuracy: Classically computed equilibrium profiles (chemical-potential sweeps) yield open-circuit voltage and differential capacity predictions in precise agreement with experiment (Yamada plateau), with the predicted $_2$9 peak within $_4$0 mV and FWHM 12.6 mV (experimental: $_4$120 mV) for LiFePO$_4$2.

Strong Numerical and Structural Claims

• Voltage Prediction Fidelity: For LiFePO$_4$6, after a single voltage anchor fixing $_4$7 at 3.45 V, the algorithm reproduces the differential capacity peak position to within 2.8 mV of experiment—inside the reproducibility band—and FWHM within the experimental uncertainty. Solid-solution shoulders and plateau extent are robust across DFT parameter uncertainty.
• Multi-Carrier Effects: For LiMn$_4$8O$_4$9 at half-filling, the multi-polaron, active-Coulomb sector captures a two-plateau structure (split 125 mV versus experimental 150 mV, $^+$0 error) derived from exact electrostatics with no fitted parameters.
• Rigorous Resource Estimates: For full-scale, fault-tolerant execution, LiFePO$^+$1 requires 376 logical qubits and $^+$2 T-gates per Trotter step. A full $^+$3 curve (with parallelization over $^+$4 measurement trajectories) totals $^+$5 T-gates, dominated by measurement cost. LiMn$^+$6O$^+$7 analogously requires 368 logical qubits and $^+$8 T-gates per two-plateau voltage curve (see Table below). These are exact emitted circuit counts.

System	Logical Qubits	T-gates/Step	T-gates/Curve
LiFePO$^+$9	376	$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$0	$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$1
LiMn$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$2O$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$3	368	$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$4	$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$5

Portability and Theoretical Implications

The same chain-mapped kernel naturally extends to other open-system Hamiltonians of the same algebraic structure. Demonstrations include:

This algebraic generality underscores the relevance of the framework to a broad class of correlated open-system problems beyond batteries—such as phonon-dressed conduction, two-band superconductivity, and quantum dissipative transitions.

Key Technical Innovations

• Chain-Mapped Open-System Embedding: Enables strong-coupling, non-Markovian, and multi-reservoir quantum simulation, overcoming the limitations of master-equation-based dilations or classical tensor network/HEOM approaches, which scale exponentially in system size, carrier number, or bath complexity.
• Locality Collapse in On-Diagonal Multiplexing: By exploiting the at-most-bilinear structure of polaron-bath and Coulomb interactions, the algorithm bypasses exponential channel-state multiplexing, achieving polynomial scaling with respect to system size.
• Reversible Quantum Arithmetic for Non-Polynomial Potentials: The only non-polynomial term, the soft-Coulomb Li–polaron interaction, is implemented with reversible $$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$7 arithmetic at T-counts competitive with conventional logic, without loss of accuracy.

Practical and Theoretical Implications

Practical

• Predictive Battery Materials Design: The approach enables first-principles, DFT-to-observable workflow for dynamical battery properties—open-circuit voltage, differential capacity, and impedance—essential for the discovery/optimization of new cathode chemistries (e.g., cobalt-free, manganese-rich).
• Resource Scaling Feasibility: The resource estimates (logical qubits, gate counts) bring materials-scale quantum simulation of open systems within the scope of anticipated future fault-tolerant hardware, at wall-clock times of order 7–15 hours per voltage curve with physically reasonable parallelization.

Theoretical

• Algorithmic Benchmark: The end-to-end resource estimate, including measurement cost and finite-temperature/dynamics, sets a milestone for future quantum hardware in simulating quantum materials far beyond what is classically tractable (tensor network, HEOM, ML-MCTDH breakdown).
• Criterion for Quantum Advantage: For this conjunctive regime—multi-carrier, strong coupling, finite temperature, three-reservoir open system—no existing classical algorithm offers full tractability; the presented approach offers conditional quantum advantage based on rigorous scaling arguments.

Outlook and Future Directions

• Hardware Realization: Execution of the full production trajectory for LiMn$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$8O$$H = -\tfrac{1}{2m_{\mathrm{Li}}} \partial_R^2 + V_{\mathrm{wb}}(R) + \sum_i \varepsilon_i n_i - t_{ij} \sum_{\langle i,j \rangle} (c_i^\dagger c_j + \mathrm{h.c.}) + \sum_i V_{\mathrm{LiP}}(R-x_i) n_i + \tfrac{1}{2} \sum_i \omega_{\mathrm{ph}}(q_i^2 + p_i^2) + g\sum_i n_i q_i + \sum_{i<j} U_{ij} n_in_j$$9 dynamical two-plateau behavior marks a critical milestone for practical quantum simulation of complex materials.
• Method Generalization: Extension to non-Ohmic baths, 3D carrier motion, and explicit solid–electrolyte interphase modeling are natural future directions.
• Portability: Direct relevance to numerous correlated materials problems—colossal magnetoresistive oxides, polaronic semiconductors, and open-system quantum coherence phenomena.

Conclusion

The manuscript delivers a comprehensive quantum algorithmic framework to simulate real, correlated battery materials in open-system, strong-coupling regimes. The approach merges algorithmic advances (chain mapping, locality collapse, reversible arithmetic) with rigorous validation and complete resource accounting, establishing a new benchmark for practical, predictive quantum simulations of material properties relevant to energy storage and beyond. The work substantiates a critical link between quantum simulation theory and realistic, high-impact applications in materials science.