Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lattice Anisotropy-Driven Reduction of Phonon Velocities in Black Phosphorus

Phonon dynamics and transport determine how heat is utilized and dissipated in materials. In 2D systems for optoelectronics and thermoelectrics, the impact of nanoscale material structure on phonon propagation is central to controlling thermal conduction. Here, we directly observe in-plane coherent acoustic phonon propagation in black phosphorus (BP) using ultrafast electron microscopy. We identify a significant reduction of the group velocities in directions intermediate to the armchair and zigzag lattice directions. Using a machine learning-based model with an >8000 atom supercell, we find that this slowing results from the mixing of in-plane transverse and longitudinal acoustic phonons and is independent of broken symmetries of edge reconstructions. In conclusion, this work demonstrates how coherent phonon transport is sensitive to propagation direction in the lattice plane.

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

Evidence for superionic H 2 O and diffusive He–H 2 O at high temperature and high pressure

In this work, we present the evidence of superionic phase formed in H 2 O and, for the first time, diffusive H 2 O–He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H 2 O and diffusive He–H 2 O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H 2 O phase forms gradually in pure H 2 O over the temperature range of 1350–1400 K at 23 GPa, but the diffusive He–H 2 O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He–H 2 O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He–H 2 O.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Multipole groups and fracton phenomena on arbitrary crystalline lattices

Multipole symmetries are of interest in multiple contexts, from the study of fracton phases, to nonergodic quantum dynamics, to the exploration of new hydrodynamic universality classes. However, prior explorations have focused on continuum systems or hypercubic lattices. In this work, we systematically explore multipole symmetries on arbitrary crystal lattices. We explain how, given a crystal structure (specified by a space group and the occupied Wyckoff positions), one may systematically construct all consistent multipole groups. We focus on two-dimensional crystal structures for simplicity, although our methods are general and extend straightforwardly to three dimensions. We classify the possible multipole groups on all two-dimensional Bravais lattices, and on the Kagome and breathing Kagome crystal structures to illustrate the procedure on general crystal lattices. Using Wyckoff positions, we provide an in-principle classification of all possible multipole groups in any space group. We explain how, given a valid multipole group, one may construct a consistent lattice Hamiltonian and a low-energy field theory. We then explore the physical consequences, beginning by generalizing certain results originally obtained on hypercubic lattices to arbitrary crystal structures. Next, we identify two apparently novel phenomena: an emergent, robust subsystem symmetry on the triangular lattice, and an exact multipolar symmetry on the breathing Kagome lattice that does not include conservation of charge (monopole), but instead conserves a vector charge. This makes clear that there is new physics to be found by exploring the consequences of multipolar symmetries on arbitrary lattices, and this work provides the map for the exploration thereof, as well as guiding the search for emergent multipolar symmetries and the attendant exotic phenomena in real materials based on nonhypercubic lattices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗
DOE OSTI2026

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences↗
DOE OSTI2024

Unitary coupled-channel three-body amplitude with pions and kaons

Three-body dynamics above threshold is required for the reliable extraction of many amplitudes and resonances from experiment and lattice QCD. The S-matrix principle of unitarity can be used to construct dynamical coupled-channel approaches in which three particles scatter off each other, rearranging two-body subsystems by particle exchange. This paper reports the development of a three-body coupled-channel, amplitude including pions and kaons. The unequal-mass amplitude contains two-body S- and P-wave subsystems (“isobars”) of all isospins, 𝐼 = 0,1/2,1,3/2,2 , and it also allows for transitions within a given isobar. The 𝑓 0 ⁡(500)⁢(𝜎),𝑓 0 ⁡(980),𝜌⁡(700),𝐾$^{*}_{0}$⁡(700)⁢(𝜅), and 𝐾*⁡(892) resonances are included, apart from repulsive isobars. Different methods to evaluate the amplitude for physical momenta are discussed. Production amplitudes for 𝑎 1 quantum numbers are shown as a proof of principle for the numerical implementation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Phenomenology of Intermediate Molecular Dynamics at Metal-Oxide Interfaces

Reaction intermediates buried within a solid-liquid interface are difficult targets for physiochemical measurements. They are inherently molecular and locally dynamic, while their surroundings are extended by a periodic lattice on one side and the solvent dielectric on the other. Challenges compound on a metal-oxide surface of varied sites and especially so at its aqueous interface of many prominent reactions. Recently, phenomenological theory coupled with optical spectroscopy has become a more prominent tool for isolating the intermediates and their molecular dynamics. The following article reviews three examples of the SrTiO3-aqueous interface subject to the oxygen evolution from water: reaction-dependent component analyses of time-resolved intermediates, a Fano resonance of a mode at the metal-oxide–water interface, and reaction isotherms of metastable intermediates. The phenomenology uses parameters to encase what is unknown at a microscopic level to then circumscribe the clear and macroscopically tuned trends seen in the spectroscopic data.

Chemistry↗
DOE OSTI2022

Samarium: from a distorted-fcc phase to melting under dynamic compression using in-situ x-ray diffraction

Abstract Lattice and electronic structure interactions for f -electrons are fundamental challenges for lanthanide equation of state development. Difficulties in first-principles calculations, such as density functional theory (DFT), emphasize the need for well-characterized experimental data. Here, we measure in-situ x-ray diffraction of shocked samarium (Sm) and temperature along the Hugoniot for the first time, providing direct evidence for phase transitions. We report direct evidence of a distorted fcc (dfcc) phase at 23 GPa. Shocked samarium melts from the dfcc phase starting at 33 GPa (1333 K), with complete melt at 40 GPa (1468 K). Previous work indicated shock melt at 27 GPa (1200 K), underscoring the significance of x-ray measurements for detecting phase transitions. Interestingly, our observed melting is in sharp contrast with the melting reported by a diamond anvil cell study. These experimental data can tightly constrain first principles calculations and serve as key touchstones for equation of state modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Algebraic Compression of Free Fermionic Quantum Circuits: Particle Creation, Arbitrary Lattices and Controlled Evolution

