def METHOD_NAME(self):
    with lib.light_speed(20.):
        pt = mp.mp2.MP2(mf.x2c())
        pt.frozen = [0,1,10,11,12]
        gscan = pt.nuc_grad_method().as_scanner().as_scanner()
        e, g1 = gscan(mol)
        ps = pt.as_scanner()
        e1 = ps([[8 , (0. , 0.     , 0.)],
                 [1 , (0. , -0.757 , 0.5871)],
                 [1 , (0. ,  0.757 , 0.587)]])
        e2 = ps([[8 , (0. , 0.     , 0.)],
                 [1 , (0. , -0.757 , 0.5869)],
                 [1 , (0. ,  0.757 , 0.587)]])
        self.assertAlmostEqual(g1[1,2], (e1-e2)/0.0002*lib.param.BOHR, 5)