Berry-phase treatment of the homogeneous electric field perturbation in insulators

A perturbation theory of the static response of insulating crystals to homogeneous electric fields, that combines the modern theory of polarization (MTP) with the variation-perturbation framework is developed, at unrestricted order of perturbation. First, we address conceptual issues related to the definition of such a perturbative approach. In particular, in our definition of an electric-field-dependent energy functional for periodic systems, the position operator appearing in the perturbation term is replaced by a Berry-phase expression, along the lines of the MTP. Moreover, due to the unbound nature of the perturbation, a regularization of the Berry-phase expression for the polarization is needed in order to define a numerically-stable variational procedure. Regularization is achieved by means of discretization, which can be performed either before or after the perturbation expansion. We compare the two possibilities and apply them to a model tight-binding Hamiltonian. Lowest-order as well as generic formulas are presented for the derivatives of the total energy, the normalization condition, the eigenequation, and the Lagrange parameters.