The polarization energy is defined using fragment electrical response functions (FERFs). For each fragment, the orbital hessian is calculated via coupled-perturbed self-consistent field (CPSCF) equations and is used to compute the multipole matrix M. The µ component of the fragment response matrix D, which couples the occupied and virtual orbitals, can be calculated from the multipole matrix and orbital Hessian.