Search NASA⌕ Search

Engineering topics

Stack, Andrew G.

Publications and source records attributed to Stack, Andrew G..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Iron Impurity Impairs the CO 2 Capture Performance of MgO: Insights from Microscopy and Machine Learning Molecular Dynamics

Magnesium oxide (MgO) is a promising sorbent for direct air capture (DAC) of carbon dioxide. Iron (Fe) is a common impurity in naturally occurring MgO and minerals used to produce MgO, yet a molecular-scale understanding of Fe-doping effects on carbonation is lacking. Here, in this study, we observed reduced carbonation performance in Fe-doped MgO experimentally. The energetics of adsorbing a (bi)carbonate ion on pristine and Fe-doped MgO(001) surfaces were further investigated using ab initio and machine learning potential molecular dynamics coupled with metadynamics simulations. Both pristine and Fe-doped surfaces exhibited a basic (OH – ) hydration layer, where the (bi)carbonate ion adsorption is thermodynamically favorable. However, the dissolution of surface Fe had smaller energy barriers and was more favorable than Mg. Leached Fe likely neutralized the near-surface basicity, yielding reduced reactivity on Fe-doped MgO. Our observations offer critical insights for material selection and emphasize the importance of evaluating the geologic origin of earth materials used for DAC.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Variation in Cation Adsorption Mechanism Controlled by Chemical and Structural Heterogeneities at the Quartz (101)–Water Interface

Mineral–water interfacial reactions are central to chemical processes that control the fate of nutrients and contaminants in natural environments. Mineral surfaces commonly have complex structures and compositions whose impact on interfacial reactivity is poorly understood. Here, in this work, we investigated the effects of surface heterogeneities on Rb + sorption at the quartz (101)–10 mM RbCl solution interface at pH 9.8 using in situ high-resolution X-ray reflectivity. Two surface locales (i.e., Spots A and B) having distinct interfacial structures were chosen: Spot A was characterized by its low defect density (≤20% topmost Si vacancies) and Rb + adsorption occurred predominantly as an inner-sphere complex. In comparison, Spot B had a higher defect density (~50% vacancies) and was covered with poorly crystalline SiO 2 . A substantially larger Rb + uptake (i.e., 7-times higher coverage) was observed on this defective surface where Rb + incorporated in the vacancy sites (confirmed by density functional tight binding-based molecular dynamics simulations) or adsorbed directly on the disordered film. These results provide a direct quantification of how surface heterogeneity influences the geochemical behavior of mineral–water interfaces, in particular highlighting the important role of chemical and structural defects on the sorbate speciation and coverage at silicate mineral surfaces.

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

Ion Correlations Decrease Particle Aggregation Rate by Increasing Hydration Forces at Interfaces

The connection between solution structure, particle forces, and emergent phenomena at solid–liquid interfaces remains ambiguous. In this case study on boehmite aggregation, we established a connection between interfacial solution structure, emerging hydration forces between two approaching particles, and the resulting structure and kinetics of particle aggregation. In contrast to expectations from continuum-based theories, we observed a nonmonotonic dependence of the aggregation rate on the concentration of sodium chloride, nitrate, or nitrite, decreasing by 15-fold in 4 molal compared to 1 molal solutions. These results are accompanied by an increase in repulsive hydration forces and interfacial oscillatory features from 0.27–0.31 nm in 0.01 molal to 0.38–0.52 nm in 2 molal. Moreover, molecular dynamics (MD) simulations indicated that these changes correspond to enhanced ion correlations near the interface and produced loosely bound aggregates that retain electrolyte between the particles. In conclusion, we anticipate that these results will enable the prediction of particle aggregation, attachment, and assembly, with broad relevance to interfacial phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of magnetic and electric fields on the free energy to form a calcium carbonate ion-pair

Electromagnetic fields are used in water treatment and desalination to regulate scale formation and extend the lifetime of membranes. External electric and magnetic fields can promote or suppress mineral nucleation and growth. However, the molecular-scale mechanisms of such processes remain unknown. Computing the free energies needed to form ion pairs under external fields provides important insights into understanding the elemental steps during the initial formation of mineral scales. In this paper, we used molecular dynamics combined with metadynamics simulations to investigate the free energies of forming the [Ca–CO 3 ] 0 ion pair, a fundamental building block of carbonate scales, under a range of magnetic (up to 10 T) and electric (up to 10 V m −1 ) fields in water. The presence of constant magnetic or electric fields favored the ion pairing reaction and lowered the free energies by up to 3% to 6%. The internal energy and entropic components of the free energy showed significant changes and exhibited non-linear behavior with increasing field strength. The [Ca–CO 3 ] 0 ion pairing is an entropy-driven process in the absence of an external field, but the mechanism shifts to an internal energy-driven process under selected external fields, suggesting possible changes in the nucleation pathways.

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

