Source author record

Jaron T. Krogel

Jaron T. Krogel 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

12works
6topics
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

12 published item(s)

preprint2022arXiv

A new generation of effective core potentials from correlated and spin-orbit calculations: selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference.

preprint2022arXiv

Correlative nanoscale imaging of strained hBN spin defects

Spin defects like the negatively charged boron vacancy color center ($V_B^-$) in hexagonal boron nitride (hBN) may enable new forms of quantum sensing with near-surface defects in layered van der Waals heterostructures. Here, we reveal the effect of strain associated with creases in hBN flakes on $V_B^-$ and $V_B$ color centers in hBN with correlative cathodoluminescence and photoluminescence microscopies. We observe strong localized enhancement and redshifting of the $V_B^-$ luminescence at creases, consistent with density functional theory calculations showing $V_B^-$ migration toward regions with moderate uniaxial compressive strain. The ability to manipulate these spin defects with highly localized strain offers intriguing possibilities for future 2D quantum sensors.

preprint2022arXiv

Magnetic measures of purity for MnBi$_2$Te$_4$

The intrinsically anti-ferromagnetic topological insulator, MnBi$_2$Te$_4$ (MBT), has garnered significant attention recently. The excitement for this layered van der Waals bonded compound stems from its potential to host numerous exotic topological quantum states. For instance, quantum anomalous Hall states are predicted for odd-layer compounds and axion insulator states for the even-layer compounds. Unfortunately, the realization of these phenomena has been hindered by experimental challenges such as the existence of negative charge carriers, i.e., electron doping, which have been linked to anti-site defects among the Mn and Bi sub-lattices. Based on high level diffusion Monte Carlo (DMC) and DMC-tuned DFT+U calculations, we provide benchmark quality results for the bulk Mn magnetization as well as for Mn$_{Bi}$ and Bi$_{Mn}$ defects. We use this information to refine and extend models that estimate the anti-site defect concentration in actual MBT samples when combined with data from magnetic susceptibility and intermediate field magnetization measurements. Our models are validated through favorable comparison with prior experimental studies that obtained both magnetic and site occupancy data. We then extend our estimates to a larger set of prior samples to identify a probable zone of low defect density that has yet to be reached in synthesis. We anticipate our theoretically based magnetic purity measures may be used as minimization targets in the cycle of refinement needed to synthesize MBT samples with low anti-site defect concentrations and more reproducible topological properties.

preprint2022arXiv

Origin of Metal-Insulator Transitions in Correlated Perovskite Metals

The mechanisms that drive metal-to-insulator transitions (MIT) in correlated solids are not fully understood. For example, the perovskite (PV) SrCoO3 is a FM metal while the oxygen-deficient (n-doped) brownmillerite (BM) SrCoO2.5 is an anti-ferromagnetic (AFM) insulator. Given the magnetic and structural transitions that accompany the MIT, the driver for such a MIT transition is unclear. We also observe that the perovskite metals LaNiO3, SrFeO3, and SrCoO3 also undergo MIT when n-doped via high-to-low valence compositional changes. Also, pressurizing the insulating BM SrCoO2.5 phase, drives a gap closing. Using DFT and correlated diffusion Monte Carlo approaches we demonstrate that the ABO3 perovskites most prone to MIT are self hole-doped materials, reminiscent of a negative charge-transfer system. Upon n-doping away from the negative-charge transfer metallic phase, an underlying charge-lattice (or e-phonon) coupling drives the system to a bond-disproportionated gapped state, thereby achieving ligand hole passivation at certain sites only, leading to charge-disproportionated states. The size of the gap opened is correlated with the size of the hole-filling at these ligand sites. This suggests that the interactions driving the gap opening to realize a MIT even in correlated metals is the charge-transfer energy, but it couples with the underlying phonons to enable the transition to the insulating phase. Other orderings (magnetic, charge, etc.) driven by weaker interactions are secondary and may assist gap openings at small dopings, but its the charge-transfer energy that predominantly determines the bandgap, with a negative energy preferring the metallic phase. This n-doping can be achieved by modulations in stoichiometry or composition or pressure. Hence, controlling the amount of the ligand-hole is key in controlling MIT. We compare our predictions to experiments where possible.

preprint2020arXiv

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

preprint2016arXiv

Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo

We have applied the diffusion quantum Monte Carlo (DMC) method to calculate the cohesive energy and the structural parameters of the binary oxides CaO, SrO, BaO, Sc2O3, Y2O3 and La2O3. The aim of our calculations is to systematically quantify the accuracy of the DMC method to study this type of metal oxides. The DMC results were compared with local, semi-local and hybrid Density Functional Theory (DFT) approximations as well as with experimental measurements. The DMC method yields cohesive energies for these oxides with a mean absolute deviation from experimental measurements of 0.18(2) eV, while with local, semi-local and hybrid DFT approximations the deviation is 3.06, 0.94 and 1.23 eV, respectively. For lattice constants, the mean absolute deviation in DMC, local, semi-local and hybrid DFT approximations, are 0.017(1), 0.07, 0.05 and 0.04 Å, respectively. DMC is highly accurate method, outperforming the DFT approximations in describing the cohesive energies and structural parameters of these binary oxides.

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.

preprint2015arXiv

Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo

We have applied the many-body ab-initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure, and the energetics of the oxygen vacancy, zinc interstitial and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O$_2$, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type, and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

preprint2014arXiv

Ab initio quantum Monte Carlo calculations of spin superexchange in cuprates: the benchmarking case of Ca$_2$CuO$_3$

In view of the continuous theoretical efforts aimed at an accurate microscopic description of the strongly correlated transition metal oxides and related materials, we show that with continuum quantum Monte Carlo (QMC) calculations it is possible to obtain the value of the spin superexchange coupling constant of a copper oxide in a quantitatively excellent agreement with experiment. The variational nature of the QMC total energy allows us to identify the best trial wave function out of the available pool of wave functions, which makes the approach essentially free from adjustable parameters and thus truly ab initio. The present results on magnetic interactions suggest that QMC is capable of accurately describing ground state properties of strongly correlated materials.

preprint2014arXiv

Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study

We develop an energy density matrix that parallels the one-body reduced density matrix (1RDM) for many-body quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For mean-field systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques implemented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences shows a quantitative connection between the occupation energies and electron addition and removal energies for the electron gas. For the oxygen atom, the association between the ground state occupation energies and particle addition and removal energies becomes only qualitative. The energy density matrix provides a new avenue for describing energetics with quantum Monte Carlo methods which have traditionally been limited to total energies.

preprint2013arXiv

The Quantum Energy Density: Improved Efficiency for Quantum Monte Carlo

We establish a physically meaningful representation of a quantum energy density for use in Quantum Monte Carlo calculations. The energy density operator, defined in terms of Hamiltonian components and density operators, returns the correct Hamiltonian when integrated over a volume containing a cluster of particles. This property is demonstrated for a helium-neon "gas," showing that atomic energies obtained from the energy density correspond to eigenvalues of isolated systems. The formation energies of defects or interfaces are typically calculated as total energy differences. Using a model of delta-doped silicon (where dopant atoms form a thin plane) we show how interfacial energies can be calculated more efficiently with the energy density, since the region of interest is small. We also demonstrate how the energy density correctly transitions to the bulk limit away from the interface where the correct energy is obtainable from a separate total energy calculation.