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.
Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application
Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.
Quantum Computation of Dynamical Quantum Phase Transitions and Entanglement Tomography in a Lattice Gauge Theory
Strongly coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute nonequal-time correlation functions and perform entanglement tomography of nonequilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently predicted dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent nonequilibrium dynamics in three ways: (i) overlap echos signaling dynamical transitions, (ii) nonequal-time correlation functions with an underlying topological nature, and (iii) the entanglement structure of nonequilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.
Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect
Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).
A new self-adaptive reconstruction method to identify defects through Wigner–Seitz approach
A new self-adaptive reconstruction method based on local atomic structure at any given molecular dynamics (MD) step has been developed in this article. The method can be used in Wigner–Seitz defect analysis approach to correctly and efficiently explore the information of both point defects and complex defect clusters (e.g. dislocation loops and voids) formed after a displacement cascade where the cascade interacts with grain boundaries and/or dislocations. The algorithm and validation are provided in detail. Results for identification of radiation defects during and after cascades interacting with a dislocation network show that the new method can well recognize all simple and complex defects and defect clusters. Thus, this new method provides a totally new way to explore the density and size of radiation defects at atomic scale after complex MD evolution processes, providing correct information to understand and predict radiation damage in materials through atomic simulations.
Direct Observation of Ultrafast Hydrogen Bond Strengthening in Liquid Water
Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network, including unusually efficient vibrational energy redistribution and relaxation. An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments. However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here we measured the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of ~1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on the femtosecond timescales. Our experiment and simulations unveiled the intermolecular character of the water vibration preceding the relaxation of the OH stretch.
Full-dimensional quantum studies of vibrational energy transfer dynamics between H 2 O and Ar: theory assessing experiment
In this work we report the first full-dimensional quantum mechanical calculations of the ro-vibrational inelastic scattering dynamics between water molecules and argon atoms on an accurate potential energy surface, using a recently developed time-independent quantum method based on the close-coupling approach. The state-to-state integral cross-sections and rate coefficients show strong observance of gap laws. The calculated thermal rate coefficients for the relaxation of the stretching fundamental states of H 2 O are in good agreement with experimental values, while those for the bending overtone state are approximately five times smaller than the values extracted through a previous kinetic modeling of fluorescence decay data. Our state-specific quantum scattering results suggest the need to reassess the kinetic modeling of the experimental data. This work advanced our understanding of the quantum dynamics of vibrationally inelastic energy transfer processes involving polyatomic molecules.
Atomic Origin of the Autocatalytic Reduction of Monoclinic CuO in a Hydrogen Atmosphere
Reducibility is key for the use of bulk metal oxides in chemical transformation involving redox reactions but probing microscopic processes of oxide reduction is challenging. This is because the insulating nature of bulk oxides restricts ion and electron spectroscopic measurements of oxide surfaces. Herein, using a combination of environmental transmission electron microscopy and atomistic modeling, we report direct in-situ atomic-scale observations of the surface and subsurface dynamics and show that the hydrogen induced CuO reduction occurs via the surface decay of Cu-O/Cu bilayer atomic steps, formation of partially reduced CuO superstructure by the self-ordering of O vacancies in the subsurface, and collapse of Cu-O layers in the bulk. In this work, all these substeps can be traced back to the progressively increased concentration and activity of O vacancies in the surface and subsurface of the oxide, thereby leading to the self-accelerated oxide reduction. These results demonstrate the microscopic details that may have broader applicability in modulating various redox processes.
In-Situ Atomic-Scale Revelation of Amorphous Metallic Iron Formation during Hydrogen-Driven Reduction of Iron Oxides
The transition to hydrogen as a green reductant in metal production is critical for decarbonizing the metallurgical industry, yet atomic-scale mechanisms governing reduction pathways and phase evolution remain unresolved. Using in-situ environmental transmission electron microscopy, we identify a hidden pathway that reveals dynamic formation of amorphous metallic iron (Fe) during the hydrogen-driven reduction of ferrous oxides of Fe 3 O 4 and FeO. Real-time imaging uncovers three coexisting transformation routes: (i) Fe 3 O 4 → FeO, (ii) Fe 3 O 4 → amorphous Fe, and (iii) FeO → amorphous Fe. The resulting amorphous Fe exhibits fluid-like mobility, enabling its rapid aggregation and crystallization into core-shell nanostructures, with a crystalline core enveloped by an amorphous shell. Complementary ab initio molecular dynamics simulations trace the amorphous Fe formation to interfacial strain at the metal/oxide interfaces, where large lattice mismatches destabilize the metal lattice during initial metallization. This interplay between thermodynamics and kinetics governs phase evolution: thermodynamics favors a self-limiting amorphous Fe overlayer, while rapid oxide reduction kinetics drives amorphous overgrowth. Our findings demonstrate that amorphous intermediates bypass rate-limiting crystalline steps, providing mechanistic insights to optimize H 2 -based processes for sustainable steelmaking. In conclusion, these insights bridge the gap between macroscopic process engineering and atomic-scale dynamics, with broader implications for catalysis and nanostructured material synthesis, where oxide reduction pathways critically shape functional phases and microstructures.
Polyamorphism in a solute-lean Al–Ce metallic glass
Polyamorphism discovered in lanthanide-rich metallic glasses (MGs) has been attributed to the electronic transition of the lanthanide element as a solvent element. In this work, we report that pressure-induced polyamorphism still exists in a Ce-poor Al 93 Ce 7 binary MG where the 4f electron element serves as a solute and solute–solute avoidance is expected. The polyamorphic transition, observed by in situ high-pressure synchrotron x-ray diffraction, is accompanied by a volume collapse of ~0.78% and occurs over a narrow pressure range from ~0.8 to ~1.8 GPa. Further synchrotron Ce L 3 -edge x-ray absorption spectroscopy measurements reveal that pressure-induced 4f electron delocalization underlies the polyamorphic transition. Molecular dynamics simulations confirm that the Ce atoms in the MG are completely isolated by the solvent Al atoms. This result demonstrates that 4f element-bearing alloys with extremely dilute concentrations can also exhibit polyamorphic states originating from electronic transitions, extending the compositional space of polyamorphism of MGs into very dilute regions. Our work suggests that tunable properties under compressive stress could be achieved in MGs by even minor doping of elements prone to electronic transitions.
Molecular level insight into non-bilayer structure formation in thylakoid membranes: a molecular dynamics study
In oxygenic photosynthetic organisms, the light reactions are performed by protein complexes embedded in the lipid bilayer of thylakoid membranes (TMs). The organization of the bulk lipid molecules into bilayer structures provide optimal conditions for the build-up of the proton motive force (pmf) and its utilization for ATP synthesis. However, the lipid composition of TMs is dominated by the non-bilayer lipid species monogalactosyl diacylglycerol (MGDG), and functional plant TMs, besides the bilayer, contain large amounts of non-bilayer lipid phases. Bulk lipids have been shown to be associated with lumenal, stromal-side and marginal-region proteins and proposed to play roles in the self-assembly and photoprotection of the photosynthetic machinery. Furthermore, it has recently been pointed out that the generation and utilization of pmf for ATP synthesis according to the ‘protet’ or protonic charge transfer model Kell (Biochim Biophys Acta Bioenerg 1865(4):149504, 2024), requires high MGDG content Garab (Physiol Plant 177(2):e70230, 2025). In this study, to gain better insight into the structural and functional roles of MGDG, we employed all atom and coarse-grained molecular dynamics simulations to explore how temperature, hydration levels and varying MGDG concentrations affect the structural and dynamic properties of bilayer membranes constituted of plant thylakoid lipids. Our findings reveal that MGDG promotes increased membrane fluidity and dynamic fluctuations in membrane thickness. MGDG-rich stacked bilayers spontaneously formed inverted hexagonal phases; these transitions were enhanced at low hydration levels and at elevated but physiologically relevant temperatures. It can thus be inferred that MGDG plays important roles in heat and drought stress mechanisms.
Thermal Scattering Law Data Development for Paraffin Wax
Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).
Phonon confinement and interface lattice dynamics of ultrathin high- k rare earth sesquioxide films: the case of Eu 2 O 3 on YSZ(001)
The spatial confinement of atoms at surfaces and interfaces significantly alters the lattice dynamics of thin films, heterostructures and multilayers. Ultrathin films with high dielectric constants (high-k) are of paramount interest for applications as gate layers in current and future integrated circuits. Here we report a lattice dynamics study of high-k Eu 2 O 3 films with thicknesses of 21.3, 2.2, 1.3, and 0.8 nm deposited on YSZ(001). The Eu-partial phonon density of states (PDOS), obtained from nuclear inelastic scattering, exhibits broadening of the phonon peaks accompanied by up to a four-fold enhancement of the number of low-energy states compared to the ab initio calculated PDOS of a perfect Eu 2 O 3 crystal. Our analysis demonstrates that while the former effect reflects the reduced phonon lifetimes observed in thin films due to scattering from lattice defects, the latter phenomenon arises from an ultrathin EuO layer formed between the thin Eu 2 O 3 film and the YSZ(001) substrate. Thus, our work uncovers another potential source of vibrational anomalies in thin films and multilayers, which has to be cautiously considered.
Efficient diffusion of superdense lithium via atomic channels for dendrite-free lithium–metal batteries
The non-uniform aggregation of fast-diffused Li on an anode surface would aggravate its tip-effect-induced nucleation/growth, leading to the notorious dendrite growth in Li metal batteries (LMBs). Tuning the Li diffusion on the anode surface has been regarded previously as a mainstream method to induce its uniform deposition, while the diffusion of Li in the anode bulk is usually ignored. In this study, conceptually different from the classic surface modification, we propose a molecular tunnelling strategy to construct atomic channels in graphite bulk, which enables the fast diffusion of superdense Li. Density functional theory calculations and ab initio molecular dynamics simulations prove that the bulk diffusion through atomic channels could become a new and dominating path. Its reversible and efficient diffusion has been further visualized by in situ transmission electron microscopy. As a result, when coupled with high-loading LiFePO 4 cathodes (20 mg cm –2 ), a high areal capacity and 100% capacity retention are achieved over 370 cycles. Through this work a new strategy is developed based on the bulk-diffusion of superdense Li for dendrite-free LMBs, which can be pervasive in other high-performance energy storage systems.
Inhibition of Amyloid β-Induced Lipid Membrane Permeation and Amyloid β Aggregation by K162
Alzheimer's disease (AD) is characterized by progressive neurodegeneration associated with amyloid β (Aβ) peptide aggregation. The aggregation of Aβ monomers (AβMs) leads to the formation of Aβ oligomers (AβOs), the neurotoxic Aβ form, capable of permeating the cell membrane. Here, we investigated the effect of a fluorene-based active drug candidate, named K162, on both Aβ aggregation and AβO toxicity toward the bilayer lipid membrane (BLM). Electrochemical impedance spectroscopy (EIS), atomic force microscopy (AFM), and molecular dynamics (MD) were employed to show that K162 inhibits AβOs-induced BLM permeation, thus preserving BLM integrity. In the presence of K162, only shallow defects on the BLM surface were formed. Apparently, K162 modifies Aβ aggregation by bypassing the formation of toxic AβOs, and only nontoxic AβMs, dimers (AβDs), and fibrils (AβFs) are produced. Unlike other Aβ toxicity inhibitors, K162 preserves neurologically beneficial AβMs. This unique K162 inhibition mechanism provides an alternative AD therapeutic strategy that could be explored in the future.
Retrieval of the molecular orientation distribution from atom-pair angular distributions
Imaging laser-induced rotational dynamics is an important and active field due to its applications in capturing reactions in the molecular frame and in molecular imaging. Experimental measurement of the molecular orientation distribution, as a function of the Euler angles, has only been demonstrated for special cases when the detectable signal is generated along the molecular symmetry axis. Here we developed the general theory that maps the probability density distribution of the molecular orientation to the atom-pair angular distributions for nonlinear molecules. With the theory, the molecular orientation distribution can be retrieved from the measured atom-pair angular distribution, which we demonstrate experimentally using ultrafast electron diffractive imaging of impulsively aligned trifluoro-iodomethane molecules. Here, the retrieved molecular orientation distribution is in good agreement with direct numerical simulations of the time-dependent Schrodinger equation using the experimental conditions. Unlike the existing retrieval methods, the retrieval method does not require solving Schrodinger equation, works for any alignment method, and is in principle applicable to asymmetric top molecules.
The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy
In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.
Neutrons in Structural Biology: Challenges and Opportunities (Workshop Report)
Gaining a thorough understanding of biological systems requires building our knowledge about biological processes from the level of atoms and electrons, and up to whole organisms. Such comprehensive knowledge will allow for a predictive understanding of complex biological systems behavior. It will guide us in the design and development of novel therapeutics and vaccines to tackle existing health threats and to prepare for future pandemics, and it will provide information necessary to create new biomaterials and bio-inspired technologies through manipulation of biological macromolecules, their assemblies, single cells and even microorganisms. Reaching these goals will require a synergistic combination of multiple experimental techniques with molecular calculations and predictive simulations, and the design and development of new techniques and capabilities that bridge current knowledge and technology gaps. Neutron scattering provides unique information about the biomacromolecular structure and function and can play a major role in achieving these goals. A workshop was held to engage the scientific community in identifying pressing challenges in biochemistry, structural biology, enzymology and structure-guided drug design not solved with the current neutron scattering technologies or utilizing other structural biology techniques such as X-ray crystallography, NMR, and cryo-EM. The workshop brought together structural biology, biochemistry and computational experts, as well as early career researchers and students, creating a forum for discussing scientific advancement and collaboration. The workshop included a one-day satellite training workshop where graduate students and postdoctoral researchers were educated in the application of neutron crystallography and small-angle scattering in structural biology. Furthermore, the Instrument Scientific Advisory Board (ISAB) for the development of a macromolecular neutron diffractometer at ORNL’s Second Target Station was introduced at the workshop. The major outcome was that neutrons can provide atomic-level understanding of biomacromolecular structure, function and dynamics which is of paramount importance for addressing the identified challenges. Neutron crystallography, in particular, can resolve long-standing biochemical issues regarding enzyme function by delineating the underlying chemistry and can have a major impact on the design of small-molecule therapeutics, especially in combination with molecular computation (quantum chemistry and molecular dynamics simulations) and the emerging artificial intelligence (AI)-assisted drug design technologies. The unique properties of neutrons, including their high sensitivity to hydrogen and their non-destructive nature, make them ideal probes of biological matter. There is a palpable need in the scientific community to expand and enhance the impact of neutron sciences on biology. Neutron crystallography is the only structural biology method capable of determining positions of all hydrogen atoms in proteins, nucleic acids and their complexes at near-physiological temperatures and of unstable species at cryogenic temperatures. Moreover, neutron analysis is non-ionizing, non-destructive and does not perturb the structure or redox chemistry of active site metal centers and clusters in proteins, which can be invaluable for studying radiation-sensitive metalloprotein complexes. Further, neutron energies used in scattering applications are similar to atomic motions, permitting neutron spectroscopies to characterize the dynamics of biomacromolecules on the picosecond to microsecond timescales. The different sensitivities of neutrons to protium (H) and deuterium (D) isotopes of hydrogen allow enhanced visibility of specific parts of biological complexes through isotopic labeling. The impact of neutrons will be most powerful when neutron scattering is combined with complementary experimental techniques that use photons and electrons, and with high-performance computing. The interconnection and mutuality of the experimental and theoretical capabilities will drive discoveries in biological and health sciences to generate more complete picture of complex biological systems. The major limitation in the field of biological neutron crystallography has been signal-to-noise, demanding large samples that are difficult to produce for the majority of biomacromolecules and limiting the applicability of this technique in biological sciences. A neutron crystallography instrument at the Second Target Station will revolutionize biological science with neutrons by engaging a large scientific community of structural biologists, enabling successful neutron diffraction experiments from radically smaller biomacromolecular crystals, resolving unanswered biochemical questions, and meaningfully contributing to rational drug design. The meeting highlighted 10 grand challenges that will be addressed with this advanced capability over the next decade and beyond, and the recommendations required to help address them are given below.