Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonequilibrium lattice models”

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 37 records · Page 2

Ultrafast simultaneous manipulation of multiple ferroic orders through nonlinear phonon excitation

Recent experimental studies have demonstrated the possibility of utilizing strong terahertz pulses to manipulate individual ferroic orders on pico- and femtosecond timescales. Here, we extend these findings and showcase the simultaneous manipulation of multiple ferroic orders in BiFeO 3 , a material that is both ferroelectric and antiferromagnetic at room temperature. We find a concurrent enhancement of ferroelectric and antiferromagnetic second-harmonic generation (SHG) following the resonant excitation of a high-frequency fully-symmetric phonon mode. Based on first-principles calculations and phenomenological modeling, we ascribe this observation to the inherent coupling of the two ferroic orders to the nonequilibrium distortions induced in the crystal lattice by nonlinearly driven phonons. Our finding highlights the potential of nonlinear phononics as a technique for manipulating multiple ferroic order parameters at once. In addition, this approach provides a promising avenue to studying the dynamical magnetic and polarization behavior, as well as their intrinsic coupling, on ultrashort timescales.

36 MATERIALS SCIENCE↗

Nonequilibrium magnons studied by inelastic neutron scattering

In this work, we explore time-resolved dynamics of non-equilibrium magnons in Rb2MnF4, a model two-dimensional square lattice Heisenberg antiferromagnet using a novel laser-neutron pump-probe technique. This technique uses a nanosecond pulsed laser to periodically excite magnons to out-of-equilibrium states and then probe the resulting dynamics using INS.

36 MATERIALS SCIENCE↗

Measurements of nonequilibrium interatomic forces using time-domain x-ray scattering

Here, we demonstrate an experimental approach to determining the excited-state interatomic forces using femtosecond x-ray pulses from an x-ray free-electron laser. We determine experimentally the excited-state interatomic forces that connect photoexcited carriers to the nonequilibrium lattice dynamics in the prototypical Peierls-distorted material, bismuth. The forces are obtained by a constrained least-squares fit of a pairwise interatomic force model to the excited-state phonon dispersion relation as measured by the time- and momentum-resolved x-ray diffuse scattering. We find that photoexcited carriers weaken predominantly the nearest-neighbor forces, which drives the measured softening of the transverse acoustic modes throughout the Brillouin zone as well as the zone-center A 1g optical mode. This demonstrates a bond-selective approach to measuring electron-phonon coupling relevant to a broad range of photoinduced phase transitions and transient light-driven states in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗

Laser-Induced Cooperative Transition in Molecular Electronic Crystal

The competing and non-equilibrium phase transitions, involving dynamic tunability of cooperative electronic and magnetic states in strongly correlated materials, show great promise in quantum sensing and information technology. To date, the stabilization of transient states is still in the preliminary stage, particularly with respect to molecular electronic solids. Here, a dynamic and cooperative phase in potassium-7,7,8,8-tetracyanoquinodimethane (K-TCNQ) with the control of pulsed electromagnetic excitation is demonstrated. Simultaneous dynamic and coherent lattice perturbation with 8 ns pulsed laser (532 nm, 15 MW cm -2 , 10 Hz) in such a molecular electronic crystal initiates a stable long-lived (over 400 days) conducting paramagnetic state (≈42 Ωcm), showing the charge–spin bistability over a broad temperature range from 2 to 360 K. Comprehensive noise spectroscopy, in situ high-pressure measurements, electron spin resonance (ESR), theoretical model, and scanning tunneling microscopy/spectroscopy (STM/STS) studies provide further evidence that such a transition is cooperative, requiring a dedicated charge–spin–lattice decoupling to activate and subsequently stabilize nonequilibrium phase. Finally, the cooperativity triggered by ultrahigh-strain-rate (above 10 6 s -1 ) pulsed excitation offers a collective control toward the generation and stabilization of strongly correlated electronic and magnetic orders in molecular electronic solids and offers unique electro-magnetic phases with technological promises.

36 MATERIALS SCIENCE↗

