Publications

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.

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

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 wave functions 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 (CPU) and graphical processing unit (GPU) systems. We detail the program’s capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org.

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.