Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic fine structure”

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

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

At least 145 records · Page 8

Structural Changes to the Gd‐DTPA Complex at Varying Ligand Protonation State

Abstract Diethylenetriaminepentaacetic acid (DTPA) is a chelating agent whose complex with the Gd 3+ ion is used in medical imaging. DTPA is also used in lanthanide‐actinide separation processes. As protonation of the DTPA ligand can facilitate dissociation of the Gd 3+ ion from the Gd‐DTPA complex, this work investigates the coordination structures of the aqueous Gd 3+ ion and its environment when chelated by DTPA in eight different DTPA protonation states. Both classical and ab initio molecular dynamics (MD) simulations are conducted to model the solvated complexes. Extended X‐ray absorption fine structure (EXAFS) measurements of the Gd 3+ aqua ion, and the Gd‐DTPA complex at pH 1 and 11, are compared to EXAFS spectra predicted from the MD simulations to verify the accuracy of the MD structures. The findings of this work provide atomic‐level details into the fluctuating Gd‐DTPA complex environment as the DTPA ligand gradually detaches from the Gd 3+ ion with increased protonation.

Chemistry↗

Structure and Oxygen Evolution Activity of β-NiOOH: Where Are the Protons?

Ni oxides and oxyhydroxides (NiO x ) have been studied for a long time as cathode materials for alkaline batteries and electrocatalysts for the oxygen evolution reaction (OER). Yet, understanding of the connection between their atomic and electronic structures and electrochemical performance or stability is still incomplete. In this work, we use first-principles density functional theory (DFT) calculations to revisit the structure, electronic properties, and OER activity of β-NiOOH, the catalytically active phase of NiO x . Following extensive DFTbased screening, we identify a hitherto overlooked structure characterized by a uniform distribution of H atoms on the NiO 2 layers. All the Ni 3+ cations in this structure exhibit an identical t g 6 e g 1 electronic configuration with an occupied 3d z2 orbital. Comparison of the calculated bond lengths with extended X-ray absorption fine structure (EXAFS) data unequivocally supports this structure relative to all other low-energy configurations. Based on this structure, we uncover and detail defect-dominated OER mechanisms on the basal β-NiOOH (001) surface, with overpotentials as low as 0.39 V. The present results should provide a valuable contribution to ongoing efforts for understanding and developing enhanced transition-metal hydroxide catalysts for the OER.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient and interpretable graph network representation for angle-dependent properties applied to optical spectroscopy

Abstract Graph neural networks are attractive for learning properties of atomic structures thanks to their intuitive graph encoding of atoms and bonds. However, conventional encoding does not include angular information, which is critical for describing atomic arrangements in disordered systems. In this work, we extend the recently proposed ALIGNN (Atomistic Line Graph Neural Network) encoding, which incorporates bond angles, to also include dihedral angles (ALIGNN-d). This simple extension leads to a memory-efficient graph representation that captures the complete geometry of atomic structures. ALIGNN-d is applied to predict the infrared optical response of dynamically disordered Cu(II) aqua complexes, leveraging the intrinsic interpretability to elucidate the relative contributions of individual structural components. Bond and dihedral angles are found to be critical contributors to the fine structure of the absorption response, with distortions that represent transitions between more common geometries exhibiting the strongest absorption intensity. Future directions for further development of ALIGNN-d are discussed.

36 MATERIALS SCIENCE↗

Computational exploration of DIII-D QH-mode with extended MHD

This work will enable improved understanding of the scientific basis for DIII-D Quiescent H-mode (QH-mode) by leveraging recent computational developments of the NIMROD code. Experimentally, it is established that the QH-mode tokamak operation regime is free of edge-localized modes (ELMs) with good energy confinement. ELMs produce an impulsive heat-flux on the divertor that is unacceptable under burning-plasma conditions. Beyond DIII-D in future high-power-density tokamaks, accommodation of a high heat flux at the divertor represents a key challenge to overcome. Proposed solutions involve operation at high divertor density to achieve divertor detachment, operation with a liquid lithium wall and/or introduction of impurities to enhance radiative losses. Additionally, the first wall material is likely to be a high-Z material such as Tungsten in future devices as materials such as graphite tiles restrict the lifetime of the first wall due to a high erosion rate. While the atomic species composition of the plasma is known to impact the edge-pedestal structure and stability, the specifics of the underlying physics basis for this interaction is not firm. Studying this interaction is the goal of this proposal without focusing specifically on the fine details of the first wall and divertor solution. It is the goal is to lay the ground work to establish the compatibility of the divertor solutions and wall material with pedestal solutions that eliminate ELMs by leveraging present DIII-D experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressure-Driven Helium Insertion for Structural Stability of CH 3 NH 3 PbBr 3 Hybrid Perovskites

