Recently, an equation of state (EoS) for the Lennard-Jones truncated and shifted (LJTS) fluid has become available. As it describes metastable and unstable states well, it is suited for predicting density profiles in vapor-liquid interfaces in combination with density gradient theory (DGT). DGT is usually applied to describe interfaces in Cartesian one-dimensional scenarios. In the present work, the perturbed LJ truncated and shifted (PeTS) EoS is implemented into a three-dimensional phase field (PF) model which can be used for studying inhomogeneous gas-liquid systems in a more general way. The results are compared with the results from molecular dynamics simulations for the LJTS fluid that are carried out in the present work and good agreement is observed. The PF model can therefore be used to overcome the scale limit of molecular simulations. A finite element approach is applied for the implementation of the PF model. This requires the first and second derivatives of the PeTS EoS which are calculated using hyper-dual numbers. Several tests and examples of applications of the new PeTS PF model are discussed.