This paper presents an efficient algorithm for energy gradients in valence bond self-consistent field(VBSCF) method with non-orthogonal orbitals.The frozen core approximation method is extended to the case of non-orthogonal orbitals.The expressions for the total energy and its gradients are presented by introducing auxiliary orbitals,where inactive orbitals are orthogonal,while active orbitals are non-orthogonal themselves but orthogonal to inactive orbitals.It is shown that our new algorithm has a low scaling of(Na+1)m4,where Na and m are the numbers of the active orbitals and basis functions,respectively,and is more efficient than the existing VBSCF algorithms.