Real-Time Atomic-Scale Structural Analysis Resolves the Amorphous to Crystalline CaCO 3 Mechanism Controversy

Amorphous calcium carbonate (ACC) occurs as a precursor to geological and biogenic calcium carbonate (CaCO 3 ), yet its transformation pathways and reaction mechanisms remain inconsistent and controversial. In this study, we investigated the transformation of ACC to calcite under both solution and dry conditions, in the presence and absence of impurity ions, utilizing operando time-resolved synchrotron X-ray diffraction (TRXRD) and reactive transport modeling. Results demonstrate that TRXRD techniques allow us to differentiate dissolution-reprecipitation versus solid-state transformation mechanisms for amorphous to crystalline phase transitions. Specifically, we observe that in environments with abundant water, ACC transforms to calcite through a dissolution-reprecipitation mechanism. This features an activation energy of 63 ± 2 kJ/mol and unit cell volume contraction during calcite crystal growth. Conversely, under water-limited conditions, ACC to calcite transformation proceeds via a solid-state transformation mechanism, with an activation energy of 210 ± 2 kJ/mol, three times greater than the dissolution-reprecipitation route, and a unit cell expansion during crystalline calcite growth. Further, to illustrate the magnitude of these effects, the rates of calcite growth were similar during dissolution-reprecipitation at 3 °C [0.00207(35) s –1 ] and solid-state transformation at 280 °C [0.00134(11) s –1 ]. Moreover, the incorporation of an impurity, strontium, significantly retards the rate of calcite growth while expanding its unit cell but whose incorporation is history dependent. Reactive transport modeling of the dissolution–precipitation kinetics suggests that ACC must be dissolving as compact aggregates. These various transformation mechanisms drive diverse geological and biological carbonate formations, impacting their use as paleoenvironmental markers and functional materials synthesis.

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

Effect of Nucleation Heterogeneity on Mineral Precipitation in Confined Environments

Abstract The formation of new mineral phases in confined environments, especially in porous media, is crucial for various geological processes like mineralization and diagenesis. The nucleation and precipitation of minerals are initiated at the microscale through fluid‐rock interaction, where dissolution of primary phases leads to supersaturated conditions and nucleation and growth of secondary ones. Previous research has focused primarily on either precipitation or nucleation, without fully exploring their combined impact. Our study introduces a computational framework that integrates classical nucleation theory with the micro‐continuum method. We validated our model by comparing with experiments, and discovered that different surface nucleation rate changes the mode of precipitation from a preferential to uniform precipitate textures. Furthermore, our study uncovered that the conventional deterministic precipitation method tends to underestimate the permeability of the porous matrix. In contrast, the new framework significantly improves model accuracy by incorporating preferential precipitation and heterogeneous nucleation.

58 GEOSCIENCES↗
DOE OSTI2024

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 NO 3 –) 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. Further, 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.

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

Variations in proton transfer pathways and energetics on pristine and defect-rich quartz surfaces in water: Insights into the bimodal acidities of quartz

Hypothesis. Understanding the mechanisms of proton transfer on quartz surfaces in water is critical for a range of processes in geochemical, environmental, and materials sciences. The wide range of surface acidities (>9 pKa units) found on the ubiquitous mineral quartz is caused by the structural variations of surface silanol groups. Molecular scale simulations provide essential tools for elucidating the origin of site-specific surface acidities. Simulations. Here, we used density-functional tight-binding-based molecular dynamics combined with rare-event metadynamics simulations to probe the mechanisms of deprotonation reactions from ten representative surface silanol groups found on both pristine and defect-rich quartz (1 0 1) surfaces with Si vacancies. Findings. The results show that deprotonation is a highly dynamic process where both the surface hydroxyls and bridging oxygen atoms serve as the proton acceptors, in addition to water. Deprotonation of embedded silanols through intrasurface proton transfer exhibited lower pKa values with less H-bond participation and higher energy barriers, suggesting a new mechanism to explain the bimodal acidity observed on quartz surface. Defect sites, recently shown to comprise a significant portion of the quartz (1 0 1) surface, diversify the coordination and local H-bonding environments of the surface silanols, changing both the deprotonation pathways and energetics, leading to a wider range of pKa values (2.4 to 11.5) than that observed on pristine quartz surface (10.4 and 12.1).

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

