In this paper, we present numerical simulations that demonstrate the effect of the particular choice of the equation of state (EoS) relating the surfactant concentration to the surface tension in surfactant-driven thin liquid films. Previous choices of the model EoS have been an ad-hoc decreasing function. Here, we instead propose an empirically motivated EoS; this provides a route to resolve some discrepancies and raises new issues to be pursued in future experiments. In addition, we test the influence of the choice of initial conditions and values for the non-dimensional groups. We demonstrate that the choice of EoS improves the agreement in surfactant distribution morphology between simulations and experiments, and influences the dynamics of the simulations. Because an empirically motivated EoS has regions with distinct gradients, future mathematical models may be improved by considering more than one timescale. We observe that the non-dimensional number controlling the relative importance of gravitational versus capillary forces has a larger influence on the dynamics than the other non-dimensional groups, but is nonetheless not a likely cause of discrepancy between simulations and experiments. Finally, we observe that the experimental approach using a ring to contain the surfactant could affect the surfactant and fluid dynamics if it disrupts the intended initial surfactant distribution. However, the fluid meniscus itself does not significantly affect the dynamics.