A nonlinear PDE model of thin film growth model, with or without slope selection, are presented in the talk. A global in time solution with Gevrey regularity is established for the one with slope selection. For the numerical simulation, an idea of convex-concave splitting of the corresponding physical energy is applied, which gives to an implicit treatment for the convex part and an explicit treatment for the concave part. That in turn leads to a numerical scheme with a non-increasing energy. A first order accurate linear splitting and a second order accurate linear iteration algorithm are also considered, with some numerical simulation results presented.