In this paper, the use of the neural network (NN) method with exponential neurons for directlȳ tting ab initio data to generate potential energy surfaces (PESs) in sum-of-product form will be discussed. The utility of the approach will be highlighted using¯ts of CS 2 , HFCO, and HONO ground state PESs based upon high-level ab initio data. Using a generic interface between the neural network PES¯tting, which is performed in MATLAB, and the Heidelberg multi-conguration time-dependent Hartree (MCTDH) software package, the PESs have been tested via comparison of vibrational energies to experimental measurements. The review demonstrates the potential of the PES¯tting method, combined with MCTDH, to tackle high-dimensional quantum dynamics problems.