In this paper we present an approach aimed at performing many-body calculations of Born-effective charges of crystalline insulators by including the electron-correlation effects. The scheme is implemented entirely in the real space, using Wannier functions as single-particle orbitals. Correlation effects are computed by including virtual excitations from the Hartree-Fock mean field, and the excitations are organized as per a BetheGoldstone-like many-body hierarchy. The results of our calculations suggest that the approach presented here is promising. The Born effective charge (BEC) of a periodic solid is an important phenomenological quantity which connects the electronic structure of the system to its phononic properties.
1Of late, in the context of ferroelectric materials and their phase transitions, BEC has generated tremendous amount of interest.2 Using BEC, one can also describe the lattice dynamics, and its coupling to infrared radiation, in a simple intuitive manner. 3 Most of the modern calculations of BECs are based upon the Berry-phase-based theory of macroscopic polarization developed by King-Smith and Vanderbilt. 4 The aforesaid formalism is based upon single-particle orbitals, and, therefore, can be implemented in a straightforward manner within ab initio density-functional theory (DFT), 2,3 or the Hartree-Fock (HF) framework.5 As far as many-body calculations of polarization properties are concerned, Martin and co-workers have proposed several approaches which, to the best of our knowledge, have not been implemented within an ab initio methodology. 6,7 Filippetti and Spaldin have recently implemented an ab initio method aimed at including correlation effects by using a self-interaction-corrected (SIC) density-functional approach. 8 Recently, we have developed a wave-function-based ab initio methodology aimed at performing electronic structure calculations on crystalline insulators. [9][10][11][12] The approach uses Wannier functions as single-particle orbitals obtained at the Hartree-Fock level, which can subsequently be used to include the electron correlation effects, if needed. The approach has been applied to calculate ground state geometries, cohesive energies, and elastic properties of crystalline insulators at the Hartree-Fock level, 9,10 as well as at the correlated level. 11,12 Moreover, within the Berry-phase formalism of King-Smith and Vanderbilt, 4 we have also used our approach to compute the BECs of several ionic insulators at the HF level. 13 The purpose behind the present work is to use our Wannier-function-based methodology to perform correlated calculations of the BEC's of insulators. Since ours is a real-space approach, we start with the following expression for the electronic contribution to the polarization per unit cell (P e ) valid for insulators: 14where is a parameter governing the state of crystal (for the present case, it represents atomic displacements), ⍀ is the volume of the unit cell, q e is the electronic charge, N͑→ϱ͒, represents the total number of unit cells in th...