Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Surface Dynamics of Selective Hydration of Rutile TiO 2 : A Kinetic Monte Carlo Approach

The utility of atomic layer deposition (ALD) continues to expand beyond conformal thin-film fabrication to include area- or site-selective ALD. We previously identified a strategy for site-selective ALD (SS-ALD) via the evaluation of thermodynamically preferred hydration of rutile TiO 2 surfaces, as deduced from electronic structure calculations. Here, we report a novel kinetic Monte Carlo (KMC) model that allows for the investigation of surface dynamics and kinetics that improves our understanding of and intuition for the selective hydration strategy. Here we demonstrate the validity of the strategy with respect to step-edge defects for the lowest energy (110) facet as well as report results for the other common facets which agree with experimental STM observations. The results here indicate that the selective hydration strategy is feasible both thermodynamically (evaluated in our previous publication) and kinetically (from the KMC model). Because diffusion has a slower rate than others, we find that any proximity effects between terrace and defect sites are unlikely to affect the selective hydration strategy for rutile TiO 2 . The KMC model further provides relevant timescales for achieving selectivity experimentally and establishes the kinetic viability of the selective hydration approach to SS-ALD.

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

Ultraviolet spectroscopy of pressurized and supercritical carbon dioxide

Abstract Carbon dioxide (CO 2 ) is prevalent in planetary atmospheres and sees use in a variety of industrial applications. Despite its ubiquitous nature, its photochemistry remains poorly understood. In this work we explore the density dependence of pressurized and supercritical CO 2 electronic absorption spectra by vacuum ultraviolet spectroscopy over the wavelength range 1455-2000 Å. We show that the lowest absorption band transition energy is unaffected by a density increase up to and beyond the thermodynamic critical point (137 bar, 308 K). However, the diffuse vibrational structure inherent to the spectrum gradually decreases in magnitude. This effect cannot be explained solely by collisional broadening and/or dimerization. We suggest that at high densities close proximity of neighboring CO 2 molecules with a variety of orientations perturbs the multiple monomer electronic state potential energy surfaces, facilitating coupling between binding and dissociative states. We estimate a critical radius of ~4.1 Å necessary to cause such perturbations.

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

Facile electrocatalytic proton reduction by a [Fe–Fe]-hydrogenase bio-inspired synthetic model bearing a terminal CN — ligand

An azadithiolate bridged CN — bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe–Fe] H 2 ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H + to H 2 between pH 0–3 at diffusion-controlled rates (10 11 M —1 s —1 ) i.e. 10 8 s —1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN— ligand increases the p K a of the cluster enabling hydrogen production from its Fe(I)–Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)–Fe(0) state. The formation of a terminal Fe–H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H 2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.

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

Enabling hydrogen chemisorption on charged graphene

Two-dimensional (2D) materials, including graphitic derivatives, have long been of interest for hydrogen storage applications, due to their high theoretical storage capacity, low weight, and other useful properties. However, poor kinetics for hydrogen adsorption and surface diffusion as part of the proposed spillover process for hydrogenation have limited their technological potential. Here, we use first-principles calculations to study electronic doping as a means to improve hydrogen chemisorption on graphene, which we use here as a proxy for graphitic derivatives more broadly. We find that positively charged graphene sheets have vastly improved kinetics for hydrogen diffusion and adsorption, while they limit unwanted hydrogen desorption. This combination of effects should favor hydrogen chemisorption via spillover. We connect these trends to the C–H bond, which introduces states near the Fermi level. These states are depopulated as electrons are removed, thereby lowering the bond energy and permitting more facile movement of hydrogen. Our results suggest that spillover mechanisms for hydrogen chemisorption should be revisited if strategies to apply a large charge to graphitic systems can be realized. Moreover, switchable application of the charge may lead to the reversible chemisorption of hydrogen. While the large magnitude of charging required suggests that graphene itself may not be suitable for reversible hydrogen chemisorption, the factors we identify and discuss could significantly boost the prospects of graphitic derivatives and other 2D or layered materials for hydrogen storage applications.

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

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2021

Alpha particle channeling in ITER

