Screening strongly correlated materials in solid oxide cells with quantum simulation
Ouns El Harzli, Jelle Bor, Linde Wester, Leonid Zhukov, Hanl Park, Jean-François Bobier, Matt Langione, Antonino Marco Saitta, Jonathan Owens, Richard Hart, Yudong Cao, François Jamet
quant-phcond-mat.mtrl-scicond-mat.str-elSubmitted Oct 22, 2026• Announced Oct 23, 2026• 0 citations
We develop an end-to-end framework for assessing a quantum-computing application in industrial materials science: the screening of strongly correlated oxygen-electrode materials for SOFCs and SOECs. We formulate the materials problem in terms of interface observables including oxygen-vacancy formation energies, vacancy migration barriers, and electronic conductivity. We describe how these observables can be obtained using a diagrammatic GW-based embedding, either GW+DMFT or self-energy embedding theory (SEET), applied to a DFT-screened shortlist. The atomistic electrode/electrolyte interface defines the electronic structure problem from which the correlated impurity models are constructed. Quantum impurity solvers would provide the finite-temperature impurity Green's functions required by the embedding. The computational challenges that classical solvers may encounter as the interacting active space grows, together with the polynomial scaling of quantum Hamiltonian simulation for efficiently represented discretized models, open an opportunity for quantum advantage. We identify the benchmarks needed to assess this opportunity through convergence of vacancy and migration energetics and comparison with the best applicable classical solvers. We provide illustrative logical qubit counts and identify the thermal state preparation, simulation, measurement, and repeated-solve costs needed for a complete resource assessment. We connect microscopic observables to macroscopic engineering metrics using a continuum model for mixed ionic--electronic conducting oxygen electrodes that couples oxygen surface exchange, solid-state oxygen transport, electrode microstructure, and finite electrode thickness. Finally, we estimate the potential worldwide business value of improved electrode materials under SOFC and SOEC deployment scenarios informed by public technical targets and electrolyser deployment data.
Comments
52 pages and 4 figures, including 36 pages of main text
We develop an end-to-end framework for assessing a quantum-computing application in industrial materials science: the screening of strongly correlated oxygen-electrode materials for SOFCs and SOECs. We formulate the materials problem in terms of interface observables including oxygen-vacancy formation energies, vacancy migration barriers, and electronic conductivity. We describe how these observables can be obtained using a diagrammatic GW-based embedding, either GW+DMFT or self-energy embedding theory (SEET), applied to a DFT-screened shortlist. The atomistic electrode/electrolyte interface defines the electronic structure problem from which the correlated impurity models are constructed. Quantum impurity solvers would provide the finite-temperature impurity Green's functions required by the embedding. The computational challenges that classical solvers may encounter as the interacting active space grows, together with the polynomial scaling of quantum Hamiltonian simulation for efficiently represented discretized models, open an opportunity for quantum advantage. We identify the benchmarks needed to assess this opportunity through convergence of vacancy and migration energetics and comparison with the best applicable classical solvers. We provide illustrative logical qubit counts and identify the thermal state preparation, simulation, measurement, and repeated-solve costs needed for a complete resource assessment. We connect microscopic observables to macroscopic engineering metrics using a continuum model for mixed ionic--electronic conducting oxygen electrodes that couples oxygen surface exchange, solid-state oxygen transport, electrode microstructure, and finite electrode thickness. Finally, we estimate the potential worldwide business value of improved electrode materials under SOFC and SOEC deployment scenarios informed by public technical targets and electrolyser deployment data.
Quantum computers (QC) could harbor the potential to significantly advance materials simulations, particularly at the atomistic scale involving strongly correlated fermionic systems, where an accurate description of quantum many-body effects scales unfavorably with size. While a full-scale treatment of condensed matter systems with currently available noisy quantum computers remains elusive, quantum embedding schemes like dynamical mean-field theory (DMFT) allow the mapping of an effective, reduced subspace Hamiltonian to available devices to improve the accuracy of ab initio calculations such as density functional theory (DFT). Here, we report on the development of a hybrid quantum-classical DFT + DMFT simulation framework which relies on a quantum impurity solver based on the Lehmann representation of the impurity Green’s function. Hardware experiments with up to 14 qubits on the IBM Quantum system are conducted, using advanced error mitigation methods and a novel calibration scheme for an improved zero-noise extrapolation to effectively reduce adverse effects from inherent noise on current quantum devices. We showcase the utility of our quantum DFT + DMFT workflow by assessing the correlation effects on the electronic structure of a real material, Ca 2 CuO 2 Cl 2 , which is mapped to an effective single-band Hubbard Hamiltonian and the subsequently derived Anderson impurity model solved with up to 6 bath sites on available quantum hardware. Further, we carefully benchmark our quantum results with respect to exact reference solutions and experimental spectroscopy measurements. While challenges remain to scale our approach to larger, multi-orbital and multi-site systems with more bath sites, the present work marks an important milestone towards achieving utility-scale quantum computation in materials simulation.
Understanding emergent many-body phenomena in correlated materials remains one of the grandest challenges in solid-state quantum chemistry calculations. The individuality of electrons diminishes as correlation increases, which leads to the emergence of new features in the system. Seeking to address such problems accurately requires the aid of non-perturbative many-body techniques as well as the corresponding high-quality basis sets. In periodic systems, a given element may be present in different spatial arrangements displaying vastly different physical and chemical properties, and an elemental basis set independent of the physical properties of materials may lead to significant simulation inaccuracies. In fact, with the rapid progress of quantum chemistry methods in condensed-phase simulations, the need for a library of reliable Gaussian basis sets explicitly designed for periodic calculations has become urgent. In the meantime, the development of parameter-free, systematic and reliable quantum chemistry methods for simultaneously treating weak and strong correlations in periodic systems has never stopped. Their ability to treat realistic solid materials needs to be tested, and the obtained results can serve as benchmarks for empirical parameter-guided solid-state calculations. In this thesis, we present numerical studies of various solid-state systems, emphasizing the design and optimization of Gaussian basis sets and assessing the band structure properties of the periodic solid systems. The work begins by introducing the theoretical framework of the electronic structure theory and presenting the fundamentals of a fully self-consistent GW approximation and the self-energy embedding theory (SEET) in Chapters 2 and 3. It continues with a description of the Gaussian basis set optimization scheme we devised for solid-state quantum chemistry calculations. The scheme we present is designed to avoid a lack of material specificity within a given basis set by simultaneously minimizing the total energy of the system and optimizing the band energies when compared to the reference plane wave calculation while accounting for the overlap matrix condition number. We compare the quality of the Gaussian basis sets generated via our method against the existing basis sets. And the optimization scheme is tested to yield improved results. We finally present a quantitative study of electronic properties of realistic materials with the treatment of Green's function based weakly correlated self-consistent GW approximation and the strongly correlated self-energy embedding theory. In the transition metal oxide BiVO3, we found the inaccurate illustration of orbitals including t_{2g} in transition metal Vanadium by DFT. Our calculation shows a systematic trend of DFT failure in illustrating both insulating and metallic solutions, demonstrating the dangers of using DFT when the experimental data is scarce.
Strongly correlated electronic systems underpin a wide range of emerging quantum technologies, including solid-state defects for quantum information processing, molecular magnets, spin qubits, and quantum sensing platforms. Predictive computational design of these materials requires electronic structure methods capable of simultaneously describing strong electron correlation, excited-state phenomena, spin-dependent interactions, and environmental effects. While multireference wave function methods provide the necessary accuracy, their computational cost severely limits their application to realistic molecular and condensed-phase systems. This dissertation addresses this challenge through the development and application of quantum embedding and fragmentation methodologies that extend multireference electronic structure theory to increasingly complex quantum materials. The first part of this dissertation develops embedding-based approaches for describing optical properties of strongly correlated defects in extended solids. Complete Active Space Density Matrix Embedding Theory (CAS-DMET) and NEVPT2-DMET are applied to oxygen vacancy defects in magnesium oxide, where a systematic double-extrapolation protocol is introduced to approach both the non-embedding and thermodynamic limits. This framework enables quantitatively reliable calculations of absorption and emission energies while significantly reducing the computational cost of conventional multireference treatments of periodic systems. The second part extends embedding methodologies to the description of spin-dependent electronic structure. Density Matrix Embedding Pair-Density Functional Theory (DME-PDFT) is developed to recover dynamic electron correlation beyond the embedded active space, providing an efficient route for studying transition metal complexes and other strongly correlated molecular systems. Building upon these developments, CAS-DMET is further extended to investigate spin-phonon relaxation in molecular magnets and molecular crystals, demonstrating that quantum embedding can accurately describe spin dynamics and decoherence processes while substantially reducing computational cost. Beyond classical electronic structure calculations, this dissertation explores fragmentation strategies for hybrid quantum-classical algorithms through the development of efficient methodologies within the Localized Active Space Unitary Selective Coupled Cluster (LAS-USCC) framework. These developments improve the scalability of quantum chemistry algorithms for large active spaces and establish connections between embedding concepts and emerging quantum computing approaches. Finally, electronic structure calculations based on time-dependent density functional theory are employed to investigate fluorescent protein chromophores as potential biocompatible quantum sensors. These studies demonstrate how first-principles calculations can provide insight into spin-dependent properties in biologically relevant systems and expand the scope of computational approaches for quantum materials beyond traditional inorganic and molecular platforms. Collectively, the work presented in this dissertation establishes embedding and fragmentation as versatile frameworks for accurately describing the optical and spin properties of strongly correlated quantum systems across multiple length scales. By combining multireference electronic structure methods with physically motivated subsystem decompositions,these approaches bridge the gap between chemical accuracy and computational scalability, providing practical computational tools for the study and design of quantum materials relevant to quantum information science, molecular spintronics, and quantum sensing.
Corrosion presents a major challenge to the longevity and reliability of products across various industries, particularly in the aerospace sector. Corrosion arises from chemical processes occurring on an atomistic scale, which lead to macroscopic degradation. Addressing this issue requires multi-scale modeling approaches, which rely on microscopic parameters that are challenging to measure experimentally or model with conventional quantum chemistry techniques. In this work, we develop and demonstrate a hybrid quantum-classical workflow tailored for atomistic simulations of corrosion processes, with a specific focus on the initial step of the oxygen reduction reaction -- a critical trigger for the corrosion of aluminum alloys widely used in modern aircraft. Using a combination of classical quantum chemistry methods and quantum computing frameworks, we identify reaction geometries characterized by multi-configurational electronic structures that are ideal for exploring with quantum algorithms. For the first time in this context, we explore both noisy intermediate-scale quantum and fault-tolerant quantum algorithms for these multi-configurational system, integrating them within a workflow designed to bridge atomistic simulations with macroscopic modeling approaches, such as finite element methods. Furthermore, we conduct a detailed quantum resource estimation to assess when and how quantum computers may play a meaningful role in tackling these problems. Our results demonstrate that significant advancements in quantum hardware but also in algorithms and error correction techniques are needed to make quantum computation practically viable for this class of problems. Nevertheless, this work establishes a critical foundation for applying quantum computation to corrosion modeling and highlights its potential to address complex, business-relevant challenges in materials science.
Quantum-chemical simulations are essential for predicting energies of chemical reactions. Accurately solving the many-body Schrödinger equation for reagent and product states of most relevant chemical processes is, however, unfeasible. Quantum computing offers a pathway for predicting energies of correlated electronic systems with localized interactions. Here, we apply a quantum embedding approach for investigating oxygen reduction reactions at the electrode surface in Lithium batteries, a representative example of energetic analysis in localized chemical reactions. We employ an active space selection method based on density difference analysis for identifying the orbitals involved in the reaction. Leveraging the Local Unitary Cluster Jastrow ansatz for state preparation, the active-space orbitals are then processed on a quantum computer. As quantum algorithms, we use Sample-based Quantum Diagonalization, SQD, and its extended version, Ext-SQD, which integrates electronic excitations into the quantum-selected electronic configuration subspace. The largest configurations are represented by quantum circuits mapped onto 80 qubits of an IBM Heron R2 quantum processing unit. For up to 12 orbitals, we are able to benchmark the quantum-computed reaction energies against results obtained with Complete Active Space Configuration Interaction. For benchmarking results in active spaces as large as 32 orbitals, we resort to Heat-Bath Configuration Interaction and Coupled Cluster Singles and Doubles calculations, respectively. At 27 orbitals, the Ext-SQD results exhibit prediction accuracy improvements with regard to the standard, quantum-chemical reference methods that remain computationally feasible at that scale. The results indicate the potential of sample-based quantum diagonalization for performing high-accuracy reaction modeling in chemistry and materials science.
Quantitative descriptions of strongly correlated materials pose a considerable challenge in condensed matter physics and chemistry. A promising approach to address this problem is quantum embedding methods. In particular, the dynamical mean-field theory (DMFT) maps the original system to an effective quantum impurity model comprising correlated orbitals embedded in an electron bath. The biggest bottleneck in DMFT calculations is numerically solving the quantum impurity model, i.e., computing the Green's function. Past studies have proposed theoretical methods to compute the Green's function of a quantum impurity model in polynomial time using a quantum computer. So far, however, efficient methods for computing the imaginary-time Green's functions have not been established despite the advantages of the imaginary-time formulation. We propose a quantum-classical hybrid algorithm for computing imaginary-time Green's functions on quantum devices with limited hardware resources by applying the variational quantum simulation. Using a quantum circuit simulator, we verified this algorithm by computing Green's functions for a dimer model as well as a four-site impurity model obtained by DMFT calculations of the single-band Hubbard model, although our method can be applied to general imaginary-time correlation functions.