Organic−inorganic metal halide perovskites (MHPs) have garnered significant attention due to their outstanding performance in optoelectronic devices, but they are prone to degradation through a variety of pathways. High-pressure studies are being used to understand the mechanical and structural stability of MHPs and investigate new methods for improving these. In this study, we map the high-pressure, low-temperature structural phase diagram of CH 3 NH 3 PbBr 3 (MAPbBr 3 ) between 15−300 K and up to 1.5 GPa. We first compare the temperature and pressure effects on the global and local structures of MAPbBr 3 using synchrotron X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS), respectively, and find evidence of PbBr 6 octahedral distortion. Our results also suggest that He inserts into the MAPbBr 3 structure. The thermodynamics of He insertion into MAPBHe x are calculated to be plausible at room temperature for x ≤ 1 and P = 1−3 GPa and shown to stabilize higher-symmetry phases. This phenomenon increases the fitted bulk modulus of MAPbBr 3 . Structural stabilization by inert atom insertion is proposed as a method to improve the mechanical robustness of MHPs and other soft hybrid materials.

Diffraction↗

Crystal Structure of Moganite and Its Anisotropic Atomic Displacement Parameters Determined by Synchrotron X-ray Diffraction and X-ray/Neutron Pair Distribution Function Analyses

The crystal structure of moganite from the Mogán formation on Gran Canaria has been re-investigated using high-resolution synchrotron X-ray diffraction (XRD) and X-ray/neutron pair distribution function (PDF) analyses. Our study for the first time reports the anisotropic atomic displacement parameters (ADPs) of a natural moganite. Rietveld analysis of synchrotron XRD data determined the crystal structure of moganite with the space group I2/a. The refined unit-cell parameters are a = 8.7363(8), b = 4.8688(5), c = 10.7203(9) Å, and β = 90.212(4)°. The ADPs of Si and O in moganite were obtained from X-ray and neutron PDF analyses. The shapes and orientations of the anisotropic ellipsoids determined from X-ray and neutron measurements are similar. The anisotropic ellipsoids for O extend along planes perpendicular to the Si-Si axis of corner-sharing SiO4 tetrahedra, suggesting precession-like movement. Neutron PDF result confirms the occurrence of OH over some of the tetrahedral sites. We postulate that moganite nanomineral is stable with respect to quartz in hypersaline water. The ADPs of moganite show a similar trend as those of quartz determined by single-crystal XRD. In short, the combined methods can provide high-quality structural parameters of moganite nanomineral, including its ADPs and extra OH position at the surface. This approach can be used as an alternative means for solving the structures of crystals that are not large enough for single-crystal XRD measurements, such as fine-grained and nanocrystalline minerals formed in various geological environments.

58 GEOSCIENCES↗

Nonadiabatic coupling of the dynamical structure to the superconductivity in YSr 2 Cu 2.75 Mo 0.25 O 7.54 and Sr 2 CuO 3.3

A crucial issue in cuprates is the extent and mechanism of the coupling of the lattice to the electrons and the superconductivity. Here we report Cu K edge extended X-ray absorption fine structure measurements elucidating the internal quantum tunneling polaron (iqtp) component of the dynamical structure in two heavily overdoped superconducting cuprate compounds, tetragonal YSr 2 Cu 2.75 Mo 0.25 O 7.54 with superconducting critical temperature, T c = 84 K and hole density p = 0.3 to 0.5 per planar Cu, and the tetragonal phase of Sr 2 CuO 3.3 with T c = 95 K and p = 0.6. In YSr 2 Cu 2.75 Mo 0.25 O 7.54 changes in the Cu-apical O two-site distribution reflect a sequential renormalization of the double-well potential of this site beginning at T c , with the energy difference between the two minima increasing by ~6 meV between T c and 52 K. Sr 2 CuO 3.3 undergoes a radically larger transformation at T c , >1-Å displacements of the apical O atoms. The principal feature of the dynamical structure underlying these transformations is the strongly anharmonic oscillation of the apical O atoms in a double-well potential that results in the observation of two distinct O sites whose Cu–O distances indicate different bonding modes and valence-charge distributions. The coupling of the superconductivity to the iqtp that originates in this nonadiabatic coupling between the electrons and lattice demonstrates an important role for the dynamical structure whereby pairing occurs even in a system where displacements of the atoms that are part of the transition are sufficiently large to alter the Fermi surface. The synchronization and dynamic coherence of the iqtps resulting from the strong interactions within a crystal would be expected to influence this process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Zr and Sc additions on coarsening- and creep resistance of AlSi10Mg fabricated by laser powder bed fusion

