Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical entropy”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A random unitary circuit model for black hole evaporation

Inspired by the Hayden-Preskill protocol for black hole evaporation, we consider the dynamics of a quantum many-body qudit system coupled to an external environment, where the time evolution is driven by the continuous limit of certain 2-local random unitary circuits. We study both cases where the unitaries are chosen with and without a conserved U(1) charge and focus on two aspects of the dynamics. First, we study analytically and numerically the growth of the entanglement entropy of the system, showing that two different time scales appear: one is intrinsic to the internal dynamics (the scrambling time), while the other depends on the system-environment coupling. In the presence of a U(1) conserved charge, we show that the entanglement follows a Page-like behavior in time: it begins to decrease in the middle stage of the “evaporation”, and decreases monotonically afterwards. Second, we study the time needed to retrieve information initially injected in the system from measurements on the environment qudits. Based on explicit numerical computations, we characterize such time both when the retriever has control over the initial configuration or not, showing that different scales appear in the two cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

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

Energy-entropy relation for asymptotically Lifshitz spacetimes with universal horizons

We numerically solve for 2 + 1 asymptotically Lifshitz universal horizon solutions in Horava-Lifshitz gravity for dynamical exponents z = 2 through z = 16 . We find that for all z there is a thermodynamical first law. Furthermore, we find that the energy-entropy relation or Smarr formula expected for a thermal state in a two-dimensional Lifshitz field theory, E = 2 z + 2 T S , is also satisfied for universal horizons, including the correct z scaling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2022

A topological kagome magnet in high entropy form

AbstractTopological kagome magnets RMn 6 Sn 6 (R = rare earth element) attract numerous interests due to their non-trivial band topology and room-temperature magnetism. Here, we report a high entropy version of kagome magnet, (Gd 0.38 Tb 0.27 Dy 0.20 Ho 0.15 )Mn 6 Sn 6 . Such a high entropy material exhibits multiple spin reorientation transitions, which is not seen in all the related parent compounds and can be understood in terms of competing magnetic interactions enabled by high entropy. Furthermore, we also observed an intrinsic anomalous Hall effect, indicating that the high entropy phase preserves the non-trivial band topology. These results suggest that high entropy may provide a route to engineer the magnetic structure and expand the horizon of topological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2025

Electronic density of states as the descriptor of elastic bond strength, ductility, and local lattice distortion in BCC refractory alloys

Although electronic density of states (DOS) is fundamental to materials properties, its general relationship to mechanical properties of alloys is not well established. In this paper, using density functional theory (DFT) calculations, we show that the electronic occupancy at the Fermi level, N(E f ), obtained from DOS is a key descriptor of alloy strength and ductility. Our comprehensive analysis of numerous body centered cubic (BCC) refractory high entropy alloys (RHEAs) shows an overwhelming correlation that low N(E f ) indicates strong bonds that have high stiffness resulting in high elastic constants. High bond stiffness indicates presence of covalent nature of bonds that are directional in nature resulting in resistance to deformation leading to high bulk (B) and shear (G) moduli. Consequently, N(E f ) provides a direct correlation to the tendency of alloy ductility evidenced in the Pugh ratio (G/B). As stiffer bonds result in lower local lattice distortion (LLD), N(E f ) are LLD are also found to be corelated which opens up a correlation to solid solution strengthening and yield strength. Thus, this work unveils fundamental correlations between N(E f ) and (1) elastic bond strength, (2) ductility, and (3) LLD. These correlations open opportunities for the design of high strength high ductile RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Entropy–Preserving and Entropy–Stable Relaxation IMEX and Multirate Time–Stepping Methods

In this work, we propose entropy-preserving and entropy-stable partitioned Runge–Kutta (RK) methods. In particular, we extend the explicit relaxation Runge–Kutta methods to IMEX–RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum entropy as a harbinger of factorizability

