Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “entropy functional”

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 775 records · Page 43

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗

Fluctuations in the entropy of Hawking radiation

We use the gravitational path integral (GPI) to compute the fluctuations of the Hawking radiation entropy around the Page curve in a two-dimensional model introduced by Penington Before the Page time, we find that δ S = e − S / 2 , where S is the black hole entropy. This result agrees with the Haar-averaged entropy fluctuations in a bipartite system. After the Page time, we find that δ S ∼ e − S , up to a prefactor that depends logarithmically on the width of the microcanonical energy window. This is not symmetric under exchange of subsystem sizes and so does not agree with the Haar average for a subsystem of fixed Hilbert space dimension. The discrepancy can be attributed to the fact that the black hole Hilbert space dimension is not fixed by the state preparation: even in a microcanonical ensemble with a top-hat smearing function, the GPI yields an additive fluctuation in the number of black hole states. This result, and the fact that the Page curve computed by the GPI is smooth, all point towards an ensemble interpretation of the GPI. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytic continuation of noisy data using Adams Bashforth residual neural network

We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). Additionally, the advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.

97 MATHEMATICS AND COMPUTING↗

Controlling the corrosion resistance of multi-principal element alloys

Multi-principal element alloys (MPEAs) offer the possibility of many degrees of freedom in the choice of alloying elements to produce either single phase solid solutions or more complex multiphase microstructures. Large ranges of material properties have been observed for MPEAs and mastery of the selection of elements and their compositions can enable novel combinations of properties not possible in traditional alloys. From the aqueous corrosion perspective, optimization of phase stability, control of heterogeneities, passive film identity and its protectiveness, as well as substrate properties such as metal-metal bond strength and activation energy associated with dissolution, can all be controlling factors governed by alloy composition and structure. These factors can mediate the electrochemical reactions controlling spontaneous corrosion. The quest for superior properties based on well-informed element choice is suggested as a path forward guiding MPEA formulations for corrosion performance. However, gaps in fundamental knowledge exist regarding (a) the specific functions of each element, (b) the behavior of elements in unusual combinations, and (c) the formation of complex protective oxides. These issues currently hold back progress in optimization of corrosion properties.

Scully, John R.↗

Linking electronic structure calculations to generalized stacking fault energies in multicomponent alloys

Abstract The generalized stacking fault energy is a key ingredient to mesoscale models of dislocations. Here we develop an approach to quantify the dependence of generalized stacking fault energies on the degree of chemical disorder in multicomponent alloys. We introduce the notion of a “configurationally-resolved planar fault” (CRPF) energy and extend the cluster expansion method from alloy theory to express the CRPF as a function of chemical occupation variables of sites surrounding the fault. We apply the approach to explore the composition and temperature dependence of the unstable stacking fault energy (USF) in binary Mo–Nb alloys. First-principles calculations are used to parameterize a formation energy and CRPF cluster expansion. Monte Carlo simulations show that the distribution of USF energies is significantly affected by chemical composition and temperature. The formalism is broadly applicable to arbitrary crystal structures and alloy chemistries and will enable the development of rigorous models for deformation mechanisms in high-entropy alloys.

Natarajan, Anirudh Raju (ORCID:0000000202529480)↗

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and here we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

36 MATERIALS SCIENCE↗

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗

Information Theoretic Approaches to Rapid Discovery of Relationships in Large Climate Data Sets

Mutual information as the asymptotic Bayesian measure of independence is an excellent starting point for investigating the existence of possible relationships among climate-relevant variables in large data sets, As mutual information is a nonlinear function of of its arguments, it is not beholden to the assumption of a linear relationship between the variables in question and can reveal features missed in linear correlation analyses. However, as mutual information is symmetric in its arguments, it only has the ability to reveal the probability that two variables are related. it provides no information as to how they are related; specifically, causal interactions or a relation based on a common cause cannot be detected. For this reason we also investigate the utility of a related quantity called the transfer entropy. The transfer entropy can be written as a difference between mutual informations and has the capability to reveal whether and how the variables are causally related. The application of these information theoretic measures is rested on some familiar examples using data from the International Satellite Cloud Climatology Project (ISCCP) to identify relation between global cloud cover and other variables, including equatorial pacific sea surface temperature (SST), over seasonal and El Nino Southern Oscillation (ENSO) cycles.

Knuth, Kevin H.↗

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong↗

Nonharmonic contributions to the high-temperature phonon thermodynamics of Cr

