Proportions of the gas and solid species as a function of temperature, calculated by thermodynamic equilibrium using the HSC software for four different experimental runs with different redox conditions (Rad-100, Rad-75L, Tan-95, Rad-50). The molecular species considered in the calculation are mostly common with Matsuno et al. (2021) with addition with following S bearing species, where (g) denotes gas and (s) denotes solids; S(s), S(g), S2(g), S₃(g), S₄(g), S₅(g), S₆(g), S₇(g), S₈(g), S₂O(g), SO(g), SO₂(g), SO₃(g), SiS(g), SiS(s), SiS₂(g), SiS₂(s), Fe_0.877S(s), FeS(s), FeS(g), FeS₂(s), FeSO₄(s), Fe₂(SO₄)₃(s), MgS(g), MgS(s),MgSO₄(s). The abundance is normalized by the amount of Si in the starting material.