Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ion coordination”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of iron redox ratio on the structures of boroaluminosilicate glasses

Iron oxide is commonly found in natural or industrial glass compositions and can exist as ferrous (Fe 2+ ) or ferric (Fe 3+ ) species, with their ratios depending on glass composition, temperature, pressure and the redox reactions during the glass forming process. The iron redox ratio plays an important role on silicate glass structures and consequently various properties. This work aims to study the effect of iron oxide, and particularly the iron redox ratio, on the structures of borosilicate and boroaluminosilicate glasses using molecular dynamics simulations with newly developed iron potential parameters that are compatible with the borosilicate potentials. The results provide detailed cation coordination states of both iron species and the effect of redox ratio on boron coordination and other structural features. Particularly, competition for charge compensating modifier cations (such as Na + ) among the fourfold-coordinated cations such as B 3+ , Al 3+ , and Fe 3+ is investigated by calculating the cation–cation pair distribution functions and coordination preferential ratios. The results show that the trivalent ferric ions, with a shorter Fe–O bond distance and better defined first coordiation shell with mainly four-fold coordination, act as a glass former whereas the divalent ferrous ions mainly play the role of glass modifier. The ferrous/ferric ratio (Fe 2+ /Fe 3+ ) was found to affect the glass chemistry and hence glass properties by regulating the amount of four-coordinated boron, the fraction of non-briding oxygen and other features. Finally, the results are compared with available experimental data to gain insights of the complex structures and charge compensation schemes of the glass system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synthesis and structural characterization of a hydrated sodium–caesium tetracosatungstate(VI), Na 5 Cs 19 [W 24 O 84 ]·21H 2 O

Crystal formation of penta­sodium nona­deca­cesium tetra­cosa­tungstate(VI) heneikosahydrate, Na 5 Cs 19 [W 24 O 84 ]·21H 2 O, was successfully achieved by the conversion of [H 2 W 12 O 42 ] 10– through the addition of excess Cs + . The crystal structure comprising the toroidal isopolyoxidometalate is presented, as well as its Raman spectrum. Na 5 Cs 19 (H 2 O) 21 W 24 O 84 crystallizes in the rhombohedral space group R$\overline{3}$ with an obverse centering. The title compound represents the addition of a new member to the isopolytungstate family with mixed alkali counter-ions and contains rarely observed five-coordinate tungsten(VI) atoms in the [W 24 O 84 ] 24– anion (site symmetry C 3i ) arising from the conversion mediated by Cs + counter-ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2021

Ionic Contraction across the Lanthanide Series Decreases the Temperature-Induced Disorder of the Water Coordination Sphere

In liquid, temperature affects the structures of lanthanide complexes in multiple ways that depend upon complex interactions between ligands, anions, and solvent molecules. The relative simplicity of lanthanide aqua ions (Ln 3+ ) make them well suited to determine how temperature induces structural changes in lanthanide complexes. In this work, we performed a combination of ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) measurements, both at 25 and 90 °C, to determine how temperature affects the first- and second-coordination spheres of three Ln 3+ (Ce 3+ , Sm 3+ , and Lu 3+ ) aqua ions. AIMD simulations show first lanthanide coordination spheres that are similar at 25 and 90 °C, more so for the Lu 3+ ion that remains as eight-coordinate than for the Ce 3+ and Sm 3+ ions that change their preferred coordination number from nine (at 25 °C) to eight (at 90 °C). The measured EXAFS spectra are very similar at 25 and 90 °C, for the Ce 3+ , Sm 3+ , and Lu 3+ ions, suggesting that the dynamical disorder of the Ln 3+ ions in liquid water is sufficient such that temperature-induced changes do not clearly manifest changes in the structure of the three ions. Both AIMD simulations and EXAFS measurements show very similar structures of the first coordination sphere of the Lu 3+ ion at 25 and 90 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Precise Modulation of CO 2 Sorption in Ti 8 Ce 2 –Oxo Clusters: Elucidating Lewis Acidity of the Ce Metal Sites and Structural Flexibility

