Search

Using ab initio based molecular dynamics and electronic structure calculations, we show that Zn impurities in hydrated amorphous calcium carbonate (ACC) have a much lower coordination number than other divalent impurities due to covalent interactions between the 3d Zn shell and the oxygen atoms of the carbonate and water groups. The local structure around Zn in ACC, including the predicted low coordination number, is confirmed by X-ray absorption spectroscopy of synthetic Zn-bearing ACC. The strong Zn-O chemical interaction leads to substantial water dissociation and slightly disrupts the hydrogen bonding network. Implications of Zn 2+ incorporation for ACC stability are discussed.

Mineral-associated soil organic matter (SOM) is critical for stabilizing organic carbon and mitigating climate change. However, mineral-SOM interactions at the molecular scale, particularly synergetic adsorption through organic-organic interaction on the mineral surface known as organic multilayering, remain poorly understood. This study investigates the impact of organic multilayering on mineral-SOM interactions, by integrating macroscale experiments and molecular-scale simulations that assess the individual and sequential adsorption of major SOM compounds-lauric acid (lipid), pentaglycine (amino acid), trehalose (carbohydrate), and lignin onto soil minerals. Ferrihydrite, Al-hydroxide, and calcite are exposed to SOM compounds to determine adsorption affinities and binding energies. Results show that lauric acid has 20-40 times higher K d than pentaglycine, following the order K d (ferrihydrite) > K d (Al-hydroxide) ≫ K d (calcite). Molecular-scale simulations confirm that lauric acid has a higher binding energy (30.8 kcal/mol) on ferrihydrite than pentaglycine (6.0 kcal/mol), attributed to lipid hydrophobicity. The lower binding energy of pentaglycine results from its hydrophilic amide groups, facilitating partitioning into water. Sequential experiments examine how the first layer of lipid or amino acid affects the adsorption of carbohydrate/lignin, which show little or no individual adsorption affinities. Macroscale results reveal that lipid and amino acid adsorption induce ferrihydrite particle repulsion increasing reactive surface area and enhancing carbohydrate/lignin adsorption independently and synergistically through organic multilayering. Molecular-scale results reveal that amino acid adsorbed on ferrihydrite interacts more readily with lignin macroaggregates (preformed in solution) than with individual lignin units, indicating organic multilayering via H-bonding. These findings reveal the molecular mechanisms of SOM-mineral interactions, crucial for enhancing soil carbon stabilization.

Amelogenin is an intrinsically disordered protein essential to tooth enamel formation in mammals. Using advanced small angle X-ray scattering (SAXS) capabilities at synchrotrons and computational models, we revisited measuring the quaternary structure of murine amelogenin as a function of pH and phosphorylation at serine-16. The SAXS data shows that at the pH extremes, amelogenin exists as an extended monomer at pH 3.0 (R g  = 38.4 Å) and nanospheres at pH 8.0 (R g  = 84.0 Å), consistent with multiple previous observations. At pH 5.0 and above there was no evidence for a significant population of monomeric species. Instead, at pH 5.0, ∼80 % of the population is a heterogenous dimeric species that increases to ∼100 % at pH 5.5. The dimer population was observed at all pH > 5 conditions in dynamic equilibrium with a species in the pentamer range at pH < 6.5 and nanospheres at pH 8.0. At pH 8.0, ∼40 % of the amelogenin remained in the dimeric state. In general, serine-16 phosphorylation of amelogenin appears to modestly stabilize the population of the dimeric species., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

