skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Importance of elastic finite-size effects: Neutral defects in ionic compounds

Abstract

Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptotic behaviour and thus easily corrected with linear elastic theory. Here we show that elastic effect are also important in the description of defects in ionic compounds and can lead to qualitatively incorrect conclusions if inadequatly small supercells are used; moreover, the spurious self-interaction does not follow the behaviour predicted by linear elastic theory. Considering the exemplar cases of metal oxides with fluorite structure, we show that numerous previous studies, employing 96-atom supercells, misidentify the ground state structure of (charge neutral) Schottky defects. We show that the error is eliminated by employing larger cells (324, 768 and 1500 atoms), and careful analysis determines that elastic effects, not electrostatic, are responsible. The spurious self-interaction was also observed in non-oxide ionic compounds and irrespective of the computational method used, thereby resolving long standing discrepancies between DFT and force-field methods, previously attributed to themore » level of theory. The surprising magnitude of the elastic effects are a cautionary tale for defect calculations in ionic materials, particularly when employing computationally expensive methods (e.g. hybrid functionals) or when modelling large defect clusters. We propose two computationally practicable methods to test the magnitude of the elastic self-interaction in any ionic system. In commonly studies oxides, where electrostatic effects would be expected to be dominant, it is the elastic effects that dictate the need for larger supercells | greater than 96 atoms.« less

Authors:
 [1];  [2]
  1. Univ. of New South Wales, Sydney, NSW (Australia)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE); USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1396125
Alternate Identifier(s):
OSTI ID: 1390645
Report Number(s):
LA-UR-17-22161
Journal ID: ISSN 2469-9950; PRBMDO; TRN: US1702843
Grant/Contract Number:  
AC52-06NA25396; AC52- 06NA25396
Resource Type:
Accepted Manuscript
Journal Name:
Physical Review B
Additional Journal Information:
Journal Volume: 96; Journal Issue: 9; Journal ID: ISSN 2469-9950
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY

Citation Formats

