Publications

Abstract not provided.

Outstanding problems in the high-pressure phase diagram of hydrogen have demonstrated the need for more accurate ab initio methods for thermodynamic sampling. One promising method that has been deployed extensively above 100 GPa is coupled electron-ion Monte Carlo (CEIMC), which treats the electronic structure with quantum Monte Carlo (QMC). However, CEIMC predictions of the deuterium principal Hugoniot disagree significantly with experiment, overshooting the experimentally determined peak compression density by 7% and lower temperature gas-gun data by well over 20%. By deriving an equation relating the predicted Hugoniot density to underlying equation of state errors, we show that QMC and many-body methods can easily spoil the error cancellation properties inherent in the Rankine-Hugoniot relation, and very likely suffer from error addition. By cross validating QMC based on systematically improvable trial functions against post-Hartree-Fock many-body methods, we find that these methods introduce errors of the right sign and magnitude to account for much of the observed discrepancy between CEIMC and experiment. We stress that this is not just a CEIMC problem, but that thermodynamic sampling based on other many-body methods is likely to experience similar difficulties.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.