We have used density functional theory calculations to investigate the interaction between methylamine (CH3NH2) and the dominant surfaces of mackinawite (FeS), where the surface and adsorption properties of mackinawite have been characterized using the DFT-D2 method of Grimme. Our calculations show that while the CH3NH2 molecule only interacts weakly with the most stable FeS(001), it adsorbs relatively strongly on the FeS(011) and FeS(100) surfaces releasing energies of 1.26 eV and 1.51 eV, respectively. Analysis of the nature of the bonding reveals that the CH3NH2 molecule interacts with the mackinawite surfaces through the lone-pair of electrons located on the N atom. The electron density built up in the bonding region between N and Fe is very much what one would expect of covalent type of bonding. We observe no significant adsorption-induced changes of the FeS surface structures, suggesting that amine capping agents would not distort the FeS nanoparticle surfaces required for active heterogeneous catalytic reactions. The vibrational frequencies and the infrared spectra of adsorbed methylamine have been calculated and assignments for vibrational modes are used to propose a kinetic model for the desorption process, yielding a simulated temperature programmed desorption with a relative desorption temperature of <140 K at the FeS(011) surface and <170 K at FeS(100) surface.