def METHOD_NAME():
    fname = "rot_Be_STO.wfnoj.xml"
    basis, mo = read_qmcpack.parse_qmc_wf(fname, ["Be"])
    wf = Wavefunction_Be_STO(basis["Be"], mo)
    VP = np.zeros(wf.param_size)
    r = np.array([[0.7, 2.0, 3.0], [1.2, 1.5, 0.5], [1.5, 1.6, 1.5], [0.7, 1.0, 1.2]])
    run_qmc(r, wf, VP, nstep=100, nblock=10)