Source author record

Lucas K. Wagner

Lucas K. Wagner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

20 published item(s)

preprint2022arXiv

Revisiting the Dark Matter Interpretation of Excess Rates in Semiconductors

In light of recent results from low-threshold dark matter detectors, we revisit the possibility of a common dark matter origin for multiple excesses across numerous direct detection experiments, with a focus on the excess rates in semiconductor detectors. We explore the interpretation of the low-threshold calorimetric excess rates above 40 eV in the silicon SuperCDMS Cryogenic Phonon Detector and above 100 eV in the germanium EDELWEISS Surface detector as arising from a common but unknown origin, and demonstrate a compatible fit for the observed energy spectra in both experiments, which follow a power law of index $α= 3.43^{+0.11}_{-0.06}$. Despite the intriguing scaling of the normalization of these two excess rates with approximately the square of the mass number $A^2$, we argue that the possibility of common origin by dark matter scattering via nuclear recoils is strongly disfavored, even allowing for exotic condensed matter effects in an as-yet unmeasured kinematic regime, due to the unphysically-large dark matter velocity required to give comparable rates in the different energy ranges of the silicon and germanium excesses. We also investigate the possibility of inelastic nuclear scattering by cosmic ray neutrons, solar neutrinos, and photons as the origin, and quantitatively disfavor all three based on known fluxes of particles.

preprint2021arXiv

Identifying materials with charge-spin physics using charge-spin susceptibility computed from first principles

The authors present a quantity termed charge-spin susceptibility, which measures the charge response to spin degrees of freedom in strongly correlated materials. This quantity is simple to evaluate using both standard density functional theory and many-body electronic structure techniques, enabling comparison between different levels of theory. A benchmark on 28 layered magnetic materials shows that large values of charge-spin susceptibility correlate with unconventional ground states such as disordered magnets and unconventional superconductivity.

preprint2020arXiv

A comparison of computed and experimental neutron diffraction intensity at large momentum for MnO and NiO

Magnetic neutron scattering measures spin-spin correlations giving information about the long-range spin order as well as the shape of the spin density in magnetic materials. Similarly, detailed first principles calculations directly compute the spin density in materials. In this work, the authors carefully compare experimentally measured magnetic neutron intensities to three levels of theory: density functional theory in two approximations, and fixed-node diffusion Monte Carlo. While each theory performs similarly for the simple antiferromagnet MnO, there are significant differences between density functional theory and diffusion Monte Carlo in NiO. For both materials, fixed-node diffusion Monte Carlo shows the lowest RMS error with respect to experiment for the form factor and the intensity. Through connection of the intensities and form factors to the real-space spin densities, it is shown that diffusion Monte Carlo spin density becomes more diffuse by spreading spin away from the bond directions. This benchmark is of importance when considering the efficacy of an ab initio calculation in capturing both the long-range magnetic correlations and the microscopic spin details of magnetic systems.

preprint2020arXiv

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, $\partial E/\partial p$, is a key ingredient in variational quantum Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve Monte Carlo estimate of this derivative suffers from an infinite variance which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of $\partial E/\partial p$ for a multi-Slater-Jastrow trial wave function on the LiH molecule. This regularized estimator is a simple and efficient estimator of $\partial E/\partial p$ for VMC optimization techniques.

preprint2020arXiv

Combining density functional theories to correctly describe the energy, lattice structure and electronic density of functional oxide perovskites

Functional oxide perovskites are the pillar of cutting-edge technological applications. Density functional theory (DFT) simulations are the theoretical methods of choice to understand and design perovskite materials. However, tests on the reliability of DFT to describe fundamental properties of oxide perovskites are scarce and mostly ill-defined due to a lack of rigorous theoretical benchmarks for solids. Here, we present a quantum Monte Carlo benchmark study of DFT on the archetypal perovskite BaTiO$_{3}$ (BTO). It shows that no DFT approximation can simultaneously reproduce the energy, structure, and electronic density of BTO. Traditional protocols to select DFT approximations are empirical and fail to detect this shortcoming. An approach combining two different non-empirical DFT schemes, "SCAN" and "HSE06", is able to holistically describe BTO with accuracy. Combined DFT approaches should thus be considered as a promising alternative to standard methods for simulating oxide perovskites.

