Predicting the equilibrium solubility of solid polycyclic aromatic hydrocarbons and dibenzothiophene using a combination of MOSCED plus molecular simulation or electronic structure calculations

<p>Methods to predict the equilibrium solubility of non-electrolyte solids are important for the design of novel separation processes. Here we demonstrate how conventional molecular simulation free energy calculations or electronic structure calculations in a continuum solvent, here SMD or SM8, can be used to predict parameters for the <u>MO</u>dified <u>S</u>eparation of <u>C</u>ohesive <u>E</u>nergy <u>D</u>ensity (MOSCED) method. The method is applied to the solutes naphthalene, anthracene, phenanthrene, pyrene and dibenzothiophene, compounds of interested to the petroleum industry and for environmental remediation. Adopting the melting point temperature and enthalpy of fusion of these compounds from experiment, we are able to predict equilibrium solubilities. Comparing to a total of 422 non-aqueous and 193 aqueous experimental solubilities, we find the proposed method is able to well correlate the data. The use of MOSCED is additionally advantageous as it is a solubility parameter-based method useful for intuitive solvent selection and formulation.</p