Source author record

Paul R. C. Kent

Paul R. C. Kent 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

13works
5topics
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

13 published item(s)

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

Doped NiO: the Mottness of a charge transfer insulator

The evolution of the electronic structures of strongly correlated insulators with doping has long been a central fundamental question in condensed matter physics; it is also of great practical relevance for applications. We have studied the evolution of NiO under hole {\em and} electron doping using high-quality thin film and a wide range of experimental and theoretical methods. The evolution is in both cases very smooth with dopant concentration. The band gap is asymmetric under electron and hole doping, consistent with a charge-transfer insulator picture, and is reduced faster under hole than electron doping. For both electron and hole doping, occupied states are introduced at the top of the valence band. The formation of deep donor levels under electron doping and the inability to pin otherwise empty states near the conduction band edge is indicative that local electron addition and removal energies are dominated by a Mott-like Hubbard $U$-interaction even though the global bandgap is predominantly a charge-transfer type gap.

preprint2018arXiv

Doping a Bad Metal: Origin of Suppression of Metal-Insulator Transition in Non-Stoichiometric VO$_2$

Rutile ($R$) phase VO$_2$ is a quintessential example of a strongly correlated bad-metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic phase below T$_{MIT} \sim 340K$. It has been experimentally shown that one can control this transition by doping VO$_2$. In particular, doping with oxygen vacancies ($V_O$) has been shown to completely suppress this MIT {\em without} any structural transition. We explain this suppression by elucidating the influence of oxygen-vacancies on the electronic-structure of the metallic $R$ phase VO$_2$, explicitly treating strong electron-electron correlations using dynamical mean-field theory (DMFT) as well as diffusion Monte Carlo (DMC) flavor of quantum Monte Carlo (QMC) techniques. We show that $V_O$'s tend to change the V-3$d$ filling away from its nominal half-filled value, with the $e_{g}^π$ orbitals competing with the otherwise dominant $a_{1g}$ orbital. Loss of this near orbital polarization of the $a_{1g}$ orbital is associated with a weakening of electron correlations, especially along the V-V dimerization direction. This removes a charge-density wave (CDW) instability along this direction above a critical doping concentration, which further suppresses the metal-insulator transition. Our study also suggests that the MIT is predominantly driven by a correlation-induced CDW instability along the V-V dimerization direction.

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.

preprint2013arXiv

Spin-Resolved Self-Doping Tunes the Intrinsic Half-Metallicity of AlN Nanoribbons

We present a first-principles theoretical study of electric field-and strain-controlled intrinsic half-metallic properties of zigzagged aluminium nitride (AlN) nanoribbons. We show that the half-metallic property of AlN ribbons can undergo a transition into fully-metallic or semiconducting behavior with application of an electric field or uniaxial strain. An external transverse electric field induces a full charge screening that renders the material semiconducting. In contrast, as uniaxial strain varies from compressive to tensile, a spin-resolved selective self-doping increases the half-metallic character of the ribbons. The relevant strain-induced changes in electronic properties arise from band structure modifications at the Fermi level as a consequence of a spin-polarized charge transfer between pi-orbitals of the N and Al edge atoms in a spin-resolved self-doping process. This band structure tunability indicates the possibility ofdesigning magnetic nanoribbons with tunable electronic structure by deriving edge states from elements with sufficiently different localization properties. Finite temperature molecular dynamics reveal a thermally stable half-metallic nanoribbon up to room temperature.

preprint2012arXiv

Understanding Controls on Interfacial Wetting at Epitaxial Graphene: Experiment and Theory

The interaction of interfacial water with graphitic carbon at the atomic scale is studied as a function of the hydrophobicity of epitaxial graphene. High resolution X-ray reflectivity shows that the graphene-water contact angle is controlled by the average graphene thickness, due to the fraction of the film surface expressed as the epitaxial buffer layer whose contact angle (contact angle θ_c = 73°) is substantially smaller than that of multilayer graphene (θ_c = 93°). Classical and ab initio molecular dynamics simulations show that the reduced contact angle of the buffer layer is due to both its epitaxy with the SiC substrate and the presence of interfacial defects. This insight clarifies the relationship between interfacial water structure and hydrophobicity, in general, and suggests new routes to control interface properties of epitaxial graphene.

