def METHOD_NAME(self):
    mol = gto.Mole()
    mol.atom = [
        ['O' , (0. , 0.     , 0.)],
        [1   , (0. , -0.757 , 0.587)],
        [1   , (0. , 0.757  , 0.587)] ]
    mol.basis = '631g'
    mol.build()
    mf = scf.UHF(mol)
    mf.conv_tol = 1e-14
    e0 = mf.scf()
    g1 = mf.Gradients().grad()