Ultrafast laser-driven many-body dynamics and Kondo coherence collapse

Ultrafast laser pulse has provided a systematic way to inspect the dynamics of electrons in condensed matter systems. In this paper, by means of time-dependent density matrix renormalization group, we study an ultrafast laser-driven Kondo lattice model, in which conduction electrons are strongly coupled with magnetically local moments. The single-particle spectral function due to strong correlation effects and photon emission in the nonequilibrium states under laser driving are calculated. We find laser field excited collective doublon-hole pairs and an associated transient melting of Kondo coherence phase, signifying the collapse of Kondo energy gap. Moreover, we show that the photon emission, induced by a strong laser field, exhibits a different intensity characteristics than in the equilibrium Kondo insulator, which could be explained by the Kondo collapse and related suppression of both intraband and interband contribution in Kondo melting liquid. This theoretical insight is accessible with time- and angle-resolved photoemission spectroscopy and high-harmonic generation spectroscopy and will stimulate the investigation of nonequilibrium dynamics and nonlinear phenomenon in heavy fermion systems.

36 MATERIALS SCIENCE↗

Concurrent probing of electron-lattice dephasing induced by photoexcitation in 1 T -TaSeTe using ultrafast electron diffraction

It has been technically challenging to concurrently probe the electrons and the lattices in materials during nonequilibrium processes, allowing their correlations to be determined. Here, in a single set of ultrafast electron diffraction patterns taken on the charge-density-wave (CDW) material 1T-TaSeTe, we discover a temporal shift in the diffraction intensity measurements as a function of scattering angle. With the help of dynamic models and theoretical calculations, we show that the ultrafast electrons probe both the valence-electron and lattice dynamic processes, resulting in the temporal shift measurements. Our results demonstrate unambiguously that the CDW is not merely a result of the periodic lattice deformation ever present in 1T-TaSeTe but has significant electronic origin. This method demonstrates an approach for studying many quantum effects that arise from electron-lattice dephasing in molecules and crystals for next-generation devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗

Quantum simulation of nonequilibrium dynamics and thermalization in the Schwinger model

In this work, we present simulations of nonequilibrium dynamics of quantum field theories on digital quantum computers. As a representative example, we consider the Schwinger model, a (1+1)-dimensional U(1) gauge theory, coupled through a Yukawa-type interaction to a thermal environment described by a scalar field theory. We use the Hamiltonian formulation of the Schwinger model discretized on a spatial lattice. With the thermal scalar fields traced out, the Schwinger model can be treated as an open quantum system and its real-time dynamics are governed by a Lindblad equation in the Markovian limit. The interaction with the environment ultimately drives the system to thermal equilibrium. In the quantum Brownian motion limit, the Lindblad equation is related to a field theoretical Caldeira-Leggett equation. By using the Stinespring dilation theorem with ancillary qubits, we perform studies of both the nonequilibrium dynamics and the preparation of a thermal state in the Schwinger model using IBM's simulator and quantum devices. The real-time dynamics of field theories as open quantum systems and the thermal state preparation studied here are relevant for a variety of applications in nuclear and particle physics, quantum information and cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamical properties of a driven dissipative dimerized $S$ =1/2 chain

We consider the dynamical properties of a gapped quantum spin system coupled to the electric field of a laser, which drives the resonant excitation of specific phonon modes that modulate the magnetic interactions. Here, we deduce the quantum master equations governing the time-evolution of both the lattice and spin sectors, by developing a Lindblad formalism with bath operators providing an explicit description of their respective phonon-mediated damping terms. We investigate the nonequilibrium steady states (NESS) of the spin system established by a continuous driving, delineating parameter regimes in driving frequency, damping, and spin-phonon coupling for the establishment of physically meaningful NESS and their related nontrivial properties. Focusing on the regime of generic weak spin-phonon coupling, we characterize the NESS by their frequency and wave-vector content, explore their transient and relaxation behavior, and discuss the energy flow, the system temperature, and the critical role of the type of bath adopted. Our study lays a foundation for the quantitative modeling of experiments currently being designed to control coherent many-body spin states in quantum magnetic materials.

