Search NASA⌕ Search

SEARCH · Search NASA

Results for “Trajectory analysis”

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

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure–Activity Relationships in Ether-Functionalized Solid-State Metal–Organic Framework Electrolytes

The structure–property relationships of metal–organic framework (MOF)-based solid-state electrolytes are not well understood. Herein, a systematic investigation of 12 Zr(IV)-based UiO-66 MOFs with varying ether-chain functional groups was carried out to elucidate the critical microscopic interactions that facilitate improved solid-state electrolyte performance. Enhanced sampling molecular dynamics (MD) simulations were employed and revealed a three-tier ion hopping mechanism: linker–linker hopping, linker-counterion hopping, and counterion-counterion hopping. Detailed structural analysis of the MD trajectories revealed that the chemistry and morphology of the linker groups affect the relative stability and population distribution of the electrolyte components, such that crown-ether-based linker groups enhance the probability of extended, low-barrier ion percolation pathways. As a result, we were able to tune the ionic conductivities by rationally manipulating the counterion distributions, linker binding strengths, and the configurational entropy (multivariability of the linkers). The resulting performance of these MOF-based solid-state electrolytes was significantly enhanced, with a methoxy-functionalized framework (UiO-66-L1 100 ) achieving high ionic conductivities of 2.32 × 10 –4 S/cm and 2.07 × 10 –3 S/cm at 30 °C and 90 °C, respectively, an order of magnitude greater than other all-solid-state MOF electrolyte systems. The electrolyte stability was evaluated with LiIn|LPSCl|MOF:LiTFSI|LPSCl|LiIn symmetric cells, showing excellent Li plating/stripping processes for over 2 months.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

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

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Capturing Ring Opening in Photoexcited Enolic Acetylacetone upon Hydrogen Bond Dissociation by Ultrafast Electron Diffraction

Photoinduced biological and chemical reactions are often based on key structural transformations of a molecule driven across multiple electronic states. Acetylacetone (AcAc) is a prototypical system for complex chemical pathways involving several conical intersections (CI) and singlet–triplet intersystem crossings (ISC) characterized by distinct geometries. In the gas phase, AcAc is predominantly in a planar ring-like enolic form stabilized by a strong intramolecular O–H···O hydrogen bond. Following excitation into the S 2 (ππ*) state at 266 nm, acetylacetone undergoes rapid internal conversion followed by intersystem crossing. Such relaxation pathways are associated with structural changes including ring opening, deplanarization, and bond elongation. In this work, ultrafast electron diffraction (UED) at the SLAC MeV-UED setup is employed as a direct structural probe with a time resolution of 160 fs. Together with trajectory surface hopping simulations, analysis of the UED data provides a new perspective on the early time nuclear dynamics in acetylacetone. Specifically, AcAc is observed to undergo ring opening, deplanarization, and bond elongation all within the first 700 fs after photoexcitation. The monitored dynamics is associated mainly with the nuclear motion on the S 1 potential energy surface, formed after very rapid transfer from S 2 to S 1 , allowing AcAc to reach the conical intersection to intersystem crossing. Such time scales of nuclear motion are contrasted with the time scales of electronic transitions in AcAc that were previously characterized with spectroscopic methods, specifically internal conversion (<100 fs) and intersystem crossing (∼1.5 ps).

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

Evaluation of Amorphous Oxide Coatings for High-Voltage Li-Ion Battery Applications Using a First-Principles Framework

