The paper is organized as follows. In section II, expressions for the electric dipole-magnetic dipole and electric dipole-electric quadrupole polarizabilities are derived that satisfy periodic boundary conditions. Section III describes the computational procedure for the simulations. In section IV, simulations on cluster models and crystalline systems are presented. Section V provides concluding remarks.