ps = np.array([list(conf.GetAtomPosition(i)) for i in range(mol.GetNumAtoms())])