Second-Derivative-Corrected FANPT for Robust Continuation of Nonlinear Wavefunction Equations

Abstract

We present a second-derivative-corrected extension of the Flexible Ansatz for N-body Perturbation Theory (FANPT) for solving nonlinear Flexible Ansatz for N-body Configuration Interaction (FANCI) wavefunction equations. The original quasilinear FANPT approximation neglects second- and higher-order derivatives of the determinant overlap with respect to wavefunction parameters. In this work, we retain the overlap Hessian and neglect only third- and higher-order parameter derivatives, thereby including the leading nonlinear response of the wavefunction ansatz while preserving the same response-matrix structure used in the original FANPT formulation. The resulting additional terms enter only through the constant vector of the response equations and are implemented for coupled-cluster wavefunctions in the FanPy/FANCI framework. The new approximation is tested on the Lithium Hydride molecule and the insertion of Beryllium into H$_2$ using seniority-restricted coupled-cluster wavefunctions in the STO-6G basis. The main advantage of the second-derivative correction is that it provides a better initial guess for solving the projected FANCI equations at the next point along the adiabatic connection. This improvement is diagnosed by comparing the FANPT-propagated parameters with the independently optimized FANCI parameters at the same value of $λ$. For nonlinear coupled-cluster ansatzes, especially when larger $λ$-steps are used, the corrected approximation reduces the same-$λ$ parameter deviations and suppresses large parameter-space excursions. These results show that the leading nonlinear overlap correction improves the reliability of FANPT as a continuation strategy for nonlinear wavefunction equations.