Experimental vapor–liquid equilibrium (VLE) data for binary systems of ethanol or water with four ionic liquids (ILs) comprising 1-butyl-3-methylimidazolium cation and different anions were measured at atmospheric pressure within a broad region of IL concentrations, IL mole fraction ranging between 0 and 0.7 or even 0.87. Siwoloboff procedure of vapor–liquid phase transition temperature measurement was applied to estimate VLE of the binary mixtures. Obtained t–x data were fitted assuming ideal vapor and non-ideal liquid phase behavior, using the original NRTL equation to calculate component activity coefficients. VLE description results are presented in form of equilibrium diagrams and as the variation of activity coefficients with mixture composition. Quality of the VLE data fitting was assessed based on root mean square deviations and mean deviations in temperature. Some experimental binary t–x data sets were compared to those available in literature (measured in much narrower concentration range) or predicted from ternary VLE data. The selection of ILs in this study was related to the separation of ethanol−water mixture by extractive distillation. Principal IL selection criterium was its ability to affect the ethanol–water relative volatility. Thermodynamic description of the VLE data of IL-containing systems within a wide concentration range is an essential building block of the mathematical model of separation devices, especially in case of regeneration equipment where concentrated IL is present.