The features and the electronic character of the states for the Ni 2p x-ray photoelectron spectroscopy (XPS) of Ni(OH)2 were analyzed. This detailed analysis is based on ab initio molecular orbital wavefunctions for a cluster model of Ni(OH)2. The theory is validated by comparison with experiment. Then, advanced methods are used to explain and contrast the properties of different groups of ionic states. An important conclusion is that in most cases, the ionic states cannot be described with a single configuration or determinant. Despite this essential many-body character of the XPS, we demonstrate that it is possible to understand the origin of the main and satellite XPS features in terms of their orbital character., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We establish the connection between the measured small angle x-ray scattering signal and the charge-charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a0, and the inverse periodicity length, Q0. Above the KTs, we find a universal scaling of a0∝c-ζ/3 and Q0 ∝ c1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Ceria nanomaterials with facile Ce III/IV redox behavior are used in sensing, catalytic, and therapeutic applications, where inclusion of Ce III has been correlated with reactivity. Understanding assembly pathways of CeO 2 nanoparticles (NC-CeO 2 ) in water has been challenged by "blind" synthesis, including rapid assembly/precipitation promoted by heat or strong base. Here, we identify a layered phase denoted Ce-I with a proposed formula Ce IV (OH) 3 (NO 3 )· x H 2 O ( x ≈ 2.5), obtained by adding electrolytes to aqueous cerium ammonium nitrate (CAN) to force precipitation. Ce-I represents intermediate hydrolysis species between dissolved CAN and NC-CeO 2 , where CAN is a commonly used Ce IV compound that exhibits unusual aqueous and organic solubility. Ce-I features Ce-(OH) 2 -Ce units, representing the first step of hydrolysis toward NC-CeO 2 formation, challenging prior assertions about Ce IV hydrolysis. Structure/composition of poorly crystalline Ce-I was corroborated by a pair distribution function, Ce-L3 XAS (X-ray absorption spectroscopy), compositional analysis, and 17 O nuclear magnetic resonance spectroscopy. Formation of Ce-I and its transformation to NC-CeO 2 is documented in solution by small-angle X-ray scattering (SAXS) and in the solid-state by transmission electron microscopy (TEM) and powder X-ray diffraction. Morphologies identified by TEM support form factor models for SAXS analysis, evidencing the incipient assembly of Ce-I . Finally, two morphologies of NC-CeO 2 are identified. Sequentially, spherical NC-CeO 2 particles coexist with Ce-I , and asymmetric NC-CeO 2 with up to 35% Ce III forms at the expense of Ce-I , suggesting direct replacement.

Nanoparticle aggregates in solution controls surface reactivity and function. Complete dispersion often requires additive sorbents to impart a net repulsive interaction between particles. Facet engineering of nanocrystals offers an alternative approach to produce monodisperse suspensions simply based on facet-specific interaction with solvent molecules. Here, we measure the dispersion/aggregation of three morphologies of hematite (α-Fe 2 O 3 ) nanoparticles in varied aqueous solutions using ex situ electron microscopy and in situ small-angle x-ray scattering. We demonstrate a unique tendency of (104) hematite nanoparticles to maintain a monodisperse state across a wide range of solution conditions not observed with (001)- and (116)-dominated particles. Density functional theory calculations reveal an inert, densely hydrogen-bonded first water layer on the (104) facet that favors interparticle dispersion. Results validate the notion that nanoparticle dispersions can be controlled through morphology for specific solvents, which may help in the development of various nanoparticle applications that rely on their interfacial area to be highly accessible in stable suspensions.

Two-dimensional (2D) materials have attracted intense interest due to their potential for applications in fields ranging from chemical sensing to catalysis, energy storage, and biomedicine. Recently, peptoids, a class of biomimetic sequence-defined polymers, have been found to self-assemble into 2D crystalline sheets that exhibit unusual properties, such as high chemical stability and the ability to self-repair. The structure of a peptoid is close to that of a peptide except that the side chains are appended to the amide nitrogen rather than the α carbon. In this study, we investigated the effect of peptoid sequence on the mechanism and kinetics of 2D assembly on mica surfaces using in situ AFM and time-resolved X-ray scattering. We explored three distinct peptoid sequences that are amphiphilic in nature with hydrophobic and hydrophilic blocks and are known to self-assemble into 2D sheets. The results show that their assembly on mica starts with deposition of aggregates that spread to establish 2D islands, which then grow by attachment of peptoids, either monomers or unresolvable small oligomers, following well-known laws of crystal step advancement. Extraction of the solubility and kinetic coefficient from the dependence of the growth rate on peptoid concentration reveals striking differences between the sequences. The sequence with the slowest growth rate in bulk and with the highest solubility shows almost no detachment; i.e., once a growth unit attaches to the island edge, there is almost no probability of detaching. Furthermore, a peptoid sequence with a hydrophobic tail conjugated to the final carboxyl residue in the hydrophilic block has enhanced hydrophobic interactions and exhibits rapid assembly both in the bulk and on mica. These assembly outcomes suggest that, while the π-π interactions between adjacent hydrophobic blocks play a major role in peptoid assembly, sequence details, particularly the location of charged groups, as well as interaction with the underlying substrate can significantly alter the thermodynamic stability and assembly kinetics.

Understanding multiple lengthscale correlations in the pair distribution functions (PDFs) of aq. electrolytes is a persistent challenge. Here, the coordination chemistry of polyoxoanions supports an ion-network of cation-coordination polyhedra in NaNO 3(aq) and NaNO 2(aq) that induce long-range solution structure. Oxygen correlations associated with Na + -coordination polyhedra have two characteristics lengthscales; 3.5-5.5 Å and 5.5-7.5 Å, the latter solely associated oligomers. The PDF contraction between 5.5-7.5 Å observed in many electrolytes is attributed to the distinct O⋯O correlation found in dimers and dimer subunits within oligomers.