Investigating the structure-property correlation in porous materials is a fundamental and consistent focus in various scientific domains, especially within sorption research. Metal oxide clusters with capping ligands, characterized by intrinsic cavities formed through specific solid-state packing, demonstrate significant potential as versatile platforms for sorption investigations due to their precisely tunable atomic structures and inherent long-range order. This study presents a series of Ti 8 Ce 2 -oxo clusters with subtle variations in coordinated linkers and explores their sorption behavior. Notably, Ti 8 Ce 2 -BA (BA denotes benzoic acid) manifests a distinctive two-step profile during CO 2 adsorption, accompanied by a hysteresis loop. This observation marks a new instance within the metal oxide cluster field. Of intrigue, the presence of unsaturated Ce(IV) sites was found to be correlated with the stepped sorption property. Moreover, the introduction of an electrophilic fluorine atom, positioned ortho or para to the benzoic acid, facilitated precise control over gate pressure and stepped sorption quantities. Advanced in-situ techniques systematically unraveled the underlying mechanism behind this unique sorption behavior. The findings elucidate that robust Lewis base-acid interactions are established between CO 2 molecules and Ce ions, consequently altering the conformation of coordinated linkers. Conversely, the F atoms primarily contribute to gate pressure variation by influencing the Lewis acidity of the Ce sites. This research advances the understanding in fabricating metal-oxo clusters with structural flexibility and provides profound insights into their host-guest interaction motifs. Furthermore, these insights hold substantial promise across diverse fields and offer valuable guidance for future adsorbent designs grounded in fundamental theories of structure-property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Alkali Metal Ions Dictate the Structure and Reactivity of an Iron(II) Imido Complex

The presence of redox innocent metal ions has been proposed to modulate the reactivity of metal ligand multiple bonds; however, insight from structure/function relationships is limited. Here, alkali metal reduction of the Fe(III) imido complex [Ph 2 B( t BuIm) 2 Fe=NDipp] (1) provides the series of structurally characterized Fe(II) imido complexes [Ph 2 B( t BuIm) 2 Fe=NDippLi(THF) 2 ] (2), [Ph 2 B( t BuIm) 2 Fe=NDippNa(THF) 3 ] (3), and [Ph 2 B( t BuIm) 2 Fe=NDippK] 2 (4), in which the alkali metal cations coordinate the imido ligand. Structural investigations demonstrate that the alkali metal ions modestly lengthen the Fe=N bond distance from that in the charge separated complex [Ph 2 B( t BuIm) 2 Fe=NDipp][K(18-C-6)THF 2 ] (5), with the longest bond observed for the smallest alkali metal ion. In contrast to 5, the imido ligands in 2–4 can be protonated and alkylated to afford Fe(II) amido complexes. Combined experimental and computational studies reveal that the alkali metal polarizes the Fe=N bond, and the basicity of imido ligand increases according to 5 < 4 ≈ 3 < 2. The basicity of the imido ligands influences the relative rates of reaction with 1,4-cyclohexadiene, specifically by gating access to complex 5, which is the species that is active for HAT. All complexes 2–4 react with benzophenone form metastable Fe(II) intermediates that subsequently eliminate the metathesis product Ph 2 C=NDipp, with relative rates dependent on the alkali metal ion. By contrast, the same reaction with 5 does not lead to the formation of Ph 2 C=NDipp. Furthermore, these results demonstrate that the coordination of alkali metal ions dictate both the structure and reactivity of the imido ligand and moreover can direct the reactivity of reaction intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2022

Binding, Release and Functionalization of Intact Pnictogen Tetrahedra Coordinated to Dicopper Complexes

