In this work, we present the second version of the Donostia natural orbital functional software, an open-source program for natural orbital functional calculations. The new release incorporates improved optimization algorithms, capabilities for excited-state computations, support for ab initio molecular dynamics, and integration with the libcint library. DoNOF 2.0 also extends its property toolbox by enabling the evaluation of nonlinear optical responses, including static polarizabilities and higher-order hyperpolarizabilities via a finite-field Romberg–Richardson scheme.
Lew-Yee et al. (2026) studied this question.