We employ the popular all-atom optimized potential for liquid simulations, OPLSAA, force-field to model 17 different alcohols in the liquid state. Using the standard simulation protocol for few hundred nanosecond time periods, we find that 1-octanol, 1-nonanol, and 1-decanol undergo spontaneous transition to a crystalline state at temperatures which are 35-55 K higher than the experimental melting temperatures. Nevertheless, the crystal structures obtained from the simulations are very similar to those determined by X-ray powder diffraction data for several n-alcohols. Although some degree of deviations from the experimental freezing points are to be expected, for 1-nonanol and 1-decanol, the elevation of the freezing temperature warrants special attention because at room temperature, these alcohols are liquids; however, if simulated by the OPLSAA force-field, they will crystallize. This behavior is likely a consequence of exaggerated attractive interactions between the alkane chains of the alcohols. To circumvent this problem, we combined the OPLSAA model with the L-OPLS force-field. We adopted the L-OPLS parameters to model the hydrocarbon tail of the alcohols, whereas the hydroxyl head group remained as in the original OPLSAA force-field. The resulting alcohols stayed in the liquid state at temperatures above their experimental melting points, thus, resolving the enhanced freezing observed with the OPLSAA force-field. In fact, the mixed-model alcohols did not exhibit any spontaneous freezing even at temperatures much lower than the experimental values. However, a series of simulations in which these mixed-OPLSAA alcohols started from a coexistence configuration of the liquid and solid phases resulted in freezing transitions at temperatures 14-25 K lower than the experimental values, confirming the validity of the proposed model. For all of the other alcohols, the mixed model yields results very similar to the OPLSAA force-field and is in good agreement with the experimental data. Thus, for simulating alcohols in the liquid phase, the mixed-OPLSAA model is necessary for large (7 carbons and above) hydrocarbon chains.