The bridging MeCN ligand in the dicopper(I) complexes [(DPFN)Cu 2 (μ,η 1 : η 1 -MeCN)][X] 2 (X=weakly coordinating anion, NTf 2 (1 a), FAl[OC 6 F 10 (C 6 F 5 )] 3 (1 b), Al[OC(CF 3 ) 3 ] 4 (1 c)) was replaced by white phosphorus (P 4 ) or yellow arsenic (As 4 ) to yield [(DPFN)Cu 2 (μ,η 2 : η 2 -E 4 )][X] 2 (E=P (2 a-c), As (3 a-c)). The molecular structures in the solid state reveal novel coordination modes for E 4 tetrahedra bonded to coinage metal ions. Experimental data and quantum chemical computations provide information concerning perturbations to the bonding in coordinated E 4 tetrahedra. Reactions with N-heterocyclic carbenes (NHCs) led to replacement of the E 4 tetrahedra with release of P 4 or As 4 and formation of [(DPFN)Cu 2 (μ,η 1 : η 1 - Me NHC)][X] 2 (4 a,b) or to an opening of one E-E bond leading to an unusual E 4 butterfly structural motif in [(DPFN)Cu 2 (μ,η 1 : η 1 -E 4 Dipp NHC)][X] 2 (E=P (5 a,b), E=As (6)). With a cyclic alkyl amino carbene ( Et CAAC), cleavage of two As-As bonds was observed to give two isomers of [(DPFN)Cu 2 (μ,η 2 : η 2 -As 4 Et CAAC)][X] 2 (7 a,b) with an unusual As 4 -triangle+1 unit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies

