Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical phase transitions”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2023

First-principles molten salt phase diagrams through thermodynamic integration

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.

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

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Bridging material models across scales: An integrated approach to equation of state and molecular dynamics modeling of copper

New uncertainty-aware equation of state (EOS) and electrical conductivity (EC) models for copper have been developed. The multiphase EOS/EC models are fit to experimental solid/liquid EC isobar measurements as well as density-functional theory molecular dynamics (DFT-MD) EC calculations in both expanded and compressed regimes (0.1–16 g/ cm 3 ⁠). The liquid and solid EOS phases were fit to available experimental data along with additional DFT-MD data over the same range as the EC. Leveraging the DFT-MD data, a corresponding machine-learned interatomic potential (MLIAP) for copper was trained using genetic-algorithm optimization. The copper MLIAP was constrained by EOS shock points at high compressions. The final EOS bounded MLIAP proves to be stable over a large density range (approximately 0.1–20 g/ cm 3 ) with good agreement to an isothermal compression curve, shock Hugoniot, and liquid speed of sound measurements at high pressures (100s of GPa).

Acoustic measurements and instrumentation↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interplay between aging and rejuvenation in spin glasses

Aging in a single-crystal spin glass (Cu 0.92 ⁢Mn 0.08 ) has been measured using ac susceptibility techniques over a temperature range of 0.3–0.8𝑇 𝑔 . In these studies, traditional aging experiments (or “quench” aging protocols) are compared to aging curves constructed from finite-cooling-rate curves. By comparing the growth rates of spin glass order between the two types of aging curves, it is determined that quantitative comparisons between protocols that are taken by quenching and protocols using a finite cooling rate are not possible without a deeper understanding of the interplay between aging and rejuvenation. Finally, we then demonstrate that the data presented indicate that rejuvenation, rather than cumulative aging, is the cause of the discrepancies between the two growth rates.

dynamical phase transitions↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗
DOE OSTI2020

Nonequilibrium renormalization group fixed points of the quantum clock chain and the quantum Potts chain

In this work, we derive an exact renormalization group recursion relation for the Loschmidt amplitude of the quantum Q-state clock model and the quantum Q-state Potts model in one dimension. The renormalization group flow is discussed in detail. The fixed points of the renormalization group flow are found to be complex in general. These fixed points control the dynamical phases of the two models, giving rise to nonanalyticities in its Loschmidt rate function, for both the pure and the disordered system. For the quench protocols studied, dynamical quantum phase transitions are found to occur in the clock model for all Qs considered, while in the Potts model, they only occur when Q < 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Instabilities of the normal state in current-biased narrow superconducting strips

We study the current-voltage characteristic of narrow superconducting strips in the gapless regime near the critical temperature in the framework of the Ginzburg-Landau model. Our focus is on its instabilities occurring at high current biases. The latter are consequences of dynamical states with periodic phase-slip events in space and time. We analyze their structure and derive the value of the reentrance current at the onset of the instability of the normal state. It is expressed in terms of the kinetic coefficient of the time-dependent Ginzburg-Landau equation and calculated numerically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Measurement of the sound velocity and Grüneisen parameter of polystyrene at inertial confinement fusion conditions

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of ICF experiments. The Grüneisen parameter increases with density below ~3.1 g/cm 3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. Finally, the discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Switching of the electron-phonon interaction in 1 T –VSe 2 assisted by hot carriers

We apply an intense infrared laser pulse in order to perturb the electronic and vibrational states in the three-dimensional charge density wave material 1T-VSe 2 . Ultrafast snapshots of the light-induced hot carrier dynamics and nonequilibrium quasiparticle spectral function are collected using time- and angle-resolved photoemission spectroscopy. Furthermore, the hot carrier temperature and time-dependent electronic self-energy are extracted from the time-dependent spectral function, revealing that incoherent electron-phonon interactions heat the lattice above the charge density wave critical temperature on a timescale of (200 ± 40) fs. Density functional perturbation theory calculations establish that the presence of hot carriers alters the overall phonon dispersion and quenches efficient low-energy acoustic phonon scattering channels, which results in a new quasiequilibrium state that is experimentally observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Skyrmion ratchet in funnel geometries

