Theory of polarization of crystalline solids

We consider the change in polarization ΔP which occurs upon making an adiabatic change in the Kohn-Sham Hamiltonian of the solid. A simple expression for ΔP is derived in terms of the valence-band wave functions of the initial and final Hamiltonians. We show that physically ΔP can be interpreted as a displacement of the center of charge of the Wannier functions. The formulation is successfully applied to compute the piezoelectric tensor of GaAs in a first-principles pseudopotential calculation.