Reaction Layer Formation on MgO in the Presence of Humidity

Mineralization by MgO is an attractive potential strategy for direct air capture (DAC) of CO 2 due to its tendency to form carbonate phases upon exposure to water and CO 2 . Hydration of MgO during this process is typically assumed to not be rate limiting, even at ambient temperatures. However, surface passivation by hydrated phases likely reduces the CO 2 capture capacity. Here, we examine the initial hydration reactions that occur on MgO(100) surfaces to determine whether they could potentially impact CO 2 uptake. We first used atomic force microscopy (AFM) to explore changes in reaction layers in water (pH = 6 and 12) and MgO-saturated solution (pH = 11) and found the reaction layers on MgO are heterogeneous and nonuniform. To determine how relative humidity (R.H.) affects reactivity, we reacted samples at room temperature in nominally dry N 2 (~11–12% R.H.) for up to 12 h, in humid (>95% R.H.) N 2 for 5, 10, and 15 min, and in air at 33 and 75% R.H. for 8 days. X-ray reflectivity and electron microscopy analysis of the samples reveal that hydrated phases form rapidly upon exposure to humid air, but the growth of the hydrated reaction layer slows after its initial formation. Reaction layer thickness is strongly correlated with R.H., with denser reaction layers forming in 75% R.H. compared with 33% R.H. or nominally dry N 2 . The reaction layers are likely amorphous or poorly crystalline based on grazing incidence X-ray diffraction measurements. After exposure to 75% R.H. in air for 8 days, the reaction layer increases in density as compared to the sample reacted in humid N 2 for 5–15 min. Further, this may represent an initial step toward the crystallization of the reaction layer. Overall, high R.H. favors the formation of a hydrated, disordered layer on MgO. Based on our results, DAC in a location with a higher R.H. will be favorable, but growth may slow significantly from initial rates even on short timescales, presumably due to surface passivation.

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

Multilevel atomic structural model for interstratified opal materials

The structure of opal has long fascinated scientists. It occurs in a number of structural states, ranging from amorphous to exhibiting features of stacking disorder. Opal-CT, where C and T signify cristobalite- and tridymite-like interstratification, represents an important link in the length scales between amorphous and crystalline states. However, details about local atomic (dis)order and arrangements extending to long-range stacking faults in opal polymorphs remain incompletely understood. Here, a multilevel modeling approach is reported that considers stacking states in correlation with the abundance of C and T segments as a high-level structural parameter (i.e. not each atom). Optimization accounting for inter-tetrahedral bond lengths and angles and the regularity of the silicate tetrahedra is included as lower levels of structural parameters. Together, a set of parameters with both coarse-grained and atomistic features for different levels of structural details is refined. Structural disorder at the ~10–100 Å distance scale is evaluated using experimental pair distribution function and diffraction datasets, comparing peak intensities, widths and asymmetry. Here this work presents a complete multilevel structural description of natural opal-CT and explains many of the unusual features observed in X-ray powder diffraction patterns. This modeling approach can be adopted generally for analyzing layered materials and their assembly into 3D structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure of the Bastnäsite (001) Surface by Crystal Truncation Rod X-ray Diffraction and Ab Initio Molecular Dynamics: Implications for Separations of a Rare Earth Ore Mineral

Bastnäsite ((Ce,La)FCO 3 ) is the primary mineral source of light rare earth elements, but its surface structure is not well understood. This presents a major challenge in improving beneficiation strategies. In this work, a synergistic combination of X-ray scattering and ab initio molecular dynamics (AIMD) was used to gain atomistic insight into the interfacial structure of bastnäsite. Surface X-ray scattering was used to measure crystal truncation rods (CTRs) of the bastnäsite (001) surface, a significant crystal face with a previously unknown termination. The best-fit atomic-scale model of the CTR data features a carbonate layer at the surface, which is stabilized by the relaxation of carbonate groups from their bulk structural positions. AIMD simulations predict similar surface relaxations, which are shown to be influenced by the protonation of oxygen atoms at the surface. Evidence of ordered water at the interface is also observed in the best-fit model and AIMD simulations. Further, the presence of a carbonate layer at this dominant crystal surface is significant for improving separation technologies because most commonly used ligands utilize anionic functional groups to chelate metal cations at particle surfaces. Without modification, anionic ligands are expected to have poor affinity for the carbonate-terminated (001) surface.

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

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Here, our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

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

