Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ion trajectory simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

Waveform-dependent air fluorescence from neutral and ionic nitrogen molecules

Laser-induced air fluorescence in the ultraviolet regime is primarily attributed to transitions between the C and B states in excited neutral N 2 molecules and between the B and X states in N$^+_2$ ions. However, the mechanism underlying the former remains contentious, as direct population to the C state by light fields is forbidden by electron spin constraints. In this work, we investigate the mechanism of air fluorescence from excited neutral N 2 molecules by carrier-envelope phase–stabilized sub–4 femtosecond pulses. Our results show that fluorescence from N$^+_2$ ions reaches a maximum with cosine-like pulses, while fluorescence from excited neutral N 2 molecules peaks with sine-like pulses. In addition, by scanning the chirp of the driving pulse, we find that ionic fluorescence is maximized with chirp-free pulses, whereas neutral fluorescence favors negatively chirped pulses. These observations, supported by classical trajectory Monte Carlo simulations, support the mechanism of intersystem crossing from excited spin-singlet states, which are populated via recollision-induced strong-field excitation.

Science & Technology - Other Topics↗
DOE OSTI2023

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗
DOE OSTI2020

3D xRAGE simulation of inertial confinement fusion implosion with imposed mode 2 laser drive asymmetry

Low-mode asymmetries represent an important obstacle to achieving high-gain inertial confinement fusion implosions. As a step in learning how to control such effects, an OMEGA experiment with imposed mode 2 laser drive asymmetries was done to study the expected signatures of this type of asymmetry [M. Gatu Johnson et al., PRE 2018]. In the present work, a 3D xRAGE simulation including the stalk mount has been brought to bear on the data from that experiment. Furthermore, comprehensive comparisons between simulated and measured observables are made. Good agreement between simulated and measured x-ray image-inferred shell trajectories, bang times and neutron emission widths are seen, showing that the hydrodynamics are well captured in the simulation. Asymmetries seen in simulated and measured time-resolved and time-integrated x-ray images and areal densities also compare well, showing impact of both stalk and mode 2. On the other hand, important differences in measured and simulated neutron emission histories, yield, and ion temperature (T ion ) asymmetries are seen, suggesting that the simulation is overestimating shock yield. Overall, the results clearly demonstrate the importance of considering all asymmetry sources when interpreting measured signatures of asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

Pulse-level noisy quantum circuits with QuTiP

The study of the impact of noise on quantum circuits is especially relevant to guide the progress of Noisy Intermediate-Scale Quantum (NISQ) computing. In this paper, we address the pulse-level simulation of noisy quantum circuits with the Quantum Toolbox in Python (QuTiP). We introduce new tools in qutip-qip, QuTiP's quantum information processing package. These tools simulate quantum circuits at the pulse level, leveraging QuTiP's quantum dynamics solvers and control optimization features. We show how quantum circuits can be compiled on simulated processors, with control pulses acting on a target Hamiltonian that describes the unitary evolution of the physical qubits. Various types of noise can be introduced based on the physical model, e.g., by simulating the Lindblad density-matrix dynamics or Monte Carlo quantum trajectories. In particular, the user can define environment-induced decoherence at the processor level and include noise simulation at the level of control pulses. We illustrate how the Deutsch-Jozsa algorithm is compiled and executed on a superconducting-qubit-based processor, on a spin-chain-based processor and using control optimization algorithms. We also show how to easily reproduce experimental results on cross-talk noise in an ion-based processor, and how a Ramsey experiment can be modeled with Lindblad dynamics. Finally, we illustrate how to integrate these features with other software frameworks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Flying focus with arbitrary directionality for spatiotemporal control of laser intensity

Flying focus techniques produce laser pulses whose focal points travel at arbitrary, controllable velocities. While this flexibility can enhance a broad range of laser-based applications, existing techniques constrain the motion of the focal point to the propagation direction of the pulse. Here, we introduce a flying focus configuration that decouples the motion of the focus from the propagation direction. Here, a chirped laser pulse focused and diffracted by a diffractive lens and grating creates a focal point that can move both along and transverse to the propagation direction. The focal length of the lens, grating period, and chirp can be tuned to control the direction and velocity of the focus. Simulations demonstrate this control for a holographic configuration suited to high-power pulses, in which two off-axis pump beams with different focal lengths encode the equivalent phase of a chromatic lens and grating in a gas or plasma. For low-power pulses, conventional solid-state or adaptive optics can be used instead. Multi-dimensional control over the focal trajectory enables new configurations for applications, including laser wakefield acceleration of ions, nonlinear Thomson scattering, and surface-plasmon emission of THz radiation.

Classical optics↗
DOE OSTI2023

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Extreme radiation emission regime for electron beams in strong focusing ion channels and undulators

A fundamental comparison between a magnetic undulator and an ion channel, or betatron, radiation from relativistic electrons is presented. While conventional theories nominally range from the undulator (𝐾 <1) to the wiggler (𝐾 >1) regime, they are only applicable for sufficiently large Lorentz factors (𝛾 0 ≫𝐾). They therefore do not account for high 𝐾/𝛾 0 cases, for which we show that particle trajectories and radiation characteristics strongly deviate from the linear predictions in both magnetic undulators and ion channels. This problem arises from the fundamental differences between a magnetostatically and electrostatically induced oscillation. A reformulation of both the ion channel betatron wavelength and amplitude, as well as the same parameters in a magnetic undulator, permits us to compare cases with equivalent oscillation period and amplitude in the two different scenarios. The notable differences in spectral features of the two radiation mechanisms can then be addressed via numerical simulations of single particle as well as full beam dynamics. Additionally, we identify and quantify a novel transverse orbit precession effect in ion channels for particles with initial angular momentum relative to the device axis. This effect, which is significant in cases of strong transverse kinetic energy oscillations, alters both the radiation divergence and the beam emittance. In this paper, we present this new theoretical framework and compare its results with numerical simulation applied to realizable experimental tests of such radiation sources.

