Current-density implementation for calculating flexoelectric coefficients
Abstract
The flexoelectric effect refers to polarization induced in an insulator when a strain gradient is applied. We have developed a first-principles methodology based on density-functional perturbation theory to calculate the elements of the bulk, clamped-ion flexoelectric tensor. In order to determine the transverse and shear components directly from a unit cell calculation, we calculate the current density induced by the adiabatic atomic displacements of a long-wavelength acoustic phonon. Previous implementations based on the charge-density response required supercells to capture these components. Our density-functional-theory implementation requires the development of an expression for the current density that is valid for the case of nonlocal pseudopotentials, and long-wavelength phonon perturbations. We benchmark our methodology on simple systems of isolated noble gas atoms, and apply it to calculate the clamped-ion flexoelectric constants for a variety of technologically important cubic oxides. We also discuss some technical issues that are associated with the definition of current density in a nonlocal pseudopotential context, and their relevance to the calculation of macroscopic response properties of crystals.
Keywords
Cite
@article{arxiv.1802.06390,
title = {Current-density implementation for calculating flexoelectric coefficients},
author = {Cyrus E. Dreyer and Massimiliano Stengel and David Vanderbilt},
journal= {arXiv preprint arXiv:1802.06390},
year = {2018}
}
Comments
23 pages, 3 figures