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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Seasonality of Formic Acid (HCOOH) in London during the ClearfLo Campaign

Abstract Following measurements in the winter of 2012, formic acid (HCOOH) and nitric acid (HNO 3 ) were measured using a chemical ionization mass spectrometer (CIMS) during the Summer Clean Air for London (ClearfLo) campaign in London, 2012. Consequently, the seasonal dependence of formic acid sources could be better understood. A mean formic acid concentration of 1.3 ppb and a maximum of 12.7 ppb was measured which is significantly greater than that measured during the winter campaign (0.63 ppb and 6.7 ppb, respectively). Daily calibrations of formic acid during the summer campaign gave sensitivities of 1.2 ion counts s −1 parts per trillion (ppt) by volume −1 and a limit of detection of 34 ppt. During the summer campaign, there was no correlation between formic acid and anthropogenic emissions such as NO x and CO or peaks associated with the rush hour as was identified in the winter. Rather, peaks in formic acid were observed that correlated with solar irradiance. Analysis using a photochemical trajectory model has been conducted to determine the source of this formic acid. The contribution of formic acid formation through ozonolysis of alkenes is important but the secondary production from biogenic VOCs could be the most dominant source of formic acid at this measurement site during the summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

NASA analysis of space mission options for the 2025 planetary defense conference hypothetical asteroid impact threat scenario

The 2025 Planetary Defense Conference (PDC) hypothetical asteroid impact threat exercise is being conducted in coordination with the United Nations-endorsed Space Mission Planning Advisory Group (SMPAG), exercising SMPAG's processes for assessing space mission options and communicating to decision makers. Here, in this paper, we describe the work performed by our NASA-led team and present the results we contributed to the SMPAG effort for the exercise. Our NASA-led team assessed mission options for asteroid reconnaissance (flyby and rendezvous), deflection, and robust disruption using several viable techniques: kinetic impactors, ion beam deflection, and nuclear explosive devices. Our simulations and analyses considered how much change-in-velocity an asteroid can tolerate before fragmentation onset, and we worked towards establishing requirements for robustly disrupting an asteroid. Heuristics informed by the simulation results are incorporated into optimization of deflection and disruption mission campaign options. Finally, we make some observations about useful generalizations from these results with potential applicability to any planetary defense scenario.

Asteroid deflection↗
DOE OSTI2022

Uncovering the role of nanoscale precipitates on martensitic transformation and superelasticity

Here we characterize the role of coherent nanoscale B2 Ni 50 Al 50 precipitates on the temperature- and stress-induced martensitic phase transformation in nanocrystalline Ni 63 Al 37 shape memory alloys using multi-million-atoms molecular dynamics (MD) simulations. We studied two types of precipitates: one with single crystal precipitates (SXP) and a second where grain boundaries cut through precipitates (PXP). Simulations reveal that the presence of B2 precipitates stabilizes the cyclic flag-shaped stress-strain response, characteristic of superelasticity, and reduces remnant strain. In contrast, single-phase nanocrystalline Ni 63 Al 37 exhibits degradation of the reverse transformation during cyclic loading and, eventually, incomplete reversible transformation within a few cycles. This is consistent with previous experimental findings of ultra-low fatigue in Ni-Ti-Cu alloys with Ti 2 Cu precipitates. The simulations reveal that the presence of precipitates significantly improves the reversibility of the transformation by acting as elastic zones that partially shield the martensitic transformation and drive the reverse transformation. A detailed analysis of the MD trajectories reveals that the martensitic transformation of the matrix induces ultra-large elastic deformation in some of the B2 precipitates (approximately 12%) to the point of resulting in a martensite-like atomic structure.

36 MATERIALS SCIENCE↗
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 OSTI2018

Adsorption mechanisms of cesium at calcium-silicate-hydrate surfaces using molecular dynamics simulations

