Los datos de equilibrio vapor–líquido (punto de ebullición) del sistema FeSO4–Fe2(SO4)3–H2SO4–H2O se midieron mediante el método ebulliométrico cuasiestático a (30, 60, 90 y 101,3) kPa. Se encontró que los puntos de ebullición aumentaban con el aumento de las concentraciones tanto de las sales como del ácido. Mediante el uso de los datos termodaynámicos (ΔGf, 298. 15K0,ΔHf, 298.15K0 y Sf,298.15K0) de iones complejos férricos (FeSO4+, Fe(SO4)2– y FeHSO42+) predichos por la ecuación de Helgeson, se desarrolló un modelo termodinámico completo para el sistema con el modelo de electrolito de disolvente mixto (MSE). Los parámetros del modelo se determinaron mediante la regresión de nuestros propios puntos de ebullición experimentales y los coeficientes osmóticos estequiométricos (ers) de la literatura del sistema Fe2(SO4)3–H2SO4–H2O. Las desviaciones absolutas medias entre los puntos de ebullición retrocedidos y los datos experimentales fueron (0,30 y 0.36) K para los sistemas Fe2(SO4)3–H2SO4–H2O y FeSO4–Fe2(SO4)3–H2SO4–H2O, respectivamente, mientras que la desviación relativa promedio para el coeficiente osmótico estequiométrico fue de solo 0,80 %. El modelo con parámetros recién obtenidos se utilizó para predecir la distribución de especies portadoras de hierro (Fe3+, FeSO4+, Fe(SO4)2-y FeHSO42+) como funciones de temperatura y concentración de hierro. Además, las gráficas del potencial redox de Fe(III)/Fe(II) obtenidas de la literatura versus ln(aFe3+/aFe2+) calculadas por los parámetros recién obtenidos son líneas rectas, lo que indica el buen rendimiento de los parámetros.