SEARCH · Engineering Papers

Results for “dynamic structure factor”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Coherences of Photoinduced Electron Spin Qubit Pair States in Photosystem I

Here, this publication presents the first comprehensive experimental study of electron spin coherences in photosynthetic reaction center proteins, specifically focusing on photosystem I (PSI). The ultrafast electron transfer in PSI generates spin-correlated radical pairs (SCRPs), which are entangled spin pairs formed in well-defined spin states (Bell states). Since their discovery in our group in the 1980s, SCRPs have been extensively used to enhance our understanding of structure–function relationships in photosynthetic proteins. More recently, SCRPs have been utilized as tools for quantum sensing. Electron spin decoherence poses a significant challenge in realizing practical applications of electron spin qubits, particularly the creation of quantum entanglement between multiple electron spins. This work is focused on the systematic characterization of decoherence in SCRPs of PSI. These decoherence times were measured as electron spin echo decay times, termed phase memory times (T M ), at various temperatures. Decoherence was recorded on both transient SCRP states P 700 + A 1 – and thermalized states. Our study reveals that T M exhibits minimal dependence on the biological species, biochemical treatment, and paramagnetic species. The analysis indicates that nuclear spin diffusion and instantaneous diffusion mechanisms alone cannot explain the observed decoherence. As a plausible explanation we discuss the assumption that the low-temperature dynamics of methyl groups in the protein surrounding the unpaired electron spin centers is the main factor governing the loss of the spin coherence in PSI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Chiral-odd generalized parton distributions in the large-𝑁 𝑐 limit of QCD: Spin-flavor structure, polynomiality, and sum rules

We study the nonperturbative properties of the nucleon’s chiral-odd generalized parton distributions (transversity GPDs) in the large-𝑁 𝑐 limit of QCD. This includes the parametric ordering of the spin-flavor components, the polynomiality property of the moments, and the sum rules connecting the GPDs with the tensor form factors. A multipole expansion in the transverse momentum transfer is used to enumerate and interpret the structures in the nucleon matrix element of the chiral-odd partonic operator, including monopole, dipole and quadrupole terms. The 1/𝑁 𝑐 expansion of the GPDs is performed using the abstract mean-field picture of baryons in the large-𝑁 𝑐 limit and its symmetries. We derive a large-𝑁 𝑐 relation between the flavor-nonsinglet GPDs 𝐸$^{𝑢−𝑑}_𝑇$ and $\tilde{𝐻}^{𝑢−𝑑}_𝑇$ and test it with recent lattice QCD results. We show that the polynomiality property and sum rules of the GPDs are fulfilled with the restricted realization of translational and rotational invariance in the mean-field picture. The results provide a basis for the phenomenological analysis of chiral-odd GPDs and hard exclusive processes in the large-𝑁 𝑐 limit, and for calculations in specific dynamical models.

generalized parton distributions↗

NASA NTRS1994

Scaling of impact damage in fiber composites

Impact damage in fiber composite structures remains of much concern, and is often the limiting factor in establishing allowable strain levels. The complexity of impact damage formation usually dictates that experiments are required, but scaling of results from small laboratory scale specimens to large structures introduces additional uncertainty into the analysis. This presentation gives the results of an analytical and experimental investigation intended to develop procedures for prediction of damage formation and subsequent strength loss, with particular emphasis on scaling of results with respect to structure size. The experimental investigation involved both drop-weight and airgun impact on carbon/epoxy plates and cylinders. Five sizes of plates ranging from 50 by 50 by 1.072 mm to 250 by 250 by 5.36 mm, and two sizes of cylinders with diameters of 96.5 and 319 mm, were employed in the experimental program. Impact tests were carried out over a range of impact conditions, and specimens were inspected for damage by C-scan and deplying. Analysis procedures were developed for both quasistatic and dynamic impacts for both the plates and cylinders. As has been reported previously, comparison of predicted structural response and measured surface strains was quite good over the entire range of sizes employed in the program. The damage formation and strength loss after impact showed a number of interesting features that are significant with respect to scaling of size. The extent of delamination was observed to increase with specimen size more than would be expected if stresses controlled the delamination extent. This was explained on the basis that delamination is controlled by energy release rates, and thus incorporates the usual dependence on the absolute size characteristic of fracture mechanics. Additionally, the experiments indicated that delamination initiated at matrix cracks and is dependent on the absolute size of the ply group thicknesses. Both the initiation and propagation of delamination are seen to be controlled by fracture mechanics parameters, and thus show specific dependence on size that must be accounted for in extrapolating results from laboratory scale tests to full size structures.

