AbstractThe Holland and Powell internally consistent data set version 5.5 has been augmented to include pyrite, troilite, trov (Fe0.875S), anhydrite, H2S, elemental S and S2 gas. Phase changes in troilite and pyrrhotite are modelled with a combination of multiple end‐members and a Landau tricritical model. Pyrrhotite is modelled as a solid solution between hypothetical end‐member troilite (trot) and Fe0.875S (trov); observed activity–composition relationships fit well to a symmetric formalism model with a value for wtrot−trov of −3.19 kJ mol−1. The hypothetical end‐member approach is required to compensate for iron distribution irregularities in compositions close to troilite. Mixing in fluids is described with the van Laar asymmetric formalism model with aij values for H2O–H2S, H2S–CH4 and H2S–CO2 of 6.5, 4.15 and 0.045 kJ mol−1 respectively. The derived data set is statistically acceptable and replicates the input data and data from experiments that were not included in the initial regression. The new data set is applied to the construction of pseudosections for the bulk composition of mafic greenschist facies rocks from the Golden Mile, Kalgoorlie, Western Australia. The sequence of mineral assemblages is replicated successfully, with observed assemblages predicted to be stable at X(CO2) increasing with increasing degree of hydrothermal alteration. Results are compatible with those of previous work. Assemblages are insensitive to the S bulk content at S contents of less than 1 wt%, which means that volatilization of S‐bearing fluids and sulphidation are unlikely to have had major effects on the stable mineral assemblage in less metasomatized rocks. The sequence of sulphide and oxide phases is predicted successfully and there is potential to use these phases qualitatively for geobarometry. Increases in X(CO2) stabilized, in turn, pyrite–magnetite, pyrite–hematite and anhydrite–pyrite. Magnetite–pyrrhotite is predicted at temperatures greater than 410 °C. The prediction of a variety of sulphide and oxide phases in a rock of fixed bulk composition as a function of changes in fluid composition and temperature is of particular interest because it has been proposed that such a variation in phase assemblage is produced by the infiltration of multiple fluids with contrasting redox state. The work presented here shows that this need not be the case.