The slip and friction behavior of n-hexadecane, confined between organic friction modifier surfactant films adsorbed on hematite surfaces, has been studied using nonequilibrium molecular dynamics simulations. The influence of the surfactant type and coverage, as well as the applied shear rate and pressure, has been investigated. A measurable slip length is only observed for surfactant films with a high surface coverage, which provide smooth interfaces between well-defined surfactant and hexadecane layers. Slip commences above a critical shear rate, beyond which the slip length first increases with increasing shear rate and then asymptotes toward a constant value. The maximum slip length increases significantly with increasing pressure. Systems and conditions which show a larger slip length typically give a lower friction coefficient. Generally, the friction coefficient increases linearly with logarithmic shear rate; however, it shows a much stronger shear rate dependency at low pressure than at high pressure. Relating slip and friction, slip only occurs above a critical shear stress, after which the slip length first increases linearly with increasing shear stress and then asymptotes. This behavior is well-described using previously proposed slip models. This study provides a more detailed understanding of the slip of alkanes on surfactant monolayers. It also suggests that high coverage surfactant films can significantly reduce friction by promoting slip, even when the surfaces are well-separated by a lubricant.