We present a quantum-classical algorithm called LAS-QKSD for multireference systems, by combining a classical localized active space (LAS) fragment-based multireference algorithm with the quantum Krylov subspace diagonalization (QKSD) method for quantum computers. The algorithm uses wave-function information from a LAS self-consistent field (LASSCF) calculation to prepare an initial state with better overlap with the target ground state than the Hartree-Fock state. This is coupled with the use of QKSD to ultimately converge to the exact energy, providing faster convergence than starting from the Hartree-Fock state. Fragmentation has the twofold benefit of fewer configurations on the classical side of the algorithm as well as fewer state preparation gates on the quantum side. First, we compare the LAS-QKSD method to the classical LASSCF method and to QKSD with a Hartree-Fock initial state. We then examine ways to load the LASSCF wave function using direct initialization and a QKSD-motivated spectral filtering approach. Finally, using a bimetallic complex, we show that the LAS-QKSD method is an efficient alternative to highly expensive complete active space self-consistent field calculations on strongly correlated systems.
We present a quantum-classical algorithm called LAS-QKSD for multireference systems, by combining a classical localized active space (LAS) fragment-based multireference algorithm with the quantum Krylov subspace diagonalization (QKSD) method for quantum computers. The algorithm uses wave-function information from a LAS self-consistent field (LASSCF) calculation to prepare an initial state with better overlap with the target ground state than the Hartree-Fock state. This is coupled with the use of QKSD to ultimately converge to the exact energy, providing faster convergence than starting from the Hartree-Fock state. Fragmentation has the twofold benefit of fewer configurations on the classical side of the algorithm as well as fewer state preparation gates on the quantum side. First, we compare the LAS-QKSD method to the classical LASSCF method and to QKSD with a Hartree-Fock initial state. We then examine ways to load the LASSCF wave function using direct initialization and a QKSD-motivated spectral filtering approach. Finally, using a bimetallic complex, we show that the LAS-QKSD method is an efficient alternative to highly expensive complete active space self-consistent field calculations on strongly correlated systems.
We develop a hybrid quantum-classical workflow combining sample-based quantum diagonalization (SQD) and the localized active space self-consistent field method (LASSCF) to solve for the ground states of transition-metal complexes, a longstanding challenge for both classical and quantum algorithms. The resulting approach, named LASSQD, integrates quantum sampling with fragment-based multireference theory to reduce the computational cost of solving strongly correlated active spaces. We test LASSQD on multiple iron-based complexes and demonstrate that it agrees with LASSCF within 1 kcal/mol, albeit at a much reduced computational cost. The cost reduction originates from the use of a sparse approximation of the exact and combinatorially large ground-state wavefunction, which also enables LASSQD to treat fragment sizes that are computationally inaccessible to LASSCF, as demonstrated by our computation of the spin gap of iron-porphyrin. These results establish that LASSQD is a scalable strategy for generating reliable multireference wave functions, providing a robust starting point for post-SCF correlation methods that recover dynamic correlation beyond the active space, and a promising pathway toward quantum-enhanced electronic structure calculations.
Accurately and efficiently describing strongly correlated electronic systems is a central challenge in quantum computational chemistry, with classical and quantum computers. The localized active space self-consistent field method (LASSCF) uses a product of fragment active spaces as a variational space, with the Schrödinger equation solved exactly in each fragment and the fragment active-space orbitals defined in a self-consistent manner. LASSCF is accurate for systems with strong intra-fragment and weak inter-fragment correlation, and its computational cost is combinatorial with respect to the size of the individual fragment active spaces, rather than their product. However, exactly solving the Schrödinger equation in each fragment remains a substantial bottleneck. Here, we address the possibility of solving the fragment active space Schrödinger equation with approximate methods, particularly sample-based quantum diagonalization (SQD). SQD is a technique that uses a quantum computer to sample configurations from a chemically motivated quantum circuit and a classical computer to mitigate errors and solve the Schrödinger equation in a subspace of the configuration space. We apply the proposed method, LASSQD, to the [Fe(H2O)4]2bpym4+ compound and the [FeIIIFeIIIFeII(μ3-O)-(HCOO)6] complex for calculating the intermediate-spin ground state energies. We observe that LASSQD can tackle fragment sizes intractable by LASSCF, achieves within 1kcal/mol agreement to LASSCF, and delivers results that are competitive with alternative classical methods to solve the Schrödinger equation, and thus can be used as a starting point for a perturbative treatment (LASSQD-PDFT) to recover correlation external to the active space.
State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.
State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in trans-butadiene, and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits, that can be used for realistic chemical applications.
We introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2, (H 2 ) 4, and trans -butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.
Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today's quantum devices. Here we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground state energy. Our algorithm, termed "local active space unitary coupled cluster" (LAS-UCC), scales linearly with system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H2)2 into two H2 molecules and by breaking the two double bonds in trans-butadiene and resources estimates are provided for linear chains of up to 20 H2 molecules.