Wannier-based definition of layer polarizations in perovskite superlattices

In insulators, the method of Marzari and Vanderbilt [Phys. Rev. BPRBMDO0163-1829 56, 12847 (1997)10.1103/PhysRevB.56.12847] can be used to generate maximally localized Wannier functions whose centers are related to the electronic polarization. In the case of layered insulators, this approach can be adapted to provide a natural definition of the local polarization associated with each layer, based on the locations of the nuclear charges and one-dimensional Wannier centers comprising each layer. Here, we use this approach to compute and analyze layer polarizations of ferroelectric perovskite superlattices, including changes in layer polarizations induced by sublattice displacements (i.e., layer-decomposed Born effective charges) and local symmetry breaking at the interfaces. The method provides a powerful tool for analyzing the polarization-related properties of complex layered oxide systems.