Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exact-Factorization-Based Surface Hopping without Velocity Adjustment

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. Furthermore, the resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

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

Deciphering Charge Transfer Processes in Transition Metal Complexes from the Perspective of Ultrafast Electronic and Nuclear Motions

Chemical transformations in charge transfer states result from the interplay between electronic dynamics and nuclear reorganization along excited-state trajectories. Here, in this study, we investigate the ultrafast structural dynamics following photoinduced electron transfer from the metal-metal-to-ligand charge transfer state of an electron donor, a Pt dimer complex, to a covalently linked electron acceptor group using ultrafast time-resolved wide-angle X-ray scattering and optical transient absorption spectroscopy methods to disentangle the interdependence of the excited-state electronic and nuclear dynamics. Following photoexcitation, Pt-Pt bond formation and contraction takes up to 1 ps, much slower than the corresponding process in analogous complexes without electron acceptor groups. Because the Pt-Pt distance change is slow with respect to excited-state electron transfer, it can affect the rate of electron transfer. These results have potential impacts on controlling electron transfer rates via structural alterations to the electron donor group, tuning the charge transfer driving force.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2025

Spin-Vibronic Effect in Photoinduced Electron Transfer

The spin-vibronic effect (SVE) accelerates quantum-mechanically forbidden electronic transitions, but experimental manifestations of this phenomenon remain limited. In this contribution, the role of SVE in photoinduced electron transfer (PET) dynamics is probed by coherent vibrational wavepacket (CVWP) motions in Pt­(II) dimer-naphthalene diimide donor–acceptor complexes. Metal–metal-to-ligand charge-transfer (MMLCT) excitation of the donor triggers ballistic PET, driving the molecule away from equilibrium. The SVE then directs the PET trajectory from the 1 MMLCT excited state to an intermediate ligand-centered triplet state, modulating the CVWP dynamics along the reaction coordinates. Herzberg–Teller type oscillations indicate that the Pt–Pt stretching vibration, arising from the formation of the excited triplet intermediate, serves as the reaction coordinate for charge separation. Finally, our experimental findings advance the understanding of the interplay between electronic, vibrational, and spin degrees of freedom in ultrafast photoactivated processes.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

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

Desolvation Processes in Channel Solvates of Niclosamide

The antiparasitic drug niclosamide (NCL) is notable for its ability to crystallize in multiple 1:1 channel solvate forms, none of which are isostructural. Here, using a combination of time-resolved synchrotron powder X-ray diffraction and thermogravimetry, the process-induced desolvation mechanisms of methanol and acetonitrile solvates are investigated. Structural changes in both solvates follow a complicated molecular-level trajectory characterized by a sudden shift in lattice parameters several degrees below the temperature where the desolvated phase first appears. Model fitting of kinetic data obtained under isothermal heating conditions suggests that the desolvation is rate-limited by the nucleation of the solvent-free product. The desolvation pathways identified in these systems stand in contrast to previous investigations of the NCL channel hydrate, where water loss by diffusion initially yields an anhydrous isomorph that converts to the thermodynamic polymorph at significantly higher temperatures. Taking the view that each solvate lattice is a unique “pre-organized” precursor, a comparison of the pathways from different starting topologies to the same final product provides the opportunity to reevaluate assumptions of how various factors (e.g., solvent binding strength, density) influence solid-state desolvation processes.

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

COVID-19: Effects of Environmental Conditions on the Propagation of Respiratory Droplets

As coronavirus disease 2019 (COVID-19) continues to spread, a detailed understanding on the transmission mechanisms is of paramount importance. The disease transmits mainly through respiratory droplets and aerosol. Although models for the evaporation and trajectory of respiratory droplets have been developed, how the environment impacts the transmission of COVID-19 is still unclear. In this study, we investigate the propagation of respiratory droplets and aerosol particles generated by speech under a wide range of temperatures (0–40 °C) and relative humidity (0–92%) conditions. We show that droplets can travel three times farther in low-temperature and high-humidity environment, whereas the number of aerosol particles increases in high-temperature and low-humidity environments. The results also underscore the importance of proper ventilation, as droplets and aerosol spread significantly farther in airstreams. This study contributes to the understanding of the environmental impact on COVID-19 transmission.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Online Bayesian State Estimation for Real-Time Monitoring of Growth Kinetics in Thin Film Synthesis