Burr, P. A., and Cooper, M. W. D. Importance of elastic finite-size effects: Neutral defects in ionic compounds. United States: N. p., 2017. Web. doi:10.1103/PhysRevB.96.094107.
Burr, P. A., & Cooper, M. W. D. Importance of elastic finite-size effects: Neutral defects in ionic compounds. United States. doi:10.1103/PhysRevB.96.094107.
Burr, P. A., and Cooper, M. W. D. Fri . "Importance of elastic finite-size effects: Neutral defects in ionic compounds". United States. doi:10.1103/PhysRevB.96.094107. https://www.osti.gov/servlets/purl/1396125.
@article{osti_1396125,
title = {Importance of elastic finite-size effects: Neutral defects in ionic compounds},
author = {Burr, P. A. and Cooper, M. W. D.},
abstractNote = {Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptotic behaviour and thus easily corrected with linear elastic theory. Here we show that elastic effect are also important in the description of defects in ionic compounds and can lead to qualitatively incorrect conclusions if inadequatly small supercells are used; moreover, the spurious self-interaction does not follow the behaviour predicted by linear elastic theory. Considering the exemplar cases of metal oxides with fluorite structure, we show that numerous previous studies, employing 96-atom supercells, misidentify the ground state structure of (charge neutral) Schottky defects. We show that the error is eliminated by employing larger cells (324, 768 and 1500 atoms), and careful analysis determines that elastic effects, not electrostatic, are responsible. The spurious self-interaction was also observed in non-oxide ionic compounds and irrespective of the computational method used, thereby resolving long standing discrepancies between DFT and force-field methods, previously attributed to the level of theory. The surprising magnitude of the elastic effects are a cautionary tale for defect calculations in ionic materials, particularly when employing computationally expensive methods (e.g. hybrid functionals) or when modelling large defect clusters. We propose two computationally practicable methods to test the magnitude of the elastic self-interaction in any ionic system. In commonly studies oxides, where electrostatic effects would be expected to be dominant, it is the elastic effects that dictate the need for larger supercells | greater than 96 atoms.},
doi = {10.1103/PhysRevB.96.094107},
journal = {Physical Review B},
number = 9,
volume = 96,
place = {United States},
year = {2017},
month = {9}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 6 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Generalized Gradient Approximation Made Simple
journal, October 1996

  • Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
  • Physical Review Letters, Vol. 77, Issue 18, p. 3865-3868
  • DOI: 10.1103/PhysRevLett.77.3865

Point defect interactions in harmonic cubic lattices
journal, February 1967


Nonlocal theory of elastic interaction between point defects
journal, February 1978


Influence of charge states on energies of point defects and clusters in uranium dioxide
journal, April 2012


Vacancies in SiC: Influence of Jahn-Teller distortions, spin effects, and crystal structure
journal, June 1999


Vacancy–Vacancy Interaction Induced Oxygen Diffusivity Enhancement in Undoped Nonstoichiometric Ceria
journal, May 2015

  • Yuan, Fenglin; Zhang, Yanwen; Weber, William J.
  • The Journal of Physical Chemistry C, Vol. 119, Issue 23
  • DOI: 10.1021/acs.jpcc.5b01317

Molecular Dynamics Simulation of Thermal Transport in UO 2 Containing Uranium, Oxygen, and Fission-product Defects
journal, October 2016


Modelling of oxygen vacancy aggregates in monoclinic HfO 2 : can they contribute to conductive filament formation?
journal, September 2015

  • Bradley, Samuel R.; Bersuker, Gennadi; Shluger, Alexander L.
  • Journal of Physics: Condensed Matter, Vol. 27, Issue 41
  • DOI: 10.1088/0953-8984/27/41/415401

Ab initio study of solution energy and diffusion of caesium in uranium dioxide
journal, March 2009


Vacancy mediated cation migration in uranium dioxide: The influence of cluster configuration
journal, November 2014


Convergence of supercell calculations for point defects in semiconductors: Vacancy in silicon
journal, July 1998


Supercell size scaling of density functional theory formation energies of charged defects
journal, January 2009


Electrostatic interactions between charged defects in supercells
journal, December 2010

  • Freysoldt, Christoph; Neugebauer, Jörg; Van de Walle, Chris G.
  • physica status solidi (b), Vol. 248, Issue 5
  • DOI: 10.1002/pssb.201046289

Computational testing of trivalent dopants in CeO 2 for improved high-κ dielectric behaviour
journal, January 2013

  • Keating, Patrick R. L.; Scanlon, David O.; Watson, Graeme W.
  • J. Mater. Chem. C, Vol. 1, Issue 6
  • DOI: 10.1039/C2TC00385F

Point defects and non-stoichiometry in thoria
journal, December 2014


Fully Ab Initio Finite-Size Corrections for Charged-Defect Supercell Calculations
journal, January 2009


Connecting point defect parameters with bulk properties to describe diffusion in solids
journal, December 2016


Melting behavior of (Th,U)O2 and (Th,Pu)O2 mixed oxides
journal, October 2016


From ultrasoft pseudopotentials to the projector augmented-wave method
journal, January 1999


Systematic Analysis of 4 f Core Photoemission Spectra in Actinide Oxides
journal, January 1991

  • Yamazaki, Takao; Kotani, Akio
  • Journal of the Physical Society of Japan, Vol. 60, Issue 1
  • DOI: 10.1143/JPSJ.60.49

Boundaries determine the formation energies of lattice defects in two-dimensional buckled materials
journal, July 2016


Fast Parallel Algorithms for Short-Range Molecular Dynamics
journal, March 1995


Periodic boundary conditions in ab initio calculations
journal, February 1995


Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study
journal, January 1998

  • Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.
  • Physical Review B, Vol. 57, Issue 3, p. 1505-1509
  • DOI: 10.1103/PhysRevB.57.1505

Anisotropic charge screening and supercell size convergence of defect formation energies
journal, March 2013


Understanding and correcting the spurious interactions in charged supercells
journal, August 2011


Accurate prediction of defect properties in density functional supercell calculations
journal, November 2009


Density-functional theory and strong interactions: Orbital ordering in Mott-Hubbard insulators
journal, August 1995


First-principles DFT modeling of nuclear fuel materials
journal, April 2012

  • Liu, X. -Y.; Andersson, D. A.; Uberuaga, B. P.
  • Journal of Materials Science, Vol. 47, Issue 21
  • DOI: 10.1007/s10853-012-6471-6

Ab initio calculations of the energetics of the neutral Si vacancy defect
journal, January 1998

  • Mercer, J. L.; Nelson, J. S.; Wright, A. F.
  • Modelling and Simulation in Materials Science and Engineering, Vol. 6, Issue 1
  • DOI: 10.1088/0965-0393/6/1/001

Charge states of point defects in uranium oxide calculated with a local hybrid functional for correlated electrons
journal, May 2011


Conduction in polar crystals. I. Electrolytic conduction in solid salts
journal, January 1938

  • Mott, N. F.; Littleton, M. J.
  • Transactions of the Faraday Society, Vol. 34
  • DOI: 10.1039/tf9383400485

First-Principles Calculation of Point Defects in Uranium Dioxide
journal, January 2006


Vacancy clustering in zirconium: An atomic-scale study
journal, October 2014


First-principles study of defects and phase transition in UO 2
journal, October 2009


Fission gas in thoria
journal, March 2017


Understanding Oxygen-Vacancy Migration in the Fluorite Oxide CeO 2 : An Ab Initio Study of Impurity-Anion Migration
journal, December 2015

  • Genreith-Schriever, Annalena R.; Hebbeker, Pascal; Hinterberg, Judith
  • The Journal of Physical Chemistry C, Vol. 119, Issue 51
  • DOI: 10.1021/acs.jpcc.5b07813

From solid solution to cluster formation of Fe and Cr in α -Zr
journal, December 2015


Light-element diffusion in Mg using first-principles calculations: Anisotropy and elastodiffusion
journal, August 2016


Modeling of Ce O 2 , Ce 2 O 3 , and Ce O 2 x in the LDA + U formalism
journal, January 2007


Point defect modeling in materials: Coupling ab initio and elasticity approaches
journal, October 2013

  • Varvenne, Céline; Bruneval, Fabien; Marinica, Mihai-Cosmin
  • Physical Review B, Vol. 88, Issue 13
  • DOI: 10.1103/PhysRevB.88.134102

Method for locating low-energy solutions within DFT + U
journal, November 2010


Temperature accelerated dynamics study of migration process of oxygen defects in UO2
journal, February 2009


Structural phase transformation through defect cluster growth in Gd-doped ceria
journal, November 2011


Point defects and clustering in uranium dioxide by LSDA + U calculations
journal, March 2008


First-principles study on defect chemistry and migration of oxide ions in ceria doped with rare-earth cations
journal, January 2009

  • Nakayama, Masanobu; Martin, Manfred
  • Physical Chemistry Chemical Physics, Vol. 11, Issue 17
  • DOI: 10.1039/b900162j

Chemical Strain and Point Defect Configurations in Reduced Ceria
journal, June 2014

  • Wang, Bu; Xi, Xiaoning; Cormack, Alastair N.
  • Chemistry of Materials, Vol. 26, Issue 12
  • DOI: 10.1021/cm500946s

The elastic dipole tensor for point defects in ionic crystals
journal, March 1984


Linking atomic and mesoscopic scales for the modelling of the transport properties of uranium dioxide under irradiation
journal, July 2015


Occupation matrix control of d- and f-electron localisations using DFT + U
journal, January 2014

  • Allen, Jeremy P.; Watson, Graeme W.
  • Phys. Chem. Chem. Phys., Vol. 16, Issue 39
  • DOI: 10.1039/C4CP01083C

Association of defects in doped non-stoichiometric ceria from first principles
journal, January 2016

  • Grieshammer, Steffen; Nakayama, Masanobu; Martin, Manfred
  • Physical Chemistry Chemical Physics, Vol. 18, Issue 5
  • DOI: 10.1039/C5CP07537H

Lanthanum energetics in cubic ZrO2 and UO2 from DFT and DFT+U studies
journal, July 2011


A combined DFT + U and Monte Carlo study on rare earth doped ceria
journal, January 2014

  • Grieshammer, Steffen; Grope, Benjamin O. H.; Koettgen, Julius
  • Physical Chemistry Chemical Physics, Vol. 16, Issue 21
  • DOI: 10.1039/c3cp54811b

First-principles calculation of defect free energies: General aspects illustrated in the case of bcc Fe
journal, August 2015


Advances in first-principles modelling of point defects in UO 2 : f electron correlations and the issue of local energy minima
journal, July 2013


A many-body potential approach to modelling the thermomechanical properties of actinide oxides
journal, February 2014


First-principles calculations of momentum distributions of annihilating electron–positron pairs in defects in UO 2
journal, November 2016


Kinetically evolving irradiation-induced point defect clusters in UO 2 by molecular dynamics simulation
journal, September 2009


Impact of uniaxial strain and doping on oxygen diffusion in CeO2
journal, August 2014

  • Rushton, M. J. D.; Chroneos, A.
  • Scientific Reports, Vol. 4, Issue 1
  • DOI: 10.1038/srep06068

The energy and elastic dipole tensor of defects in ionic crystals calculated by the supercell method
journal, February 1985


First-principles calculation of the elastic dipole tensor of a point defect: Application to hydrogen in α -zirconium
journal, December 2016


Stability of oxygen point defects in UO 2 by first-principles DFT + U calculations: Occupation matrix control and Jahn-Teller distortion
journal, July 2010


Comparison of interatomic potentials for UO2. Part I: Static calculations
journal, June 2007


On the ionicity of carbon in Be 2 C. A theoretical investigation of its electronic structure
journal, January 1985

  • Herzig, P.; Redinger, J.
  • The Journal of Chemical Physics, Vol. 82, Issue 1
  • DOI: 10.1063/1.448755

First-principles study of noble gas impurities and defects in UO 2
journal, October 2011


Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set
journal, October 1996


First-principles calculations for point defects in solids
journal, March 2014

  • Freysoldt, Christoph; Grabowski, Blazej; Hickel, Tilmann
  • Reviews of Modern Physics, Vol. 86, Issue 1
  • DOI: 10.1103/RevModPhys.86.253

Defect cluster formation in M2O3-doped cubic ZrO2
journal, February 2000


Analysis of Intrinsic Defects in CeO 2 Using a Koopmans-Like GGA+ U Approach
journal, November 2011

  • Keating, Patrick R. L.; Scanlon, David O.; Morgan, Benjamin J.
  • The Journal of Physical Chemistry C, Vol. 116, Issue 3
  • DOI: 10.1021/jp2080034

Orbital Ordering, Jahn-Teller Distortion, and Anomalous X-Ray Scattering in Manganates
journal, May 1999


Development of Xe and Kr empirical potentials for CeO 2 , ThO 2 , UO 2 and PuO 2 , combining DFT with high temperature MD
journal, August 2016


Systematic Analysis of Core Photoemission Spectra for Actinide Di-Oxides and Rare-Earth Sesqui-Oxides
journal, January 1992

  • Kotani, Akio; Yamazaki, Takao
  • Progress of Theoretical Physics Supplement, Vol. 108
  • DOI: 10.1143/PTPS.108.117

DFT + U investigation of charged point defects and clusters in UO 2
journal, July 2014


Electrostatics-based finite-size corrections for first-principles point defect calculations
journal, May 2014


The elastic interaction of point defects
journal, September 1955


The elastic strengths of point defects
journal, November 1968


Experimental electronic structure of Be 2 C
journal, September 1998


Ab initio analysis of the defect structure of ceria
journal, April 2013


Trapping and diffusion of fission products in ThO2 and CeO2
journal, July 2011


What Is the Actual Local Crystalline Structure of Uranium Dioxide, UO 2 ? A New Perspective for the Most Used Nuclear Fuel
journal, December 2016


First-principles calculation of intrinsic defect formation volumes in silicon
journal, November 2005


Ordering of oxygen vacancies and excess charge localization in bulk ceria: A DFT + U study
journal, September 2014

  • Murgida, G. E.; Ferrari, V.; Ganduglia-Pirovano, M. Verónica
  • Physical Review B, Vol. 90, Issue 11
  • DOI: 10.1103/PhysRevB.90.115120

First-principles study of noble gas stability in ThO 2
journal, July 2017


First-order phase transition in UO 2 : 235 U and 17 O NMR study
journal, February 2001


Probing oxygen vacancy concentration and homogeneity in solid-oxide fuel-cell cathode materials on the subunit-cell level
journal, August 2012

  • Kim, Young-Min; He, Jun; Biegalski, Michael D.
  • Nature Materials, Vol. 11, Issue 10
  • DOI: 10.1038/nmat3393