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 145 records · Page 8

Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model

In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.

97 MATHEMATICS AND COMPUTING↗

Universal tripartite entanglement signature of ungappable edge states

Gapped two-dimensional topological phases can feature ungappable edge states which are robust even in the absence of protecting symmetries. In this Letter we argue that a multipartite entanglement measure recently proposed in the context of holography, the Markov gap, provides a universal diagnostic of ungappable edge states. Defined as a difference of the reflected entropy and mutual information h(A:B)=S R (A:B)-I(A:B) between two parties, we argue that for A,B being adjacent subregions in the bulk h=$\frac{c+}{3}$ ln 2, where c + is the minimal total central charge of the boundary theory. As evidence, we prove that h=0 for string-net models and numerically verify that h=$\frac{|C|}{3}$ ln 2 for a Chern-C insulator. Our Letter establishes a unique bulk entanglement criteria for the presence of a conformal field theory on the boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modulating effect of evanescent waves on thin film growth

Atomic-scale smooth thin films are keys to successful integration and proper function of many multilayer-structured devices. However, the intrinsic islandlike growth mode prevents human beings from realizing such ultrasmooth films of many important functional materials. To solve this problem, we propose a negative entropy-infusing method that employs evanescent waves to enhance the downward interlayer diffusion of adatoms and thus transform the islandlike growth mode to the layer-by-layer growth mode. In this work, the formulas of the optical force and the lowered diffusion barrier were derived, and the application of this theory on a simplified example demonstrates significantly improved surface morphologies through numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations

We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.

Bailo, Rafael (ORCID:0000000180183799)↗

Probing Postmeasurement Entanglement without Postselection

We study the problem of observing quantum collective phenomena emerging from large numbers of measurements. These phenomena are difficult to observe in conventional experiments because, in order to distinguish the effects of measurement from dephasing, it is necessary to postselect on sets of measurement outcomes with Born probabilities that are exponentially small in the number of measurements performed. An unconventional approach, which avoids this exponential “postselection problem”, is to construct cross-correlations between experimental data and the results of simulations on classical computers. However, these cross-correlations generally have no definite relation to physical quantities. We first show how to incorporate classical shadows into this framework, thereby allowing for the construction of quantum information-theoretic cross-correlations. We then identify cross-correlations that both upper and lower bound the measurement-averaged von Neumann entanglement entropy, as well as cross-correlations that lower bound the measurement-averaged purity and entanglement negativity. These bounds show that experiments can be performed to constrain postmeasurement entanglement without the need for postselection. To illustrate our technique, we consider how it could be used to observe the measurement-induced entanglement transition in Haar-random quantum circuits. We use exact numerical calculations as proxies for quantum simulations and, to highlight the fundamental limitations of classical memory, we construct cross-correlations with tensor-network calculations at finite bond dimension. Our results reveal a signature of measurement-induced criticality that can be observed using a quantum simulator in polynomial time and with polynomial classical memory. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗

Sample-efficient verification of continuously-parameterized quantum gates for small quantum processors

Most near-term quantum information processing devices will not be capable of implementing quantum error correction and the associated logical quantum gate set. Instead, quantum circuits will be implemented directly using the physical native gate set of the device. These native gates often have a parameterization (e.g., rotation angles) which provide the ability to perform a continuous range of operations. Verification of the correct operation of these gates across the allowable range of parameters is important for gaining confidence in the reliability of these devices. In this work, we demonstrate a procedure for sample-efficient verification of continuously-parameterized quantum gates for small quantum processors of up to approximately 10 qubits. This procedure involves generating random sequences of randomly-parameterized layers of gates chosen from the native gate set of the device, and then stochastically compiling an approximate inverse to this sequence such that executing the full sequence on the device should leave the system near its initial state. We show that fidelity estimates made via this technique have a lower variance than fidelity estimates made via cross-entropy benchmarking. This provides an experimentally-relevant advantage in sample efficiency when estimating the fidelity loss to some desired precision. We describe the experimental realization of this technique using continuously-parameterized quantum gate sets on a trapped-ion quantum processor from Sandia QSCOUT and a superconducting quantum processor from IBM Q, and we demonstrate the sample efficiency advantage of this technique both numerically and experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Universal signatures of Dirac fermions in entanglement and charge fluctuations

