Variational charge relaxation in ionic crystals: An efficient treatment of statics and dynamics

Andrew Chizmeshya, Frank M. Zimmermann, Randall A. LaViolette, George Wolf

Research output: Contribution to journalArticle

18 Citations (Scopus)

Abstract

An ab initio variationally induced breathing (VIB) description of the energetics and dynamics of ionic crystals is developed within the spherical-ion-pair approximation and the Gordon-Kim ansatz. In the VIB method, the total crystal charge density is given by the overlap of Kohn-Sham ionic charge densities. Spherical charge deformation of the ions is accomplished by adding an effective many-body crystal site potential to the atomic potential used in the Kohn-Sham procedure. Parameters defining this effective site potential are treated variationally such that the total crystal electronic energy is a minimum for a given structural configuration. We have explored a number of forms of the effective potential, including Watson-sphere types where both the charge and radius of the sphere are variationally determined. In addition, we have investigated anion-cation charge transfer by incorporating an additional variational parameter. We also present a formulation of the lattice dynamics of ionic crystals based on the VIB prescription in which the electronic variational parameters are treated as dynamical variables within the adiabatic approximation. With the computed phonon density of states, the structural parameters of a crystal can then be determined at any temperature and pressure (or stress condition) by mini- mizing the quasiharmonic Gibbs free energy with respect to the structural parameters along the electronic adiabatic surface. The complete pressure- and temperature-dependent thermal and elastic properties of a crystal can also be determined. We have applied the VIB procedure to calculate the equations of state, elastic properties, phonon dispersion relations, and phase stability of the alkaline earth oxides (MgO, CaO, and SrO). The calculated properties are generally in quite good agreement with experiment. Spherical charge deformation within the ionic description is responsible for a considerable improvement over rigid-ion elastic constants and bulk moduli, but some of the remaining discrepancies cannot in principle be eliminated within a spherical-ion model. In particular, the effect of spherical breathing on the optical vibrational frequencies is not dramatic and further improvement (e.g., reduction of LO-TO splitting) will require the inclusion of nonspherical charge relaxation (polarizability).

Original languageEnglish (US)
Pages (from-to)15559-15574
Number of pages16
JournalPhysical Review B
Volume50
Issue number21
DOIs
StatePublished - 1994

Fingerprint

ionic crystals
breathing
Crystals
elastic properties
crystals
Ions
ions
electronics
Charge density
alkaline earth oxides
Gibbs free energy
approximation
bulk modulus
Lattice vibrations
Phase stability
equations of state
thermodynamic properties
Elastic constants
Vibrational spectra
charge transfer

ASJC Scopus subject areas

  • Condensed Matter Physics

Cite this

Variational charge relaxation in ionic crystals : An efficient treatment of statics and dynamics. / Chizmeshya, Andrew; Zimmermann, Frank M.; LaViolette, Randall A.; Wolf, George.

In: Physical Review B, Vol. 50, No. 21, 1994, p. 15559-15574.

Research output: Contribution to journalArticle

Chizmeshya, Andrew ; Zimmermann, Frank M. ; LaViolette, Randall A. ; Wolf, George. / Variational charge relaxation in ionic crystals : An efficient treatment of statics and dynamics. In: Physical Review B. 1994 ; Vol. 50, No. 21. pp. 15559-15574.
@article{15bd10ba76744050b418089450a5f99d,
title = "Variational charge relaxation in ionic crystals: An efficient treatment of statics and dynamics",
abstract = "An ab initio variationally induced breathing (VIB) description of the energetics and dynamics of ionic crystals is developed within the spherical-ion-pair approximation and the Gordon-Kim ansatz. In the VIB method, the total crystal charge density is given by the overlap of Kohn-Sham ionic charge densities. Spherical charge deformation of the ions is accomplished by adding an effective many-body crystal site potential to the atomic potential used in the Kohn-Sham procedure. Parameters defining this effective site potential are treated variationally such that the total crystal electronic energy is a minimum for a given structural configuration. We have explored a number of forms of the effective potential, including Watson-sphere types where both the charge and radius of the sphere are variationally determined. In addition, we have investigated anion-cation charge transfer by incorporating an additional variational parameter. We also present a formulation of the lattice dynamics of ionic crystals based on the VIB prescription in which the electronic variational parameters are treated as dynamical variables within the adiabatic approximation. With the computed phonon density of states, the structural parameters of a crystal can then be determined at any temperature and pressure (or stress condition) by mini- mizing the quasiharmonic Gibbs free energy with respect to the structural parameters along the electronic adiabatic surface. The complete pressure- and temperature-dependent thermal and elastic properties of a crystal can also be determined. We have applied the VIB procedure to calculate the equations of state, elastic properties, phonon dispersion relations, and phase stability of the alkaline earth oxides (MgO, CaO, and SrO). The calculated properties are generally in quite good agreement with experiment. Spherical charge deformation within the ionic description is responsible for a considerable improvement over rigid-ion elastic constants and bulk moduli, but some of the remaining discrepancies cannot in principle be eliminated within a spherical-ion model. In particular, the effect of spherical breathing on the optical vibrational frequencies is not dramatic and further improvement (e.g., reduction of LO-TO splitting) will require the inclusion of nonspherical charge relaxation (polarizability).",
author = "Andrew Chizmeshya and Zimmermann, {Frank M.} and LaViolette, {Randall A.} and George Wolf",
year = "1994",
doi = "10.1103/PhysRevB.50.15559",
language = "English (US)",
volume = "50",
pages = "15559--15574",
journal = "Physical Review B-Condensed Matter",
issn = "0163-1829",
publisher = "American Institute of Physics Publising LLC",
number = "21",

}