Microstructure and creep properties are studied in a eutectic AlSi10Mg alloy modified with Zr and Sc additions (Al-9.8Si-0.32Mg-0.70Zr-0.23Sc, wt%) manufactured through laser powder-bed fusion (L-PBF). Three types of Zr/Sc-bearing powders - elemental, master-alloy, and pre-alloyed - are employed in the fabrication process, with the pre-alloyed powders providing the highest incorporation of Zr and Sc in the Al matrix. Here, the as-printed alloy exhibits a fine cellular Al-Si eutectic structure which fragments and coarsens into micron-sized Si particles during aging at 300 ºC, leading to a steady drop in alloy microhardness between 0.1 and 1000 h. Coarsening of the eutectic Si phase during aging is not measurably affected by Zr and Sc in solid solution, which precipitate during aging and increases strength, compensating weakening from Si coarsening between 1 and 200 h at 300 ºC. Atom-probe analysis in the peak-aged condition (96 h/300 °C) confirms the presence of Al 3 (Sc,Zr) secondary nano-precipitates with an average radius of 1.0 nm and some Si solubility. Micron-sized grains are present throughout the alloy, with ultra-fine-grained regions at the melt pool boundaries, neither of which coarsen during long-term aging at 300 ºC. Under creep conditions at 300 °C, the Zr/Sc-bearing alloy with Al 3 (Sc,Zr) nano-precipitates exhibits power-law behavior, with a high apparent stress exponent (n a = 9) and a high threshold stress (σ th = 43 MPa), exhibiting nearly double the strength of a Zr/Sc-free AlSi10Mg control alloy (σ th = 22 MPa).

36 MATERIALS SCIENCE↗

Observing 3-hydroxyanthranilate-3,4-dioxygenase in action through a crystalline lens

The synthesis of quinolinic acid from tryptophan is a critical step in the de novo biosynthesis of nicotinamide adenine dinucleotide (NAD + ) in mammals. Herein, the nonheme iron-based 3-hydroxyanthranilate-3,4-dioxygenase responsible for quinolinic acid production was studied by performing time-resolved in crystallo reactions monitored by UV-vis microspectroscopy, electron paramagnetic resonance (EPR) spectroscopy, and X-ray crystallography. Seven catalytic intermediates were kinetically and structurally resolved in the crystalline state, and each accompanies protein conformational changes at the active site. Among them, a monooxygenated, seven-membered lactone intermediate as a monodentate ligand of the iron center at 1.59-Å resolution was captured, which presumably corresponds to a substrate-based radical species observed by EPR using a slurry of small-sized single crystals. Other structural snapshots determined at around 2.0-Å resolution include monodentate and subsequently bidentate coordinated substrate, superoxo, alkylperoxo, and two metal-bound enol tautomers of the unstable dioxygenase product. These results reveal a detailed stepwise O-atom transfer dioxygenase mechanism along with potential isomerization activity that fine-tunes product profiling and affects the production of quinolinic acid at a junction of the metabolic pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning Structural and Electronic Configuration of FeN 4 via External S for Enhanced Oxygen Reduction Reaction

The Fe–N–C material represents an attractive oxygen reduction reaction electrocatalyst, and the FeN 4 moiety has been identified as a very competitive catalytic active site. Fine tuning of the coordination structure of FeN 4 has an essential impact on the catalytic performance. Herein, we construct a sulfur‐modified Fe–N–C catalyst with controllable local coordination environment, where the Fe is coordinated with four in‐plane N and an axial external S. The external S atom affects not only the electron distribution but also the spin state of Fe in the FeN 4 active site. The appearance of higher valence states and spin states for Fe demonstrates the increase in unpaired electrons. With the above characteristics, the adsorption and desorption of the reactants at FeN 4 active sites are optimized, thus promoting the oxygen reduction reaction activity. This work explores the key point in electronic configuration and coordination environment tuning of FeN 4 through S doping and provides new insight into the construction of M–N–C‐based oxygen reduction reaction catalysts.

25 ENERGY STORAGE↗

Application of Single-Site Catalysts in the Hydrogen Economy

Catalysts development is an evolutionary process. The generation of new single-metal-site catalysts, with isolated metal atoms dispersed on the surface of substrates, offers maximum atom-utilization efficiency, high specific activity, and novel properties. The unique structure of single-site catalysts also enables their potential to bridge the gap between homogeneous and heterogeneous catalysis. This field is at an early stage, characterized by important discoveries and novel demonstrations of catalytic properties associated with the fine structure of single metal sites; however, specific applications, especially catalytic hydrogenation reactions, are still lacking. Here, we focus on applications of single-site catalysts in the hydrogen economy based on our studies of hydrogen generation/storage via single-metal-site catalysts. We hope this work will inspire the development of single-site catalysts for other applications.

