Gradients in valence bond theory

A gradient method for very general valence bond (VB) wavefunctions is presented. This method introduces the electronic energy as a Lagrange multiplier, and evaluates the contributions of the derivatives of the normalisation and of the first- and second-order cofactors present in the VB energy expression. The correctness of the method is illustrated with classic and breathing-orbital VB calculations on the HF molecule.