In this work [1], we extend our recently introduced algebraic circuit compression algorithms [2], [3] that can compress time evolution circuits of free fermionic Hamiltonians on an n-site 1D chain, equation H(t)=∑i=1n-1 (hi(t)cici+1+pi(t)cici+1)+h.c., 1 equation in three significant ways: (1) we allow for compression of free fermionic Hamiltonians on arbitrary lattices, (2) we incorporate particle creation/annihilation operators into the compression schemes, and (3) we extend the compression scheme to controlled time-evolution operators. We illustrate the effectiveness of our approach by simulating the dynamics of a fermion on a 4× 4 2D square lattice on ibmq_washington, both in the presence and absence of disorder. Our quantum simulations show a remarkably high fidelity which is enabled through the compressed circuits.

Kökcü, Efekan↗
DOE OSTI2021

Two-dimensional turbulence in a massless fluid with a relativistic lattice Boltzmann model

We investigate a relativistic adaptation of the lattice Boltzmann method that reproduces the equations of motion for a turbulent, two-dimensional, massless hydrodynamic system. The classical lattice Boltzmann method and its extension to relativistic fluid dynamics are described. The numeric formulation is evaluated using a zero-averaged stirring force introduced into the numerics to induce turbulence, and the flow characteristics produced are compared to properties of a classical turbulent hydrodynamic flow. The model can reasonably be expected to offer quantitative simulations of electron fluid flows in graphene or Kagome lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗
DOE OSTI2020

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Probing Dopant Redistribution, Phase Propagation, and Local Chemical Changes in the Synthesis of Layered Oxide Battery Cathodes

Achieving the targeted control of layered oxide properties calls for more fundamental studies to mechanistically probe their evolution during their synthesis. Herein, dopant distribution, phase propagation, and local chemical changes as well as their interplay in multielement-doped LiNiO 2 materials are investigated using spectroscopic, imaging, and scattering techniques. It is shown that dopants undergo dynamic redistribution in the Ni(OH) 2 host lattice at the early stage of calcination (below 300 °C). Such redistribution behavior exhibits strong dopant-dependent characteristics, allowing for targeted surface and bulk doping control. The Ni oxidation process exhibits depth-dependent characteristics and the most rapid Ni oxidation takes place between 300 and 700 °C. Using Ni oxidation state as the proxy for the phase transformation, the buildup of heterogenous phase propagation in the early stage of calcination is shown, especially along the radial direction of secondary particles. Furthermore, the radial heterogenous phase distribution gradually decreases upon completing the calcination. However, a high degree of mosaic-like heterogeneity may still be present in the final product, departing from the perfect layered oxide. Overall, the present study offers fundamental insights into manipulating multiscale materials properties during calcination for obtaining stable, high-energy layered oxide cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic Structure and Safety Insights into Prussian Blue Analog Cathode Behavior at Elevated Temperatures in Sodium-Ion Batteries

Prussian blue analogs (PBAs) represent promising cathode materials for sodium-ion batteries (SIBs) due to their high theoretical capacity, open framework structure, and use of earth-abundant elements. However, the high-temperature structural evolution, water content effects, and thermal safety of PBAs, particularly in charged states, remain poorly understood, hindering their practical deployment. Here, we investigate Na 2 Fe[Fe(CN) 6 ]·2H 2 O using thermogravimetric analysis (TGA), ex situ and in situ temperature-dependent X-ray absorption spectroscopy (XAS), and accelerated rate calorimetry (ARC). TGA and ex situ XAS confirm water loss between 150 and 200 °C, resulting in Fe 2+ oxidation, enhanced local symmetry, and uniform redox behavior that improves electrochemical performance. In situ XAS reveals irreversible structural changes above 240 °C, including ligand loss, Fe site distortion, and increased disorder, while ARC on charged electrodes shows minimal self-heating rates (<0.1 °C/min) up to 300 °C, indicating exceptional thermal stability without lattice oxygen release. These insights elucidate PBA thermal dynamics, demonstrating improved electrochemical performance of water-deficient PBAs and informing future material design and safety assessment for SIB applications.

batteries↗
DOE OSTI2020

Moving beyond the Solvent-Tip Approximation to Determine Site-Specific Variations of Interfacial Water Structure through 3D Force Microscopy

Although interfacial solution structure impacts environmental, biological, and technological phenomena, including colloidal stability, protein assembly, heterogeneous nucleation, and water desalination, its molecular details remain poorly understood. Here, we visualize the three-dimensional (3D) hydration structure at the boehmite(010)-water interface using fast force mapping (FFM). Using a self-consistent scheme to decouple long-range tip-sample interactions from short-range solvation forces, we obtain the solution structure with lattice resolution. The results are benchmarked against molecular dynamics simulations that explicitly include the effects of the tip with different levels of approximation and systematically account for tip size, chemistry, and confinement effects. We find four laterally structured water layers within 1 nm of the surface, with the highest water densities at sites adjacent to hydroxyl groups. The key features beyond the first two layers can only be predicted using a full-scale simulation of the boehmite-water-silica system. Our findings further reveal a complex relationship between site-specific chemistry, water density, and long-range particle interactions; and present important advances toward quantitative data interpretation in 3D FFM.

36 MATERIALS SCIENCE↗