This release modifies the _eval_xc_sf() to add the third-order derivatives of functionals, which supports the sftda_grad module in pyscf-forge package in the future. The module sftda_grad aims to calculate the analytical gradients of the spin-flip excited states.