Alpha particle channeling, the attempt to extract energy from high-energy fusion alpha particles using a combination of high- and low-frequency waves introduced in the plasma with antennae, is studied using a guiding center code. We examine the effect of the combination of cooling caused by a cyclotron mode, particle diffusion due to Alfvén modes introduced for ash removal, the inclusion of microturbulent diffusion, and the slowing down of the alpha particles on electrons. Induced cooling is achieved with a cyclotron resonance, producing additional losses of particles above 3.5 MeV in an amount about equal to the prompt loss. The removal of the low-energy cooled particles by additional modes or by microturbulent diffusion is the hardest challenge for alpha channeling. The cooled alpha particles remain deep in the plasma center, so for this removal, the method used must be capable of providing transport from the plasma center to the edge without also causing losses of high-energy particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Atomic layer deposition of nanofilms on porous polymer substrates: Strategies for success

Atomic layer deposition (ALD) is a versatile technique for engineering the surfaces of porous polymers, imbuing the flexible, high-surface-area substrates with inorganic and hybrid material properties. Previously reported enhancements include fouling resistance, electrical conductance, thermal stability, photocatalytic activity, hydrophilicity, and oleophilicity. However, there are many poorly understood phenomena that introduce challenges in applying ALD to porous polymers. In this paper, we address five common challenges and ways to overcome them: (1) entrapped precursor, (2) embrittlement, (3) film fracture, (4) deformation, and (5) pore collapse. These challenges are often interrelated and can exacerbate one another. To investigate these phenomena, we applied various ALD chemistries to porous polymers including polyethersulfone, polysulfone, polyvinylidene fluoride, and polycarbonate track-etched membranes. Reaction-diffusion modeling revealed why certain precursors and processing conditions result in embrittling subsurface material growth, entrapment of unreacted precursors, and nongrowth. We quantify the limits of ALD processing temperatures that are dictated by thermal expansion mismatch and can lead to fractured ALD films. The results herein allow us to make recommendations to avoid, mitigate, or overcome the difficulties encountered when performing ALD and plasma-enhanced ALD on porous polymers. We intend this article to serve as a “lessons learned” guide informed by previous experience to provide a better understanding of the difficulties and limitations of ALD on porous polymers and knowledge-based guidelines for successful depositions. This knowledge can accelerate future research and help experimentalists navigate and troubleshoot as they expose porous polymers to reactive precursor vapors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predictive Integrated Modeling of Heated Tokamak Plasmas (Final Technical Report)

The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analysis of Tempered Fractional Operators

Tempered fractional operators provide an improved predictive capability for modeling anomalous effects that cannot be captured by standard partial differential equations. These effects include subdiffusion and superdiffusion (i.e. the mean square displacement in a diffusion process is proportional to a fractional power of the time), that often occur in, e.g., geoscience and hydrology. We analyze tempered fractional operators within the nonlocal vector calculus framework in order to assimilate them to the rigorous mathematical structure developed for nonlocal models. First, we show they are special instances of generalized nonlocal operators by means of a proper choice of the nonlocal kernel. Then, we present a plan for showing tempered fractional operators are equivalent to truncated fractional operators. These truncated operators are useful because they are less computationally intensive than the tempered operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Chromium Tolerant, Highly Active and Stable Electrocatalytic Internal Surface Coating for Cathode of Commercial SOFCs (Final Report)

