Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ion fluxes and memory effects in an Ar–O 2 modulated radiofrequency-driven atmospheric pressure plasma jet

In this work we used molecular beam mass spectrometry to characterize ion fluxes produced by a modulated radiofrequency (RF)-driven atmospheric pressure plasma jet operating in a homogeneous gas environment (Ar + 1% O 2 ). The influence of the RF modulation frequency (100 Hz to 20 kHz) on the ion fluxes was investigated by time-resolved measurements, and the lifetimes of the dominant ions, O 2 + , NO + , O − , O 2 − and O 3 − , were found to be 28 ± 2, 117 ± 8, 7.3 ± 0.4, 17 ± 1 and 23 ± 2 μs, respectively. The absolute ion densities in the near afterglow region were found to be of the order of 1011 cm−3. Significant differences in the dynamics of the positive and negative ions were found; these are explained by the presence of electrons at similar densities in the afterglow produced by electron detachment reactions from negative ions due to the large concentrations of atomic oxygen and singlet delta oxygen. Transitions in ion flux dynamics for different modulation frequencies and at the startup of the plasma were analyzed together with intensified charge coupled device images recording the plasma propagation in order to assess the dynamics of plasma plume propagation and how it is impacted by ‘memory effects’. Quantitative measurements of the ion densities enabling these memory effects are reported. The results in this paper highlight the tremendous impact of memory effects on plasma propagation and the corresponding pre-ionization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2022

Simulations of graphite boronization: A molecular dynamics study of amorphization resulting from bombardment

The molecular dynamics code LAMMPS was used to simulate the bombardment of a graphite structure by atomic boron with impact energies ranging from 50–250 eV. The transient structural evolution, penetration depth, and amorphous layer thickness were analyzed. Simulations show that larger impact energies lead to a greater volume of amorphization and penetration of boron, but that the growth rate of the amorphous layer decreases with increasing fluence. Furthermore, the change in surface chemistry of the amorphized structures was studied using the ReaxFF formalism, which found that the amorphization process introduces dangling bonds thus increasing reactivity in the amorphous region.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Plasma-assisted atomic layer etching of single-crystal diamond

Applications of near-surface nitrogen-vacancy (NV) centers in diamond are often limited by surface defects created during processing. Understanding and controlling plasma-induced surface damage is important for preserving the optical and spin properties of diamond NV centers. We report molecular dynamics simulations of a novel form of plasma-aided atomic layer etching of diamond. In this proposed scheme, the initial surface modification step consists of Ar + ion bombardment. This creates an amorphous layer at the surface, the thickness of which is controlled by the ion energy. Amorphization is known to help smooth the surface, at least locally, and would also serve to sputter clean an initially contaminated surface. A second step impacts the amorphous carbon layer with O + between about 1 and 5 eV. Simulations show that this energy range will remove the amorphous carbon but will not etch the underlying diamond. Though this energy range is not trivial to achieve in a conventional low-temperature plasma, potential methods for creating these low-energy O + ions are discussed. In addition to a-C etching, the O + impacts are predicted to remove any isolated (100) diamond terraces by selectively attacking the edges of the terraces. The proposed atomic layer etching (ALE) approach reverses the conventional ALE sequence in which the surface modification step is usually chemical modification (oxidation in this case), followed by a removal step using Ar + impacts. This plasma ALE procedure is predicted to create a diamond surface that is atomically flat and defect free.