Here, phonon densities of states (DOSs) of body-centered cubic chromium were measured by time-of-flight inelastic neutron scattering at temperatures up to 1493 K. Density functional theory calculations with both quasiharmonic (QH) and anharmonic (AH) methods were performed at temperatures above the Néel temperature. Features in the phonon DOSs decrease in energy (soften) substantially with temperature. A Born–von Kármán analysis using fits to the experimental DOSs reveals a softening of almost 17% of the high-transverse phonon branch between 330 and 1493 K. The low-transverse branch changes by approximately half this amount. The AH calculations capture the observed behavior of the two transverse phonon branches, but the QH calculations give some inverted trends. Vibrational entropies from phonons and electrons are obtained, and their sum is in excellent agreement with the entropy of chromium obtained by calorimetry, indicating that above 330 K, no explicit temperature-dependent magnetic contributions are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the lapse contour in the gravitational path integral

The gravitational path integral is usually implemented with a covariant action by analogy with other gauge field theories, but the gravitational case is different in important ways. A key difference is that the integrand has an essential singularity, which occurs at zero lapse where the spacetime metric degenerates. The lapse integration contour required to impose the local time reparametrization constraints must run from − ∞ to + ∞ , yet must not pass through zero. This raises the question: for an application—such as a partition function—where the constraints should be imposed, what is the correct integration contour, and why? We study that question by starting with the reduced phase space path integral, which involves no essential singularity. We observe that if the momenta are to be integrated before the lapse, to obtain a configuration space path integral, the lapse contour should pass below the origin in the complex lapse plane. This contour is also consistent with the requirement that quantum field fluctuation amplitudes have the usual short distance vacuum form, and with obtaining the Bekenstein-Hawking horizon entropy from a Lorentzian path integral. Published by the American Physical Society 2025

Banihashemi, Batoul (ORCID:0000000228679209)↗

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗

Thickness dependent OER electrocatalysis of epitaxial thin film of high entropy oxide

High entropy oxides (HEOs), which contain multiple elements in the same crystallographic site, are a promising platform for electrocatalysis in oxygen evolution reaction (OER). Investigating these materials in epitaxial thin film form expands the possibility of tuning OER activity by several means, which are not realizable in polycrystalline samples. To date, very few such studies have been reported. In this work, the OER activity of single-crystalline thin films of (La 0.2 Pr 0.2 Nd 0.2 Sm 0.2 Eu 0.2 )NiO 3 , grown on NdGaO 3 substrates have been investigated in 0.1 M KOH electrolyte as a function of film thickness. The OER activity increases with the thickness of the film. X-ray absorption spectroscopy measurements find an increase in Ni $d$-O $p$ covalency and a decrease in charge transfer energy with the increase in film thickness. These facilitate higher charge transfer between Ni and surface adsorbates, increasing OER activity. However, the OER process leads to excessive leaching of thicker films and the OER activity of a 75 unit cell thick film is found to be optimal in the present study. In conclusion, this work demonstrates that the thickness of perovskite oxides can be used as a parameter to enhance OER activity.

36 MATERIALS SCIENCE↗

The Application B3LYP to Large Systems

The application of density functional theory (DFT), using the B3LYP functional, to a series of chemical problems is described. The first involves the calculation of silica-adsorbate bond energies, including both chemical bonds and weak hydrogen bonding. The calculation of vibrational frequencies for large organic systems is discussed. For the closed shell neutral systems, the B3LYP results are similar to the self- consistent- field results, however, for the positive ions, only the B3LYP level of theory is accurate and sufficiently inexpensive to allow the study of large systems. The final application involves the calculation of successive metal-ligand bond energies. The B3LYP bond energies and entropies are shown to be in good agreement with experiment.

Bauschlicher, Charles W.↗

Towards Pareto optimal high entropy hydrides via data-driven materials discovery

The ability to rapidly screen material performance in the vast space of high entropy alloys is of critical importance to efficiently identify optimal hydride candidates for various use cases. Given the prohibitive complexity of first principles simulations and large-scale sampling required to rigorously predict hydrogen equilibrium in these systems, we turn to compositional machine learning models as the most feasible approach to screen on the order of tens of thousands of candidate equimolar high entropy alloys (HEAs). Critically, we show that machine learning models can predict hydride thermodynamics and capacities with reasonable accuracy (e.g. a mean absolute error in desorption enthalpy prediction of ~5 kJ mol H 2 –1 ) and that explainability analyses capture the competing trade-offs that arise from feature interdependence. We can therefore elucidate the multi-dimensional Pareto optimal set of materials, i.e., where two or more competing objective properties can't be simultaneously improved by another material. This provides rapid and efficient down-selection of the highest priority candidates for more time-consuming density functional theory investigations and experimental validation. Various targets were selected from the predicted Pareto front (with saturation capacities approaching two hydrogen per metal and desorption enthalpy less than 60 kJ mol H 2 –1 ) and were experimentally synthesized, characterized, and tested amongst an international collaboration group to validate the proposed novel hydrides. Finally, additional top-predicted candidates are suggested to the community for future synthesis efforts, and we conclude with an outlook on improving the current approach for the next generation of computational HEA hydride discovery efforts.