We study the entanglement entropy (EE) and charge fluctuations in models where the low-energy physics is governed by massless Dirac fermions. We focus on the response to flux insertion which, for the EE, is widely assumed to be universal, i.e., independent of the microscopic details. We provide an analytical derivation of the EE and charge fluctuations for the seminal example of graphene, using the dimensional reduction of its tight-binding model to the one-dimensional Su-Schrieffer-Heeger model. Our asymptotic expression for the EE matches the conformal field theory prediction. We show that the charge variance has the same asymptotic behavior, up to a constant prefactor. To check the validity of universality arguments, we numerically consider several models, with different geometries and numbers of Dirac cones, and either for strictly two-dimensional models or for a gapless surface mode of three-dimensional topological insulators. We also show that the flux response does not depend on the entangling surface geometry as long as it encloses the flux. Finally we consider the universal corner contributions to the EE. We show that in the presence of corners, the Kitaev-Preskill subtraction scheme provides nonuniversal, geometry-dependent results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Feasibility of measuring the speed of sound of the quark-gluon plasma from the multiplicity and mean 𝑝 𝑇 of ultracentral heavy-ion collisions

The mean transverse momentum ⟨𝑝 𝑇 ⟩ of hadrons has been observed experimentally and in numerical simulations to have a power-law dependence on the hadronic multiplicity 𝑁 in ultracentral relativistic heavy-ion collisions: ⟨𝑝 𝑇 ⟩∝𝑁 𝑏 UC . It has been put forward that this exponent 𝑏 UC is the speed of sound of quark-gluon plasma measured at a temperature determined from ⟨𝑝 𝑇 ⟩. We study step by step the connection between (i) the energy and entropy of hydrodynamic simulations and (ii) experimentally measurable observables. We show that an argument based on energy and entropy should yield an exponent equal to the pressure over energy density 𝑃/ɛ, rather than the speed of sound 𝑐$_s^2$; however, we also observe that ⟨𝑝 𝑇 ⟩ and 𝑁 are not sufficiently accurate proxies for the energy and entropy to make this possible in practice. From simulations, we find that the exponent 𝑏 UC is significantly different whether the “effective volume” is strictly constant or not, a condition that cannot be enforced experimentally. Additional tests using a modified equation of state find that the exponent 𝑏 UC exhibits a variable degree of correlations with the speed of sound and with 𝑃/ɛ, but is not an accurate measurement of either quantity in general.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗

Deformation and failure of the CrCoNi medium-entropy alloy subjected to extreme shock loading

The extraordinary work hardening ability and fracture toughness of the face-centered cubic (fcc) high-entropy alloys render them ideal candidates for many structural applications. Here, the deformation and failure mechanisms of an equiatomic CrCoNi medium-entropyalloy (MEA) were investigated by powerful laser-driven shock experiments. Multiscale characterization demonstrates that profuse planar defects including stacking faults, nanotwins, and hexagonal nanolamella were generated during shock compression, forming a three-dimensional network. During shock release, the MEA fractured by strong tensile deformation and numerous voids was observed in the vicinity of the fracture plane. High defect populations, nanorecrystallization, and amorphization were found adjacent to these areas of localized deformation. Molecular dynamics simulations corroborate the experimental results and suggest that deformation-induced defects formed before void nucleation govern the geometry of void growth and delay their coalescence. Our results indicate that the CrCoNi-based alloys are impact resistant, damage tolerant, and potentially suitable in applications under extreme conditions.

36 MATERIALS SCIENCE↗

Towards absolutely stable ergodicity breaking in two and three dimensions

We propose physically reasonable systems capable of avoiding ergodicity at infinite time in the thermodynamic limit, even with generic perturbations and when coupled to a heat bath. In two dimensions, the rainbow loop soup has (stretched) exponentially numerous absolutely stable nonergodic states with diverging energy but vanishing energy density. In three dimensions the rainbow membrane soup has (stretched) exponentially numerous nonergodic states with diverging energy barriers, leading to infinite-time robust ergodicity breaking that even survives coupling to a nonzero temperature heat bath. We describe our results in the language of exact emergent symmetries and demonstrate how the systems avoid common instabilities. Furthermore, our construction naturally connects to quantum dimer models, topologically ordered systems, the group word construction, and Hamiltonians whose low-energy eigenstates exhibit anomalous entanglement entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Elaboration of Metallic Materials by SPS: Processing, Microstructures, Properties, and Shaping