1-dimensional spin chains↗

Elastic forces drive nonequilibrium pattern formation in a model of nanocrystal ion exchange

The widely used process of nanocrystal ion exchange operates out of thermodynamic equilibrium and can require mixing components of varying sizes. Here we use theory and computer simulation to study a simple model which captures these two basic features of ion exchange reactions. We show that a strong driving force for exchange among different-sized species creates nonequilibrium patterns within model nanocrystals. We further demonstrate that such patterns can be thermodynamically stable within core/shell nanocrystals. These results help us understand the heterostructures formed in ion-exchanged nanocrystals and suggest strategies for leveraging elasticity to design patterns in nanoscale materials. Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal’s surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition—from modulated to unstructured—as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally patterned nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cold equation of state in a strong magnetic field - Effects of inverse beta-decay

The influence of a high magnetic field (B is greater than 10 exp 12 G) on the degenerate matter equation of state appropriate to a neutron star is studied. The regime dominated by relativistic electrons up to the neutron drip density is highlighted. The equilibrium matter composition and equation of state, allowing for inverse beta-decay. Two different equilibrium models are determined: an ideal neutron-proton-electron (npe) gas and the more realistic model of Baym, Pethick, and Sutherland (1971) consisting of a Coulomb lattice of heavy nuclei embedded in an electron gas. For a sufficiently high field strength, the magnetic field has an appreciable effect, changing the adiabatic index of the matter and the nuclear transition densities. The influence of a strong field on some simple nonequilibrium processes, including beta-decay and inverse beta-decay (electron capture) is also considered. The effects produced by the magnetic field are mainly due to the changes in the transverse electron quantum orbits and the allowed electron phase space induced by the field.

Lai, Dong↗

Heat current fluctuations and anomalous transport in low-dimensional carbon lattices

Molecular dynamics simulations and methods of importance sampling are used to study the heat transport of low-dimensional carbon lattices. For both carbon nanotubes and graphene sheets, heat transport is found to be anomalous, violating Fourier’s law of conduction with a system size dependent thermal conductivity and concomitant nonlinear temperature profiles. For carbon nanotubes, the thermal conductivity is found to increase as the square root of the length of the nanotube, while for graphene sheets the thermal conductivity is found to increase as the logarithm of the length of the sheet over the system sizes considered. Furthermore, the particular length dependence and nonlinear temperature profiles place carbon lattices into a universality class with nonlinear lattice models, and suggest that heat transport through carbon nanostructures is better described by a Levy walk rather than simple diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonequilibrium Crystal Growth During Laser Annealing of Silicon

The Systematic studies of dopant (Group III, IV, or V) incorporation into silicon as a result of pulsed-laser annealing are reviewed. At low concentrations, dopant profiles measured after laser annealing can be compared to model calculations for redistribution by liquid-phase diffusion in order to determine the interfacial distribution coefficient K('). At high concentrations there is a limit (C sub s) to dopant incorporation into substitutional lattice sites. The mechanisms which limit dopant incorporation are discussed as well as the dependence of K(') and (max) on crystal orientation. Finally, measurements of interface temperatures, temperature profiles, and thermal gradients during rapid solidification are summarized.

White, C. W.↗

Random Quantum Circuits

Quantum circuits—built from local unitary gates and local measurements—are a new playground for quantum many-body physics and a tractable setting to explore universal collective phenomena far from equilibrium. These models have shed light on longstanding questions about thermalization and chaos, and on the underlying universal dynamics of quantum information and entanglement. In addition, such models generate new sets of questions and give rise to phenomena with no traditional analog, such as dynamical phase transitions in quantum systems that are monitored by an external observer. Quantum circuit dynamics is also topical in view of experimental progress in building digital quantum simulators that allow control of precisely these ingredients. Randomness in the circuit elements allows a high level of theoretical control, with a key theme being mappings between real-time quantum dynamics and effective classical lattice models or dynamical processes. Many of the universal phenomena that can be identified in this tractable setting apply to much wider classes of more structured many-body dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