Cathode surface coatings are widely used industrially as a means to suppress degradation and improve electrochemical performance of lithium-ion batteries. However, developing an optimal coating is challenging, as different coating materials may enhance one aspect of performance while hindering another. Furthermore, to elucidate the fundamental thermodynamic and transport properties of amorphous cathode coating materials, here, we present a framework for calculating and analyzing the Li + and O 2- transport and the stability against delithiation in such materials. Our framework includes systematic workflows of ab-initio molecular dynamics calculations to obtain amorphous structures and diffusion trajectories coupled with an analysis of critical changes of the active-ion local environment during diffusion. Based on these data, we provide an estimate of room-temperature diffusivities, including statistical error bars, and the evaluation of the coating suitability in terms of its ability to facilitate Li + transport while blocking O 2- transport. Finally, we add the thermodynamic stability analysis of the coating chemistry within the operating voltage of common Li-ion cathodes. We apply this framework to two commonly used amorphous coating materials, Al 2 O 3 and ZnO. We find that (1) in general, a higher Li + content increases both Li + and O 2- diffusivities in both Al 2 O 3 and ZnO. Also, Li + and O 2- diffuse much faster in ZnO than in Al 2 O 3 . (2) However, neither Al 2 O 3 nor ZnO is expected to retain a significant concentration of Li + at high charge. (3) ZnO performs much more poorly in terms of O 2- blocking, and hence, Al 2 O 3 is preferred for high-voltage cathode applications. These results will help to quantitatively evaluate amorphous materials, such as metal oxides and fluorides, for different performance metrics and facilitate the development of optimal cathode coatings.

25 ENERGY STORAGE↗
DOE OSTI2021

Redox Mediated Control of Electrochemical Potential in Liquid Cell Electron Microscopy

Liquid cell electron microscopy enables the study of nanoscale transformations in solvents with high spatial and temporal resolution, but for the technique to achieve its potential requires a new level of control over the reactivity caused by radical generation under electron beam irradiation. An understanding of how to control electron-solvent interactions is needed to further advance the study of structural dynamics for complex materials at the nanoscale. Additionally, we developed an approach that scavenges radicals with redox species that form well-defined redox couples and control the electrochemical potential in situ . This approach enables the observation of electrochemical structural dynamics at near-atomic resolution with precise control of the liquid environment. Analysis of nanocrystal etching trajectories indicates that this approach can be generalized to several chemical systems. The ability to simultaneously observe heterogeneous reactions at near-atomic resolution and precisely control the electrochemical potential enables the fundamental study of complex nanoscale dynamics with unprecedented detail.

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

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗
DOE OSTI2024

Community Public Mobility Using On-Demand, Low-Speed Electric Vehicles: A Case Study in Downtown St. Louis, Missouri

Legacy fixed route transit systems designed to serve commuters struggle to provide efficient and effective service for short neighborhood trips and for population groups unable to access and egress transit stops using active modes (e.g., elderly, disabled). Neighborhood on-demand transit (ODT) services using low-speed electric vehicles (LSEV) are an innovative technological solution that can help fill this gap in service (e.g., short, high-frequency trips) for diverse populations and trip types. This study evaluated user characteristics and travel behavior for a neighborhood ODT service (using LSEVs) in downtown St. Louis, Missouri using responses from a community survey (n=244), ridership data, and vehicle trajectory information. A comparative analysis between neighborhood ODT, fixed route transit, and transportation network companies (TNC) was also conducted from the perspectives of total travel time, cost, and greenhouse gas emissions. Ultimately, the goal of the analysis was to motivate and inform holistic public mobility systems where different services are optimized to meet specific community needs. Findings indicate that the neighborhood ODT was effective at reaching diverse populations (elderly (20%), lower income (27%), and households with limited access to private vehicles (34%)). ODT reduced total travel time by 32% compared to fixed route transit, produced 2.4 - 4.3 times less greenhouse gas emissions per passenger mile (compared to transit and TNCs), and was more affordable (free to users) than alternative options ($1 for transit, $10-12 for TNCs). Overall satisfaction rates were high, with 80% of respondents rating the service a 4 or 5 out of 5.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Western Interconnection Baseline Study

