We extend correlated sampling to quantum-classical auxiliary-field quantum Monte Carlo (QC-AFQMC), enabling accurate nuclear force evaluation in strongly correlated systems. Computing forces via finite differences typically incurs prohibitive statistical noise in stochastic methods. We suppress this noise by maximizing correlation between geometries through synchronized random streams, orbital alignment, deterministic integral decomposition, and consistent classical shadow measurements. Crucially, a single shadow ensemble defined at the reference geometry suffices for all displaced structures, eliminating additional quantum measurements. This approach substantially reduces force variance while preserving accuracy. We validate the method on hydrogen chains across varying correlation regimes and demonstrate accurate forces for N 2 dissociation and stretched H 4 in strongly correlated regions where restricted coupled cluster methods fail qualitatively. Application to the MEA-CO 2 carbon capture reaction, integrating quantum information metrics for active space selection and matchgate shadows for overlap estimation, demonstrates that QC-AFQMC delivers accurate forces for complex reaction pathways in strongly correlated systems where conventional methods are unreliable or prohibitively expensive.
We extend correlated sampling to quantum-classical auxiliary-field quantum Monte Carlo (QC-AFQMC), enabling accurate nuclear force evaluation in strongly correlated systems. Computing forces via finite differences typically incurs prohibitive statistical noise in stochastic methods. We suppress this noise by maximizing correlation between geometries through synchronized random streams, orbital alignment, deterministic integral decomposition, and consistent classical shadow measurements. Crucially, a single shadow ensemble defined at the reference geometry suffices for all displaced structures, eliminating additional quantum measurements. This approach substantially reduces force variance while preserving accuracy. We validate the method on hydrogen chains across varying correlation regimes and demonstrate accurate forces for N 2 dissociation and stretched H 4 in strongly correlated regions where restricted coupled cluster methods fail qualitatively. Application to the MEA-CO 2 carbon capture reaction, integrating quantum information metrics for active space selection and matchgate shadows for overlap estimation, demonstrates that QC-AFQMC delivers accurate forces for complex reaction pathways in strongly correlated systems where conventional methods are unreliable or prohibitively expensive.
Using trial wave functions prepared on quantum devices to reduce the bias of auxiliary-field quantum Monte Carlo (QC-AFQMC) has established itself as a promising hybrid approach to the simulation of strongly correlated many body systems. Here, we further reduce the required quantum resources by decomposing the trial wave function into classical and quantum parts, respectively treated by classical and quantum devices, within the contextual subspace projection formalism. Importantly, we show that our algorithm is compatible with the recently developed matchgate shadow protocol for efficient overlap calculation in QC-AFQMC. Investigating the nitrogen dimer and the reductive decomposition of ethylene carbonate in lithium-based batteries, we observe that our method outperforms a number of established algorithm for ground state energy computations, while reaching chemical precision with less than half of the original number of qubits.
to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.
To make useful connections with experimental measurements, correlated electronic structure theories must accurately predict chemical properties in addition to energies. We present a finite-difference based algorithm to compute first-order response properties with phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) that relies on a branching correlated sampling approach. Focusing on electric dipole moments, we show that mean-field trial wave functions are sufficient to obtain high accuracy relative to CCSD(T) and experimental measurements for a set of 21 molecules, ranging in size from 2 to 18 atoms in their equilibrium geometries. As with energies, the quality of predicted dipole moments can be systematically improved with the use of correlated trial wave functions, even (or especially) in strongly correlated regimes. We show that the challenges faced by low-order perturbation theories in predicting the dipole moment of hydrogen fluoride across its dissociation coordinate are overcome with ph-AFQMC when using relatively simple trials. The key advantage of our approach over those previously reported for ph-AFQMC is its scalability to large system sizes with a phaseless bias no worse than that of a typical ground-state energy calculation; we routinely converge dipole moments for systems with more than 1000 basis functions.
High Resolution Image Download MS PowerPoint Slide Quantum Monte Carlo (QMC) methods are powerful approaches for solving electronic structure problems. Although they often provide high-accuracy solutions, the precision of most QMC methods is ultimately limited by the trial wave function that must be used. Recently, an approach has been demonstrated to allow the use of trial wave functions prepared on a quantum computer [Huggins et al., Unbiasing fermionic quantum Monte Carlo with a quantum computer. Nature 2022, 603, 416] in the auxiliary-field QMC (AFQMC) method using classical shadows to estimate the required overlaps. However, this approach has an exponential post-processing step to construct these overlaps when performing classical shadows obtained using random Clifford circuits. Here, we study an approach to avoid this exponential scaling step by using a fixed-node Monte Carlo method based on full configuration interaction quantum Monte Carlo. This method is applied to the local unitary cluster Jastrow ansatz. We consider H 4, ferrocene, and benzene molecules using up to 12 qubits as examples. Circuits are compiled to native gates for typical near-term architectures, and we assess the impact of circuit-level depolarizing noise on the method. We also provide a comparison of AFQMC and fixed-node approaches, demonstrating that AFQMC is more robust to errors, although extrapolations of the fixed-node energy reduce this discrepancy. Although the method can be used to reach chemical accuracy, the sampling cost to achieve this is high even for small active spaces, suggesting caution about the prospect of outperforming conventional QMC approaches.
We introduce QC-CBT-AFQMC, a hybrid algorithm that incorporates computational basis tomography (CBT) into the quantum-classical auxiliary field quantum Monte Carlo (QC-AFQMC) method proposed by Huggins [], replacing the use of classical shadows. While the original QC-AFQMC showed high accuracy for quantum chemistry calculations, it required exponentially costly postprocessing. Subsequent work using Matchgate shadows [Wan ] improved scalability, but still suffers from prohibitive computational requirements that limit practical applications. Our QC-CBT-AFQMC approach uses shallow Clifford circuits with a quadratic reduction of two-qubit gates over the original algorithm, significantly reducing computational requirements and enabling accurate calculations under limited measurement budgets. We demonstrate its effectiveness on the hydroxyl radical, ethylene, and nitrogen molecule, producing potential energy curves that closely match established benchmarks. We also examine the influence of CBT measurement counts on accuracy, showing that subtracting the active space AFQMC energy mitigates measurement-induced errors. Furthermore, we apply QC-CBT-AFQMC to estimate reaction barriers in [ 3 + 2 ] cycloaddition reactions, achieving agreement with high-level references and successfully incorporating complete basis set extrapolation techniques. These results highlight QC-CBT-AFQMC as a practical quantum-classical hybrid method that bridges the capabilities of quantum devices and accurate chemical simulations.
Quantum-classical auxiliary-field quantum Monte Carlo (QC-AFQMC) is a promising route to using quantum computers for strongly correlated electronic structure problems, but its scalability is limited by a classical post-processing bottleneck: the measured quantum trial state must be used to evaluate overlaps, force biases, and local energies for many walkers throughout the AFQMC propagation. Existing matchgate-shadow approaches avoid the exponential cost of generic classical shadows, but their post-processing remains expensive for large systems. In this work, we show that particle-number-preserving fermionic shadows, implemented as orbital-rotation shadows, are able to drastically reduce QC-AFQMC post-processing costs. We derive efficient estimators for overlaps and their derivatives, giving O(N3) scaling for overlaps and force biases and O(N4) scaling for local energies per shadow and walker and with O(N) Cholesky vector. For fixed normalized target determinants, we also prove that the overlap estimator is unbiased with variance upper bounded by O(η), where η is the particle number, yielding a sufficient budget of O(η/ε2) ideal independent measurements for additive overlap error ε. We combine these estimators with a novel virtual-correlation-energy treatment and a GPU implementation and demonstrate QC-AFQMC using orbital-rotation shadows. Performance projections for π-conjugated molecules up to pentacene with a (22e, 22o) active space and 356 total orbitals show that the full classical post-processing workflow can be completed in slightly more than two days on 8×A100 graphics processing units (GPUs). These results significantly reduce classical post-processing costs for QC-AFQMC and motivate further improvements in preparation and measurement of the trial state on quantum hardware.