Hugh G. A. Burton

It has recently been shown that configuration state functions (CSFs) with local orbitals can provide a compact reference state for low-spin open-shell electronic structures, such as antiferromagnetic states. However, optimizing a low-spin configuration using self-consistent field (SCF) theory has been a long-standing challenge since each orbital must be an eigenfunction of a different Fock operator. We introduce a low-spin restricted open-shell Hartree-Fock (ROHF) algorithm to optimize any CSF at mean-field cost. This algorithm employs quasi-Newton Riemannian optimization on the orbital constraint manifold to provide robust convergence, extending the geometric direct minimization approach to open-shell electronic structures with arbitrary genealogical spin coupling. Numerical calculations on transition metal aquo complexes show improved convergence over existing methodology, while the possibility of local CSF energy minima is demonstrated for iron-sulfur complexes. Finally, open-shell CSFs with different spin coupling patterns are used to qualitatively study the singlet ground state in polyacenes, revealing the onset of polyradical character as the chain length increases.

We present a new application of the generator coordinate method (GCM) as an electronic structure method for strong electron correlation in molecular systems. We identify spin fluctuations as an important generator coordinate responsible for strong static electron correlation that is associated with bond-breaking processes. Spin-constrained unrestricted HF (c-UHF) states are used to define a manifold of basis states for the Hill–Wheeler equations, which are discretized and solved as a non-orthogonal configuration interaction expansion. The method was tested on two-electron systems that are dominated by static and/or dynamic correlations. In a minimal basis set for H2, the resulting GCM quickly captures the ground-state full configuration interaction energy with just a few c-UHF states, whereas second-order perturbation theory on top of the GCM is needed to recover over 90% of the correlation energy in the cc-pVDZ basis set.

The performance of quantum algorithms for eigenvalue problems, such as computing Hamiltonian spectra, depends strongly on the overlap of the initial wave function and the target eigenvector. In a basis of Slater determinants, the representation of energy eigenstates of systems with N strongly correlated electrons requires a number of determinants that scales exponentially with N . On classical processors, this restricts simulations to systems where N is small. Here, we show that quantum computers can efficiently simulate strongly correlated molecular systems by directly encoding the dominant entanglement structure in the form of spin-coupled initial states. This avoids resorting to expensive classical or quantum state preparation heuristics and instead exploits symmetries in the wave function. We provide quantum circuits for deterministic preparation of a family of spin eigenfunctions with ( N N / 2 ) Slater determinants with depth O ( N ) and O ( N 2 ) local gates. Their use as highly entangled initial states in quantum algorithms reduces the total runtime of quantum phase estimation and related fault-tolerant methods by orders of magnitude. Furthermore, we assess the application of spin-coupled wave functions as initial states for a range of heuristic quantum algorithms, namely, the variational quantum eigensolver, adiabatic state preparation, and different versions of quantum subspace diagonalization (QSD) including QSD based on real-time-evolved states. We also propose a QSD algorithm that exploits states obtained through adaptive quantum eigensolvers. For all algorithms, we demonstrate that using spin-coupled initial states drastically reduces the quantum resources required to simulate strongly correlated ground and excited states. Our work provides a crucial component for enabling scalable quantum simulation of classically challenging electronic systems.

Density functional approximations can reduce the spin symmetry breaking observed for self-consistent field (SCF) solutions compared to Hartree-Fock theory, but the amount of exact Hartree-Fock (HF) exchange appears to be a key determinant in broken Ŝ² symmetry. To elucidate the precise role of exact exchange, we investigate the energy landscape of unrestricted Hartree-Fock and Kohn-Sham density functional theory for benzene and square cyclobutadiene, which provide paradigmatic examples of closed-shell and open-shell electronic structures, respectively. We find that increasing the amount of exact exchange leads to more local SCF minima, which can be characterized as combinatorial arrangements of unpaired electrons in the carbon π system. Furthermore, we studied the pathways connecting local minima to understand the relationships between different solutions. Our analysis reveals a subtle balance between one- and two-body interactions in determining SCF symmetry breaking, shedding new light on the physical driving forces for spin-symmetry-broken solutions in SCF approaches.

Approximating the electronic wave function for strongly correlated systems remains a major theoretical challenge. Emerging quantum computers can enable new types of wave-function ansatz to be considered, with the potential to overcome the exponential memory storage for strong correlation. I have recently introduced the tiled Unitary Product States (tUPS) ansatz, which successfully combines the preservation of particle-number and spin symmetry with shallow quantum circuits and local qubit connectivity [H. G. A. Burton, Phys. Rev. Res., 2024, 6, 023300]. In this contribution, I investigate the accuracy of this tUPS hierarchy for strongly-correlated Hamiltonians. I consider the picket-fence pairing Hamiltonian and the two-dimensional Hubbard lattice, which collectively describe a range of strong correlation mechanisms found in molecules. Numerical results demonstrate that highly accurate energies can be achieved with a compact approximation for both weak and strong correlation in the Hubbard model, and the repulsive pairing regime. These data provide valuable insights into the applicability of the tUPS hierarchy for strong electron correlation.

The ability of quantum computers to overcome the exponential memory scaling of many-body problems is expected to transform quantum chemistry. Quantum algorithms require accurate representations of electronic states on a quantum device, but current approximations struggle to combine chemical accuracy and gate efficiency while preserving physical symmetries, and rely on measurement-intensive adaptive methods that tailor the wave function Ansatz to each molecule. In this contribution, we present a symmetry-preserving and gate-efficient Ansatz that provides chemically accurate molecular energies with a well-defined circuit structure. Our approach exploits local qubit connectivity, orbital optimization, and connections with generalized valence bond theory to maximize the accuracy that is obtained with shallow quantum circuits. Numerical simulations for molecules with weak and strong electron correlation, including benzene, water, and the singlet-triplet gap in tetramethyleneethane, demonstrate that chemically accurate energies are achieved with as much as 84% fewer two-qubit gates compared to state-of-the-art adaptive Ansatz techniques.

State-specific complete active space self-consistent field (SS-CASSCF) theory has emerged as a promising route to accurately predict electronically excited energy surfaces away from molecular equilibria. However, its accuracy and practicality for chemical systems of photochemical interest have yet to be fully determined. We investigate the performance of the SS-CASSCF theory for the low-lying ground and excited states in the double bond rotation of ethylene. We show that state-specific approximations with a minimal (2e,2o) active space provide comparable accuracy to state-averaged calculations with much larger active spaces, while optimizing the orbitals for each excited state significantly improves the spatial diffusivity of the wave function. However, the incorrect ordering of state-specific solutions causes excited state solutions to coalesce and disappear, creating unphysical discontinuities in the potential energy surface. Our findings highlight the theoretical challenges that must be overcome to realize practical applications of state-specific electronic structure theory for computational photochemistry.