Ab initio analytical infrared intensities for periodic systems through a coupled perturbed Hartree-Fock/Kohn-Sham method.

A fully analytical method for calculating Born charges and, hence, infrared intensities of periodic systems, is formulated and implemented in the CRYSTAL program, which uses a local gaussian type basis set. Our efficient formalism combines integral gradients with first-order coupled perturbed Hartree-Fock/Kohn Sham electronic response to an electric field. It avoids numerical differentiation with respect to wave vectors, as in some Berry phase approaches, and with respect to atomic coordinates. No perturbation equations for the atomic displacements need to be solved. Several tests are carried out to verify numerical stability, consistency in one, two, and three dimensions, and applicability to large unit cells. Future extensions to piezoelectricity and Raman intensities are noted.