Incorporation of a Ba impurity in amorphous calcium carbonate (ACC) is shown with ab initio molecular dynamics simulations to have a long-range effect on its atomic-level structure and to be energetically favoured relative to incorporation in crystalline calcium carbonate polymorphs. The ability of carbonate ions to rotate and of ACC to undergo local density changes explain ACC's propensity for incorporating divalent metal impurities with a wide range of ionic radii. These findings provide an atomic-level basis for understanding the significant effects of low concentrations of impurities on the structure of ACC.

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 · n H 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

Amelogenin, the dominant organic component (>90%) in the early stages of amelogenesis, orchestrates the mineralization of apatite crystals into enamel. The self-association properties of amelogenin as a function of pH and protein concentration have been suggested to play a central role in this process; however, the large molecular weight of the self-assembled quaternary structures has largely prevented structural studies of the protein in solution under physiological conditions using conventional approaches. Here, using perdeuterated murine amelogenin (0.25 mM, 5 mg/mL) and TROSY-based NMR experiments to improve spectral resolution, we assigned the 1 H- 15 N spectra of murine amelogenin over a pH range (5.5 to 8.0) where amelogenin is reported to exist as oligomers (pH ≤∼6.8) and nanospheres (pH ≥∼7.2). The disappearance or intensity reduction of amide resonances in the 1 H- 15 N HSQC spectra was interpreted to reflect changes in dynamics (intermediate millisecond-to-microsecond motion) and/or heterogenous interfaces of amide nuclei at protein-protein interfaces. The intermolecular interfaces were concentrated toward the N-terminus of amelogenin (L3-G8, V19-G38, L46-Q49, and Q57-L70) at pH 6.6 (oligomers) and at pH 7.2 (nanospheres) including the entire N-terminus up to Q76 and regions distributed through the central hydrophobic region (Q82-Q101, S125-Q139, and F151-Q154). At all pH levels, the C-terminus appeared disordered, highly mobile, and not involved in self-assembly, suggesting nanosphere structures with solvent-exposed C-termini. These findings present unique, residue-specific insights into the intermolecular protein-protein interfaces driving amelogenin quaternary structure formation and suggest that nanospheres in solution predominantly contain disordered, solvent-exposed C-termini.

Characterizing the chemical state and physical disposition of uranium that has persisted over geologic time scales is key for modeling the long-term geologic sequestration of nuclear waste, accurate uranium-lead dating, and the use of uranium isotopes as paleo redox proxies. X-ray absorption spectroscopy coupled with molecular dynamics modeling demonstrated that pentavalent uranium is incorporated in the structure of 1.6 billion year old hematite (α-Fe 2 O 3 ), attesting to the robustness of Fe oxides as waste forms and revealing the reason for the great success in using hematite for petrogenic dating. The extreme antiquity of this specimen suggests that the pentavalent state of uranium, considered a transient, is stable when incorporated into hematite, a ubiquitous phase that spans the crustal continuum. Thus, it would appear overly simplistic to assume that only the tetravalent and hexavalent states are relevant when interpreting the uranium isotopic record from ancient crust and contained ore systems.

