Two models are presented based on the phase-field methodology to simulate thin film growth during physical vapor deposition (PVD), including subsurface microstructure evolution, for isotropic single-phase polycrystalline materials. The first model couples previous phase-field modeling efforts on ballistic deposition of single-phase materials and grain orientation evolution in polycrystalline materials in a sequential simulation algorithm. The second model incorporates both PVD and grain evolution dynamics into a single free energy functional for use in a phase-field model. To illustrate the capability of the proposed models in capturing combined thin film growth and subsurface grain evolution, PVD simulations of a generic single-phase polycrystalline metal are performed on substrates with different grain sizes. In both models, when the initial substrate grain sizes are smaller than the expected surface features, the thin film grains coarsen via grain boundary (GB) migration until the GBs become aligned with the valleys between the columnar surface features. Thus, each columnar feature is associated with a distinct subsurface grain, in qualitative agreement with experimental observations. Differences between the models arise when initial substrate grain sizes are larger than the surface columnar features. For example, when using the single free energy functional approach, grains contain noticeable internal low-angle variations, which are not captured using the coupled model.