Swanson, Stephen R.↗

DOE OSTI2023

The kinetics of the PuO 2 to Pu 2 O 3 conversion

Here in an oxidizing environment, the oxide formed on plutonium (Pu) metal is composed of a plutonium dioxide (PuO 2 ) top layer and a thin cubic plutonium sesquioxide (Pu 2 O 3 ) middle layer. In a reducing environment, the PuO 2 layer auto-reduces to cubic Pu 2 O 3 . The speed and extent of this conversion depend on the combination of temperature and time. While PuO 2 provides a strong diffusion barrier against unwanted Pu corrosion by gaseous species (like hydrogen), Pu 2 O 3 does not, since its crystal structure has chains of oxygen vacancies. The kinetics of the PuO 2 reduction are, therefore, of fundamental interest and enable researchers to better protect Pu from corrosion. In this report, the oxygen-diffusion-limited kinetics of the dioxide to sesquioxide conversion were obtained by dynamically heating a PuO 2 -covered Pu sample from 294 to 418 K in a high-vacuum vessel equipped with an in situ spectroscopic ellipsometer. The physical/chemical constraints in the conversion process were combined with the ellipsometry method of multi-sample analysis to track the percentage of PuO 2 and to compute the extent of Pu 2 O 3 formation. The resulting diffusion coefficients were compared against and then combined with complementary literature data to produce a comprehensive set of kinetic parameters for reliably modeling oxide conversion over a larger temperature range than spanned by prior studies. The extracted thermal activation energy barrier (43.7 kJ/mol) and pre-exponential factor (5.0 × 10 -10 cm 2 /s) for the oxygen-diffusion-limited process can be used to accurately model the PuO 2 to Pu 2 O 3 transformation in vacuum and/or inert gas applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Development of Composite Photocatalyst Materials that are Highly Selective for Solar Hydrogen Production and their Evaluation in Z-Scheme Reactor Designs

The key technology gap preventing a vertically stacked dual-bed particle suspension reactor from achieving the DOE MYRD&D ultimate cost target for H 2 production remains the lack of materials in particle form factor that exhibit ≥10% solar-to-H 2 energy conversion (STH) efficiency as a suspension. Therefore, our project goals centered around strategies to increase the STH efficiency by enhancing photophysical properties of perovskite oxide particles including increased visible-light absorption, increased selectivity for electrocatalysis of the H 2 evolution reaction (HER) and the O 2 evolution reaction (OER) through development of ultrathin oxide coatings, correlating composition and structure to function, and improving understanding of multiscale transport and kinetic processes.

08 HYDROGEN↗

NASA NTRS2014

High-Fidelity Buckling Analysis of Composite Cylinders Using the STAGS Finite Element Code

Results from previous shell buckling studies are presented that illustrate some of the unique and powerful capabilities in the STAGS finite element analysis code that have made it an indispensable tool in structures research at NASA over the past few decades. In particular, prototypical results from the development and validation of high-fidelity buckling simulations are presented for several unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells along with a discussion on the specific methods and user-defined subroutines in STAGS that are used to carry out the high-fidelity simulations. These simulations accurately account for the effects of geometric shell-wall imperfections, shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform applied end loads, and elastic boundary conditions. The analysis procedure uses a combination of nonlinear quasi-static and transient dynamic solution algorithms to predict the prebuckling and unstable collapse response characteristics of the cylinders. Finally, the use of high-fidelity models in the development of analysis-based shell-buckling knockdown (design) factors is demonstrated.

Hilburger, Mark W.↗

NASA NTRS2022

Understanding Model-Observation Discrepancies in Satellite Retrievals of Atmospheric Temperature Using GISS ModelE

