Molecular dynamics simulations of a slab exposing a dominant crystal facet of ibuprofen, generated with the aid of the functionalities implemented in VESTA and solvated using the insert-molecules utility as implemented in Gromacs 5.1.4 in each of the ten different solvents previously mentioned, were performed. The slab thickness was set up so as the volume occupied by the crystal is half of the volume occupied by the solvent in order to prevent surface-surface interactions through periodic boundary conditions. The Gen-eralised Amber Force Field (GAFF) was used to represent system properties. For all systems considered in this work GAFF is able to reproduce properties consistent with experimental data. In support of this we report solvent densities in the Supplementary Material. Force field parameters for solvent molecules were obtained from the Virtual Chemistry solvent database or parametrised following the standard Amber procedure with antechamber. A standard cut-off distance of 1.0 nm for the non-bonded interactions was chosen, along with including long-range intermolecular interactions using the Particle Mesh Ewald (PME) approach. For computational efficiency, a time step of 2 fs was used. Temperature and pressure control have been implemented through the use of the Bussi-Donadio-Parrinello thermostat, Berendsen barostat 28 and Parrinello-Rahman barostat.