attracted considerable experimental [4,5] and theoretical attention [6][7][8][9][10].Like pyrite [11], mackinawite (FeS) is an exceptional scavenger for heavy metals in water [12,13], while the similarity between its structural unit [14] and reactive biological clusters [15] has suggested that FeS structures could have catalysed the formation of the first prebiotic molecules in deep sea hydrothermal vents [16,17]. However, contrary to FeS 2 , the interface between water and FeS has not been widely investigated [18,19], in part because of the lack of force field parameters for the system.In previous work, we have derived a force field, based on pairwise interactions, for molecular mechanics (MM) simulations of mackinawite [20]. In this model, each interaction can be written as a function of the interatomic distances r ij :where A, ρ, and C are the coefficients of Buckingham potentials, and q i the electrostatic charges of the atoms. In addition, to take into account polarisability effects, each sulphur atom is represented by a core and a massless shell (S-S shell ), which interact through a harmonic potential [21]:where k s is the shell force constant, and r ij the core-shell distance. The parameters, reported in Table 1, have highlighted the stability of the (001) surface compared to all other surfaces, in agreement with ab initio calculations [22] and experimental findings [23].Here, we expand the FeS interatomic potential model to include its interaction with water, the latter described byAbstract We introduce a force field for the description of the mackinawite/water interface, which we derive by refining, and consistently merging with the SPC/Fw model of water, a set of existing interatomic potentials for the mineral. The thermal behaviour predicted for bulk mackinawite is in good agreement with experiment. The adsorption of water on the low-index surfaces of mackinawite reproduces results from density functional theory calculations at different water coverages, while the behaviour of water intercalated into the mineral is also remarkably similar to that from ab initio results. The force field is therefore suitable to model bulk mackinawite and its surfaces in an aqueous environment.