We revisit the numerical solutions of vibrational eigenstates of weakly bound homonuclear and heteronuclear noble gas pairs by applying a Fortran program based on the Numerov method. The harmonic, Lennard-Jones (LJ), and Improved Lennard- Jones (ILJ) potential models have been implemented to represent the potential energy curves (PECs). The obtained vibrational energies spectrum was tested on the experimental data and accurateab initiocalculations at CCSD(T)/CBS level. The vibrational eigenvalues and eigenfunctions can be reproduced accurately within the ILJ potential model. Moreover, considering from the calculated lifetime of van der Waals (vdW) complexes, the implementation of ILJ rather than standard LJ potential model has a significant impact on the systems dynamics by providing more representative atomic trajectories when the function is incorporated in force fields for molecular dynamics (MD) simulations.Overall, the ILJ function is the best suited potential model for the representation of vibrational motions and the determination of vibrational energy levels of weakly bound systems, both at equilibrium and non-equilibrium conditions.