Progress in understanding crystallization pathways depends on the ability to unravel relationships between intermediates and final crystalline products at the nanoscale, which is a particular challenge at elevated pressure and temperature. Here we exploit a high-pressure atomic force microscope to directly visualize brucite carbonation in water-bearing supercritical carbon dioxide (scCO 2 ) at 90 bar and 50 °C. On introduction of water-saturated scCO 2 , in situ visualization revealed initial dissolution followed by nanoparticle nucleation consistent with amorphous magnesium carbonate (AMC) on the surface. This is followed by growth of nesquehonite (MgCO 3 ·3H 2 O) crystallites. In situ imaging provided direct evidence that the AMC intermediate acts as a seed for crystallization of nesquehonite. In situ infrared and thermogravimetric-mass spectrometry indicate that the stoichiometry of AMC is MgCO 3 ·xH 2 O (x = 0.5-1.0), while its structure is indicated to be hydromagnesite-like according to density functional theory and X-ray pair distribution function analysis. Our findings thus provide insight for understanding the stability, lifetime and role of amorphous intermediates in natural and synthetic systems., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Gibbsite [α-Al(OH) 3 ] is the solubility limiting phase for aluminum across a wide pH range, and it is a common mineral phase with many industrial applications. The growth mechanism of this layered-structure material, however, remains incompletely understood. Synthesis of gibbsite at low to circumneutral pH yields nanoplates with substantial interlayer disorder. Here we examine defects in this material in detail, and the effects of recrystallization in highly alkaline sodium hydroxide solution at 80 °C. We employed a multimodal approach, including scanning electron microscopy, magic-angle spinning nuclear magnetic resonance (MAS-NMR), Raman and infrared spectroscopies, X-ray diffraction (XRD), and X-ray total scattering pair distribution function (XPDF) analysis to characterize the ageing of the nanoplates over several days. XRD and XPDF indicate that gibbsite nanoplates precipitated at circumneutral pH contain dense, truncated sheets imparting a local difference in interlayer distance. These interlayer defects appear well described by flat Al 13 aluminum hydroxide nanoclusters nearly isostructural with gibbsite sheets present under synthesis conditions and trapped as interlayer inclusions during growth. Ageing at elevated temperature in alkaline solutions gradually improves crystallinity, showing a gradual increase in H-bonding between interlayer OH groups. Between 7 to 8 vol% of the initial gibbsite nanoparticles exhibit this defect, with the majority of differences disappearing after 2-4 hours of recrystallization in alkaline solution. The results not only identify the source of disorder in gibbsite formed under acidic/neutral conditions but also point to a possible cluster-mediated growth mechanism evident through inclusion of relict oligomers with gibbsite-like topology trapped in the interlayer spaces.