We report understanding the mechanisms of interactions and immobilization of radionuclides in cement systems at the molecular level is key to developing safe and novel solutions for the long-term storage of nuclear wastes. The molecular level adsorption mechanisms of cesium (Cs + ) ions onto tobermorite 9 Å, tobermorite 14 Å, and jennite – three crystalline structural analogues for calcium-silicate-hydrates, the main components in cement-based materials, were investigated using molecular dynamics (MD) simulations. Convergence monitoring of the simulations, using the root mean square displacement with average correlation analysis, indicated that MD trajectories of 10–15 ns were needed given the size of the selected computational cell to properly capture the dynamic process of adsorption. Cs + ions adsorbed as inner-sphere complexes at the tobermorite surfaces, while they showed lower affinity for jennite. The charge and structure of the surface layer strongly influenced the adsorption mechanisms of Cs + ions. Strong association with the hexagonal cavities of the silica sites were seen at the tetrahedral SiO 4 surface of the tobermorite structures. Cs + ion adsorption was attributed to co-ion adsorption at the octahedral CaO 6 surface of tobermorite 9 Å, cation exchange at that of tobermorite 14 Å, and electrostatic interaction with the surface hydroxyl groups for jennite. Distribution coefficients of Cs + ions onto cement pastes derived from molecular dynamics by taking the surface structure and dynamics of Cs + ions directly into account showed a range of overlap with experimental results.

42 ENGINEERING↗
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 OSTI2018

Machine-Learning Enabled New Insights into the Coil-to-Globule Transition of Thermosensitive Polymers Using a Coarse-Grained Model

In this paper, we present a computational framework that integrates coarse-grained (CG) molecular dynamics (MD) simulations and a data-driven machine-learning (ML) method to gain insights into the conformations of polymers in solutions. We employ this framework to study conformational transition of a model thermosensitive polymer, poly(N-isopropylacrylamide) (PNIPAM). Here, we have developed the first of its kind, a temperature-independent CG model of PNIPAM that can accurately predict its experimental lower critical solution temperature (LCST) while retaining the tacticity in the presence of an explicit water model. The CG model was extensively validated by performing CG MD simulations with different initial conformations, varying the radius of gyration of chain, the chain length, and the angle between the adjacent monomers of the initial configuration of PNIPAM (total simulation time = 90 μs). Moreover, for the first time, we utilize the nonmetric multidimensional scaling (NMDS) method, a data-driven ML approach, to gain further insights into the mechanisms and pathways of this coil-to-globule transition by analyzing CG MD simulation trajectories. In conclusion, NMDS analysis provides entirely new insights and shows multiple metastable states of PNIPAM during its coil-to-globule transition above the LCST.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Nanoscale Chemical Reaction Imaging at the Solid-Liquid Interface via TERS

Plasmon-enhanced chemical transformations at the solid-liquid interface can be imaged with high sensitivity, chemical selectivity, and nanoscale precision through tip-enhanced Raman scattering (TERS). We demonstrate the latter for the first time through measurements aimed at (i) locating plasmonic hotspots at the solid-liquid interface at which chemical transformations take place, (ii) monitoring the evolution from reactants to products through their distinct Raman spectra, and (iii) 2D correlation analysis of Raman time trajectories to unambiguously extract the spectral components that mark chemical transformation and to understand the correlations between the product and parent signatures. For our proof-of-principle study, we select a model plasmon-enhanced chemical process, namely, the dimerization of p-nitrothiophenol to dimercaptoazobenzene, but now, at the solid-liquid interface. Our plasmonic construct otherwise consists of chemically functionalized gold microplates in aqueous solution, which we image using a gold-coated TERS probe irradiated at 633-nm. Overall, we demonstrate chemical reaction imaging in aqueous solution via TERS for the first time, herein, at a pixel-limited lateral spatial resolution of 10 nm.

Bhattarai, Ashish↗
DOE OSTI2017

Comparative Study of the Effect of Defects on Selective Adsorption of Butanol from Butanol/Water Binary Vapor Mixtures in Silicalite-1 Films

