Gold(I)-bisulfite complexation in hydrothermal nanodroplets: A molecular dynamics study

WCH Hui and KH Lemke, GEOCHIMICA ET COSMOCHIMICA ACTA, 406, 317-325 (2025).

DOI: 10.1016/j.gca.2025.01.015

Water nanodroplets present a unique environment for gold hydrothermal transport, with fluid properties in aqueous nanodroplets distinct from bulk liquid and vapor phases. By performing classical and ab initio molecular dynamics simulations, we have probed the stability of water nanodroplets (H2O)(n) (n = 100, 1000) at 25 degrees C and 100 degrees C. The solvation and complexation of gold(I)-bisulfite AuHSO3 in nanodroplet environments were also examined, with a particular focus on surface and interior solvation. Classical TIP4P/2005 molecular dynamics simulations reveal extreme densities in the interior of (H2O)(100) and (H2O)(1000) nanodroplets compared to droplet surface regions. At 25 degrees C, the interior region of (H2O)(100) exhibits fluctuating densities at 1.016-1.079 g/cm(3), with two maxima at 1.079 g/cm(3) and 1.074 g/cm(3), corresponding to pressures of similar to 2.23 kbar and 2.07 kbar, respectively; Reduced densities are predicted for the larger (H2O)(1000) systems, these being 1.013 g/cm(3) (25 degrees C, 370 bar) and 0.968 g/cm(3) (100 degrees C, 220 bar). The outer regions, on the other hand, featured densities intermediate between saturated liquid and vapor conditions, as part of a transition from liquid to vapor-like densities at the edge of the droplet. Born-Oppenheimer molecular dynamics simulations at 100 degrees C show that the gold(I)-bisulfite complex H2O-AuHSO3 maintains a near linear solvation structure (theta(O-Au-S) = 172 degrees-174 degrees) in bulk aqueous fluids and at surface and interior sites of (H2O)(100) nanodroplets. Distance constrained simulations reveal that, upon extension of the gold(I)-bisulfite Au-S contact (equilibrium r(Au-S) = 2.3 & Aring;), HSO3- is displaced by a water molecule, forming a two-water solvation shell around Au+. Thermodynamic integration gives gold(I)-bisulfite dissociation energies (Delta G) of 17.65 +/- 0.37 kcal/mol (bulk), 20.22 +/- 0.38 kcal/mol (nanodroplet surface), and 18.31 +/- 0.31 kcal/mol (nanodroplet interior). Our ab initio molecular dynamics results demonstrate that water nanodroplets are stable at hydrothermal conditions and would play an important role in the speciation and transport of gold in volcanic and hydrothermal vapors.

Return to Publications page