36 MATERIALS SCIENCE↗

Free fermion entanglement with a semitransparent interface: the effect of graybody factors on entanglement islands

We study the entanglement entropy of free fermions in 2d in the presence of a partially transmitting interface that splits Minkowski space into two half-spaces. We focus on the case of a single interval that straddles the defect, and compute its entanglement entropy in three limits: Perturbing away from the fully transmitting and fully reflecting cases, and perturbing in the amount of asymmetry of the interval about the defect. Using these results within the setup of the Poincaré patch of AdS _2 2 statically coupled to a zero temperature flat space bath, we calculate the effect of a partially transmitting AdS _2 2 boundary on the location of the entanglement island region. The partially transmitting boundary is a toy model for black hole graybody factors. Our results indicate that the entanglement island region behaves in a monotonic fashion as a function of the transmission/reflection coefficient at the interface.

Kruthoff, Jorrit↗

Stability of Pseudobrookite-Type Titanium Oxides

Orthorhombic, (Bbmm), (Al, Fe, Cr, Ti)(sub 2) TiO5-(Mg, Fe)Ti2O5 solid solutions (pseudobrookites, s.l.) are found either as an oxidation product of ilmenite and/or spinel or a primary crystallizing phase in igneous and metamorphic rocks on Earth (e.g., basalt flows, crustal and mantle xenoliths, hornfels), and basaltic rocks on the Moon. Moreover, orthorhombic oxides are often part of the crystalline matrix in glass/ceramics with useful applications, and play a major role in the industrial production of TiO2. To fully exploit the potential of these compounds as petrogenetic indicators and/or useful materials we need to quantitatively understand the factors controlling their properties and stability, and thus, to extrapolate beyond the calibrating experiments. For that purpose, we need to combine thermochemistry, phase equilibrium, and in situ P-V-T-cation disorder experimental data that presently either are incomplete or lacking. Perhaps, the most complete data set is that for MgTi2O5 (karrooite) which allows the calibration of models for the Gibbs free energy of the MgTi2O5 as a function of pressure, temperature, and the Mg2+-Ti4+ distribution between the two nonequivalent octahedral sites. Consequently, the effect of cation disorder on MgTi2O5 stability, and the phase relations among MgTi2O5, other titanium oxides, and silicate minerals can be examined. Calculated phase relations in the Mg-Ti-Si-O system and phase equilibrium experiments in Fe-bearing compositions suggest that pseudobrookite-type oxides may be a more common in rocks than previously realized. However, homogeneous and heterogeneous equilibria, and crystallization paths likely affect their stability. For example, isobaric increases in temperature favor disordering and thus entropy-stabilization, in contrast, isothermal increases in pressure have the opposite effect. Although, currently, the potential effect of composition to cation disorder cannot be fully explored, it appears that enrichment in trivalent cations probably enhances entropy-stabilization and thus may increase the stability of (Al, Fe, Cr, Ti)-rich pseudobrookites relative to that of (Mg, Fe)-rich ones. In addition, high-temperature, nearly isothermal, decompression paths of olivine+orthopyroxene+oxide assemblages may favor pseudobrookites (s.l.) over rutile and/or ilmenite, in contrast, cooling at low pressures seems to favor ilmenite and/or rutile. In the case of crustal and mantle xenoliths, the presence or absence of orthorhombic oxides is probably controlled by reactions with olivine, orthopyroxene, ilmenite, and rutile. In oceanic mantle xenoliths such reactions may also involve a TiO2-enriched but not SiO2-enriched melt/fluid, because pseudobrookites (s.l.) would react with the SiO2-enriched melt/fluid to form orthopyroxene and rutile. Parenthetically, experiments and model calculations in the Mg-Ti-Si-O system suggest that low degree partial melting of low-TiO2 bulk compositions may produce Ti-enriched liquids in equilibrium with olivine, orthopyroxen ad=nd MgTi2O5, rutile or ilmenite.

Xirouchakis, Dimistrios↗