Frazzitta, Andrea [University of Rome “La Sapienza↗
DOE OSTI2025

Investigation of post-breakup Coulomb acceleration using a trajectory model

Intermediate mass fragments ejected during the deexcitation of excited projectilelike fragments may promptly decay following ejection; the daughter particles that are subsequently produced are subject to interactions with the residual nucleus that affect final-state observables, a process herein referred to as post-breakup Coulomb acceleration. A simple classical Coulomb interaction model was used to study modification of 8 Be (2 + ), 5 Li (3/2 – ), 7 Li (7/2 – ), 7 Be (7/2 – ), and states in 12 B emitted following heavy-ion collisions of 28 Si + 12 C at 35 MeV/nucleon. Here, in contrast to previous work studying 8 Be (2 + ), excellent agreement between simulation and experiment was obtained using only Coulomb forces when either a Lorentzian or R-matrix line shape was used to describe the initial relative energy rather than a Gaussian. In consideration of the obtained results, improvements to the model and evaluation of experimental data are discussed as future directions, but it was concluded that the effects observed in the present data can be accurately described using only elements of classical mechanics and that the process is largely understood for a wide range of state lifetimes and mass and charge (a)symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2025

String-based axial and helicity-flip GPDs: A comparison to lattice QCD

We construct an analytic, string-based representation of the nucleon’s axial and helicity-flip conformal moments of generalized parton distributions (GPDs) that holds for skewness and for both the quark and gluon channels. The starting point is the Mellin-Barnes resummation of the conformal partial-wave expansion, where the moments are parametrized by open- (Reggeon) and closed-string (Pomeron) trajectories with slopes determined by experimental form factors and meson/glueball spectroscopy. The forward limits are fixed by the empirical unpolarized and polarized parton distributions. Polynomiality, crossing symmetry, and support are satisfied by construction. After next-to-leading order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi/Efremov-Radyushkin-Brodsky-Lepage evolution to μ = 2 GeV our analytic framework (i) reproduces some of the currently available lattice moments of E and H˜ in the nonsinglet sector, (ii) predicts sea-quark and gluon polarized moments that will be testable by forthcoming simulations and experiments at Jefferson Lab and the future Electron-Ion Collider, and (iii) yields axial and helicity-flip GPDs in x-space in reasonable agreement with lattice QCD.

Gauge-gravity dualities↗
DOE OSTI2025

Light control of intramolecular nuclear dynamics by vortex electron localization

In strong-field ionization of molecules, intense light pulses are thought to have a negligible direct influence on atomic nuclei. Molecular dissociation is thus expected to be determined by the geometrical configuration of the molecular ion at the ionization instant. Contrary to this picture, we observe a counterintuitive electron-proton angular correlation and the formation of proton vortices following strong-field ionization of H 2 molecules by bicircularly polarized two-color laser fields. We explain this phenomenon by the pathway interference and localization of the residual H 2 ⁡+ electron in different angular-momentum states formed in the tail of the driving laser pulse. We validate this interpretation by combining a quantum-mechanical numerical simulation of the field-driven coupled electronic-nuclear dynamics and a semiclassical-trajectory model for the phase accumulation of the laser-driven electronic-nuclear wave packet. Our joint experimental and theoretical study reveals a general picture of vortex electron localization which can be used for controlling molecular-bond breaking with circularly polarized laser fields.

Atomic & molecular processes in external fields↗
DOE OSTI2021

Strong-field ionization of water. II. Electronic and nuclear dynamics en route to double ionization

We investigate the role of nuclear motion and strong-field-induced electronic couplings during the double ionization of deuterated water using momentum-resolved coincidence spectroscopy. By examining the three-body dicationic dissociation channel, D$^{+}$/D$^{+}$/O, for both few- and multi-cycle laser pulses, strong evidence for intra-pulse dynamics is observed. The extracted angle- and energy-resolved double ionization yields are compared to classical trajectory simulations of the dissociation dynamics occurring from different electronic states of the dication. In contrast with measurements of single photon double ionization, pronounced departure from the expectations for vertical ionization is observed, even for pulses as short as 10~fs in duration. We outline numerous mechanisms by which the strong laser field can modify the nuclear wavefunction en-route to final states of the dication where molecular fragmentation occurs. Specifically, we consider the possibility of a coordinate-dependence to the strong-field ionization rate, intermediate nuclear motion in monocation states prior to double ionization, and near-resonant laser-induced dipole couplings in the ion. These results highlight the fact that, for small and light molecules such as D$_2$O, a vertical-transition treatment of the ionization dynamics is not sufficient to reproduce the features seen experimentally in the strong field coincidence double-ionization data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Decoding the Pair Distribution Function of Uranium in Molten Fluoride Salts from X-Ray Absorption Spectroscopy Data by Machine Learning

Thermal properties of actinides in molten salts are linked to the strongly disordered local environment of actinide ions. Here, we illustrate both the limitations of the commonly used fitting method for analysis of extended X-ray absorption fine structure (EXAFS) spectra in molten UF 4 and a possible solution using an "objective neural network - EXAFS" (ONNE) method. ONNE provides both extraction of the pair distribution function, as validated by its application to the EXAFS spectra calculated on molecular dynamics trajectory, and the EXAFS data reconstruction. The ONNE analysis of the molten UF4 has revealed reduction of the first nearest neighbor U-F coordination number, expansion of the U-F bond length and smaller contribution to the second shell compared to its crystalline counterpart. This method is therefore an attractive alternative to conventional EXAFS analysis and molecular dynamics simulations for studies of disordered environment of actinides in molten salts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