Cody Melton

Spin-orbit interactions in electronic structure quantum Monte Carlo methods

We develop generalization of the fixed-phase diffusion Monte Carlo method for Hamiltonians which explicitly depend on particle spins such as for spin-orbit interactions. The method is formulated in zero variance manner and is similar to treatment of nonlocal operators in commonly used static- spin calculations. Tests on atomic and molecular systems show that it is very accurate, on par with the fixed-node method. This opens electronic structure quantum Monte Carlo methods to a vast research area of quantum phenomena in which spin-related interactions play an important role.

QMCPACK: An open source ab initio quantum Monte Carlo package for the electronic structure of atoms, molecules and solids

QMCPACK is an open source quantum Monte Carlo package for ab initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wavefunctions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary-field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit and graphical processing unit systems. We detail the programs capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://qmcpack.org.

A new generation of effective core potentials for correlated calculations

We outline ideas on desired properties for a new generation of effective core potentials (ECPs) that will allow valence-only calculations to reach the full potential offered by recent advances in many-body wave function methods. The key improvements include consistent use of correlated methods throughout ECP constructions and improved transferability as required for an accurate description of molecular systems over a range of geometries. The guiding principle is the isospectrality of all-electron and ECP Hamiltonians for a subset of valence states. We illustrate these concepts on a few first- and second-row atoms (B, C, N, O, S), and we obtain higher accuracy in transferability than previous constructions while using semi-local ECPs with a small number of parameters. In addition, the constructed ECPs enable many-body calculations of valence properties with higher (or same) accuracy than their all-electron counterparts with uncorrelated cores. This implies that the ECPs include also some of the impacts of core-core and core-valence correlations on valence properties. The results open further prospects for ECP improvements and refinements.

Quantum Monte Carlo with variable spins

We investigate the inclusion of variable spins in electronic structure quantum Monte Carlo, with a focus on diffusion Monte Carlo with Hamiltonians that include spin-orbit interactions. Following our previous introduction of fixed-phase spin-orbit diffusion Monte Carlo, we thoroughly discuss the details of the method and elaborate upon its technicalities. We present a proof for an upper-bound property for complex nonlocal operators, which allows for the implementation of T-moves to ensure the variational property. We discuss the time step biases associated with our particular choice of spin representation. Applications of the method are also presented for atomic and molecular systems. We calculate the binding energies and geometry of the PbH and Sn2 molecules, as well as the electron affinities of the 6p row elements in close agreement with experiments.

A new generation of effective core potentials from correlated calculations: 2nd row elements

Very recently, we have introduced correlation consistent effective core potentials (ccECPs) derived from many-body approaches with the main target being their use in explicitly correlated methods, while still usable in mainstream approaches. The ccECPs are based on reproducing excitation energies for a subset of valence states, namely, achieving near-isospectrality between the original and pseudo Hamiltonians. In addition, binding curves of dimer molecules were used for refinement and overall improvement of transferability over a range of bond lengths. Here we apply similar ideas to the 2nd row elements and study several aspects of the constructions in order to find the high accuracy solutions within the chosen ccECP forms with 3s, 3p valence space (Ne-core). Our new constructions exhibit accurate low-lying atomic excitations and equilibrium molecular bonds (on average within ≈0.03 eV and 3 mÅ); however, the errors for Al and Si oxide molecules at short bond lengths are notably larger for both ours and existing effective core potentials. Assuming this limitation, our ccECPs show a systematic balance between the criteria of atomic spectra accuracy and transferability for molecular bonds. In order to provide another option with much higher uniform accuracy, we also construct He-core ccECPs for the whole 2nd row with typical discrepancies of ≈0.01 eV or smaller.

A new generation of effective core potentials from correlated calculations: 3d transition metal series

Recently, we have introduced a new generation of effective core potentials (ECPs) designed for accurate correlated calculations but equally useful for a broad variety of approaches. The guiding principle has been the isospectrality of all-electron and ECP Hamiltonians for a subset of valence many-body states using correlated, nearly-exact calculations. Here we present such ECPs for the 3d transition series Sc to Zn with Ne-core, i.e., with semi-core 3s and 3p electrons in the valence space. Besides genuine many-body accuracy, the operators are simple, being represented by a few gaussians per symmetry channel with resulting potentials that are bounded everywhere. The transferability is checked on selected molecular systems over a range of geometries. The ECPs show a high overall accuracy with valence spectral discrepancies typically ≈0.01-0.02 eV or better. They also reproduce binding curves of hydride and oxide molecules typically within 0.02-0.03 eV deviations over the full non-dissociation range of interatomic distances.

Fixed-node and fixed-phase approximations and their relationship to variable spins in quantum Monte Carlo

We compare the fixed-phase approximation with the better known, but closely related fixed-node approximation on several testing examples. We found that both approximations behave very similarly with the fixed-phase results being very close to the fixed-node method whenever nodes/phase were of high and comparable accuracy. The fixed-phase exhibited larger biases when the trial wave functions errors in the nodes/phase were intentionally driven to unrealistically large values. We also present a formalism that enables to describe wave functions with the full antisymmetry in spin-spatial degrees of freedom using our recently developed method for systems with spins as fully quantum variables. This opens new possibilities for simulations of fermionic systems in the fixed-phase approximation formalism.