Engineering Papers⌕ Search

Engineering topics

Sassi, Michel

Publications and source records attributed to Sassi, Michel.

Formation and transformation of iron oxy-hydroxide precursor clusters to ferrihydrite

Iron (Fe) oxy-hydroxide minerals such as ferrihydrite (Fh) are ubiquitous in Earth-surface environments and important in biogeochemical element cycling. Recent research has suggested that their formation is preceded by the precipitation of ultrasmall (~1 nm) Keggin-like Fe oxy-hydroxide clusters. However, relatively little is understood about the structure of the precursor clusters and the impacts of pH and time on their growth and transformation to more stable phases. Here, we used a new method that involves mixed flow reactors (MFR) to synthesize these Fe oxy-hydroxide precursor clusters at pH 1.0, 1.5, 2.5, and 4.5. In situ and ex situ synchrotron scattering measurements and laboratory small-angle X-ray scattering (SAXS) were used to study the structure and size of Fe oxy-hydroxide clusters and their transformation products, respectively. Results show that with increasing pH, the particle size and structural order of samples increase, forming solids that resemble 2-line Fh at pH 4.5. The experimental data were compared with X-ray pair distribution functions (PDF) calculated for a range of Fe(III) oxyhydroxide clusters, including Fe 13 Keggin isomers computed previously using density functional theory (DFT), which yielded at best only partial agreement at short range (<5 Å). Aging of the clusters synthesized at pH 1.5 and 2.5 results in growth and transformation via Ostwald ripening to mixtures of goethite (Gt) and lepidocrocite (Lp). This process was inhibited by immediately reacting the early-formed clusters with phosphate (PO 4 3– ), suggesting that oxyanion surface complexes can stabilize the initial clusters by preventing growth and crystallization to more stable phases.

58 GEOSCIENCES↗

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗

Crystal dissolution by particle detachment

Crystal dissolution, which is a fundamental process in both natural and technological settings, has been predominately viewed as a process of ion-by-ion detachment into a surrounding solvent. Here we report a mechanism of dissolution by particle detachment (DPD) that dominates in mesocrystals formed via crystallization by particle attachment (CPA). Using liquid phase electron microscopy to directly observe dissolution of hematite crystals — both compact rhombohedra and mesocrystals of coaligned nanoparticles — we find that the mesocrystals evolve into branched structures, which disintegrate as individual sub-particles detach. The resulting dissolution rates far exceed those for equivalent masses of compact single crystals. Applying a numerical generalization of the Gibbs-Thomson effect, we show that the physical drivers of DPD are curvature and strain inherently tied to the original CPA process. Based on the generality of the model, we anticipate that DPD is widespread for both natural minerals and synthetic crystals formed via CPA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the Impact of Electric Fields on the Local Structure and Migration of Potassium Ions at the Orthoclase (001) Surface

Developing an understanding of the response of mineral/water interfaces to applied electric fields is central to detecting and interpreting signatures of interfacial processes in the subsurface. Here we focus density functional theory calculations on understanding K + cation binding and diffusion across the (001) surface of orthoclase feldspar under various applied electric fields with and without surface hydration. The calculations reveal how water ligands labilize surface K + cations for diffusion while also increasing their sensitivity to electric field effects on the binding energy at different surface sites. The calculations also show how the direction and strength of the electric field systematically affect surface cation mobility, sorption, and hydration behavior. Specifically, electric fields directed toward the surface increase the potassium binding energies, facilitate K + diffusion across multiple binding sites, and increase the propensity of surface K + ions to lose their hydration shells. Furthermore, the findings help fill a major knowledge gap into the impact of electric fields on mineral/water interface structure and dynamics more generally, featuring in this case a commonly found type of feldspar involved in a multitude of atmospheric and geochemical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the mechanisms of anisotropic dissolution in metal oxides by applying radiolysis simulations to liquid-phase TEM

