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.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The scale and complexity of the quantum system to which real-space quantum Monte Carlo (QMC) can be applied in part depends on the representation and memory usage of the trial wavefunction. B-splines, the computationally most efficient basis set, can have memory requirements exceeding the capacity of a single computational node. This situation has traditionally forced a difficult choice of either using slow internode communication or a potentially less accurate but smaller basis set such as Gaussians. Here, we introduce a hybrid representation of the single particle orbitals that combine a localized atomic basis set around atomic cores and B-splines in the interstitial regions to reduce the memory usage while retaining the high speed of evaluation and either retaining or increasing overall accuracy. We present a benchmark calculation for NiO demonstrating a superior accuracy while using only one eighth of the memory required for conventional B-splines. The hybrid orbital representation therefore expands the overall range of systems that can be practically studied with QMC.

In this work, we study the interlayer interactions between sheets of blue phosphorus with quantum Monte Carlo (QMC) methods. We find that as previously observed in black phosphorus, interlayer binding of blue phosphorus cannot be described by van der Waals (vdW) interactions alone within the density functional theory framework. Specifically, while some vdW density functionals produced reasonable binding curves, none of them could provide a correct, even qualitatively, description of charge redistribution due to interlayer binding. We also show that small systematic errors in common practice QMC calculations, such as the choice of optimized geometry and finite-size corrections, are non-negligible given the energy and length scales of this problem. We mitigate some of the major sources of error and report QMC-optimized lattice constant, stacking, and interlayer binding energy for blue phosphorus. It is strongly suggested that these considerations are important and quite general in the modeling of two-dimensional phosphorus allotropes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Forsterite (Mg2SiO4) single crystals were shock compressed to pressures between 200 and 950 GPa using independent plate-impact steady shocks and laser-driven decaying shock compression experiments. Additionally, we performed density functional theory-based molecular dynamics to aid interpretation of the experimental data and to investigate possible phase transformations and phase separations along the Hugoniot. We show that the experimentally obtained Hugoniot cannot distinguish between a pure liquid Mg2SiO4 and an assemblage of solid MgO plus liquid magnesium silicate. The measured reflectivity is nonzero and increases with pressure, which implies that the liquid is a poor electrical conductor at low pressures and that the conductivity increases with pressure.

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.

Abstract not provided.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.