Source author record

James Shee

James Shee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.str-el physics.class-ph quant-ph

Catalog footprint

What is connected

5works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo

Phaseless Auxiliary-Field Quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium to large-sized molecules, due to its accuracy and favorable polynomial scaling with system size. Unfortunately the memory footprint of standard energy evaluation algorithms are non-trivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low rank structure of the Coulombic integrals have been successful, but are significantly limited by high prefactors, rendering the utility limited to very large systems. Here, we present a complementary, cubic scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to phaseless AFQMC. We show that the error due to this approximation, which we term Localized Orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable, and is computationally favorable even for small systems. We present results demonstrating a robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized orbital based methods, including the singlet-triplet gaps of polyacenes benzene through pentacene, the heats of formation for a set of platonic hydrocarbon cages, and the total energy of ferrocene (Fe(Cp)$_2$). Finally, we reproduce our previous result of the gas phase ionization energy of Ni(Cp)$_2$, agreeing with full ph-AFQMC to within statistical error while using less than a fifteenth of the computer time.

preprint2022arXiv

Say NO to Optimization: A Non-Orthogonal Quantum Eigensolver

A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures a significant portion of the exact wavefunction in a highly compact manner, and that allows computation of the resulting energies and wavefunctions at polynomial cost with a quantum computer. This enables a significant improvement over the corresponding classical non-orthogonal solver, which incurs an exponential cost when evaluating off-diagonal matrix elements between the ansatz states, and is therefore intractable. We implement the non-orthogonal quantum eigensolver (NOQE) here with an efficient ansatz parameterization inspired by classical quantum chemistry methods that succeed in capturing significant amounts of electronic correlation accurately. By taking advantage of classical methods for chemistry, NOQE provides a flexible, compact, and rigorous description of both static and dynamic electronic correlation, making it an attractive method for the calculation of electronic states of a wide range of molecular systems.

preprint2020arXiv

Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo

Transition metal complexes are ubiquitous in biology and chemical catalysis, yet they remain difficult to accurately describe with ab initio methods due to the presence of a large degree of dynamic electron correlation, and, in some cases, strong static correlation which results from a manifold of low-lying states. Progress has been hindered by a scarcity of high quality gas-phase experimental data, while exact ab initio predictions are usually computationally unaffordable due to the large size of the systems. In this work, we present a data set of 34 3d metal-containing complexes with gas-phase ligand-dissociation energies that have reported uncertainties of $\leq$ 2 kcal/mol. We perform all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multi-determinant trial wavefunctions selected by a blackbox procedure. We compare the results with those from DFT with various functionals, and DLPNO-CCSD(T). We find MAE of 1.09 $\pm$ 0.28 kcal/mol for our best ph-AFQMC method, vs 2.89 kcal/mol for DLPNO-CCSD(T) and 1.57 - 3.87 kcal/mol for DFT. We find maximum errors of 2.96 $\pm$ 1.71 kcal/mol for our best ph-AFQMC method, vs 9.15 kcal/mol for DLPNO-CCSD(T) and 5.98 - 13.69 kcal/mol for DFT. The reasonable performance of several functionals is in stark contrast to the much poorer accuracy previously demonstrated for diatomics, suggesting a moderation in electron correlation due to ligand coordination. However, the unpredictably large errors for a small subset of cases with both DFT and DLPNO-CCSD(T) leave cause for concern, especially due to the unreliability of common multi-reference indicators. In contrast, the robust and, in principle, systematically improvable results of ph-AFQMC for these realistic complexes establish it as a useful tool for elucidating the electronic structure of transition metal-containing complexes and predicting their gas-phase properties.

preprint2020arXiv

The Role of High-Order Electron Correlation Effects in a Model System for Non-valence Correlation-bound Anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H$_2$O)$_4$ cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree-Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger electron binding energies are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution, and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

preprint2016arXiv

Classical Mechanics in Hilbert Space: Path Integral Formulation, and a Quantum Correction

While it is well-known that quantum mechanics can be reformulated in terms of a path integral representation, it will be shown that such a formulation is also possible in the case of classical mechanics. From Koopman-von Neumann theory, which recasts classical mechanics in terms of a Hilbert space wherein the Liouville operator acts as the generator of motion, we derive a path integral representation of the classical propagator and suggest an efficient numerical implementation using fast fourier transform techniques. We then include a first quantum correction to derive a revealing expression for the semi-classical path integral, which augments the classical picture of a single trajectory through phase space with additional wave-like spreading.