The present paper proposes to use the group contribution (GC) polar perturbed-chain-statistical associating fluid theory (GC-PPC-SAFT) equation of state (EoS), that has already been used with success on various organic mixtures, and extend it to model simultaneously the liquidÀliquid equilibrium (LLE) and vaporÀliquid equilibrium (VLE) of hydrocarbons þ water systems, in wide ranges of pressure and temperature. Mixtures of water with aliphatics, aromatics, alcohols, carbon dioxide, and hydrogen sulfide have been investigated. Pure water is assumed associative (according to the 4C association scheme) and dipolar; the aromatic compounds are quadrupolar. Alcohols are autoassociative with a 3B association scheme. A cross-association between water and alcohols or H 2 S is taken into account. Cross association between water and other polar molecules (CO 2 or aromatic molecules) was also taken into account explicitly. Only one set of cross association parameters ε cross /k and k cross values were used for all the water þ aromatic mixtures considered here. ε cross /k was adjusted on experimental data, whereas k cross is set to the value found for pure water. For each system, the same binary interaction parameter k ij was used for simultaneous modeling LLE and VLE. This parameter was correlated to pseudo-ionization energy parameters for pure compounds through London's dispersion force theory, and reused from previous works [