In order to include proper orbital relaxation effects on the first-order IIP corrections, the resolvent operators R 10 and R 01 in Eqs. ( ) and ( ) can be replaced respectively by the operators C A and C B that generate the solutions of the coupled-perturbed Hartree-Fock (CPHF) or coupled-perturbed Kohn-Sham (CPKS) equations,