This project is aimed to develop a chromium (Cr) tolerant, highly active, and stable coating layer on the internal surfaces of the porous composite cathode from commercially available SOFCs. Such coating layer was developed using the additive manufacturing process of Atomic Layer Deposition (ALD) and has been applied on the cathode consisting of either an electronic conductor of LaxSr 1-x MnyO 3-δ (LSM) or mixed ionic and electronic conducting La x Sr 1-x Co y Fe 1-y O 3-δ (LSCF). PI's work has demonstrated that the internal surface of cathode from the commercial cells, can be further tailored using ALD coating to dramatically enhance the cell performance. For instance, ALD layer consisting heterostructured nano composite of nano-Pt and nano-(Mn 0.8 Co 0.2 ) 3 O 4 oxide on the internal surface of porous LSM/YSZ cathode from SOFCs, has resulted in the large reduction of the cell polarizations resistance by up to 53%, and enormous increase of power density over 370%. For the cells with LSCF/Sm 2 O 3 doped CeO 2 (SDC) cathode, the conformal layer of nano-composite consisting of superjacent CoOx and subjacent minimum amount of Pt nano-grains has resulted in the power density enhancement by 126% for the large scale industry tubular cells at 750°C, and both the performance enhancement and nanostructure of the ALD layer are stable over ~ 2000 h continuous operation performed at industry test station. In the meanwhile, those ALD coating layer developed by PI's work is also inherently Cr-tolerant, and could act as physical barrier for preventing Cr diffusion into the cathode backbone, so as to mitigate the Cr poisoning effect on the cathode. In this project, the impact of Cr on the performance of those ALD coated commercial cells has been evaluated. Based on evolution of the cell performance, the ALD coating layer chemistry and ALD coating layer thickness has been optimized to maximize the overall Cr tolerance, cell power density and cell longevity. Different ALD coating has been applied onto the internal surface of LSM/YSZ and LSCF/SDC backbone respectively. The architecture/scaffold structures on the internal surface of different cathode, designed by this project, was catalogued and analyzed using High Resolution Transmission Electron Microscopy (HRTEM), and cell power/durability performance are assured via comprehensive electrochemical performance testing in the industry operation relevant conditions. The impact of the electrochemical operation current density, the water humidity, the cell operation temperature, and cell operation duration on the Cr tolerance of ALD coated cells has been systematically investigated. There is completely different nanostructure degradation mechanisms between LSM and LSCF cells induced by Cr contamination. For the LSCF/SDC baseline cell, With the Cr source, there is no apparent Sr surface segregation phase even for the baseline cell operated for 3000 h at 750 °C. With the Cr source, there is significant amorphous (SrCr)Ox phase accumulated in the original pore region. For the commercial baseline cells, Cr contaminants on the LSM electrode severely impacted the entire cell's electrochemical performance and nanostructure degradation. Those degradations include (1). Peak power density loss of 64 % after 109 h of operation. The dramatic increase in Rp (2). They are cracking at LSM/SSZ interface, LSM grains. SSZ remains intact but with (CrMn)Ox. By contrast, ALD coating (MnCo)Ox/Pt dramatically improves the Cr resistance, as follows (1). ALD-coated cell with a power density is 280-380 % of the baseline cell, depending on the ALD layer thickness. (2). For a cell with a 20 nm thick ALD layer, there is a large performance enhancement (> 200 % power density) induced by ALD coating of Cr-tolerant Mn 0.8 Co 0.2 Ox. (3). For a cell with a 20 nm thick ALD layer, after 168 h at 750 °C power density of the ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination for 109 h. The ALD coating on the internal surface of cathode developed by this project integrated multi-functions. Those multi-functions include (1). Dramatically improving the cell power density for the commercial cells; (2). Dramatically improving contamination resistance of the cathode, for being an excellent protection coating layer sealing off Cr contamination. (3). Dramatically increasing the cell longevity by potentially preventing the microstructure evolution and grain coarsening of the cathode. Overall, this project will provide a simple solution to simultaneously enhance power density and increase the reliability, robustness, and endurance of commercial SOFCs, over the entire operating temperature range of 650-800 °C. For the inherently functional SOFC, the ALD coating of LSM based cathode mitigate the Cr-contamination. Power density of ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination. In addition to SOFCs, the novel on-demand design approach and creation of multifunctional heterogeneous architecture on the electrode surface presented in this work opens further research for their application in other types of fuel cells, batteries, and sensors for which electrochemical reactions on the surface are similarly critical.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of irradiation temperature on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents the experimental findings obtained from November 2023 to September 2024. The study aims to investigate the effects of temperature on the microstructure, deuterium (D) retention, and lithium (Li) loss in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation to a high dose. Sequential ion irradiation has been a key method in our previous studies to simulate the behavior of γ-LiAlO 2 under neutron irradiation. In addition to dose and dose rate, irradiation temperature is a critical factor influencing microstructural and compositional changes. LiAl 5 O 8 precipitates have been observed in γ-LiAlO 2 pellets irradiated with reactor neutrons at ~300°C. These precipitates also form during ion irradiation at an elevated temperature of 500°C, but not at 300°C. In our ion irradiation experiments, the dose rate is typically three orders of magnitude higher than that of neutron irradiation, leading to a more rapid damage production. To better emulate the microstructural features in neutron-irradiated pellets using ion irradiation, a higher irradiation temperature is needed to accelerate the diffusion of point defects and enhance defect recovery rates, thereby compensating for the effects of the higher dose rate. The microstructural changes observed are the result of competing processes occurring during ion irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ion diffusion coefficients in poly(3-alkylthiophenes) for energy conversion and biosensing: role of side-chain length and microstructure

