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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin–orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3 (A) leading to the NH(X 3 Σ – , a 1 Δ) + H 2 channel

Here we simulate the photodissociation of NH 3 originating from its first excited singlet state S 1 into the NH 2 + H (radical) and NH + H 2 (molecular) channels. The states considered are the ground singlet state S 0 , the first excited singlet state S 1 and the lowest-lying triplet state T 1 , which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S 0 , S 1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T 1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H 2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH 2 + H channel. The NH(X 3 Σ_) radical produced through the intersystem crossing from S 0 to T 1 is rare (~0.2%) compared to NH(a 1 Δ) due to the process being spin forbidden.

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

Stable excited dication: trapping on the S 1 state of formaldehyde dication after strong field ionization

Combined theoretical and experimental work examines the dynamics of dication formaldehyde produced by strong field ionization. Trajectory surface hopping dynamics on the first several singlet electronic states of the formaldehyde dication are used to examine the relaxation pathways and dissociation channels, while kinetic energy distributions after strong field ionization of formaldehyde and deuterated formaldehyde are used to confirm the theoretical predictions. Here, we find that the first excited state of the formaldehyde dication is stable, neither decays to the ground state nor dissociates, even though the ground state and higher lying states are directly dissociative. The stability of the first excited state is explained by its symmetry which does not allow for radiative or nonradiative transitions to the ground state and by large barriers to dissociate on the excited state surface.

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

Predicting third-body collision efficiencies for water and other polyatomic baths

Low-pressure-limit microcanonical (collisional activation) and thermal rate constants are predicted using a combination of automated ab initio potential energy surface construction, classical trajectories, transition state theory, and a detailed energy- and angular-momentum-resolved collision kernel. Several systems are considered, including CH 4 (+M) and HO 2 (+M), with an emphasis on systems where experimental information is available for comparison. The a priori approach involves no adjustable parameters, and we show that the predicted thermal rate constants are in excellent agreement with experiments, with average deviations of less than 25%. Notably, the a priori approach is shown to perform equally well for atomic, diatomic, and polyatomic baths, including M = H 2 O, CO 2 , and “fuel” baths like M = CH 4 and NH 3 . Additionally, the utility of microcanonical rate constants for interpreting trends and inferring mechanistic details in the thermal kinetics is demonstrated.

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

On the use of vibronic coherence to identify reaction coordinates for ultrafast excited-state dynamics of transition metal-based chromophores

In this study, the question of whether one can use information from quantum coherence as a means of identifying vibrational degrees of freedom that are active along an excited-state reaction coordinate is discussed. Specifically, we are exploring the notion of whether quantum oscillations observed in single-wavelength kinetics data exhibiting coherence dephasing times that are intermediate between that expected for either pure electronic or pure vibrational dephasing are vibronic in nature and therefore may be coupled to electronic state-to-state evolution. In the case of a previously published Fe(ιι) polypyridyl complex, coherences observed subsequent to 1 A 1 → 1 MLCT excitation were linked to large-amplitude motion of a portion of the ligand framework; dephasing times on the order of 200–300 fs suggested that these degrees of freedom could be associated with ultrafast (~100 fs) conversion from the initially formed MLCT excited state to lower-energy, metal-centered ligand-field excited state(s) of the compound. Incorporation of an electronically benign but sterically restrictive Cu(ι) ion into the superstructure designed to interfere with this motion yielded a compound exhibiting a ~25-fold increase in the compound’s MLCT lifetime, a result that was interpreted as confirmation of the initial hypothesis. However, new data acquired on a different chemical system – Cr(acac') 3 (where acac' represents various derivatives of acetylacetonate) – yielded results that call into question this same hypothesis. Coherences observed subsequent to 4 A 2 → 4 T 2 ligand-field excitation on a series of molecules implicated similar vibrational degrees of freedom across the series, but exhibited dephasing times ranging from 340 fs to 2.5 ps without any clear correlation to the dynamics of excited-state evolution in the system. Taken together, the results obtained on both of these chemical platforms suggest that while identification of coherences can indeed point to degrees of freedom that should be considered as candidate modes for defining reaction trajectories, our understanding of the factors that determine the interplay across coherences, dephasing times, and electronic and geometric structure is insufficient at the present time to view this parameter as a robust metric for differentiating active versus spectator modes for ultrafast dynamics.

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

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗
DOE OSTI2022

Deep learning potential of mean force between polymer grafted nanoparticles