A promising route for sustainable 1-butanol (butanol) production is ABE (acetone, butanol, ethanol) fermentation. However, recovery of the products is challenging because of the low concentrations obtained in the aqueous solution, thus hampering large-scale production of biobutanol. Membrane and adsorbent-based technologies using hydrophobic zeolites are interesting alternatives to traditional separation techniques (e.g., distillation) for energy-efficient separation of butanol from aqueous mixtures. To maximize the butanol over water selectivity of the material, it is important to reduce the number of hydrophilic adsorption sites. This can, for instance, be achieved by reducing the density of lattice defect sites where polar silanol groups are found. The density of silanol defects can be reduced by preparing the zeolite at neutral pH instead of using traditional synthesis solutions with high pH. In this work, binary adsorption of butanol and water in two silicalite-1 films was studied using in situ attenuated total reflectance–Fourier transform infrared (ATR-FTIR) spectroscopy under equal experimental conditions. One of the films was prepared in fluoride medium, whereas the other one was prepared at high pH using traditional synthesis conditions. The amounts of water and butanol adsorbed from binary vapor mixtures of varying composition were determined at 35 and 50 °C, and the corresponding adsorption selectivities were also obtained. Both samples showed very high selectivities (100–23 000) toward butanol under the conditions studied. The sample having low density of defects, in general, showed ca. a factor 10 times higher butanol selectivity than the sample having a higher density of defects at the same experimental conditions. This difference was due to a much lower adsorption of water in the sample with low density of internal defects. Analysis of molecular simulation trajectories provides insights on the local selectivities in the zeolite channel network and at the film surface.

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

Rotation and Negative Torque in Electrodynamically Bound Nanoparticle Dimers

Here, we examine the formation and concomitant rotation of electrodynamically bound dimers (EBD) of 150nm diameter Ag nanoparticles trapped in circularly polarized focused Gaussian beams. The rotation frequency of an EBD increases linearly with the incident beam power, reaching high mean values of ~ 4kHz for a relatively low incident power of 14mW. Using a coupled-dipole/effective polarizability model, we reveal that retardation of the scattered fields and electrodynamic interactions can lead to a “negative torque” causing rotation of the EBD in the direction opposite to that of the circular polarization. This intriguing opposite-handed rotation due to negative torque is clearly demonstrated using electrodynamics-Langevin dynamics simulations by changing particle separations and thus varying the retardation effects. Finally, negative torque is also demonstrated in experiments from statistical analysis of the EBD trajectories. These results demonstrate novel rotational dynamics of nanoparticles in optical matter using circular polarization and open a new avenue to control orientational dynamics through coupling to interparticle separation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2015

Stratification of discharge in noble gases from the viewpoint of the discrete dynamics

Based on the analysis of electron phase trajectories in sinusoidal electric fields, a new point of view on discharge stratification is proposed. It is shown that the positive column can be considered as a spatial resonator in which electric fields with a fundamental period length L{sub S} or higher mode length q/p L{sub S} establish, where p and q are integers and p > q. The fundamental mode length L{sub S} is equivalent to the distance on which electrons gain energy equal to the lowest excitation threshold. This distance determines a length of the S-striation. Unlike kinetic theory, in the presented model resonance properties of the discharge column are not connected with elastic collision energy losses. A point map is used to obtain the resonance trajectories of electrons in the phase plane. Stable points for the positions of inelastic collisions in the resonance trajectories have been found at the positions of field maxima in the case of integer ratios p/q . For non-integer ratios p/q , multiple resonance trajectories arise according to a more complex stability criterion. From this point of view, S-, P-, and R-striations in noble gas discharges can be explained. Due to energy losses in elastic collisions, initial electron energy distribution functions converge to the resonance trajectories (the so-called “bunch effect”). The findings of the discrete model agree with results of kinetic theory and experiment. The new approach avoids difficulties of the kinetic theory in the case of exceptionally large relaxation lengths which can even exceed the positive column length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