Abstract Accelerating analytical energy derivative calculations for large systems remains challenging within wave function-based quantum chemistry frameworks. The pair-natural orbital (PNO) method provides a powerful approach to reduce computational costs while maintaining accuracy. The authors recently reported a fully analytical nuclear gradient theory for the second-order Møller–Plesset perturbation theory using orthogonal localized virtual molecular orbital-based PNOs. This paper provides an overview of the theoretical framework and describes its extension to electronic response property calculations.
Hayashi et al. (Thu,) studied this question.