SA . They are approximate Hessians and, more precisely, are exact firstorder derivatives of the approximate electronic gradient. The matrix A is no longer symmetric (and maybe non-positive definite), so the conventional preconditioned conjugate gradients method cannot be used for improved convergence. In this study, the conjugate gradients squared (CGS) method is employed. Preconditioning is also implemented according to Algorithm 4 in Ref.