Grafting polymer chains on the surfaces of nanoparticles is a well-known route to control their self-assembly and distribution in a polymer matrix. A wide variety of self-assembled structures are achieved by changing the grafting patterns on the surface of an individual nanoparticle. However, an accurate estimation of the effective potential of mean force between a pair of grafted nanoparticles that determines their assembly and distribution in a polymer matrix is an outstanding challenge in nanoscience. We address this problem via deep learning. As a proof of concept, here we report a deep learning framework that learns the interaction between a pair of single-chain grafted spherical nanoparticles from their molecular dynamics trajectory. Subsequently, we carry out the deep learning potential of mean force-based molecular simulation that predicts the self-assembly of a large number of single-chain grafted nanoparticles into various anisotropic superstructures, including percolating networks and bilayers depending on the nanoparticle concentration in three-dimensions. The deep learning potential of mean force-predicted self-assembled superstructures are consistent with the actual superstructures of single-chain polymer grafted spherical nanoparticles. This deep learning framework is very generic and extensible to more complex systems including multiple-chain grafted nanoparticles. Finally, we expect that this deep learning approach will accelerate the characterization and prediction of the self-assembly and phase behaviour of polymer-grafted and unfunctionalized nanoparticles in free space or a polymer matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mixed quantum/classical theory for rotational energy exchange in symmetric-top-rotor + linear-rotor collisions and a case study of the ND 3 + D 2 system

The extension of mixed quantum/classical theory (MQCT) to describe collisional energy transfer is developed for a symmetric-top-rotor + linear-rotor system and is applied to ND 3 + D 2 . State-to-state transition cross sections are computed in a broad energy range for all possible processes: when both ND 3 and D 2 molecules are excited or both are quenched, when one is excited while the other is quenched and vice versa, when the ND 3 state changes its parity while D 2 is excited or quenched, and when ND 3 is excited or quenched while D 2 remains in the same state, ground or excited. In all these processes the results of MQCT are found to approximately satisfy the principle of microscopic reversibility. For a set of sixteen state-to-state transitions available from the literature for a collision energy of 800 cm –1 the values of cross sections predicted by MQCT are within 8% of accurate full-quantum results. A useful time-dependent insight is obtained by monitoring the evolution of state populations along MQCT trajectories. It is shown that, if before the collision, D 2 is in its ground state, the excitation of ND 3 rotational states proceeds through a two-step mechanism in which the kinetic energy of molecule–molecule collision is first used to excite D 2 and only then is transferred to the excited rotational states of ND 3 . Finally, it is found that both potential coupling and Coriolis coupling play important roles in ND 3 + D 2 collisions.

inelastic cross-sections↗
DOE OSTI2023

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

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

Ultrafast structural dynamics of UV photoexcited cis , cis -1,3-cyclooctadiene observed with time-resolved electron diffraction

Conjugated diene molecules are highly reactive upon photoexcitation and can relax through multiple reaction channels that depend on the position of the double bonds and the degree of molecular rigidity. Understanding the photoinduced dynamics of these molecules is crucial for establishing general rules governing the relaxation and product formation. Here, in this study, we investigate the femtosecond time-resolved photoinduced excited-state structural dynamics of cis,cis-1,3-cyclooctadiene, a large-flexible cyclic conjugated diene molecule, upon excitation with 200 nm using mega-electron-volt ultrafast electron diffraction and trajectory surface hopping dynamics simulations. We tracked the photoinduced structural changes from the Franck–Condon region through the conical intersection seam to the ground state. Our findings revealed a novel primary reaction coordinate involving ring distortion, where the ring stretches along one axis and compresses along the perpendicular axis. The nuclear wavepacket remains compact along this reaction coordinate until it reaches the conical intersection seam, and it rapidly spreads as it approaches the ground state, where multiple products are formed.

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

Challenges and strategies for probing the composite interface of PEM electrolyzers and fuel cells using operando AP-XPS

Left: cross-section schematic of a membrane electrode assembly, the working electrode changing state with applied potential. Center: the operando cell design that enables snapshot data acquisition during trajectory movement. Right: resulting spectra. Understanding the surface chemistry of electrocatalyst systems under operando conditions is central to revealing the electrocatalytic cell's working mechanisms. Determination of these catalytic processes on a molecular scale and the involved components is fundamental to streamlining material design for energy conversion and storage applications. X-ray photoelectron spectroscopy (XPS) is an established technique used to study the chemical and electronic states of materials. While the surface sensitivity of XPS is typically high, use of tender X-ray energies and technical advancements have allowed for the direct probing of solid–vapor and solid–liquid interfaces. However, protocols and documentation of experimental considerations for operando XPS probing of working electrolyzers and fuel cells remain scarce. Herein, we report an approach for the study of working polymer electrolyte membrane (PEM) electrolysis cells using ambient pressure X-ray photoelectron spectroscopy (AP-XPS). This approach directly probes the composite electrode surface on the membrane electrode assembly (MEA) in 100% relative humidity to establish a meaningful liquid layer for electrocatalysis. We carry out a systematic investigation from the cell constituent components to a fully assembled working operando electrolytic system and establish a method for AP-XPS study of the complex composite MEA, providing recommendations for data acquisition and component analysis.

Hamlyn, Rebecca↗
DOE OSTI2024

Influence of crossover on capacity fade of symmetric redox flow cells

