Exploiting the Hessian for a Better Convergence of the SCF-RDMFT Procedure

One-body reduced density matrix functional theory provides an alternative to density functional theory, which is able to treat static correlation while keeping a relatively low computation scaling. Its disadvantageous cost comes mainly from a slow convergence of the self-consistent energy optimization. To improve on that problem, we propose in this work the use of the Hessian of the energy, including the coupling term. We show that using the exact Hessian is very effective at reducing the number of iterations. However, since the exact Hessian is too expensive to use in practice, we propose an approximation based on an inexpensive exact part and BFGS updates.