The purpose of the baseline study is to evaluate the degree to which current industry planning processes meet the national 2035 decarbonization goals for the Western Interconnection. This analysis serves as a comparative baseline for the scenario analysis conducted in the NTP Study using a Western Interconnection dataset that is readily available to industry. This baseline analysis differs from the production cost modeling analysis and power flow analysis in the main NTP Study report (forthcoming). In particular, the analysis presented in this report reflects a business-as-usual future with an optimistic build out of specific planned transmission projects and foreseeable generation. In contrast, the NTP Study models a future generation and transmission expansion based on optimization from a capacity expansion model. The analysis presented herein also reflects a 2030 timeframe, whereas the main NTP Study production cost modeling analysis and power flow analysis reflect a 2035 time frame. This baseline analysis utilizes industry’s most reliable data to account for future transmission projects across various stages of development, with a particular focus on those in the permitting stage. Additionally, it incorporates projections for changes in generation capacity (both additions and retirements). This baseline analysis outlines a probable trajectory, given current process and practice, for the future of the bulk power system with a horizon extending to 2030.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Community Public Mobility Using On-Demand, Low-Speed Electric Vehicles: A Case Study in Downtown St. Louis, Missouri

Legacy fixed route transit systems designed to serve commuters struggle to provide efficient and effective service for short neighborhood trips and for population groups unable to access and egress transit stops using active modes (e.g., elderly, disabled). Neighborhood on-demand transit (ODT) services using low-speed electric vehicles (LSEV) are an innovative technological solution that can help fill this gap in service (e.g., short, high-frequency trips) for diverse populations and trip types. This study evaluated user characteristics and travel behavior for a neighborhood ODT service (using LSEVs) in downtown St. Louis, Missouri using responses from a community survey (n=244), ridership data, and vehicle trajectory information. A comparative analysis between neighborhood ODT, fixed route transit, and transportation network companies (TNC) was also conducted from the perspectives of total travel time, cost, and greenhouse gas emissions. Ultimately, the goal of the analysis was to motivate and inform holistic public mobility systems where different services are optimized to meet specific community needs. Findings indicate that the neighborhood ODT was effective at reaching diverse populations (elderly (20%), lower income (27%), and households with limited access to private vehicles (34%)). ODT reduced total travel time by 32% compared to fixed route transit, produced 2.4 - 4.3 times less greenhouse gas emissions per passenger mile (compared to transit and TNCs), and was more affordable (free to users) than alternative options ($1 for transit, $10-12 for TNCs). Overall satisfaction rates were high, with 80% of respondents rating the service a 4 or 5 out of 5.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Asc-Seurat: analytical single-cell Seurat-based web application

Abstract Background Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of transcriptomes, arising as a powerful tool for discovering and characterizing cell types and their developmental trajectories. However, scRNA-seq analysis is complex, requiring a continuous, iterative process to refine the data and uncover relevant biological information. A diversity of tools has been developed to address the multiple aspects of scRNA-seq data analysis. However, an easy-to-use web application capable of conducting all critical steps of scRNA-seq data analysis is still lacking. Summary We present Asc-Seurat, a feature-rich workbench, providing an user-friendly and easy-to-install web application encapsulating tools for an all-encompassing and fluid scRNA-seq data analysis. Asc-Seurat implements functions from the Seurat package for quality control, clustering, and genes differential expression. In addition, Asc-Seurat provides a pseudotime module containing dozens of models for the trajectory inference and a functional annotation module that allows recovering gene annotation and detecting gene ontology enriched terms. We showcase Asc-Seurat’s capabilities by analyzing a peripheral blood mononuclear cell dataset. Conclusions Asc-Seurat is a comprehensive workbench providing an accessible graphical interface for scRNA-seq analysis by biologists. Asc-Seurat significantly reduces the time and effort required to analyze and interpret the information in scRNA-seq datasets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Material coherence from trajectories via Burau eigenanalysis of braids

In this paper, we provide a numerical tool to study a material’s coherence from a set of 2D Lagrangian trajectories sampling a dynamical system, i.e., from the motion of passive tracers. We show that eigenvectors of the Burau representation of a topological braid derived from the trajectories have levelsets corresponding to components of the Nielsen–Thurston decomposition of the dynamical system. One can thus detect and identify clusters of space–time trajectories corresponding to coherent regions of the dynamical system by solving an eigenvalue problem. Unlike previous methods, the scalable computational complexity of our braid-based approach allows the analysis of large amounts of trajectories.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