After a few decades of increasing interest, spark plasma sintering (SPS) has now become a mature powder metallurgy technique, which allows assessing its performances toward fabricating enhanced materials. Here, the case of metals and alloys will be presented. The main advantage of SPS lies in its rapid heating capability enabled by the application of high intensity electric currents to a metallic powder. This presents numerous advantages balanced by some limitations that will be addressed in this review. The first section will be devoted to sintering issues, with an emphasis on the effect of the electric current on the densification mechanisms. Then, typical as-SPS microstructures and properties will be presented. In some cases, they will be compared with that of materials processed by conventional techniques. As such, examples of nanostructured materials, intermetallics, metallic glasses, and high entropy alloys, will be presented. Finally, the implementation of SPS as a technique to manufacture complex, near-net shape industrial parts will be discussed.

36 MATERIALS SCIENCE↗

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗

A simple introduction to the SiMPL method for density-based topology optimization

We introduce a novel method for solving density-based topology optimization problems: Sigmoidal Mirror descent with a Projected Latent variable (SiMPL). The SiMPL method (pronounced as “the simple method”) optimizes a design using only first-order derivative information of the objective function. The bound constraints on the density field are enforced with the help of the (negative) Fermi–Dirac entropy, which is also used to define a non-symmetric distance function called a Bregman divergence on the set of admissible designs. This Bregman divergence leads to a simple update rule that is further simplified with the help of a so-called latent variable. Because the SiMPL method involves discretizing the latent variable, it produces a sequence of pointwise-feasible iterates, even when high-order finite elements are used in the discretization. Numerical experiments demonstrate that the method outperforms other popular first-order optimization algorithms. In conclusion, to outline the general applicability of the technique, we include examples with (self-load) compliance minimization and compliant mechanism optimization problems.

Calculus of Variations and Optimization↗

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗

Finite-temperature many-body perturbation theory in the grand canonical ensemble

A finite-temperature many-body perturbation theory is presented, which expands in power series the electronic grand potential, chemical potential, internal energy, and entropy on an equal footing. Sum-over-states and sum-over-orbitals analytical formulas for the second-order perturbation corrections to these thermodynamic properties are obtained in a time-independent, nondiagrammatic, algebraic derivation, relying on the sum rules of the Hirschfelder–Certain degenerate perturbation energies in a degenerate subspace as well as nine algebraic identities for the zeroth-order thermal averages of one- through four-indexed quantities and products thereof. Lastly, they reproduce numerically exactly the benchmark data obtained as the numerical derivatives of the thermal-full-configuration-interaction results for a wide range of temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flux Ropes, Turbulence, and Collisionless Perpendicular Shock Waves: High Plasma Beta Case

With the onset of solar maximum and the expected increased prevalence of interplanetary shock waves, Parker Solar Probe is likely to observe numerous shocks in the next few years. An outstanding question that has received surprisingly little attention has been how turbulence interacts with collisionless shock waves. Turbulence in the supersonic solar wind is described frequently as a superposition of a majority 2D and a minority slab component. Here we formulate a collisional perpendicular shock-turbulence transmission problem in a way that enables investigation of the interaction and transmission of quasi-perpendicular fluctuations such as magnetic flux ropes/islands and vortices as well as entropy and acoustic modes in the large plasma beta regime. We focus on the transmission of an upstream spectrum of these modes, finding that the downstream spectral amplitude is typically increased significantly (a factor of 10 or more), and that the upstream spectral index of the inertial range, and indeed the general spectral shape, is unchanged for the downstream magnetic variance, kinetic energy, and density variance. A comparison of the theoretically predicted downstream magnetic variance, kinetic energy, and density variance spectra with those observed at 1, 5, and 84 au by Wind, Ulysses, and Voyager 2 shows excellent agreement. The overall theoretically predicted characteristics of the transmission of turbulence across shocks observed in the solar wind appear to be largely consistent with recent observational studies by Pitňa et al. and Borovsky.

79 ASTRONOMY AND ASTROPHYSICS↗