The system Fe 2 O 3 -SO 3 -H 2 O contains the most important minerals of acid mine drainage (AMD), iron oxides, and iron sulfates. For geochemical modeling of the AMD systems, reliable thermodynamic data for these phases are needed. In this work, we have determined thermodynamic data for the most acidic sulfates rhomboclase [(H 5 O 2 )Fe(SO 4 ) 2 ·2H 2 O or (H 3 O)Fe(SO 4 ) 2 ·3H 2 O] and the phase (H 3 O)Fe(SO 4 ) 2 . The actual compositions of the studied phases are (H 3 O) 1.34 Fe(SO 4 ) 2.17 (H 2 O) 3.06 (molecular mass of 344.919 g/mol) and (H 3 O) 1.34 Fe(SO 4 ) 2.17 (289.792 g/mol). Structural details for both phases were refined from synchrotron powder X-ray diffraction data. Enthalpies of formation were determined by acid-solution calorimetry. Low-temperature heat capacity was measured for rhomboclase by relaxation calorimetry but a critical analysis of entropies for several oxysalts showed that these data are too high. Entropies for both phases were estimated from a Kopp-rule algorithm. The enthalpies of formation and entropies were combined with previously published temperature-relative humidity brackets to generate an internally consistent thermodynamic data set for rhomboclase: Δ f H ° = –3202.03 kJ/mol, S ° = 378.7 J/(mol∙K); and for (H 3 O) 1.34 Fe(SO 4 ) 2.17 : Δ f H ° = –2276.25 kJ/mol, S ° = 253.2 J/(mol∙K). Solubility experiments at room temperature and at T = 4 °C agree well with previously reported data in the system Fe 2 O 3 -SO 3 -H 2 O. An inspection of the extended Pitzer model for Fe 3+ -SO 4 solutions shows that this model reproduces the general topology of the phase diagram, but the position of the calculated solubility curves deviates substantially from the experimental data. Solid-state 2 H MAS NMR spectra on deuterated rhomboclase show two isotropic chemical shifts δ iso ( 2 H) = of 8 ± 1 and 228 ± 1 ppm, assigned to D 5 O 2 + and Fe-OD 2 groups, respectively. Canonical ensemble (NVT) molecular dynamics simulations for (H 3 O)Fe(SO 4 ) 2 at T = 300 K showed that the H 3 O groups maintain their trigonal pyramidal geometry and perform different types of motion.