We examine multiple factors in the representation of satellite-retrieved atmospheric temperature diagnostics in historical simulations of climate change during the satellite era (specifically 1979–2021) using GISS ModelE contributions to the Coupled Model Intercomparison Project (Phase 6) (CMIP6). The tropospheric and stratospheric trends in these diagnostics are affected by greenhouse gases (notably carbon dioxide and ozone), coupling with the ocean, volcanic aerosols, solar activity and compositional and dynamic feedbacks. We explore the impacts of internal variability, changing forcing specifications, composition interactivity, the quality of the stratospheric circulation, vertical resolution, and possible impacts of the mis-specification of volcanic aerosol optical depths. Overall temperature trends throughout the satellite period are well captured, but discrepancies at all levels exist and have multiple distinct causes. We find that stratospheric comparisons (using Stratospheric Sounding Unit (SSU) retrievals and successor instruments) are most affected by variations in the representation of ozone depletion and feedbacks, followed by the volcanic signals. Tropospheric skill (using the Microwave Sounding Unit (MSU) retrievals) is affected by the trends in ocean heat uptake and tropospheric aerosols, but also by the representation of stratospheric processes through the impact of the Brewer-Dobson circulation on the height of the tropical tropopause. We demonstrate that no single factor is the dominant cause of the discrepancies and that almost all observations lie within the broad envelope of structural uncertainty.

satellite retrievals of atmospheric temperature↗

DOE OSTI2022

Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Real-time correlation of crystallization and segmental order in conjugated polymers

Modulating the segmental order in the morphology of conjugated polymers is widely recognized as a crucial factor for achieving optimal electronic properties and mechanical deformability. However, it is worth noting that the segmental order is typically associated with the crystallization process, which can result in rigid and brittle long-range ordered crystalline domains. To precisely control the morphology, a comprehensive understanding of how highly anisotropic conjugated polymers form segmentally ordered structures with ongoing crystallization is essential, yet currently elusive. To fill this knowledge gap, we developed a novel approach with a combination of stage-type fast scanning calorimetry and micro-Raman spectroscopy to capture the series of specimens with a continuum in the polymer percent crystallinity and detect the segmental order in real-time. Through the investigation of conjugated polymers with different backbones and side-chain structures, we observed a generally existing phenomenon that the degree of segmental order saturates before the maximum crystallinity is achieved. This disparity allows the conjugated polymers to achieve good charge carrier mobility while retaining good segmental dynamic mobility through the tailored treatment. Moreover, the crystallization temperature to obtain optimal segmental order can be predicted based on T g and T m of conjugated polymers. In conclusion, this in-depth characterization study provides fundamental insights into the evolution of segmental order during crystallization, which can aid in designing and controlling the optoelectronic and mechanical properties of conjugated polymers.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗

DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Toward large-scale simulation of residual stress and distortion in wire and arc additive manufacturing

In this study, we aim to advance the structural analysis of wire and arc additive manufacturing (WAAM) by considering the thermomechanical features inherent in direct energy deposition. Simulation approaches including the iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit finite element method (FEM) were developed for accelerating additive manufacturing stress analysis that is very time consuming by conventional numerical methods. ISM and DMRM take advantage of the strong nonlinearity phenomenon near the moving heat source by reducing the global iterations and refining the local mesh, respectively. In addition to the use of GPUs, the explicit FEM is accelerated by a time scaling technique based on the inherent strain concept. The residual stress and distortion of two large builds were analyzed, showing very consistent numerical results and good agreement with experiments. Compared with the commercial software Abaqus, the novel approaches reduced the computational cost substantially without compromising accuracy. The GPU code showed the highest computational efficiency (30~70×), while DMRM and ISM had acceleration factors of 9× and 3×, respectively. Finally, such high-fidelity modeling approaches will be very useful for building up a digital twin of WAAM to reduce development time and cost.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗

DOE OSTI2022

The TPR domain of PgaA is a multifunctional scaffold that binds PNAG and modulates PgaB-dependent polymer processing

The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗

DOE OSTI2024

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.

54 ENVIRONMENTAL SCIENCES↗