Researcher profile

Luke Shulenburger

Luke Shulenburger contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

Benchmarking Adiabatic Quantum Optimization for Complex Network Analysis

We lay the foundation for a benchmarking methodology for assessing current and future quantum computers. We pose and begin addressing fundamental questions about how to fairly compare computational devices at vastly different stages of technological maturity. We critically evaluate and offer our own contributions to current quantum benchmarking efforts, in particular those involving adiabatic quantum computation and the Adiabatic Quantum Optimizers produced by D-Wave Systems, Inc. We find that the performance of D-Wave's Adiabatic Quantum Optimizers scales roughly on par with classical approaches for some hard combinatorial optimization problems; however, architectural limitations of D-Wave devices present a significant hurdle in evaluating real-world applications. In addition to identifying and isolating such limitations, we develop algorithmic tools for circumventing these limitations on future D-Wave devices, assuming they continue to grow and mature at an exponential rate for the next several years.

preprint2016arXiv

Phase Stability of TiO$_2$ Polymorphs from Diffusion Quantum Monte Carlo

Titanium dioxide, TiO$_2$, 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$_2$ 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.

preprint2016arXiv

Quantum Monte Carlo analysis of a charge ordered insulating antiferromagnet: the Ti$_4$O$_7$ Magnéli phase

The Magnéli phase Ti$_4$O$_7$ is an important transition metal oxide with a wide range of applications because of its interplay between charge, spin, and lattice degrees of freedom. At low temperatures, it has non-trivial magnetic states very close in energy, driven by electronic exchange and correlation interactions. We have examined three low-lying states, one ferromagnetic and two antiferromagnetic, and calculated their energies as well as Ti spin moment distributions using highly accurate Quantum Monte Carlo methods. We compare our results to those obtained from density functional theory-based methods that include approximate corrections for exchange and correlation. Our results confirm the nature of the states and their ordering in energy, as compared with density-functional theory methods. However, the energy differences and spin distributions differ. A detailed analysis suggests that non-local exchange-correlation functionals, in addition to other approximations such as LDA+U to account for correlations, are needed to simultaneously obtain better estimates for spin moments, distributions, energy differences and energy gaps.

preprint2016arXiv

X-Ray Thomson scattering without the Chihara decomposition

X-Ray Thomson Scattering (XRTS) 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 (DSF). In most models, this is decomposed into three terms [Chihara, J. Phys. F: Metal Phys. {\bf 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.

preprint2015arXiv

Beyond chemical accuracy: The pseudopotential approximation in diffusion Monte Carlo calculations of the HCP to BCC phase transition in beryllium

Motivated by the disagreement between recent diffusion Monte Carlo calculations and experiments on the phase transition pressure between the ambient and beta-Sn phases of silicon, we present a study of the HCP to BCC phase transition in beryllium. This lighter element provides an oppor- tunity for directly testing many of the approximations required for calculations on silicon and may suggest a path towards increasing the practical accuracy of diffusion Monte Carlo calculations of solids in general. We demonstrate that the single largest approximation in these calculations is the pseudopotential approximation. After removing this we find excellent agreement with experiment for the ambient HCP phase and results similar to careful calculations using density functional theory for the phase transition pressure.

preprint2015arXiv

Shock Response and Phase Transitions of MgO at Planetary Impact Conditions

The moon-forming impact and the subsequent evolution of the proto-Earth is strongly dependent on the properties of materials at the extreme conditions generated by this violent collision. We examine the high pressure behavior of MgO, one of the dominant constituents in the earth's mantle, using high-precision, plate impact shock compression experiments performed on Sandia National Laboratories Z-Machine and extensive quantum simulations using Density Functional Theory (DFT) and quantum Monte Carlo (QMC). The combined data span from ambient conditions to 1.2 TPa and 42,000 K, showing solid-solid and solid-liquid phase boundaries. Furthermore our results indicate under impact that the solid and liquid phases coexist for more than 100 GPa, pushing complete melting to pressures in excess of 600 GPa. The high pressure required for complete shock melting places a lower bound on the relative velocities required for the moon forming impact.

preprint2015arXiv

The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus

An outstanding challenge of theoretical electronic structure is the description of van der Waals (vdW) interactions in molecules and solids. Renewed interest in resolving this is in part motivated by the technological promise of layered systems including graphite, transition metal dichalcogenides, and more recently, black phosphorus, in which the interlayer interaction is widely believed to be dominated by these types of forces. We report a series of quantum Monte Carlo (QMC) calculations for bulk black phosphorus and related few-layer phosphorene, which elucidate the nature of the forces that bind these systems and provide benchmark data for the energetics of these systems. We find a significant charge redistribution due to the interaction between electrons on adjacent layers. Comparison to density functional theory (DFT) calculations indicate not only wide variability even among different vdW corrected functionals, but the failure of these functionals to capture the trend of reorganization predicted by QMC. The delicate interplay of steric and dispersive forces between layers indicate that few-layer phosphorene presents an unexpected challenge for the development of vdW corrected DFT.

preprint2014arXiv

Equations of state and stability of MgSiO$_3$ perovskite and post-perovskite phases from quantum Monte Carlo simulations

We have performed quantum Monte Carlo (QMC) simulations and density functional theory (DFT) calculations to study the equations of state of MgSiO$_3$ perovskite (Pv) and post-perovskite (PPv), up to the pressure and temperature conditions of the base of Earth's lower mantle. The ground state energies were derived using QMC and the temperature dependent Helmholtz free energies were calculated within the quasi-harmonic approximation and density functional perturbation theory. The equations of state for both phases of MgSiO$_3$ agree well with experiments, and better than those from generalized gradient approximation (GGA) calculations. The Pv-PPv phase boundary calculated from our QMC equations of states is also consistent with experiments, and better than previous LDA calculations. We discuss the implications for double crossing of the Pv-PPv boundary in the Earth.

preprint2013arXiv

Quantum Monte Carlo applied to solids

We apply diffusion quantum Monte Carlo (DMC) to a broad set of solids, benchmarking the method by comparing bulk structural properties (equilibrium volume and bulk modulus) to experiment and 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 density functional theory (DFT) when using the new generation of functionals, including one hybrid functional and two dispersion corrected functionals. The excellent performance of quantum Monte Carlo (QMC) 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 pseudo-potential approximations. This test set allows for any improvements in these methods to be judged in a systematic way.

preprint2013arXiv

Theory of melting at high pressures: Amending Density Functional Theory with Quantum Monte Carlo

We present an improved first-principles description of melting under pressure based on thermodynamic integration comparing Density Functional Theory (DFT) and quantum Monte Carlo (QMC) treatments of the system. The method is applied to address the longstanding discrepancy between density functional theory (DFT) calculations and diamond anvil cell (DAC) experiments on the melting curve of xenon, a noble gas solid where van der Waals binding is challenging for traditional DFT methods. The calculations show excellent agreement with data below 20 GPa and that the high-pressure melt curve is well described by a Lindemann behavior up to at least 80 GPa, a finding in stark contrast to DAC data.

preprint2009arXiv

Relative importance of crystal field versus bandwidth to the high pressure spin transition in transition metal monoxides

The crystal field splitting and d bandwidth of the 3d transition metal monoxides MnO, FeO, CoO and NiO are analyzed as a function of pressure within density functional theory. In all four cases the 3d bandwidth is significantly larger than the crystal field splitting over a wide range of compressions. The bandwidth actually increases more as pressure is increased than the crystal field splitting. Therefore the role of increasing bandwidth must be considered in any explanation of a possible spin collapse that these materials may exhibit under pressure.

preprint2008arXiv

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

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 fulfills exact low and high density limits.