We have conducted ab initio molecular dynamics (MD) simulations to calculate the speciation of Zn(II)–HS complexes from ambient to hydrothermal-magmatic conditions (25–600 °C, up to 2000 bar). These theoretical calculations were complemented by X-ray absorption spectroscopy (XAS) measurements of Zn(II) in HS−-rich solutions at 200–500 °C and 600–1000 bar. The speciation and geometrical properties predicted by the ab initio MD simulations and the in situ XAS data are in excellent agreement. Upon heating from room temperature to 250 °C, Zn(II) speciation in HS−-rich solutions shows a transition from the sixfold octahedral hexaaquo complex [Zn(H2O)6]2+ to fourfold tetrahedral [Zn(HS)n(H2O)4−n]2−n complexes (n = 1–4). Ab initio MD simulations also show that at temperatures > 250 °C, the threefold trigonal-planar [Zn(HS)3]− complex becomes increasingly stable, and predominates in S-rich solutions; in contrast, chloro-complexes display a tetrahedral geometry at 25–500 °C, while trigonal planar ZnCl3− predominates at temperatures > 500 °C. The stability constants of Zn(II)–HS complexes were calculated by thermodynamic integration of constrained ab initio MD simulations at 200, 350 and 600 °C. The stability constants generated from this study predict that zinc can be transported by HS− at high temperature in reduced, neutral to alkaline solutions, while Zn chloride complexes dominate in acidic fluids.