Solubility experiments for selenium were performed at high ionic strengths (1.0 and 2.0 mol dm À3 NaCl) under anoxic conditions (p O2 < 10 À6 atm) from undersaturation and oversaturation directions. Equilibrium was attained in 40 days. The dominant aqueous selenium species were identified as hydrogen selenide, HSe À , at pH between 5 and 8, and polyselenide, Se 4 2À , at pH between 9 and 13 by UV-Vis absorption spectrometry. A solubility-limiting solid was determined to be amorphous selenium at pH between 5 and 8, and crystalline selenium at pH between 9 and 13. The equilibrium constants oflog K 0 ¼ À16:67 AE 0:03 ð9 < pH < 13Þand ion interaction coefficients for HSe À and Se 4 2À versus Na þ , "(HSe À , Na þ ) = À0:01 AE 0:10 and "(Se 4 2À , Na þ ) = À0:03 AE 0:02, were determined with the specific ion interaction theory. The standard molar free energies of formation of amorphous selenium and polyselenide species were determined to be Á f G m 0 (Se(am)) = 6:0 AE 2:2 kJ mol À1 , Á f G m 0 (Se 2 2À ) = 111:64 AE 1:61 kJ mol À1 , Á f G m 0 (Se 3 2À ) = 98:68 AE 1:02 kJ mol À1 , and Á f G m 0 (Se 4 2À ) = 95:14 AE 0:17 kJ mol À1 .