Quantum-mechanical interatomic potentials with electron temperature for strong-coupling transition metals.

In narrow d-band transition metals, electron temperature T(el) can impact the underlying electronic structure for temperatures near and above melt, strongly coupling the ion- and electron-thermal degrees of freedom and producing T(el)-dependent interatomic forces. Starting from the Mermin formulation of density functional theory, we have extended first-principles generalized pseudopotential theory to finite electron temperature and then developed efficient T(el)-dependent model generalized pseudopotential theory interatomic potentials for a Mo prototype. Unlike potentials based on the T(el)=0 electronic structure, the T(el)-dependent model generalized pseudopotential theory potentials yield a high-pressure Mo melt curve consistent with density functional theory quantum simulations, as well as with dynamic experiments, and also support a rich polymorphism in the high-(T,P) phase diagram.

Prospects for release-node quantum Monte Carlo.

We perform release-node quantum Monte Carlo simulations on the first row diatomic molecules in order to assess how accurately their ground-state energies can be obtained. An analysis of the fermion-boson energy difference is shown to be strongly dependent on the nuclear charge, Z, which in turn determines the growth of variance of the release-node energy. It is possible to use maximum entropy analysis to extrapolate to ground-state energies only for the low Z elements. For the higher Z dimers beyond boron, the error growth is too large to allow accurate data for long enough imaginary times. Within the limit of our statistics we were able to estimate, in atomic units, the ground-state energy of Li(2) (-14.9947(1)), Be(2) (-29.3367(7)), and B(2)(-49.410(2)).

Systematic reduction of sign errors in many-body calculations of atoms and molecules.

The self-healing diffusion Monte Carlo algorithm (SHDMC) is shown to be an accurate and robust method for calculating the ground state of atoms and molecules. By direct comparison with accurate configuration interaction results for the oxygen atom, we show that SHDMC converges systematically towards the ground-state wave function. We present results for the challenging N2 molecule, where the binding energies obtained via both energy minimization and SHDMC are near chemical accuracy (1 kcal/mol). Moreover, we demonstrate that SHDMC is robust enough to find the nodal surface for systems at least as large as C20 starting from random coefficients. SHDMC is a linear-scaling method, in the degrees of freedom of the nodes, that systematically reduces the fermion sign problem.

Neutral and charged excitations in carbon fullerenes from first-principles many-body theories.

We investigate the accuracy of first-principles many-body theories at the nanoscale by comparing the low-energy excitations of the carbon fullerenes C(20), C(24), C(50), C(60), C(70), and C(80) with experiment. Properties are calculated via the GW-Bethe-Salpeter equation and diffusion quantum Monte Carlo methods. We critically compare these theories and assess their accuracy against available photoabsorption and photoelectron spectroscopy data. The first ionization potentials are consistently well reproduced and are similar for all the fullerenes and methods studied. The electron affinities and first triplet excitation energies show substantial method and geometry dependence. These results establish the validity of many-body theories as viable alternative to density-functional theory in describing electronic properties of confined carbon nanostructures. We find a correlation between energy gap and stability of fullerenes. We also find that the electron affinity of fullerenes is very high and size independent, which explains their tendency to form compounds with electron-donor cations.

Insulator to metal transition in fluid deuterium.

We have investigated the insulator to metal transition in fluid deuterium using first principles simulations. Both density functional and quantum Monte Carlo calculations indicate that the electronic energy gap of the liquid vanishes at about ninefold compression and 3000 K. At these conditions the computed conductivity values are characteristic of a poor metal. These findings are consistent with those of recent shock wave experiments but the computed conductivity is larger than the measured value. From our ab initio results we conclude that the transition is driven by molecular dissociation rather than disorder and that both temperature and pressure play a key role in determining structural changes in the fluid.

Quantum Monte Carlo study of the optical and diffusive properties of the vacancy defect in diamond.

Fixed-node diffusion quantum Monte Carlo (DMC) calculations of the ground and excited state energetics of the neutral vacancy defect in diamond are reported. The multiplet structure of the defect is modeled using guiding wave functions of the Slater-Jastrow type with symmetrized multideterminant Slater parts. For the ground state we obtain the 1E state in agreement with experiment. The calculated energy of the lowest dipole allowed transition is consistent with the experimentally observed GR1 band, which has long been identified with the neutral vacancy in diamond, although no previous first-principles ab initio calculation of this transition exists. The calculated multiplet splitting of over 2 eV indicates the importance of a proper treatment of electron exchange and correlation in this system. DMC calculations of the formation and migration energy of the vacancy defect are presented.

Quantum Monte Carlo calculations of nanostructure optical gaps: application to silicon quantum dots.

Quantum Monte Carlo (QMC) calculations of the optical gaps of silicon quantum dots ranging in size from 0 to 1.5 nm are presented. These QMC results are used to examine the accuracy of density functional (DFT) and empirical pseudopotential based calculations. The GW approximation combined with a solution of the Bethe-Salpeter equation performs well but is limited by its scaling with system size. Optical gaps predicted by DFT vary by 1-2 eV depending on choice of functional. Corrections introduced by the time dependent formalism are found to be minimal in these systems.