Local-density-functional total energy gradients in the linear combination of Gaussian-type orbitals method.