Draney, J. S. [Princeton University, NJ (United St↗
DOE OSTI2025

Tracing Long-Lived Atomic Coherences Generated via Molecular Conical Intersections

Accessing coherences is key to fully understand and control ultrafast dynamics of complex quantum systems like molecules. Most photochemical processes are mediated by conical intersections, which generate coherences between electronic states in molecules. We show with accurate calculations performed on gas-phase methyl iodide that electronic coherences of spin-orbit-split states persist in atomic iodine after dissociation. Our simulation predicts a maximum magnitude of vibronic coherence in the molecular regime of 0.75% of the initially photoexcited state population. Upon dissociation, one-third of this coherence magnitude is transferred to a long-lived atomic coherence where vibrational decoherence can no longer occur. To trace these dynamics, we propose a tabletop experimental approach—heterodyned attosecond four-wave-mixing spectroscopy. This technique can temporally resolve small electronic coherence magnitudes and reconstruct the full complex coherence function via phase cycling. Hence, heterodyned attosecond four-wave-mixing spectroscopy leads the way to a complete understanding and optimal control of spin-orbit-coupled electronic states in photochemistry.

Chemical Physics↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2025

Resolving electrochemically triggered topological defect dynamics and structural degradation in layered oxides

Understanding topological defects-controlled structural degradation of layered oxides—a key cathode material for high-performance lithium-ion batteries—plays a critical role in developing next-generation cathode materials. Here, by constructing a nanobattery in an electron microscope enabling atomic-scale monitoring of electrochemcial reactions, we captured the electrochemically driven atomistic dynamics and evolution of dislocations—a most important topological defect in material. We deciphered how dislocations nucleate, move, and annihilate within layered cathodes at the atomic scale. Specifically, we found two types of dislocation configurations, i.e., single dislocations and dislocation dipoles. Both pure dislocation glide/climb and mixed motions were captured, and the dislocation glide and climb velocities were first experimentally measured. Moreover, dislocation activity-mediated structural degradation such as crack nucleation, phase transformation, and lattice reorientation was unraveled. Our work provides deep insights into the atomistic dynamics of electrochemically driven dislocation activities in layered oxides.

Wang, Chunyang (ORCID:0000000184613952)↗
DOE OSTI2024

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO 2

Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Peptoid-based hierarchically-structured biomimetic nanomaterials: Synthesis, characterization and applications

Peptoids (or poly-N-substituted glycines) is a promising class of bioinspired sequence-defined polymers due to their highly efficient synthesis, high chemical stability, enzyme hydrolysis resistance, and biocompatibility. By tuning the side chain chemistry of peptoids, it allows for a precise control over sequences and achieving a large side-chain diversity. Due to these unique features, in the last several years, many amphiphilic peptoids were designed as highly tunable building blocks for the preparation of biomimetic nanomaterials with well-defined hierarchical structures and desired functionalities. In this paper, we provide an overview of the recent achievements in this area by dividing them into the following three aspects. First, mica- and silica-templated peptoid self-assembly was summarized. Here, the presence of inorganic substrates provides the guarantee of investigating their self-assembly mechanisms and interactions between peptoids and substrates using nanoscale characterization techniques, particularly in situ atomic force spectroscopy (AFM) and AFM-based dynamic force spectroscopy (AFM-DFS). Second, solution-phase self-assembly of peptoids into nanotubes and nanosheets was presented, as well as their self-repair properties. Third, the applications of peptoid based nanomaterials were outlined, including the construction of catalytic nanomaterials as a template and cytosolic delivery as cargoes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Collapse of helium-filled voids in extreme deformation: Dislocation mechanisms

The mechanisms responsible for the collapse of helium-filled bubbles during the passage of shock waves in monocrystalline copper are revealed. Both internal pressure (caused by pre-existing helium atoms) and bubble size are varied in molecular dynamics simulations to understand the atomistic scale deformation as they are subjected to shock compression at pressures of 48, 123, and 170 GPa, corresponding to particle velocities of 1.0, 2.0, and 2.5 km/s. Both empty and helium filled bubbles serve as dislocation sources, generating intense, localized plastic regions. There are distinct differences in the collapse of empty voids compared to He-filled bubbles, the former requiring less stress and generating a greater density of dislocations for a given shock strength. Additionally, a generalized model for dislocation emission is proposed, where the inclusion of shear stress generated by the helium bubble increases the critical stress to generate dislocations at the defect surface, demonstrating the change in plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unconventional Pathways to Carbide Phase Synthesis via Thermal Decomposition of UI 4 (1,4-dioxane) 2

UI4(1,4-dioxane) 2 was subjected to laser-based heating-a method that enables localized, fast heating (T > 2000 °C) and rapid cooling under controlled conditions (scan rate, power, atmosphere, etc.)-to understand its thermal decomposition. A predictive computational thermodynamic technique estimated the decomposition temperature of UI 4 (1,4-dioxane) 2 to uranium (U) metal to be 2236 °C, a temperature achievable under laser irradiation. Dictated by the presence of reactive, gaseous byproducts, the thermal decomposition of UI 4 (1,4-dioxane) 2 under furnace conditions up to 600 °C revealed the formation of UO 2 , UI x , and U(C 1–x O x ) y , while under laser irradiation, UI 4 (1,4-dioxane) 2 decomposed to UO 2 , U(C 1–x O x ) y , UC 2–z O z , and UC. Despite the fast dynamics associated with laser irradiation, the central uranium atom reacted with the thermal decomposition products of the ligand (1,4-dioxane = C 4 H 8 O 2 ) instead of producing pure U metal. In conclusion, the results highlight the potential to co-develop uranium precursors with specific irradiation procedures to advance nuclear materials research by finding new pathways to produce uranium carbide.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probing Interactions between Adsorbed Polymers and γ-Al 2 O 3 Using Solid-State NMR

The interactions between polymers and metal oxide surfaces affect the physical properties of polymer composites and the selectivity of catalysts used for plastic conversion. Gauging the strengths of these interactions and how they depend on adsorption conformations and polymer structure is nevertheless quite challenging. Herein, we evaluate whether a recently introduced dynamic nuclear polarization-based surface-to-atom distance measurement method can be used to gauge the specific interactions between five vinyl polymers and an oxide support (alumina). The adsorption conformations of polymer monolayers were studied using site-specific 13 C{ 27 Al} rotational-echo saturation-pulse double-resonance (RESPDOR) experiments. The method is effective in revealing the average distance between the adsorbed polymer and the surface, as well as the fraction of the monomers that are adsorbed. Polymers that possess larger side chains are more weakly adsorbed, as evidenced by their formation of more numerous loop and tail structures. In conclusion, the measured RESPDOR dephasing data further reveals that polymer adsorption conformations are largely disordered and involve few hydrogen bonding interactions.

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

Molecular Mechanisms Underlying Surfactant-Based Plastics De-Inking

Surfactant-mediated binder removal is critical for de-inking pretreatment in plastics recycling, yet the molecular mechanisms governing surfactant performance remain poorly understood. We used all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations alongside alkaline surfactant washing experiments to investigate interactions between a series of surfactants and a polyether urethane (PEU) binder in solution and on a polyethylene (PE) surface. Experiments reveal a range of de-inking efficiencies, ranging from <25 to >95% depending upon surfactant headgroup charge and tail length. AA simulations reveal that charged surfactants reach stable levels of surfactant coverage, while nonionic surfactants aggregate on the binder. CG umbrella sampling calculations quantify the thermodynamics of binder desorption in water. In ∼0.25 M surfactant solutions, up to a 52% reduction in the free energy barrier is computed, with trends in good agreement (R 2 = 0.92, Pearson’s r = –0.96, Spearman’s ρ = –0.83) with experimental de-inking efficiencies. We find that charged surfactants are more effective than nonionic surfactants for de-inking and propose three regimes of surfactant de-inking processes: good de-inking occurs in surfactants that promote PEU desorption with a low radius of gyration (R g ); moderate de-inking occurs when surfactants stabilize PEU but increase R g ; poor de-inking occurs in surfactants that aggregate on the binder and promote extensive anchoring to the surface. Overall, these molecular-level insights have the potential to guide the design of surfactant formulations for plastics recycling applications.

ink removal↗
DOE OSTI2026

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2021

Temperature-Dependent Properties of Molten Li 2 BeF 4 Salt Using Ab Initio Molecular Dynamics

Molten lithium tetrafluoroberyllate (Li 2 BeF 4 ) salt, also known as FLiBe, with a 2:1 mixture of LiF and BeF 2 is being proposed as a coolant and solvent in advanced nuclear reactor designs, such as the molten salt reactor or the fluoride salt cooled high-temperature reactor. We present the results on the structure and properties of FLiBe over a wide range of temperatures, 0–2000 K, from high-throughput ab initio molecular dynamics simulation using a supercell model of 504 atoms. The variations in the local structures of solid and liquid FLiBe with temperature are discussed in terms of a pair distribution function, coordination number, and bond angle distribution. The temperature-dependent electronic structure and optical and mechanical properties of FLiBe are calculated. The optical and mechanical property results are reported for the first time. The results above and below the melting temperature (~732 K) are compared with the experimental data and with data for crystalline FLiBe. The electronic structure and interatomic bonding results are discussed in correlation with the mechanical strength. A novel concept of total bond order density (TBOD), an important quantum mechanical parameter, is used to characterize the internal cohesion and strength in the simulated models. The results show a variation in the rate of change in properties in solid and liquid phases with anomalous behavior across the melting region. The observed trend is the decrease in mechanical strength, band gap, and TBOD in a nonlinear fashion as a function of temperature. The refractive index shows a surprising minimum at 850 K, among the tested temperatures, which lies above the melting point. These findings provide a new platform to understand the interplay between the temperature-dependent structures and properties of FLiBe salt.

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

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