Sorption of Arsenate, Selenate, and Molybdate on the Barite (001) Surface

Interactions of ions with ionically bonded minerals such as barite (BaSO 4 ) influence the fate and transport of the ions, while the factors that control the sorption of toxic oxyanions on barite remain elusive. In this study, the sorption of arsenate, selenate, and molybdate on the barite (001) surface was examined at pH ~5 using in situ crystal truncation rod analysis, resonant anomalous X-ray reflectivity, and atomic force microscopy. The results show that arsenate and selenate mainly incorporate into the top monolayer of barite, while molybdate primarily adsorbs above the surface. The sorption coverage of arsenate is greater (by ~100%) than that of selenate but similar to that of molybdate. The different incorporation coverages between arsenate and selenate can be explained by their different protonation states at pH 5. The incorporated arsenate may be stabilized by hydrogen bonds between arsenate and oxygen atoms of neighboring sulfate compared to selenate, which exists predominantly in the deprotonated state. The adsorption of molybdate above the surface probably stems from a surface-induced oligomerization, as the anion and the oligomer may be too large for incorporation. Our observation of these different sorption mechanisms demonstrates how the physicochemical properties of the anions control the selective uptake of the toxic metals on the dominant surface of the ionically bonded mineral barite.

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

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.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Grain boundary widening controls siderite (FeCO3) replacement of limestone (CaCO3)

Abstract The microstructure of minerals and rocks can significantly alter reaction rates. This study focuses on identifying transport paths in low porosity rocks based on the hypothesis that grain boundary widening accelerates reactions in which one mineral is replaced by another (replacement reaction). We conducted a time series of replacement experiments of three limestones (CaCO 3 ) of different microstructures and solid impurity contents using FeCl 2 . Reacted solids were analyzed using chemical imaging, small angle X-ray and neutron scattering and Raman spectroscopy. In high porosity limestones replacement is reaction controlled and complete replacement was observed within 2 days. In low porosity limestones that contain 1–2% dolomite impurities and are dominated by grain boundaries, a reaction rim was observed whose width did not change with reaction time. Siderite (FeCO 3 ) nucleation was observed in all parts of the rock cores indicating the percolation of the solution throughout the complete core. Dolomite impurities were identified to act as nucleation sites leading to growth of crystals that exert force on the CaCO 3 grains. Widening of grain boundaries beyond what is expected based on dissolution and thermal grain expansion was observed in the low porosity marble containing dolomite impurities. This leads to a self-perpetuating cycle of grain boundary widening and reaction acceleration instead of reaction front propagation.

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

Solution Structure and Hydration Forces between Mica and Hydrophilic Versus Hydrophobic Surfaces

Solid-liquid interfaces are central to a range of interesting phenomena including catalysis, heterogeneous nucleation, water desalination, and biomolecular assembly. While three-dimensional Fast Force Mapping (3D FFM) has emerged as a technique for resolving interfacial solution structure at the molecular scale, key challenges for data interpretation persist, most notably regarding the influence of the probe on the measured structure. Using the mica-water system as a case study, we investigate the effect of hydrophilic and hydrophobic probes on interfacial solution structure measured by 3D FFM. Data from hydrophilic silicon-based probes are in good agreement with molecular dynamics simulations, wherein the innermost water molecules adsorb preferentially at the surface ditrigonal cavity sites followed by two subsequent hydration layers. In contrast, the hydrophobic carbon-based probes detect vertical oscillatory features, but do not show lateral patterning that matches the underlying mica lattice. At high ionic strength, up to six of these oscillatory features are observed extending 2 nm into the solution phase with an average spacing of 0.29 ± (0.04) nm. Further, we also determine that the repulsive hydration force between mica and the hydrophilic probe depends on the nature and concentration of ions in solution. Specifically, solutions with stronger ion-water and ion-ion interactions increase the interfacial viscosity, which results in a stronger repulsive hydration force as the probe approaches the surface. Based on these observations, we present a scheme for controlling the outcomes of particle attachment and aggregation by varying the solution conditions to tune the hydration force.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