Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation.

Classical theories of particle aggregation, such as Derjaguin-Landau-Verwey-Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na+, K+, and NO3-) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (∼3-4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal-basal surfaces is more favorable than between (001) edge-edge surfaces (∼5-6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (∼5-10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. The results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

Assessing an aqueous flow cell designed for in situ crystal growth under X-ray nanotomography and effects of radiolysis products.

Nucleation and growth of minerals has broad implications in the geological, environmental and materials sciences. Recent developments in fast X-ray nanotomography have enabled imaging of crystal growth in solutions in situ with a resolution of tens of nanometres, far surpassing optical microscopy. Here, a low-cost, custom-designed aqueous flow cell dedicated to the study of heterogeneous nucleation and growth of minerals in aqueous environments is shown. To gauge the effects of radiation damage from the imaging process on growth reactions, radiation-induced morphological changes of barite crystals (hundreds of nanometres to ∼1 µm) that were pre-deposited on the wall of the flow cell were investigated. Under flowing solution, minor to major crystal dissolution was observed when the tomography scan frequency was increased from every 30 min to every 5 min (with a 1 min scan duration). The production of reactive radicals from X-ray induced water radiolysis and decrease of pH close to the surface of barite are likely responsible for the observed dissolution. The flow cell shown here can possibly be adopted to study a wide range of other chemical reactions in solutions beyond crystal nucleation and growth where the combination of fast flow and fast scan can be used to mitigate the radiation effects.

Armoring of MgO by a Passivation Layer Impedes Direct Air Capture of CO2.

It has been proposed to use magnesium oxide (MgO) to separate carbon dioxide directly from the atmosphere at the gigaton level. We show experimental results on MgO single crystals reacting with the atmosphere for longer (decades) and shorter (days to months) periods with the goal of gauging reaction rates. Here, we find a substantial slowdown of an initially fast reaction as a result of mineral armoring by reaction products (surface passivation). In short-term experiments, we observe fast hydroxylation, carbonation, and formation of amorphous hydrated magnesium carbonate at early stages, leading to the formation of crystalline hydrated Mg carbonates. The preferential location of Mg carbonates along the atomic steps on the crystal surface of MgO indicates the importance of the reactive site density for carbonation kinetics. The analysis of 27-year-old single-crystal MgO samples demonstrates that the thickness of the reacted layer is limited to ∼1.5 µm on average, which is thinner than expected and indicates surface passivation. Thus, if MgO is to be employed for direct air capture of CO2, surface passivation must be circumvented.

The evolution of solvation symmetry and composition in Zn halide aqueous solutions from dilute to extreme concentrations.

The emergence of cation-anion species, or contact ion pairs, is fundamental to understanding the physical properties of aqueous solutions when moving from the ideal, low-concentration limit to the manifestly non-ideal limits of very high solute concentration or constituent ion activity. We focus here on Zn halide solutions both as a model system and also as an exemplar of the applications spanning from (i) electrical energy storage via the paradigm of water in salt electrolyte (WiSE) to (ii) the physical chemistry of brines in geochemistry to (iii) the long-standing problem of nucleation. Using a combination of experimental and theoretical approaches we quantify the halide coordination number and changing coordination geometry without embedded use of theoretical equilibrium constants. These results and the associated methods, notably including the use of valence-to-core X-ray emission spectroscopy, provide new insights into the Zn halide system and new research directions in the physical chemistry of concentrated electrolytes.

Improving Rare-Earth Mineral Separation with Insights from Molecular Recognition: Functionalized Hydroxamic Acid Adsorption onto Bastnäsite and Calcite.

Enhancing the separation of rare-earth elements (REEs) from gangue materials in mined ores requires an understanding of the fundamental interactions driving the adsorption of collector ligands onto mineral interfaces. In this work, we examine five functionalized hydroxamic acid ligands as potential collectors for the REE-containing bastnäsite mineral in froth flotation using density functional theory calculations and a suite of surface-sensitive analytical spectroscopies. These include vibrational sum frequency generation, attenuated total reflectance Fourier transform infrared, Raman, and X-ray photoelectron spectroscopies. Differences in the chemical makeup of these ligands on well-defined bastnäsite and calcite surfaces allow for a systematic relationship connecting the structure to adsorption activity to be framed in the context of interfacial molecular recognition. We show how the intramolecular hydrogen bonding of adsorbed ligands requires the inclusion of explicit water solvent molecules to correctly map energetic and structural trends measured by experiments. We anticipate that the results and insights from this work will motivate and inform the design of improved flotation collectors for REE ores.

Selenite and Selenate Sequestration during Coprecipitation with Barite: Insights from Mineralization Processes of Adsorption, Nucleation, and Growth.

Coprecipitation of selenium oxyanions with barite is a facile way to sequester Se in the environments. However, the chemical composition of Se-barite coprecipitates usually deviates from that predicted from thermodynamic calculations. This discrepancy was resolved by considering variations in nucleation and growth rates controlled by ion-mineral interactions, solubility, and interfacial energy. For homogeneous precipitation, ∼10% of sulfate, higher than thermodynamic predictions (<0.3%), was substituted by Se(IV) or Se(VI) oxyanion, which was attributed to adsorption-induced entrapment during crystal growth. For heterogeneous precipitation, thiol- and carboxylic-based organic films, utilized as model interfaces to mimic the natural organic-abundant environments, further enhanced the sequestration of Se(VI) oxyanions (up to 41-92%) with barite. Such enhancement was kinetically driven by increased nucleation rates of selenate-rich barite having a lower interfacial energy than pure barite. In contrast, only small amounts of Se(IV) oxyanions (∼1%) were detected in heterogeneous coprecipitates mainly due to a lower saturation index of BaSeO3 and deprotonation degree of Se(IV) oxyanion at pH 5.6. These roles of nanoscale mineralization mechanisms observed during composition selection of Se-barite could mark important steps toward the remediation of contaminants through coprecipitation.

Synergistic Enhancement of Lead and Selenate Uptake at the Barite (001)-Water Interface.

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. We also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion-ion interactions at mineral-water interfaces that control the fate and transport of contaminants in the environment.

Radiolysis and Radiation-Driven Dynamics of Boehmite Dissolution Observed by In Situ Liquid-Phase TEM.

Over the last several decades, there have been several studies examining the radiation stability of boehmite and other aluminum oxyhydroxides, yet less is known about the impact of radiation on boehmite dissolution. Here, we investigate radiation effects on the dissolution behavior of boehmite by employing liquid-phase transmission electron microscopy (LPTEM) and varying the electron flux on the samples consisting of either single nanoplatelets or aggregated stacks. We show that boehmite nanoplatelets projected along the [010] direction exhibit uniform dissolution with a strong dependence on the electron dose rate. For nanoplatelets that have undergone oriented aggregation, we show that the dissolution occurs preferentially at the particles at the ends of the stacks that are more accessible to bulk solution than at the others inside the aggregate. In addition, at higher dose rates, electrostatic repulsion and knock-on damage from the electron beam causes delamination of the stacks and dissolution at the interfaces between particles in the aggregate, indicating that there is a threshold dose rate for electron-beam enhancement of dissolution of boehmite aggregates.

Nanopore facilitated monohydrocalcitic amorphous calcium carbonate precipitation.

Predicting the precipitation of solids is important in both natural systems and subsurface energy applications. The factors controlling reaction mechanisms, phase selection and conversion between phases are particularly important. In this contribution the precipitation and growth of an amorphous calcium carbonate species from flowing aqueous solution in a nanoporous controlled pore glass is followed in situ with differential X-ray pair distribution function analysis. It is discovered that the local atomic structure of this phase indicates monohydrocalcite-like pair-pair correlations, yet is functionally amorphous because it lacks long-range structure. The unexpected occurrence of synthetic proto-monohydrocalcite amorphous calcium carbonate, precipitated from a solution undersaturated with respect to published solubilities, suggests that nanopore confinement facilitates formation of an amorphous phase at the expense of more favorable crystalline ones. This result illustrates that confinement and interface effects are physical factors exerting control on mineral nucleation behavior in natural and geological systems.

Concentration dependent interfacial chemistry of the NaOH(aq): gibbsite interface.

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, γ-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)-particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the γ-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5-10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH- accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH- ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion-surface interactions. The molecular-scale changes to surface speciation and competition between ion-surface vs. ion-ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

The "good," the "bad," and the "hidden" in neutron scattering and molecular dynamics of ionic aqueous solutions.

We characterize a concentrated 7.3 m CaCl2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca2+-Cl- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca2+-Ow and Ca2+-Hw contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

Organic-mineral interfacial chemistry drives heterogeneous nucleation of Sr-rich (Ba x , Sr1-x )SO4 from undersaturated solution.

Sr-bearing marine barite [(Ba x , Sr1-x )SO4] cycling has been widely used to reconstruct geochemical evolutions of paleoenvironments. However, an understanding of barite precipitation in the ocean, which is globally undersaturated with respect to barite, is missing. Moreover, the reason for the occurrence of higher Sr content in marine barites than expected for classical crystal growth processes remains unknown. Field data analyses suggested that organic molecules may regulate the formation and composition of marine barites; however, the specific organic-mineral interactions are unclear. Using in situ grazing incidence small-angle X-ray scattering (GISAXS), size and total volume evolutions of barite precipitates on organic films were characterized. The results show that barite forms on organic films from undersaturated solutions. Moreover, from a single supersaturated solution with respect to barite, Sr-rich barite nanoparticles formed on organics, while micrometer-size Sr-poor barites formed in bulk solutions. Ion adsorption experiments showed that organic films can enrich cation concentrations in the adjacent solution, thus increasing the local supersaturation and promoting barite nucleation on organic films, even when the bulk solution was undersaturated. The Sr enrichment in barites formed on organic films was found to be controlled by solid-solution nucleation rates; instead, the Sr-poor barite formation in bulk solution was found to be controlled by solid-solution growth rates. This study provides a mechanistic explanation for Sr-rich marine barite formation and offers insights for understanding and controlling the compositions of solid solutions by separately tuning their nucleation and growth rates via the unique chemistry of solution-organic interfaces.

Hydration structure and water exchange kinetics at xenotime-water interfaces: implications for rare earth minerals separation.

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral-water interfaces. Here, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5-180 picoseconds), dictated by ion-water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

Resolving local configurational contributions to X-ray and neutron radial distribution functions within solutions of concentrated electrolytes - a case study of concentrated NaOH.

Extreme conditions of complex materials often lead to a manifold of local environments that challenge characterization and require new advances at the intersection of modern experimental and theoretical techniques. In this contribution, highly caustic and viscous aqueous NaOD solutions were characterized with a combination of X-ray and neutron radial distribution function (RDF) analyses, molecular dynamics simulations and sub-ensemble analysis. While this system has been the topic of some study, the current work expands upon the state of knowledge regarding the extent to which water is perturbed within this chemically extreme solution. Further, we introduce analyses that goes beyond merely identifying the different local environments (ion solvation and coordination environments) that are present, but toward understanding their relative contributions to the ensemble solution RDF. This integrated approach yields unique insight into the experimental sensitivity of RDFs to changes in local geometries, the composition of solvation environments about ions, and the challenge of experimentally differentiating the ensemble of all superimposed local environments-a feature of increasing importance within the extreme condition of high ionic strength.

In Situ 27Al NMR Spectroscopy of Aluminate in Sodium Hydroxide Solutions above and below Saturation with Respect to Gibbsite.

Aluminum hydroxide (Al(OH)3, gibbsite) dissolution and precipitation processes in alkaline environments play a commanding role in aluminum refining and nuclear waste processing, yet mechanistic aspects underlying sluggish kinetics during crystallization have remained obscured due to a lack of in situ probes capable of isolating incipient ion pairs. At a molecular level Al is cycling between tetrahedral ( T d) coordination in solution to octahedral ( O h) in the solid. We explored dissolution of Al(OH)3 that was used to produce variably saturated aluminate (Al(OH)4-)-containing solutions under alkaline conditions (pH >13) with in situ 27Al magic angle spinning (MAS)-nuclear magnetic resonance (NMR) spectroscopy, and interrogated the results with ab initio molecular dynamics (AIMD) simulations complemented with chemical shift calculations. The collective results highlight the overall stability of the solvation structure for T d Al in the Al(OH)4- oxyanion as a function of both temperature and Al concentration. The observed chemical shift did not change significantly even when the Al concentration in solution became supersaturated upon cooling and limited precipitation of the octahedral Al(OH)3 phase occurred. However, subtle changes in Al(OH)4- speciation correlated with the dissolution/precipitation reaction were found. AIMD-informed chemical shift calculations indicate that measurable perturbations should begin when the Al(OH)4-···Na+ distance is less than 6 Å, increasing dramatically at shorter distances, coinciding with appreciable changes to the electrostatic interaction and reorganization of the Al(OH)4- solvation shell. The integrated findings thus suggest that, under conditions incipient to and concurrent with gibbsite crystallization, nominally expected contact ion pairs are insignificant and instead medium-range (4-6 Å) solvent-separated Al(OH)4-···Na+ pairs predominate. Moreover, the fact that these medium-range interactions bear directly on resulting gibbsite characteristics was demonstrated by detailed microscopic and X-ray diffraction analysis and by progressive changes in the fwhm of the O h resonance, as measured by in situ NMR. Sluggish gibbsite crystallization may arise from the activation energy associated with disrupting this robust medium-range ion pair interaction.

Wellbore Cement Porosity Evolution in Response to Mineral Alteration during CO2 Flooding.

Mineral reactions during CO2 sequestration will change the pore-size distribution and pore surface characteristics, complicating permeability and storage security predictions. In this paper, we report a small/wide angle scattering study of wellbore cement that has been exposed to carbon dioxide for three decades. We have constructed detailed contour maps that describe local porosity distributions and the mineralogy of the sample and relate these quantities to the carbon dioxide reaction front on the cement. We find that the initial bimodal distribution of pores in the cement, 1-2 and 10-20 nm, is affected differently during the course of carbonation reactions. Initial dissolution of cement phases occurs in the 10-20 nm pores and leads to the development of new pore spaces that are eventually sealed by CaCO3 precipitation, leading to a loss of gel and capillary nanopores, smoother pore surfaces, and reduced porosity. This suggests that during extensive carbonation of wellbore cement, the cement becomes less permeable because of carbonate mineral precipitation within the pore space. Additionally, the loss of gel and capillary nanoporosities will reduce the reactivity of cement with CO2 due to reactive surface area loss. This work demonstrates the importance of understanding not only changes in total porosity but also how the distribution of porosity evolves with reaction that affects permeability.

A comparative study of surface energies and water adsorption on Ce-bastnäsite, La-bastnäsite, and calcite via density functional theory and water adsorption calorimetry.

Bastnäsite, a fluoro-carbonate mineral, is the single largest mineral source of light rare earth elements (REE), La, Ce and Nd. Enhancing the efficiency of separation of the mineral from gangue through froth flotation is the first step towards meeting an ever increasing demand for REE. To design and evaluate collector molecules that selectively bind to bastnäsite, a fundamental understanding of the structure and surface properties of bastnäsite is essential. In our earlier work (J. Phys. Chem. C, 2016, 120, 16767), we carried out an extensive study of the structure, surface stability and water adsorption energies of La-bastnäsite. In this work, we make a comparative study of the surface properties of Ce-bastnäsite, La-bastnäsite, and calcite using a combination of density functional theory (DFT) and water adsorption calorimetry. Spin polarized DFT+U calculations show that the exchange interaction between the electrons in Ce 4f orbitals is negligible and that these orbitals do not participate in bonding with the oxygen atom of the adsorbed water molecule. In agreement with calorimetry, DFT calculations predict larger surface energies and stronger water adsorption energies on Ce-bastnäsite than on La-bastnäsite. The order of stabilities for stoichiometric surfaces is as follows: [101[combining macron]0] > [101[combining macron]1] > [101[combining macron]2] > [0001] > [112[combining macron]2] > [101[combining macron]4] and the most favorable adsorption sites for water molecules are the same as for La-bastnäsite. In agreement with water adsorption calorimetry, at low coverage water molecules are strongly stabilized via coordination to the surface Ce3+ ions, whereas at higher coverage they are adsorbed less strongly via hydrogen bonding interaction with the surface anions. Due to similar water adsorption energies on bastnäsite [101[combining macron]1] and calcite [101[combining macron]4] surfaces, the design of collector molecules that selectively bind to bastnäsite over calcite must exploit the structural differences in the predominantly exposed facets of these minerals.

Uncovering the Atomistic Mechanism for Calcite Step Growth.

Determining a complete atomic-level picture of how minerals grow from aqueous solution remains a challenge as macroscopic rates can be a convolution of many reactions. For the case of calcite (CaCO3 ) in water, computer simulations have been used to map the complex thermodynamic landscape leading to growth of the two distinct steps, acute and obtuse, on the basal surface. The carbonate ion is found to preferentially adsorb at the upper edge of acute steps while Ca2+ only adsorbs after CO32- . In contrast to the conventional picture, ion pairs prefer to bind at the upper edge of the step with only one ion, at most, coordinated to the step and lower terrace. Migration of the first carbonate ion to a growth site is found to be rate-limiting for kink nucleation, with this process having a lower activation energy on the obtuse step.

Heterogeneous Nucleation and Growth of Barium Sulfate at Organic-Water Interfaces: Interplay between Surface Hydrophobicity and Ba(2+) Adsorption.

Barium sulfate (BaSO4) is a common scale-forming mineral in natural and engineered systems, yet the rates and mechanisms of heterogeneous BaSO4 nucleation are not understood. To address these, we created idealized interfaces on which to study heterogeneous nucleation rates and mechanisms, which also are good models for organic-water interfaces: self-assembled thin films terminated with different functional groups (i.e., -COOH, -SH, or mixed -SH & COOH) coated on glass slides. BaSO4 precipitation on coatings from Barite-supersaturated solutions (saturation index, SI, = 1.1) was investigated using grazing-incidence small-angle X-ray scattering. After reaction for 1 h, a little amount of BaSO4 formed on hydrophilic bare and -COOH coated glasses. Meanwhile, BaSO4 nucleation was significantly promoted on hydrophobic -SH and mixed -SH & COOH coatings. This is because substrate hydrophobicity likely affected the interfacial energy and hence thermodynamic favorability of heterogeneous nucleation. The heterogeneous BaSO4 nucleation and growth kinetics were found to be affected by the amount of Ba(2+) adsorption onto the substrate and incipient BaSO4 nuclei. The importance of Ba(2+) adsorption was further corroborated by the finding that precipitation rate increased under [Ba(2+)]/[SO4(2-)] concentration ratios >1. These observations suggest that thermodynamic favorability for nucleation is governed by substrate-water interfacial energy, while given favorable thermodynamics, the rate is governed by ion attachment to substrates and incipient nuclei.

Internal Domains of Natural Porous Media Revealed: Critical Locations for Transport, Storage, and Chemical Reaction.

Internal pore domains exist within rocks, lithic fragments, subsurface sediments, and soil aggregates. These domains, termed internal domains in porous media (IDPM), represent a subset of a material's porosity, contain a significant fraction of their porosity as nanopores, dominate the reactive surface area of diverse media types, and are important locations for chemical reactivity and fluid storage. IDPM are key features controlling hydrocarbon release from shales in hydraulic fracture systems, organic matter decomposition in soil, weathering and soil formation, and contaminant behavior in the vadose zone and groundwater. Traditionally difficult to interrogate, advances in instrumentation and imaging methods are providing new insights on the physical structures and chemical attributes of IDPM, and their contributions to system behaviors. Here we discuss analytical methods to characterize IDPM, evaluate information on their size distributions, connectivity, and extended structures; determine whether they exhibit unique chemical reactivity; and assess the potential for their inclusion in reactive transport models. Ongoing developments in measurement technologies and sensitivity, and computer-assisted interpretation will improve understanding of these critical features in the future. Impactful research opportunities exist to advance understanding of IDPM, and to incorporate their effects in reactive transport models for improved environmental simulation and prediction.