We have developed a force field to describe the interaction of alkanethiols HS(CH2)n-1 CH3 (Cn for short) by fitting a set of ∼220 interaction energies for dimers of Cn (with n = 1,2,⋯6) and CH4 molecules obtained from second-order Møller-Plesset perturbation theory calculations. The derived force field, based on a sum of so-called exp-6 pairwise potentials and effective Coulomb interaction potential between the HS- heads, predicts very well the interaction energies for dimers and trimers of alkanethiols not included in the input database for the fit. Also the calculated minimum energy tilt angle of the alkyl chains for a hexagonal arrangement of alkanethiolates with a nearest neighbor distance of 5 Å is in good agreement with the available experimental data for a √3 × √3 self-assembled monolayer (SAM) on Au(111). Thus, the derived force field might be suitable for large-scale molecular dynamics and/or Monte Carlo simulations to predict the structure, stability and/or kinetics of SAMs on other surfaces. © 2010 the Owner Societies.