The FeOs-AD builds on the implementation of phase equilibrium calculations in FeOs to provide implicit automatic differentiation of properties and phase equilibria based on Helmholtz energy equations of state. Derivatives can be determined for any inputs, like temperature or pressure, but also model parameters.