Deeply inelastic scattering (DIS) is a powerful probe for investigating the QCD structure of hadronic matter and testing the standard model (SM). DIS can be described through QCD factorization theorems which separate contributions to the scattering interaction arising from disparate scales — e.g ., with nonperturbative matrix elements associated with long distances and a perturbative hard scattering kernel applying to short-distance parton-level interactions. The fundamental underpinnings of factorization may be recast in the quantum-theoretic terms of entanglement, (de)coherence, and system localization in a fashion which sheds complementary light on the dynamics at work in DIS from QCD bound states. In this Letter, we propose and quantitatively test such a quantum-information theoretic approach for dissecting factorization in DIS and its domain of validity; we employ metrics associated with quantum entanglement such as a differential quantum entropy and associated Kullback-Leibler (KL) divergences in numerical tests. We deploy these methods on an archetypal quark-spectator model of the proton, for which we monitor quantum decoherence in DIS as underlying model parameters are varied. On this basis, we demonstrate quantitatively how factorization-breaking effects may be imprinted on quantum entropies in a kinematic regime where leading-twist factorization increasingly receives large corrections from finite- Q 2 effects; our findings suggest potential applications of quantum simulation to QCD systems and their interactions.

Deep inelastic scattering↗
DOE OSTI2022

Fluctuations of subsystem entropies at late times

We study the fluctuations of subsystem entropies in closed quantum many-body systems after thermalization. Using a combination of analytics and numerics for both random quantum circuits and Hamiltonian dynamics, we find that the statistics of such entropy fluctuations is drastically different than in the classical setting. For instance, shortly after a system thermalizes, the probability of entropy fluctuations for a subregion is suppressed in the dimension of the Hilbert space of the complementary subregion. This suppression becomes increasingly stringent as a function of time, ultimately depending on the exponential of the Hilbert space dimension, until extremely late times when the amount of suppression saturates. We also use our results to estimate the total number of rare fluctuations at large timescales. We find that the “Boltzmann brain” paradox is largely ameliorated in quantum many-body systems, in contrast with the classical setting.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Topological obstructed atomic limit insulators by annihilating Dirac fermions

We show that annihilating a pair of Dirac fermions implies a topological transition from the critical semimetallic phase to an obstructed atomic limit insulator phase instead of a trivial insulator. This is shown to happen because of branch cuts in the phase of the wave functions, leading to nontrivial Zak phase along certain directions. To this end, we study their ${\mathbb{Z}}_{2}$ invariant and also study the phase transition using entanglement entropy. We use low-energy Hamiltonians and numerical result from model systems to show this effect. These transitions are observed in realistic materials, including strained graphene and buckled honeycomb group V (Sb/As).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-driven Minimum Entropy Control for Stochastic Nonlinear Systems using the Cumulant-Generating Function

Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.

42 ENGINEERING↗
DOE OSTI2024

Mixed Enthalpy–Entropy Descriptor for the Rational Design of Synthesizable High-Entropy Materials Over Vast Chemical Spaces

The practically unlimited high-dimensional composition space of high-entropy materials (HEMs) has emerged as an exciting platform for functional material design and discovery. However, the identification of stable and synthesizable HEMs and robust design rules remains a daunting challenge. Here, we propose a mixed enthalpy–entropy descriptor (MEED) that enables highly efficient, robust, high-throughput prediction of synthesizable HEMs across vast chemical spaces from first-principles. The MEED is based on two parameters: the relative formation enthalpy with respect to the most stable competing compound and the spread of the point-defect formation energy spectrum. The former measures the relative synthesizability of an HEM to its most stable competing phase, going beyond the conventional thermodynamic understanding. Further, the latter gauges the relative entropy forming ability of an HEM, entailing no sampling over numerous alloy configurations. By applying the MEED to two structurally distinct representative material systems (i.e., 3D rocksalt carbides and 2D layered sulfides), we not only successfully identify all experimentally reported HEMs within these systems but also reveal a cutoff criterion for assessing their relative synthesizability within each system. By the MEED, tens of new high-entropy carbides and 2D high-entropy sulfides are also predicted, which have the potential for a wide variety of applications such as coating in aerospace devices, energy conversion and storage, and flexible electronics.

36 MATERIALS SCIENCE↗