Anew Helmholtz energy model of mixing for random copolymer solutions based on a close-packed lattice has been developed. The model contains three terms: the contribution of the athermal mixing of polymer chain and solvent, the Helmoltz energy of mixing in a multi-component Ising lattice where the interactions between segments is accounted for, and the contribution of the dissociation of the polymer and the association of monomers. The Guggenheim model, Yang et al.'s model and the sticky-point model of Cummings, Zhou and Stell are used respectively, for the above three contributions. Comparisons between Monte Carlo simulated coexistence curves with those predicted by various theories for random copolymer solutions with various chain lengths, chain compositions and inter-segment interaction parameters show that the agreement between simulations and the predictions of this work is nearly perfect. The model can be used satisfactorily to correlate the liquid-liquid equilibria of practical random copolymer solutions. (C) 2008 Elsevier B.V. All rights reserved.