Injecting supercritical CO 2 (scCO 2 ) into basalt formations for long-term storage is a promising strategy for mitigating CO 2 emissions. Mineral carbonation can result in permanent entrapment of CO 2 ; however, carbonation kinetics in thin H 2 O films in humidified scCO 2 is not well understood. We investigated forsterite (Mg 2 SiO 4 ) carbonation to magnesite (MgCO 3 ) via amorphous magnesium carbonate (AMC; MgCO 3 · x H 2 O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H 2 O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO 2 with a H 2 O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO 2 as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by 2 to 3 orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H 2 O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO 2 fate and transport.

Hematite (α-Fe 2 O 3 ) exerts a strong control over the transport of minor but critical metals in the environment and is used in multiple industrial applications; the photocatalysis community has explored the properties of hematite nanoparticles over a wide range of transition metal dopants. Nonetheless, simplistic assumptions are used to rationalize the local coordination environment of impurities in hematite. Here, we use ab initio molecular dynamics (AIMD)-guided structural analysis to model the extended X-ray absorption fine structure (EXAFS) of Cu 2+ - and Zn 2+ -doped hematite nanoparticles. Specific defect-impurity associations were identified, and the local coordination environments of Cu and Zn both displayed considerable configurational disorder that, in aggregate, approached Jahn-Teller-like distortion for Cu but, in contrast, maintained hematite-like symmetry for Zn. This study highlights the role of defects in accommodating impurities in a nominally low-entropy phase and the limits to traditional shell-by-shell fitting of EXAFS for dopants/impurities in unprecedented bonding environments.

Ettringite is a naturally occurring mineral found in cementitious matrices that is known for its ability to incorporate environmentally mobile oxyanion contaminants. To better assess this immobilization mechanism for contaminants within cementitious waste forms intended for nuclear waste storage, this work explores how mixed oxyanion contaminants compete for ettringite incorporation and influence the evolving mineralogy. Ettringite was precipitated in the presence of TcO 4 - , IO 3 - , and/or CrO 4 2- , known contaminants of concern to nuclear waste treatment, over pre-determined precipitation periods. Solution analyses quantified contaminant removal, and the collected solid was characterized using bulk and microprobe X-ray diffraction coupled with pair distribution function and microprobe X-ray fluorescence analyses. Results suggest that ≥96% IO 3 - is removed from solution, regardless of ettringite precipitation time or the presence of TcO 4 - or CrO 4 2- . However, TcO 4 - removal remained <20%, was not significantly improved with longer ettringite precipitation times, and decreased to zero in the presence of IO 3 - . When IO 3 - is co-mingled with CrO 4 2- , calcite and gypsum are formed as secondary mineral phases, which allows for oxyanion partitioning, e.g., IO 3 - incorporation into ettringite, and CrO 4 2- incorporation into calcite. Results from this work exemplify the importance of competitive immobilization when assessing waste form performance and environmental risk of contaminant release.

Nitrite (NO2-) is a prevalent nitrogen oxyanion in environmental and industrial processes, but its behavior in solution, including ion pair formation, is complex. This solution phase complexity impacts industries such as nuclear waste treatment, where NO2- significantly affects the solubility of other constituents present in sodium hydroxide (NaOH)-rich nuclear waste. This work provides molecular scale information into sodium nitrite (NaNO2) and NaOH ion-pairing processes to provide a physical basis for later development of thermodynamic models. Solubility isotherms of NaNO2 in aqueous mixtures with NaOH and total alkalinity were also measured. Spectroscopic characterization of these solutions utilized high-field nuclear magnetic resonance spectroscopy (NMR) and Raman spectroscopy, with additional solution structure detailed by X-ray total scattering pairwise distribution function analysis (X-ray PDF). Despite the NO2- deformation Raman band's insensitivity to added NaOH in saturated NaNO2 solutions, 23Na and 15N NMR studies indicated the Na+ and NO2- chemical environments change likely due to ion pairing. The ion pairing correlates with a decrease in diffusion coefficient of solution species as measured by pulsed field gradient 23Na and 1H NMR. Two-dimensional correlation analyses of the 2800-4000 cm-1 Raman region and X-ray PDF indicated that saturated NaNO2 and NaOH mixtures disrupt the hydrogen network of water into a new structure where the length of the OO correlations is contracted relative to the typical H2O structure. Beyond describing the solubility of NaNO2 in a multicomponent electrolyte mixture, these results also indicate that nitrite exhibits greater ion pairing in mixtures of concentrated NaNO2 and NaOH than in comparable solutions with only NaNO2.

The bulk behavior of materials is often controlled by minor impurities that create nonperiodic localized defect structures due to ionic size, symmetry, and charge balance mismatches. Here, we used transmission electron microscopy (TEM) of atom-resolved dynamics to directly map the topology of Fe vacancy clusters surrounding structurally incorporated U 6+ in nanohematite (α-Fe 2 O 3 ). Ab initio molecular dynamic simulations provided additional independent constraints on coupled U, Fe, and vacancy mobility in the solid. A clearer understanding of how such an apparently incompatible element can be accommodated by hematite emerged. The results were readily interpretable without the need for sophisticated data reconstruction methods, model structures, or ultrathin samples, and with the proliferation of aberration-corrected TEM facilities, the approach is accessible. Given sufficient z -contrast, the ability to observe impurity-vacancy structures by means of atom hopping can be used to directly probe the association of impurities and such defects in other materials, with promising applications across a broad range of disciplines.

Technetium-99 immobilization in low-temperature nuclear waste forms often relies on additives that reduce environmentally mobile pertechnetate (TcO 4 - ) to insoluble Tc(IV) species. However, this is a short-lived solution unless reducing conditions are maintained over the hazardous life cycle of radioactive wastes (some ∼10,000 years). Considering recent experimental observations, this work explores how rapid formation of ettringite [Ca 6 Al 2 (SO 4 ) 3 (OH) 12 ·26(H 2 O)], a common mineral formed in cementitious waste forms, may be used to directly immobilize TcO 4 - . Results from ab initio molecular dynamics (AIMD) simulations and solid-phase characterization techniques, including synchrotron X-ray absorption, fluorescence, and diffraction methods, support successful incorporation of TcO 4 - into the ettringite crystal structure via sulfate substitution when synthesized by aqueous precipitation methods. One sulfate and one water are replaced with one TcO 4 - and one OH - during substitution, where Ca 2+ -coordinated water near the substitution site is deprotonated to form OH - for charge compensation upon TcO 4 - substitution. Furthermore, AIMD calculations support favorable TcO 4 - substitution at the SO 4 2- site in ettringite rather than gypsum (CaSO 4 ·2H 2 O, formed as a secondary mineral phase) by at least 0.76 eV at 298 K. These results are the first of their kind to suggest that ettringite may contribute to TcO 4 - immobilization and the overall lifetime performance of cementitious waste forms.

Forsterite carbonated in thin H2O films to magnesite via amorphous magnesium carbonate during reaction with H2O-bearing liquid CO2 at 25 °C. This novel reaction pathway contrasts with previous studies that were carried out at higher H2O activity and temperature, where more highly hydrated nesquehonite was the metastable intermediate.

Zn is an essential micronutrient that is often limited in tropical, lateritic soils in part because it is sequestered in nominally refractory iron oxide phases. Stable phases such as goethite and hematite, however, can undergo reductive recrystallization without a phase change under circumneutral pH conditions and release metal impurities such as Zn into aqueous solutions. Further, the process appears to be driven by Fe vacancies. In this contribution, we used ab initio molecular dynamics informed extended X-ray absorption fine structure spectra to show that Zn incorporated in the structure of hematite is associated with coupled O-Fe and protonated Fe vacancies, providing a potential link between crystal chemistry and the bioavailability of Zn.