An understanding of the fundamentals of the reaction between CuO with trace amounts of H 2 S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu 2-x S compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298–383 K) and from CuO sorbents with different crystallite sizes (2.8–40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu–S bonds and higher delocalization of the positive charge of the Cu center into (S 1– )2. This distorted CuS layer dominates the XAS signal at lower temperatures (298–323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S 1– )2 on both CuO (111) and ($\bar{1}11$) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Conformational Isomerization of Imide Anions Governs Solvation and Transport in Water-in-Salt Electrolytes

The behavior of highly concentrated electrolytes departs radically from the dilute-solution theory, yet the molecular origin of this transformation remains unresolved. Here, we identify the conformational isomerization of molecular ions as a decisive, previously unrecognized control parameter governing structure and transport in crowded aqueous electrolytes. Across a series of fluorosulfonimide anions, we show that increasing concentration drives a collective shift from extended transoid to compact cisoid conformers, revealed by small-angle X-ray scattering, vibrational spectroscopy, pulsed-field gradient NMR, and molecular dynamics simulations. This conformational transition triggers a collapse of the hydrogen-bonded water network and the emergence of densely packed ionic domains with confined water, producing a qualitative change in Li+ transport from solvent-mediated diffusion to network-confined hopping. Anion size and asymmetry systematically tune the onset of this transition, demonstrating that molecular geometry dictates mesoscale organization and dynamics in the ion-rich regime. Our results establish ion conformation, not merely composition or coordination, as a fundamental thermodynamic variable in concentrated solutions, providing a chemical framework that unifies solvation structure and transport in water-in-salt electrolytes and suggesting new principles for designing dense ionic media.

Nguyen, Huong TD↗
DOE OSTI2024

Molecular insights into the in situ early-stage assembly of metal–organic frameworks on cellulose nanofibrils

The integration of metal–organic frameworks (MOFs) into renewable matrices, including those derived from cellulose, to create cellulose/MOF hybrids has attracted significant interest due to the synergistic combination of cellulosic biopolymer properties and MOFs' multifunctional features. The interfacial interactions between cellulose and MOF not only ensure stable bonding, but ultimately also determine the bulk physical properties of the macroscopic composites. However, the mechanistic understanding of the in situ assembly between the two materials remains unclear. In this study, we employ a combination study of synthesis, experimental characterization, and molecular dynamics simulations to explore the early-stage assembly of a cellulose/MOF hybrid. Our results revealed that the growth of MOF clusters on the 2,2,6,6-tetramethyl-1-piperidinyloxy oxidized cellulose nanofibrils (TOCNF) follows an inhomogeneous sequential transformation pathway. The carboxylates of TOCNF form coordination-like bonds with the metal ions, while the hydroxyl groups of TOCNF form hydrogen bonds with MOF ligands. These interactions provide the initial nucleation sites that mediate the growth of MOF clusters and also guide the assembly of larger MOF clusters onto the TOCNF substrate. In conclusion, the fundamental insights into the in situ assembly of MOF nanoparticles on cellulosic substrates are essential for the rational design of high-performance materials with tailored morphology and optimized properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

High performance FPGA embedded system for machine learning based tracking and trigger in sPhenix and EIC

We present a comprehensive end-to-end pipeline to classify triggers versus background events in this paper. This pipeline makes online decisions to select signal data and enables the intelligent trigger system for efficient data collection in the Data Acquisition System (DAQ) of the upcoming sPHENIX and future EIC (Electron-Ion Collider) experiments. Starting from the coordinates of pixel hits that are lightened by passing particles in the detector, the pipeline applies three-stage of event processing (hits clustering, track reconstruction, and trigger detection) and labels all processed events with the binary tag of trigger versus background events. The pipeline consists of deterministic algorithms such as clustering pixels to reduce event size, tracking reconstruction to predict candidate edges, and advanced graph neural network-based models for recognizing the entire jet pattern. In particular, we apply the message-passing graph neural network to predict links between hits and reconstruct tracks and a hierarchical pooling algorithm (DiffPool) to make the graph-level trigger detection. We obtain an impressive performance (≥70% accuracy) for trigger detection with only 3200 neuron weights in the end-to-end pipeline. We deploy the end-to-end pipeline into a field-programmable gate array (FPGA) and accelerate the three stages with speedup factors of 1152, 280, and 21, respectively.

Instruments & Instrumentation↗
DOE OSTI2021

Organic-inorganic hybrid perovskite compounds

Photoactive materials comprising organic-inorganic hybrid halide perovskite compounds are provided. Photovoltaic cells and light-emitting devices incorporating the photoactive materials into their light-absorbing and light-emitting layers, respectively, are also provided. The halide perovskites have an amAMX3 perovskite crystal structure, wherein am is an alkyl diamine cation, an aromatic diamine cation, an aromatic azole cation, a cyclic alkyl diamine cation or a hydrazinediium cation; A is a monovalent alkylammonium cation or an alkali metal cation; X is a halide ion or a combination of halide ions; and M is an octahedrally coordinated bivalent metal atom.

Ke, Weijun↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Resolving metal binding properties within subunits of a multimeric enzyme Mnx by surface induced dissociation and native ion mobility mass spectrometry

Multi-subunit enzymes function as coordinated assemblies. Yet most enzymatic assays measure the summed output of all populations in solution and cannot easily differentiate contributions of individual subunits. Native mass spectrometry detects intact protein complexes in the gas phase. Surface induced dissociation further releases subunits from protein complexes while retaining compact conformations and bound ligands. Combined with ion mobility, the released subunits can then be carefully monitored for more in-depth structural analysis. Mnx is a unique bacterial multicopper oxidase complex that oxidizes Mn(II) to form MnO 2 minerals, and is composed of three subunits: MnxG, a multicopper oxidase containing the active site, and two accessory proteins, MnxE and MnxF which also bind copper ions. Other known multicopper oxidases do not require accessory proteins, therefore the functions of MnxE and MnxF are not well understood. Here, we use native mass spectrometry with surface induced dissociation and ion mobility to characterize the metal binding properties of Mnx with two metals, catalytic Cu(II) and Mn(II) substrate. We demonstrate our assay can detect subtle structural changes within each subunit, which are presumably related to the allosteric mechanism. We also noticed that ionic strength and solution composition can impact metal binding and must be carefully investigated for such experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Advanced Interactive 3D Visualization Tool for Customizable Analyses of Tomography Datasets in Material Science

Current methods for visualizing and analyzing 3D tomography datasets in materials science often lack the interactivity and depth required for detailed structural insights. This limitation restricts a researchers' ability to accurately interpret complex data, which is critical for advancing material innovations and understanding structural properties. To address this issue, we have developed a novel, web-based interactive 3D visualization and analysis tool from the Trame framework that offers customizable features to enhance data interpretability. The tool allows users to adjust parameters such as visible range, slice planes, data rotation, and layering, providing a more detailed and dynamic view of complex structures. Its user-friendly web interface increases the accessibility and ease of use for both novice and experienced researchers, to visualize large volumetric datasets. The tool supports a diverse range of data formats, making it versatile for various research applications. Unique capabilities include real-time data manipulation, automated feature detection, context-sensitive feedback, and real-time volume calculations and distributions per sliced region or layer, alongside the ability to quickly generate high-quality screenshots and videos for presentations and reports. These advancements offer a comprehensive solution for enhanced 3D data exploration, significantly improving the analysis process and communication of results in materials science.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

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 SiO 2 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.

54 ENVIRONMENTAL SCIENCES↗