Here, using a particle-based model, we simulate the behavior of a skyrmion under the influence of asymmetric funnel geometries and ac driving at zero temperature. We specifically investigate possibilities for controlling the skyrmion motion by harnessing a ratchet effect. Our results show that as the amplitude of a unidirectional ac drive is increased, the skyrmion can be set into motion along either the easy or hard direction of the funnel depending on the ac driving direction. When the ac drive is parallel to the funnel axis, the skyrmion flows in the easy direction and its average velocity is quantized. In contrast, when the ac drive is perpendicular to the funnel axis, a Magnus-induced ratchet effect occurs, and the skyrmion moves along the hard direction with a constant average velocity. For biharmonic ac driving of equal amplitude along both the parallel and perpendicular directions, we observe a reentrant pinning phase where the skyrmion ratchet vanishes. For asymmetric biharmonic ac drives, the skyrmion exhibits a combination of effects and can move in either the easy or hard direction depending on the configuration of the ac drives. These results indicate that it is possible to achieve controlled skyrmion motion using funnel geometries, and we discuss ways in which this could be harnessed to perform data transfer operations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamics and nonmonotonic drag for individually driven skyrmions

Here, we examine the motion of an individual skyrmion driven through an assembly of other skyrmions in the absence of quenched disorder. The skyrmion behavior is determined by the ratio of the damping and Magnus terms, known as the intrinsic skyrmion Hall angle. For fixed drive in the damping dominated regime, the effective viscosity decreases monotonically with increasing skyrmion density. In contrast, in the Magnus dominated regime, the velocity varies nonmonotonically with density, and there is a regime in which the skyrmion moves faster with increasing density, as well as a pronounced speed-up effect in which a skyrmion traveling through a dense medium moves more rapidly than for low densities or the single-particle limit. The velocity-force curves in the Magnus-dominated regime exhibit an overshoot effect and other differences from the damping-dominated regime. We find a finite threshold force for skyrmion motion which increases with density as well as a drive-dependent skyrmion Hall angle. We map dynamic phase diagrams showing the threshold for motion, nonlinear flow, speed-up, and saturation regimes. In some cases, increasing the density can reduce the skyrmion Hall angle while producing a velocity boost, which could be valuable for applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural and vibrational properties of methane up to 71 GPa

Single-crystal synchrotron x-ray diffraction, Raman spectroscopy, and first principles calculations have been used to identify the structure of the high-pressure (HP) phase of molecular methane above 20 GPa up to 71 GPa at room temperature. The structure of the HP phase is trigonal R3, which can be represented as a distortion of the cubic phase B, previously documented at 7–15 GPa and confirmed here. The positions of hydrogen atoms in the HP phase have been obtained from first principles calculations, which also demonstrated the stability of this structure above 260 K at 25 GPa. The molecules occupy four different crystallographic sites in phase B and 11 sites in the HP phase, which result in splitting of molecular stretching modes detected in Raman spectroscopy and assigned here based on a good agreement with the Raman spectra calculated from the first principles. Furthermore, our study points out to a single HP phase unlike up to three HP phases proposed previously based on the Raman spectroscopy results only.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Reversible motions and disordered structure of soft particles in amorphous solids

In amorphous solids, soft vibrational modes derived from normal mode analysis can be utilized to identify the soft particles that are prone to irreversible rearrangements. However, the normal mode analysis approach cannot explain why the spatial distributions of clustered soft particles do not change over time. We define a softness parameter based on the vibrational density of states calculated directly from molecular dynamics simulations with both the harmonic vibrations and anharmonic relaxations being properly captured at finite temperatures. This parameter spontaneously correlates with the real space atomic motions and the dynamics heterogeneity. Using the softness parameter, we show that the softest particles are confined within rigid cages. These particles keep rearranging reversibly without long-range diffusion. The moderately soft particles rearrange irreversibly, and the hard particles mainly participate in vibrations without rearrangement. Here, we also show that the soft particles form locally disordered structures, while the hard particles present strong ordering. These findings confirm the defective nature of soft particles, and provide insights on the nature of softness as the ability to rearrange, but not necessarily irreversibly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonequilibrium DMFT + CPA for correlated disordered systems

We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. Finally, the real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