preprint2020arXiv

Electric field effects on the band gap and edge states of monolayer 1T'-WTe2

Monolayer 1T'-WTe2 is a quantum spin Hall insulator with a gapped bulk and gapless helical edge states persisting to temperatures around 100 K. Recent studies have revealed a topological-to-trivial phase transition as well the emergence of an unconventional, potentially topological superconducting state upon tuning the carrier concentration with gating. However, despite extensive studies, the effects of gating on the band structure and the helical edge states have not yet been established. In this work we present a combined low-temperature STM and first principles study of back-gated monolayer 1T'-WTe2 films grown on graphene. Consistent with a quantum spin Hall system, the films show well-defined bulk gaps and clear edge states that span the gap. By directly measuring the density of states with STM spectroscopy, we show that the bulk band gap magnitude shows substantial changes with applied gate voltage, which is contrary to the naïve expectation that a gate would rigidly shift the bands relative to the Fermi level. To explain our data, we carry out density functional theory and model Hamiltonian calculations which show that a gate electric field causes doping and inversion symmetry breaking which polarizes and spin-splits the bulk bands. Interestingly, the calculated spin splitting from the effective Rashba-like spin-orbit coupling can be in the tens of meV for the electric fields in the experiment, which may be useful for spintronics applications. Our work reveals the strong effect of electric fields on the bulk band structure of monolayer 1T'-WTe2, which will play a critical role in our understanding of gate-induced phenomena in this system.

preprint2019arXiv

Direct comparison of many-body methods for realistic electronic Hamiltonians

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons.

preprint2018arXiv

Effective spin-orbit models using correlated first-principles wave functions

Diffusion Monte Carlo is one of the most accurate scalable many-body methods for solid state systems. However, to date, spin-orbit interactions have not been incorporated into these calcualtions at a first-principles level; only having been applied to small systems. In this technique, we use explicitly correlated first-principles quantum Monte Carlo calculations to derive an effective spin-orbit model Hamiltonian. The simplified model Hamiltonian is then solved to obtain the energetics of the system. To demonstrate this method, benchmark studies are performed in main-group atoms and monolayer tungsten disulfide, where high accuracy is obtained.

preprint2016arXiv

Competing collinear magnetic structures in superconducting FeSe by first principles quantum Monte Carlo calculations

Resolving the interplay between magnetic interactions and structural properties in strongly correlated materials through a quantitatively accurate approach has been a major challenge in condensed matter physics. Here we apply highly accurate first principles quantum Monte Carlo (QMC) techniques to obtain structural and magnetic properties of the iron selenide (FeSe) superconductor under pressure. Where comparable, the computed properties are very close to the experimental values. Of potential ordered magnetic configurations, collinear spin configurations are the most energetically favorable over the explored pressure range. They become nearly degenerate in energy with bicollinear spin orderings at around 7 GPa, when the experimental critical temperature $T_c$ is the highest. On the other hand, ferromagnetic, checkerboard, and staggered dimer configurations become relatively higher in energy as the pressure increases. The behavior under pressure is explained by an accurate analysis of the charge compressibility and the orbital occupation as described by the QMC many-body wave function, which reveals how spin, charge and orbital degrees of freedom are strongly coupled in this compound. This remarkable pressure evolution suggests that stripe-like magnetic fluctuations may be responsible for the enhanced $T_c$ in FeSe and that higher T$_c$ is associated with nearness to a crossover between collinear and bicollinear ordering.

preprint2016arXiv

Computational and experimental investigation of unreported transition metal selenides and sulphides

