On the calculation of natural orbitals by perturbation theory