preprint2011arXiv

Simple Impurity Embedded in a Spherical Jellium: Approximations of Density Functional Theory compared to Quantum Monte Carlo Benchmarks

We study the electronic structure of a spherical jellium in the presence of a central Gaussian impurity. We test how well the resulting inhomogeneity effects beyond spherical jellium are reproduced by several approximations of density functional theory (DFT). Four rungs of Perdew's ladder of DFT functionals, namely local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA and orbital-dependent hybrid functionals are compared against our quantum Monte Carlo (QMC) benchmarks. We identify several distinct transitions in the ground state of the system as the electronic occupation changes between delocalized and localized states. We examine the parameter space of realistic densities ($1 \le r_s\le 5$) and moderate depths of the Gaussian impurity ($Z<7$). The selected 18 electron system (with closed-shell ground state) presents $1d \to 2s$ transitions while the 30 electron system (with open-shell ground state) exhibits $1f \to 2p$ transitions. For the former system, the accuracy for the transitions is clearly improving with increasing sophistication of functionals with meta-GGA and hybrid functionals having only small deviations from QMC. However, for the latter system, we find much larger differences for the underlying transitions between our pool of DFT functionals and QMC. We attribute this failure to treatment of the exact exchange within these functionals. Additionally, we amplify the inhomogeneity effects by creating the system with spherical shell which leads to even larger errors in DFT approximations.

preprint2011arXiv

The Role of Polytetrahedral Structures in the Elongation and Rupture of Gold Nanowires

We report comprehensive high-accuracy molecular dynamics simulations using the ReaxFF forcefield to explore the structural changes that occur as Au nanowires are elongated, establishing trends as a function of both temperature and nanowire diameter. Our simulations and subsequent quantitative structural analysis reveal that polytetrahedral structures (e.g., icosahedra) form within the "amorphous" neck regions, most prominently for systems with small diameter at high temperature. We demonstrate that the formation of polytetrahedra diminishes the conductance quantization as compared to systems without this structural motif. We demonstrate that use of the ReaxFF forcefield, fitted to high-accuracy first principles calculations of Au, combines the accuracy of quantum calculations with the speed of semi-empirical methods.

preprint2010arXiv

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

The self-healing diffusion Monte Carlo algorithm (SHDMC) [Phys. Rev. B {\bf 79}, 195117 (2009), {\it ibid.} {\bf 80}, 125110 (2009)] 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 N$_2$ 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 C$_{20}$ 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.

preprint2009arXiv

Self-healing diffusion quantum Monte Carlo algorithms: methods for direct reduction of the fermion sign error in electronic structure calculations

We develop a formalism and present an algorithm for optimization of the trial wave-function used in fixed-node diffusion quantum Monte Carlo (DMC) methods. We take advantage of a basic property of the walker configuration distribution generated in a DMC calculation, to (i) project-out a multi-determinant expansion of the fixed-node ground-state wave function and (ii) to define a cost function that relates the fixed-node ground-state and the non-interacting trial wave functions. We show that (a) locally smoothing out the kink of the fixed-node ground-state wave function at the node generates a new trial wave-function with better nodal structure and (b) we argue that the noise in the fixed-node wave-function resulting from finite sampling plays a beneficial role, allowing the nodes to adjust towards the ones of the exact many-body ground state in a simulated annealing-like process. We propose a method to improve both single determinant and multi-determinant expansions of the trial wave-function. We test the method in a model system where benchmark configuration interaction calculations can be performed. Comparing the DMC calculations with the exact solutions, we find that the trial wave-function is systematically improved. The overlap of the optimized trial wave function and the exact ground state converges to 100% even starting from wave-functions orthogonal to the exact ground state. In the optimization process we find an optimal non-interacting nodal potential of density-functional-like form whose existence was predicted earlier[Phys.Rev. B {\bf 77}, 245110 (2008)]. We obtain the exact Kohn-Sham effective potential from the DMC data.