Expanding the library of known inorganic materials with functional electronic or magnetic behavior is a longstanding goal in condensed matter physics and materials science. Recently, the transition metal chalchogenides including selenium and sulfur have been of interest because of their correlated-electron properties, as seen in the iron based superconductors and the transition metal dichalcogenides. However, the chalcogenide chemical space is less explored than that of oxides, and there is an open question of whether there may be new materials heretofore undiscovered. We perform a systematic combined theoretical and experimental search over ternary phase diagrams that are empty in the Inorganic Crystal Structure Database containing cations, transition metals, and one of selenium or sulfur. In these 27 ternary systems, we use a probabilistic model to reduce the likelihood of false negative predictions, which results in a list of 24 candidate materials. We then conduct a variety of synthesis experiments to check the candidate materials for stability. While the prediction method did obtain previously unknown compositions that are predicted stable within density functional theory, none of the candidate materials formed in our experiments. We come to the conclusion that these phase diagrams are "empty" in the case of bulk synthesis, but it remains a possibility that alternate synthesis routes may produce some of these phases.

preprint2016arXiv

Discovering correlated fermions using quantum Monte Carlo

It has become increasingly feasible to use quantum Monte Carlo (QMC) methods to study correlated fermion systems for realistic Hamiltonians. We give a summary of these techniques targeted at researchers in the field of correlated electrons, focusing on the fundamentals, capabilities, and current status of this technique. The QMC methods often offer the highest accuracy solutions available for systems in the continuum, and, since they address the many-body problem directly, the simulations can be analyzed to obtain insight into the nature of correlated quantum behavior.

preprint2016arXiv

Using Local Operator Fluctuations to Identify Wave Function Improvements

A method is developed that allows analysis of quantum Monte Carlo simulations to identify errors in trial wave functions. The purpose of this method is to allow for the systematic improvement of variational wave functions by identifying degrees of freedom that are not well-described by an initial trial state. We provide proof of concept implementations of this method by identifying the need for a Jastrow correlation factor, and implementing a selected multi-determinant wave function algorithm for small dimers that systematically decreases the variational energy. Selection of the two-particle excitations is done using quantum Monte Carlo within the presence of a Jastrow correlation factor, and without the need to explicitly construct the determinants. We also show how this technique can be used to design compact wave functions for transition metal systems. This method may provide a route to analyze and systematically improve descriptions of complex quantum systems in a scalable way.

preprint2015arXiv

Density-matrix based determination of low-energy model Hamiltonians from ab initio wavefunctions

We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio Quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding (AIDMD). For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard $U^{*}/t$ to be $1.3 \pm 0.2$, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large scale calculations using techniques designed for lattice models.

preprint2015arXiv

Ground state of doped cuprates from first principles quantum Monte Carlo calculations

The author reports on new high-fidelity simulations of charge carriers in the high-T$_c$ cuprate materials using quantum Monte Carlo techniques applied to the first principles Hamiltonian. With this high accuracy technique, the doped ground state is found to be a spin polaron, in which charge is localized through a strong interaction with the spin. This spin polaron has calculated properties largely similar to the phenomenology of the cuprates, and may be the object which forms the Fermi surface and charge inhomogeneity in these materials. The spin polaron has some unique features that should be visible in X-ray, EELS, and neutron experiments. The results contained in this paper comprise an accurate first principles derived paradigm from which to study superconductivity in the cuprates.

preprint2015arXiv

Phase Stability and Properties of Manganese Oxide Polymorphs: Assessment and Insights from Diffusion Monte Carlo

We present an analysis of the polymorphic energy ordering and properties of the rock salt and zincblende structures of manganese oxide using fixed node diffusion Monte Carlo (DMC). Manganese oxide is a correlated, antiferromagnetic material that has proven to be challenging to model from first principles across a variety of approaches. Unlike conventional density functional theory and some hybrid functionals, fixed node diffusion Monte Carlo finds the rock salt structure to be more stable than the zincblende structure, and thus recovers the correct energy ordering. Analysis of the site-resolved charge fluctuations of the wave functions according to DMC and other electronic structure descriptions give insights into elements that are missing in other theories. While the calculated band gaps within DMC are in agreement with predictions that the zincblende polymorph has a lower band gap, the gaps themselves overestimate reported experimental values.