Conductive polymers are promising materials as active elements for energy storage and conversion devices due to mixed ion–electron conduction. The ion diffusion coefficient is a relative measure of the efficacy of ion transport, allowing for comparison between materials and electrochemical conditions. In this work, diffusion coefficients of hexafluorophosphate (PF6–) counterions in poly(3-alkylthiophene) (P3AT) materials are measured as a function of both side-chain length and microstructure using electrochemical impedance spectroscopy (EIS). For semi-crystalline films, the diffusion coefficient is found to be anomalous and nearly independent of applied electrochemical potential. Here, the anomalous behavior of diffusion indicates that spin casting yields compact films with an enthalpic barrier to ion transport, attributed to ionic trapping. Diffusion coefficient values ~10 –11 cm 2 s –1 were measured for all films, indicating interchain spacing, in the absence of strong intermolecular interactions with the electrolyte, is not a viable design strategy to control ion transport. For the prototypical system of poly(3-hexylthiophene), we observe almost no potential dependence in ion transport for regioregular and regiorandom films of comparable molecular weight, with both exhibiting anomalous diffusion. Alternatively, changing the microstructure of poly(3-hexylthiophene) to a mostly amorphous, ion-imprinted structure yields ~500× increase in the diffusion coefficient to ~2 × 10 –8 cm 2 s –1 at 0.8 V vs. Ag/Ag + with behavior closer to ordinary diffusion. Collectively, these results indicate new insight into ion transport in conductive polymers, where ionic trapping effects can be mitigated through electrodeposition protocols over post-synthesis processing (i.e. spin coating).

14 SOLAR ENERGY↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2023

Rubidium and potassium isotopic variations in chondrites and Mars: Accretion signatures and planetary overprints

As moderately volatile elements, isotopes of Rb and K can trace volatilization processes in planetary bodies. Rubidium isotopic data are however very scarce, especially for non-carbonaceous meteorites. Here, in this study, we report combined Rb and K isotopic data (δ 87/85 Rb and δ 41/39 Κ) for 7 ordinary, 6 enstatite, and 4 Martian meteorite falls to understand the causes for the variations in volatile abundances and isotopic compositions. Bulk Rb and K isotopic compositions of planetary bodies are estimated to be (Table 1): Mars +0.10 ± 0.03 ‰ for Rb and -0.26 ± 0.05 ‰ for K, bulk OCs $-0.12^{+0.15}_{-0.24}$ ‰ for Rb and $-0.72^{+0.28}_{-0.41}$ ‰ for K, bulk ECs $-0.02^{+0.29}_{-0.26}$ ‰ for Rb and $-0.33^{+0.37}_{-0.23}$ ‰ for K. The bulk K isotopic compositions of subgroup OCs are estimated to be $-0.72^{+0.26}_{-0.55}$ ‰ for H chondrites, $-0.71^{+0.23}_{-0.39}$ ‰ for L chondrites, and $-0.77^{+0.63}_{-0.30}$ ‰ for LL chondrites. A broad correlation between the Rb and K isotopic compositions of planetary bodies is observed. The correlation follows a slope that is consistent with kinetic evaporation and condensation processes, suggesting volatility-controlled mass-dependent isotope fractionation (as opposed to nucleosynthetic anomalies). Individual ordinary and enstatite chondrites show large Rb and K isotopic variations (-1.02 to +0.29 ‰ for Rb and -0.91 to -0.15 ‰ for K). Samples of lower metamorphic grades display correlated elemental and isotopic fractionations between Rb and K, while samples of higher metamorphic grades show great scatter, suggesting that chondrite parent-body processes have decoupled the two elements and their isotopes at the sample scale. Several processes could have contributed to the observed isotopic variations of Rb and K, including (i) chondrule “nugget effect”, (ii) volatilization during parent-body thermal metamorphism (heat-induced vaporization and gas transport within parent bodies), (iii) thermal diffusion during parent-body metamorphism, and (iv) impact/shock heating. Quantitative modeling of the first two processes suggests that neither of them could produce isotopic variations large enough to explain the observed isotopic variations. Volatilization during parent-body thermal metamorphism [the scenario (ii)], which has been commonly invoked to explain the isotopic variations of volatile elements, is gas transport-limited and its effect on isotopic fractionations of moderately volatile elements should be negligible. Modeling of diffusion processes suggests that (iii) could produce K isotopic variation comparable to the observed variation. The large isotopic variations in non-carbonaceous meteorites are thus most likely due to diffusive redistribution of K and Rb during metamorphism and/or shock-induced heating and vaporization.

58 GEOSCIENCES↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