Luke Shulenburger

Hybrid algorithms in quantum Monte Carlo

With advances in algorithms and growing computing powers, quantum Monte Carlo (QMC) methods have become a leading contender for high accuracy calculations for the electronic structure of realistic systems. The performance gain on recent HPC systems is largely driven by increasing parallelism: the number of compute cores of a SMP and the number of SMPs have been going up, as the Top500 list attests. However, the available memory as well as the communication and memory bandwidth per element has not kept pace with the increasing parallelism. This severely limits the applicability of QMC and the problem size it can handle. OpenMP/MPI hybrid programming provides applications with simple but effective solutions to overcome efficiency and scalability bottlenecks on large-scale clusters based on multi/many-core SMPs. We discuss the design and implementation of hybrid methods in QMCPACK and analyze its performance on current HPC platforms characterized by various memory and communication hierarchies.

Quantum Monte Carlo applied to solids

We apply diffusion quantum Monte Carlo to a broad set of solids, benchmarking the method by comparing bulk structural properties (equilibrium volume and bulk modulus) to experiment and density functional theory (DFT) based theories. The test set includes materials with many different types of binding including ionic, metallic, covalent, and van der Waals. We show that, on average, the accuracy is comparable to or better than that of DFT when using the new generation of functionals, including one hybrid functional and two dispersion corrected functionals. The excellent performance of quantum Monte Carlo on solids is promising for its application to heterogeneous systems and high-pressure/high-density conditions. Important to the results here is the application of a consistent procedure with regards to the several approximations that are made, such as finite-size corrections and pseudopotential approximations. This test set allows for any improvements in these methods to be judged in a systematic way.

Application of Diffusion Monte Carlo to Materials Dominated by van der Waals Interactions.

van der Waals forces are notoriously difficult to account for from first principles. We have performed extensive calculations to assess the usefulness and validity of diffusion quantum Monte Carlo when predicting van der Waals forces. We present converged results for noble gas solids and clusters, archetypical van der Waals dominated systems, as well as the highly relevant π-π stacking supramolecular complex: DNA + intercalating anticancer drug ellipticine. Analysis of the calculated binding energies underscores the existence of significant interatomic many-body contributions.

X-ray Thomson Scattering in Warm Dense Matter without the Chihara Decomposition.

X-ray Thomson scattering is an important experimental technique used to measure the temperature, ionization state, structure, and density of warm dense matter (WDM). The fundamental property probed in these experiments is the electronic dynamic structure factor. In most models, this is decomposed into three terms [J. Chihara, J. Phys. F 17, 295 (1987)] representing the response of tightly bound, loosely bound, and free electrons. Accompanying this decomposition is the classification of electrons as either bound or free, which is useful for gapped and cold systems but becomes increasingly questionable as temperatures and pressures increase into the WDM regime. In this work we provide unambiguous first principles calculations of the dynamic structure factor of warm dense beryllium, independent of the Chihara form, by treating bound and free states under a single formalism. The computational approach is real-time finite-temperature time-dependent density functional theory (TDDFT) being applied here for the first time to WDM. We compare results from TDDFT to Chihara-based calculations for experimentally relevant conditions in shock-compressed beryllium.

Phase stability of TiO2 polymorphs from diffusion Quantum Monte Carlo

Titanium dioxide, TiO

Spin resolved energy parametrization of a quasi-one-dimensional electron gas

_2

Mechanical and optical response of [100] lithium fluoride to multi-megabar dynamic pressures

, has multiple applications in catalysis, energy conversion and memristive devices because of its electronic structure. Most of these applications utilize the naturally existing phases: rutile, anatase and brookite. Despite the simple form of TiO

Theory of melting at high pressures: Amending density functional theory with quantum Monte Carlo

_2

Equations of state and stability of MgSiO3 perovskite and post-perovskite phases from quantum Monte Carlo simulations

and its wide uses, there is long-standing disagreement between theory and experiment on the energetic ordering of these phases that has never been resolved. We present the first analysis of phase stability at zero temperature using the highly accurate many-body fixed node diffusion Quantum Monte Carlo (QMC) method. We also include the effects of temperature by calculating the Helmholtz free energy including both internal energy and vibrational contributions from density functional perturbation theory based quasi harmonic phonon calculations. Our QMC calculations find that anatase is the most stable phase at zero temperature, consistent with many previous mean-field calculations. However, at elevated temperatures, rutile becomes the most stable phase. For all finite temperatures, brookite is always the least stable phase.

Quantum Monte Carlo analysis of a charge ordered insulating antiferromagnet: The Ti4O7 Magneli phase

By carrying out extensive lattice regularized diffusion Monte Carlo calculations, we study the spin and density dependence of the ground-state energy for a quasi-one-dimensional electron gas, with harmonic transverse confinement and long-range 1/r interactions. We present a parametrization of the exchange–correlation energy suitable for spin density functional calculations, which fulfils exact low and high density limits.