Volumetrically unbalanced compositionally symmetric cell cycling with potentiostatic (CV) or galvanostatic-with-potential-hold (CCCV) protocols is a rigorous technique for evaluating the calendar lifetime of reactants for redox flow batteries. Here, we evaluate the influence of reactant crossover through the membrane on symmetric cell cycling behavior. We tested symmetric cells of anthraquinone disulfonic acid (AQDS) with Nafion membranes of varied thickness and manufacture (NR211, NR212, N115, and N117, ranging 25–183 μm). Membranes were tested both as-received and pretreated with a common procedure of soaking in water at elevated temperature and then in dilute hydrogen peroxide. We found no significant difference in capacity fade rates of symmetric cells with any of the membranes as-received, indicating a negligible influence of crossover. However, we observed increased capacity fade with increased permeability through pretreated membranes. Supported by zero-dimensional modeling and operando UV-vis spectrophotometry, we propose a mechanism for net crossover in AQDS symmetric cells based on a higher time-averaged concentration of quinhydrone dimers in the non-capacity limiting side (NCLS) compared to the capacity limiting side (CLS), driving net crossover of AQDS reactants out of the CLS. Further, we illustrate other hypothetical scenarios of net crossover using the zero-dimensional model. Overall, many membrane–electrolyte systems used in symmetric cell studies have sufficiently low crossover flux as to avoid the influence of crossover on capacity fade, but under conditions of higher crossover flux, complex interactions of crossover and chemical reactions may result in diverse capacity fade trajectories, the mechanisms of which may be untangled with operando characterization and modeling.

25 ENERGY STORAGE↗
DOE OSTI2025

Real-space Kohn–Sham density functional theory for complex energy applications

Real-space Kohn-Sham density functional theory (real-space KS-DFT) enables large-scale electronic structure simulations that is particularly well-suited for the modern high-performance computing (HPC) architectures. This feature article reviews its theoretical foundations, highlights the algorithmic advances and recent developments, and showcases applications in complex nano systems. We aim to provide a perspective on the trajectory of real-space KS-DFT as an emerging tool for computational chemistry and materials science in the exascale era.

Zhang, Zeyi↗
DOE OSTI2025

Active learning path-dependent properties using a cloud-based materials acceleration platform

Solid state materials are central to many modern technologies in which a given material may be exposed to a variety of environments. The material properties often vary with the sequence of environments in an irreversible manner, resulting in a quintessential path-dependency in experimental observables. While sequential learning techniques have been effectively deployed for accelerating learning of state properties of materials, they often use a consistent environment path in all experiments. To elevate such techniques for making optimal decisions in experimental investigations of path-dependent properties, we introduce an iterated expected information gain acquisition function that optimizes over entire experimental trajectories. This approach is implemented within a cloud-based Materials Acceleration Platform architecture utilizing an event-driven stateful broker coupled with remote HELAO (Hierarchical Experimental Laboratory Automation and Orchestration) instances and an AI science manager. The platform's efficacy was demonstrated through a case study optimizing multi-step spectro-electrochemical experiments to identify optically stable potential windows in (Co–Ni–Sb)O z metal oxides. The system successfully integrated AI-driven experiment design, remote laboratory automation, and cloud-based data infrastructure, validating the platform's capability for managing complex, adaptive, path-dependent workflows in materials discovery.

Guevarra, Dan [California Institute of Technology ↗
DOE OSTI2026

Jahn–Teller distortion controls electron transfer in photoexcited Cu( I ) donor–acceptor systems

The Jahn–Teller distortion (JTD) is a defining structural response to electronic excitation in Cu( I )-based transition metal complexes, yet its role in photoinduced electron transfer (PET) remains largely unexplored. Here, we demonstrate that the JTD governs charge separation (CS) through a vibronically controlled conical intersection in heteroleptic Cu( I ) bisphenanthroline–naphthalene diimide (CuHETPHEN-NDI) donor–acceptor dyads. Using 20-fs broadband transient absorption spectroscopy combined with coherent vibrational wavepacket (CVWP) analysis and quantum chemical calculations, we directly track nuclear motions that steer the system from the metal-to-ligand charge-transfer ( 1 MLCT) state to the CS state. Steric bulkiness of the pendant groups at the 2,9-positions of the phenanthroline ligand systematically slows both JTD and CS. Short-time Fourier transformation and Fourier filtering analyses identify two key vibrational signatures: a low-frequency breathing mode (∼100 cm −1 ) via a bond distance change between Cu and ligated Ns (Cu–N) that modulates the NDI anion absorption and acts as a vibronic coupling coordinate, and a higher-frequency mode (∼313 cm −1 ) that evolves along the PET trajectory. Normal mode analysis and potential energy surface calculations show that the JTD brings the 1 MLCT and CS states into degeneracy, while the Cu–N breathing motion dynamically modulates donor–acceptor electronic coupling to enable ultrafast nonadiabatic electron transfer. Steric hindrance exerted by the groups at the 2,9 positions of the phenanthroline ligands suppresses this vibronic coupling, leading to faster CVWP decoherence for the 313 cm −1 mode and slower CS. These findings unravel JTD-controlled vibronic coupling at conical intersection as a governing factor for CS and provide insight into designing Cu-based photosensitizers by harnessing structural dynamics to control PET.

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

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