Iron-based redox-active minerals are ubiquitous in soils, sediments, and aquatic systems. Their dissolution is of great importance for microbial impacts on carbon cycling and the biogeochemistry of the lithosphere and hydrosphere. Despite its widespread significance and extensive prior study, the atomic-to-nanoscale mechanisms of dissolution remain poorly understood, particularly the interplay between acidic and reductive processes. Here, we use in situ liquid-phase-transmission electron microscopy (LP-TEM) and simulations of radiolysis to probe and control acidic versus reductive dissolution of akaganeite (β–FeOOH) nanorods. Informed by crystal structure and surface chemistry, the balance between acidic dissolution at rod tips and reductive dissolution at rod sides was systematically varied using pH buffers, background chloride anions, and electron beam dose. We find that buffers, such as bis-tris, effectively inhibited dissolution by consuming radiolytic acidic and reducing species such as superoxides and aqueous electrons. In contrast, chloride anions simultaneously suppressed dissolution at rod tips by stabilizing structural elements while promoting dissolution at rod sides through surface complexation. Dissolution behaviors were systematically varied by shifting the balance between acidic and reductive attacks. The findings show LP-TEM combined with simulations of radiolysis effects can provide a unique and versatile platform for quantitatively investigating dissolution mechanisms, with implications for understanding metal cycling in natural environments and the development of tailored nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Simulations of Tritium Diffusion in Al-Rich Iron Aluminide Coating Phases

Surface coatings of steels used in extreme conditions and corrosive environments generally aim to provide protection and increased durability. In the case of tritium-producing burnable absorber rods (TPBARs) used in nuclear reactors, a 316 stainless steel has been coated with Al. Scanning transmission electron microscopy (STEM) characterization of the coating found three Al-rich (>60 atom % Al) iron aluminide alloys identified as hexagonal FeNiAl 5 , monoclinic Fe 4 Al 13 , and orthorhombic Fe 2 Al 5 . Density functional theory simulations using nudged elastic band have been performed to investigate the diffusion of interstitial tritium in each Al-rich iron aluminide phase. While FeNiAl 5 and Fe 4 Al 13 can be viewed as the stacking of two layers, the structural peculiarity of Fe 2 Al 5 is that channels of variable Al vacancy content are present along the c-axis. Therefore, three stoichiometries for Fe 2 Al x phase, namely, Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 , have been considered to evaluate the impact of Al vacancy concentration on tritium diffusion behavior. Altogether, we found that at 600 K, tritium diffusion decreases from a faster rate in the channels of Fe 2 Al x phases (D T ≤ 10 –11 m 2 ·s –1 ) to Fe 4 Al 13 (D T ≈ 10–12 m2·s–1), and finally in FeNiAl5 (D T ≈ 10 –13 m 2 ·s –1 ). Here, we also find that interstitial tritium generally diffuses faster in Fe–Al coating phases than in the tritium breeding material γ-LiAlO 2 (D T ≈ 10 –14 m 2 ·s –1 ) but slightly slower than in 316 stainless steel (D T ≈ 10 –10 m 2 ·s –1 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A density functional theory analysis of the adsorption and surface chemistry of inorganic iodine species on graphitea

In the event of a nuclear accident, fission products may be released into the environment. The release of 131 I is of particular concern to human health. Iodine can be captured using a number of materials and frequently, this is accomplished with activated carbon impregnated with organic bases. Previous studies have used DFT and the graphite (0001) surface as a surrogate for adsorption, those studies focus on the species I • , I 2 , and CH 3 I. In this work we perform an ab initio study of the adsorption onto the surface of a graphite sheet of I 2 , CH 3 I, and inorganic acidic iodine species (HI, HOI, HIO 2 , and HIO 3 ), which were selected to examine the possible effect of oxidation state on adsorption. The PBE exchange-correlation functional with D3 dispersion was employed. It was found that for molecular iodine, the iodine atoms tended to either situate above the center of a hexagonal site on the graphite or directly atop a carbon atom with the lighter components resting closer to the graphite. For each species the relative binding energies spanned the range of 21–33 kJ mol -1 and graphite-iodine distance was in the range of 3.52–3.93 Å. In all cases we found no significant charge transfer between the iodine species and the graphite, thus we conclude that all the iodine species studied undergo strong physisorption to the graphite.

Ritzmann, Andrew M.↗

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