preprint2015arXiv

Towards a systematic assessment of errors in diffusion Monte Carlo calculations of semiconductors: case study of zinc selenide and zinc oxide

The fixed node diffusion Monte Carlo (DMC) method has attracted interest in recent years as a way to calculate properties of solid materials with high accuracy. However, the framework for the calculation of properties such as total energies, atomization energies, and excited state energies is not yet fully established. Several outstanding questions remain as to the effect of pseudopotentials, the magnitude of the fixed node error, and the size of supercell finite size effects. Here, we consider in detail the semiconductors ZnSe and ZnO and carry out systematic studies to assess the magnitude of the energy differences arising from controlled and uncontrolled approximations in DMC. The former include time step errors and supercell finite size effects for ground and optically excited states, and the latter include pseudopotentials, the pseudopotential localization approximation, and the fixed node approximation. We find that for these compounds, the errors can be controlled to good precision using modern computational resources, and that quantum Monte Carlo calculations using Dirac-Fock pseudopotentials can offer good estimates of both cohesive energy and the gap of these systems. We do however observe differences in calculated optical gaps that arise when different pseudopotentials are used.

preprint2014arXiv

The effect of electron correlation on the electronic structure and spin-lattice coupling of the high-T$_c$ cuprates: quantum Monte Carlo calculations

Electron correlation effects are particularly strong in the high temperature superconducting materials. Devising an accurate description of these materials has long been a challenge, with these strong correlation effects historically being considered impossible or impractical to simulate computationally. Using quantum Monte Carlo techniques, we have explicitly simulated electron correlations in several cuprate materials from first principles. These simulations accurately reproduce many important physical quantities about these materials, including the interaction-induced gap and the superexchange coupling between copper spins, with no additional parameters beyond fundamental constants. We further investigate the dimensionless spin-lattice coupling parameter in the parent materials, showing that it varies dramatically between 0.1 and 1.0 depending on the interlayer. This result indicates that the lattice and magnetic degrees of freedom are not independent in these systems, which may have ramifications for the origin of superconductivity.

preprint2012arXiv

Point-Defect Optical Transitions and Thermal Ionization Energies from Quantum Monte Carlo Methods: Application to F-center Defect in MgO

We present an approach to calculation of point defect optical and thermal ionization energies based on the highly accurate quantum Monte Carlo methods. The use of an inherently many-body theory that directly treats electron correlation offers many improvements over the typically-employed density functional theory Kohn-Sham description. In particular, the use of quantum Monte Carlo methods can help overcome the band gap problem and obviate the need for ad-hoc corrections. We demonstrate our approach to the calculation of the optical and thermal ionization energies of the F-center defect in magnesium oxide, and obtain excellent agreement with experimental and/or other high-accuracy computational results.

preprint2012arXiv

Types of single particle symmetry breaking in transition metal oxides due to electron correlation

Very accurate wave functions are calculated for small transition metal oxide molecules. These wave functions are decomposed using reduced density matrices to study the underlying correlation of electrons. The correlation is primarily of left-right type between the transition metals and the oxygen atoms, which is mediated by excitations from the nominal single Slater ground state into antibonding and d-type orbitals. In a localized representation, this correlation manifests itself in a 2-electron hopping term that is off-diagonal. This term is of similar magnitude to the commonly considered Hubbard-type on-site interaction.

preprint2010arXiv

Quantum Monte Carlo for minimum energy structures

We present an efficient method to find minimum energy structures using energy estimates from accurate quantum Monte Carlo calculations. This method involves a stochastic process formed from the stochastic energy estimates from Monte Carlo that can be averaged to find precise structural minima while using inexpensive calculations with moderate statistical uncertainty. We demonstrate the applicability of the algorithm by minimizing the energy of the H2O-OH- complex and showing that the structural minima from quantum Monte Carlo calculations affect the qualitative behavior of the potential energy surface substantially.