Rapid validation of newly predicted materials through autonomous synthesis requires real-time adaptive control methods that exploit physics knowledge, a capability that is lacking in most systems. Here, in this study, we demonstrate an approach to enable real-time control of thin film synthesis by combining in situ optical diagnostics with a Bayesian state estimation method. We developed a physical model for film growth and applied the direct filter (DF) method for real-time estimation of nucleation and growth rates during pulsed laser deposition (PLD). We validated the approach using simulated and experimental reflectivity data for WSe 2 growth and ultimately deployed the algorithm on an autonomous PLD system during the growth of 1T'-MoTe 2 . The DF robustly estimates growth parameters in real time at early stages of growth, down to 15% monolayer area coverage. This fusion of in situ diagnostics, data assimilation, and physical modeling opens new opportunities in adaptive control of synthesis trajectories toward desired material states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Li-Ion Diffusion Correlations in LiAlGeO 4 : Quasielastic Neutron Scattering and Ab Initio Simulation

Here, we investigated the impact of Li stoichiometry and host flexibility on Li + diffusion processes in LiAlGeO 4 at the microscopic level using quasielastic neutron scattering (QENS) and ab initio molecular dynamics (AIMD) simulations. Using sufficiently long AIMD trajectories, we could simulate the observed QENS signal and identify the localized dynamics of Li in crystalline LiAlGeO 4 . Such information is vital to identify the bottleneck of diffusion processes and design materials for battery application. Our AIMD simulations in LiAlGeO 4 reveal that the Li + conductivity can be significantly improved by manipulating the Li stoichiometry and/or host flexibility via amorphization. We determined that excess Li stoichiometry enhances the Coulomb repulsion of neighboring Li sites and softens the host structure to enable faster Li + diffusion along the hexagonal c-axis. In the amorphous structure, random orientations of AlO 4 and GeO 4 polyhedral units create a wide distribution of intersite distances and significantly soften the host structure, greatly enhancing the Li + diffusion. The simulations are used to understand the nature of diffusion, especially the role of the host structure, the possible hopping pathways, and the diffusion behavior in various structures of LiAlGeO 4 .

localized dynamics↗
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

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

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

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

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

In Situ Atmospheric Plume Thermometry via Carbon Monoxide Spectral Profile: Laboratory and Field Validation

For small molecules with large rotational constants, knowing the relative intensities of the ro-vibrational transitions can be used to determine the temperature within a gas plume. We demonstrate the use of carbon monoxide (CO) as an in situ spectroscopic probe of gas plume temperature by application of both laboratory and standoff Fourier transform infrared spectroscopy to monitor the CO spectral response at different temperatures. Here, the measured CO rotational contours were analyzed using a simple Boltzmann model to deduce the population distribution of the J-levels, from which the in-plume temperature is deduced. The method was vetted by comparing deduced temperatures in both static laboratory measurements of known temperatures, as well as field measurements using a simulated smokestack release. For the smokestack experiments, spectroscopically deduced temperatures were compared to readings from a series of thermocouples placed at strategically sampled distances along the plume trajectory. Both the spectroscopically-derived and thermocouple-measured temperatures revealed an expansion-induced (mixing) rapid cooling of the plume, with the infrared thermometry values displaying greater temperature values which are believed to better represent the actual plume temperatures.

Analysis of rotational structure↗
DOE OSTI2024

Chemical Insights into the Molecular Composition of Organic Aerosols in the Urban Region of Houston, Texas

Molecular functional groups, such as organosulfates (CHOS) and organonitrates (CHNO) are important tracers for field observations of secondary organic aerosols (SOA). While CHOS and CHNO are prevalent in the atmosphere, there is a lack of knowledge regarding daily and day- and night-time variations in these species in the urban atmosphere. Meteorological factors such as wind speed/direction, relative humidity (RH), and temperature can influence the formation of CHOS/CHNO. To investigate these trends, we utilized multimodal chemical imaging and advanced high resolution mass spectrometry techniques to acquire particle speciation and molecular formulas (MFs) associated with day and night sampling periods. Back trajectory analyses revealed the oceanic influence of southern wind airmasses in later June sampling periods with organic fractions <10%. Conversely, northern winds in early June sampling periods contributed to the episodic emergence of extremely low volatile organics (ELVOCs) and organic factions up to 41%. The observed unique MFs to June 3 (223 MFs) and to June 4 (144 MFs) were largely found to be of biogenic rather than anthropogenic origin. Finally, our findings reveal episodic prevalence and temporal distribution of SOA constituents across the urban region of Houston, Texas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Photochemical and Cloud and Aerosol Aqueous Contributions to Regionally-Emitted Shipping and Biogenic Non-Sea-Salt Sulfate Aerosol in Coastal California

