An efficient method to predict modified separation of cohesive energy density model (MOSCED) parameters for nonelectrolyte solids using electronic structure calculations in SMD and SM8 continuum solvents is proposed and applied to acetanilide, acetaminophen, and phenacetin. The resulting parameters are ultimately used to predict the equilibrium solubility in a range of solvents over a range of temperatures. By combining MOSCED with SMD and SM8, we are able to leverage the strengths of both methods while eliminating shortcomings that would prevent their use alone for solvent selection in design processes involving nonelectrolyte solid solutes. Comparing to 77 non-aqueous experimental solubilities of acetaminophen over the range 10-30 degrees C, the proposed method has an average absolute error of 0.03 and 0.04 mol fracs for SMD and SM8 regressed parameters, respectively. Aqueous solubilities of acetaminophen over this temperature range are predicted with an average error of 0.030 and 0.0023 mol fracs, respectively. (c) 2016 American Institute of Chemical Engineers AIChE J, 63: 781-791, 2017