TY - JOUR

T1 - Variational charge relaxation in ionic crystals

T2 - An efficient treatment of statics and dynamics

AU - Chizmeshya, Andrew

AU - Zimmermann, Frank M.

AU - LaViolette, Randall A.

AU - Wolf, George

PY - 1994

Y1 - 1994

N2 - An ab initio variationally induced breathing (VIB) description of the energetics and dynamics of ionic crystals is developed within the spherical-ion-pair approximation and the Gordon-Kim ansatz. In the VIB method, the total crystal charge density is given by the overlap of Kohn-Sham ionic charge densities. Spherical charge deformation of the ions is accomplished by adding an effective many-body crystal site potential to the atomic potential used in the Kohn-Sham procedure. Parameters defining this effective site potential are treated variationally such that the total crystal electronic energy is a minimum for a given structural configuration. We have explored a number of forms of the effective potential, including Watson-sphere types where both the charge and radius of the sphere are variationally determined. In addition, we have investigated anion-cation charge transfer by incorporating an additional variational parameter. We also present a formulation of the lattice dynamics of ionic crystals based on the VIB prescription in which the electronic variational parameters are treated as dynamical variables within the adiabatic approximation. With the computed phonon density of states, the structural parameters of a crystal can then be determined at any temperature and pressure (or stress condition) by mini- mizing the quasiharmonic Gibbs free energy with respect to the structural parameters along the electronic adiabatic surface. The complete pressure- and temperature-dependent thermal and elastic properties of a crystal can also be determined. We have applied the VIB procedure to calculate the equations of state, elastic properties, phonon dispersion relations, and phase stability of the alkaline earth oxides (MgO, CaO, and SrO). The calculated properties are generally in quite good agreement with experiment. Spherical charge deformation within the ionic description is responsible for a considerable improvement over rigid-ion elastic constants and bulk moduli, but some of the remaining discrepancies cannot in principle be eliminated within a spherical-ion model. In particular, the effect of spherical breathing on the optical vibrational frequencies is not dramatic and further improvement (e.g., reduction of LO-TO splitting) will require the inclusion of nonspherical charge relaxation (polarizability).

AB - An ab initio variationally induced breathing (VIB) description of the energetics and dynamics of ionic crystals is developed within the spherical-ion-pair approximation and the Gordon-Kim ansatz. In the VIB method, the total crystal charge density is given by the overlap of Kohn-Sham ionic charge densities. Spherical charge deformation of the ions is accomplished by adding an effective many-body crystal site potential to the atomic potential used in the Kohn-Sham procedure. Parameters defining this effective site potential are treated variationally such that the total crystal electronic energy is a minimum for a given structural configuration. We have explored a number of forms of the effective potential, including Watson-sphere types where both the charge and radius of the sphere are variationally determined. In addition, we have investigated anion-cation charge transfer by incorporating an additional variational parameter. We also present a formulation of the lattice dynamics of ionic crystals based on the VIB prescription in which the electronic variational parameters are treated as dynamical variables within the adiabatic approximation. With the computed phonon density of states, the structural parameters of a crystal can then be determined at any temperature and pressure (or stress condition) by mini- mizing the quasiharmonic Gibbs free energy with respect to the structural parameters along the electronic adiabatic surface. The complete pressure- and temperature-dependent thermal and elastic properties of a crystal can also be determined. We have applied the VIB procedure to calculate the equations of state, elastic properties, phonon dispersion relations, and phase stability of the alkaline earth oxides (MgO, CaO, and SrO). The calculated properties are generally in quite good agreement with experiment. Spherical charge deformation within the ionic description is responsible for a considerable improvement over rigid-ion elastic constants and bulk moduli, but some of the remaining discrepancies cannot in principle be eliminated within a spherical-ion model. In particular, the effect of spherical breathing on the optical vibrational frequencies is not dramatic and further improvement (e.g., reduction of LO-TO splitting) will require the inclusion of nonspherical charge relaxation (polarizability).

UR - http://www.scopus.com/inward/record.url?scp=0001733678&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001733678&partnerID=8YFLogxK

U2 - 10.1103/PhysRevB.50.15559

DO - 10.1103/PhysRevB.50.15559

M3 - Article

AN - SCOPUS:0001733678

VL - 50

SP - 15559

EP - 15574

JO - Physical Review B-Condensed Matter

JF - Physical Review B-Condensed Matter

SN - 0163-1829

IS - 21

ER -