Effects of nanoconfinement and surface charge on iron adsorption on mesoporous silica

Here, we present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(II) adsorption onto negatively charged SiO2 surfaces and Fe(III) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(II) at deprotonated surface sites (SiO-), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.