08 HYDROGEN↗

Atomically Dispersed Platinum in Surface and Subsurface Sites on MgO Have Contrasting Catalytic Properties for CO Oxidation

Atomically dispersed metals on metal oxide supports are a rapidly growing class of catalysts. Developing an understanding of where and how the metals are bonded to the supports is challenging because support surfaces are heterogeneous, and most reports lack a detailed consideration of these points. In this work, we report two atomically dispersed CO oxidation catalysts having markedly different metal–support interactions: platinum in the first layer of crystalline MgO powder and platinum in the second layer of this support. Structural models have been determined on the basis of data and computations, including those determined by extended X-ray absorption fine structure and X-ray absorption near edge structure spectroscopies, infrared spectroscopy of adsorbed CO, and scanning transmission electron microscopy. The data demonstrate the transformation of surface to subsurface platinum as the temperature of sample calcination increased. Catalyst performance data demonstrate the lower activity but greater stability of the subsurface platinum than of the surface platinum.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic insight into copper cation exchange in cadmium selenide semiconductor nanocrystals using X-ray absorption spectroscopy

In terms of producing new advances in sustainable nanomaterials, cation exchange (CE) of post-processed colloidal nanocrystals (NCs) has opened new avenues towards producing non-toxic energy materials via simple chemical techniques. The main processes governing CE can be explained by considering hard/soft acid/base theory, but the detailed mechanism of CE, however, has been debated and has been attributed to both diffusion and vacancy processes. In this work, we have performed in situ x-ray absorption spectroscopy to further understand the mechanism of the CE of copper in solution phase CdSe NCs. The x-ray data indicates clear isosbestic points, suggestive of cooperative behavior as previously observed via optical spectroscopy. Examination of the extended x-ray absorption fine structure data points to the observation of interstitial impurities during the initial stages of CE, suggesting the diffusion process is the fundamental mechanism of CE in this system.

74 ATOMIC AND MOLECULAR PHYSICS↗

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Porous Semiconducting K–Sn–Mo–S Aerogel: Synthesis, Local Structure, and Ion-Exchange Properties

Chalcogenide-based aerogels are emerging porous semiconducting nanomaterials that appeal to applications in clean energy and the environment. Here, we report a novel gel, potassium–tin–molybdenum–sulfides (KTMS), that integrates the electrostatically bound K + ions in the covalent network of Sn–Mo–S. Its gelation requires a concurrent reduction of Mo 6+ → Mo 4+/5+ and the oxidation of S 2– → Sn – (n ≈ 1) and Sn 2+ → Sn 4+ . KTMS is an amorphous semiconductor showing quantum confinement effects on band gap energies, 2.1 → 1.4 → 0.9 eV for its wet- → aero- → xerogels. Synchrotron X-ray pair distribution function (PDF) and extended X-ray absorption fine structure (EXAFS) revealed a complex local structure of KTMS consisting of molecular Mo 2 (S 2 ) 6 and Mo 3 S(S 2 ) 6 clusters. In addition, the Sn–S coordination is related to crystalline Na4Sn3S8 and SnS2. KTMS also demonstrated the removal of the radionuclides of Cs + , Sr 2+ , and UO 2 2+ from ppm to ppb levels with distribution constants (Kd) up to ≥104 mL/g. Notably, despite the lack of atomic periodicity in the amorphous KTMS, the K+ ion is ion-exchangeable with chemically diverse Sr 2+ , Cs + , and UO 2 2+ in aqueous solutions; especially the ion-exchange properties of Sr 2+ and UO 2 2+ ≡(O=U=O) 2+ is not known to any chalcogels known to date. The sequestration of Cs + and Sr 2+ was achieved by the exchange of K + in the amorphous KTMS, and the removal of [O=U 6+ =O] 2+ synergistically involves surface sorption via -S····U 6+ =O 2 2+ covalent interactions and ion-exchange via the hard–soft Lewis acid–base paradigm. Overall, cooperative roles played by the diverse bonding motifs, surface-exposed Lewis basic frameworks, and polarizability of the (poly)sulfides make it an exceptional adsorbent for chemically diverse radioactive species. This finding will guide the design of superior sorbents for chemically distinct metal ion separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidating the Structure of the Eu‐EDTA Complex in Solution at Various Protonation States

Abstract Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu 3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu‐EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X‐ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu 3+ aqua ion and Eu‐EDTA complexes at pH 3 and 11. This work details how Eu‐EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu 3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation.

Chemistry↗

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