Electric polarization in inhomogeneous crystals

Abstract

We derive the charge density up to second order in spatial gradient in inhomogeneous crystals using the semiclassical coarse graining procedure based on the wave packet method. It can be recast as divergence of polarization, whose first-order contribution consists of three parts, a perturbative correction to the original Berry connection expression, a topological part that can be written as an integral of the Chern-Simons 3-form, and a previously-unknown, quadrupolelike contribution. The topological part can be related to the quantized fractional charge carried by a vortex in two-dimensional systems. We then generalize our results to the multiband case and show that the quadrupolelike contribution plays an important role, as it makes the total polarization gauge-independent. Finally, we verify our theory in several model systems.