Aerosol nonsea-salt sulfate (NSS sulfate) forms in the atmosphere by secondary reactions of emissions from marine phytoplankton and shipping, with gas-phase as well as cloud and aerosol aqueous reactions controlling production. Twelve months of Atmospheric Radiation Measurements (ARM) during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) at Scripps Pier in La Jolla, California, showed the highest NSS sulfate mass concentrations occurred for the northwesterly back-trajectories over 64% of the year, with an average of 0.90 μg/m 3 that contributed 76% of annual NSS sulfate concentration. Multiple Linear Regression (MLR) and a refractory black carbon tracer method attributed 76–80% of the regionally emitted sulfur dioxide (SO 2 ) sources of submicron NSS sulfate to marine biogenic emissions and 20–24% to shipping emissions. MLR for oxidation processes explained 21% of the variability with Downwelling Shortwave Radiation (DSW) driving photochemical reactions to account for 34% of annual regional sulfate production, Upwind Cloud Vertical Fraction (UCVF) controlling cloud-associated oxidation to account for 29%, and relative humidity (RH) describing aerosol-phase oxidation to account for 36%. NSS sulfate was correlated moderately to UCVF during April-June and August but to RH in October-January. These findings show the apportionment of SO 2 emissions to biogenic and shipping sources and provide observational constraints for the mechanisms for sulfate production from SO 2 in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Single-Particle Studies Reveal a Nanoscale Mechanism for Elastic, Bright, and Repeatable ZnS:Mn Mechanoluminescence in a Low-Pressure Regime

Mechanoluminescent materials, which emit light in response to elastic deformation, are demanded for use as in situ stress sensors. ZnS doped with Mn is known to exhibit one of the lowest reported thresholds for appearance of mechanoluminescence, with repeatable light emission under contact pressure <10 MPa. Here, the physical basis for such behavior remains as yet unclear. Here, reliable microscopic detection of mechanoluminescence of single ZnS:Mn microparticles, in combination with nanoscale structural characterization, provides evidence that the mechanoluminescent properties of these particles result from interplay between a non-centrosymmetric crystal lattice and its defects, viz., dislocations and stacking faults. Statistical analysis of the distributions of mechanoluminescence energy release trajectories reveals two distinct mechanisms of excitation: one attributable to a piezo-phototronic effect and the other due to dislocation motion. At pressures below 8.1 MPa, both mechanisms contribute to mechanoluminescent output, with a dominant contribution from the piezo-phototronic mechanism. In contrast, above 8.1 MPa, dislocation motion is the primary excitation source. For the piezo-phototronic mechanism, we propose a specific model that accounts for elastic ZnS:Mn mechanoluminescence under very low pressure. The charged interfaces in stacking faults lead to the presence of filled traps, which otherwise would be empty in the absence of the built-in electric field. Upon application of external stress, local enhancement of the piezoelectric field at the stacking faults’ interfaces facilitates release of the trapped carriers and subsequent luminescence. This field enhancement explains how <10 MPa pressure produces thousands of photons.

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

Disentangling Rotational Dynamics and Ordering Transitions in a System of Self-Organizing Protein Nanorods via Rotationally Invariant Latent Representations

The dynamics of complex ordering systems with active rotational degrees of freedom exemplified by protein self-assembly is explored using a machine learning workflow that combines deep learning-based semantic segmentation and rotationally invariant variational autoencoder-based analysis of orientation and shape evolution. The latter allows for disentanglement of the particle orientation from other degrees of freedom and compensates for lateral shifts. The disentangled representations in the latent space encode the rich spectrum of local transitions that can now be visualized and explored via continuous variables. The time dependence of ensemble averages allows insight into the time dynamics of the system and, in particular, illustrates the presence of the potential ordering transition. Finally, analysis of the latent variables along the single-particle trajectory allows tracing these parameters on a single-particle level. The proposed approach is expected to be universally applicable for the description of the imaging data in optical, scanning probe, and electron microscopy seeking to understand the dynamics of complex systems where rotations are a significant part of the process.

representation learning↗