Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum entropy method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Low-energy enhancement of the magnetic dipole radiation in odd-mass lanthanides

We compute the magnetic dipole (M1) $\gamma$-ray strength functions ($\gamma$SF) for the odd-mass lanthanides $^{\textrm{143-151}}$Nd and $^{\textrm{147-153}}$Sm using the shell-model Monte Carlo method in combination with the static-path approximation and the maximum-entropy method. In particular, we quantify the statistical uncertainties in the calculated M1 $\gamma$SFs and show that they are under control for the excitation energies relevant to the experiments despite a Monte Carlo sign problem that originates in the projection onto an odd number of neutrons. We identify a low-energy enhancement (LEE) in the M1 $\gamma$SFs of these odd-mass lanthanides, which was recently observed experimentally in some of them. We also find a scissors mode resonance (SR) in the strongly deformed isotopes. We observe that the decrease in the LEE strength with neutron number along an isotopic chain is compensated for by an increase in the SR strength in the deformed nuclei. Furthermore, we compare our results with recent experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Floating zone crystal growth, structure, and properties of a cubic Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 garnet-type lithium-ion conductor

As a promising candidate for solid-state electrolytes in Li-ion batteries, the garnet-type Li-ion conductor series Li 5+x La 3 Nb 2-x Zr x O 12 (LLNZO) (0 ≤ x ≤ 2) exhibits high ionic conductivity at room temperature. However, no previous single-crystal growth or characterization has been reported for LLNZO compositions 0 ≤ x ≤ 1. To obtain a complete understanding of the trend in the structure–property relationship in this class of materials, we used the floating zone (FZ) method to grow a single crystal of Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 that was 4 mm in diameter and 10 mm in length. Using Laue neutron single-crystal diffraction, two distinct Li sites were observed: tetrahedral 24d and octahedral 96h sites. The maximum entropy method (MEM) based on neutron single-crystal diffraction data was used to map Li nuclear density and estimate that the bottleneck of Li transport exists between neighboring tetrahedral and octahedral sites, and that Li is delocalized between split octahedral sites. Room-temperature Li-ion conductivity in Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 measured with electrochemical impedance spectroscopy (EIS) was 1.37 × 10 -4 S cm -1 . The Li migration activation energy was estimated to be 0.50 eV from EIS and 0.47 eV from dielectric relaxation measurements. The Li-ion jump attempt rate was estimated to be 1.47 × 10 12 Hz while the time scale of successful migration is 10 -7 to 10 -6 s.

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

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Walking on the dark side: Anthropogenic factors limit suitable habitat for gray wolf ( Canis lupus ) in a large natural area covering Belarus and Ukraine

Due to successful conservation initiatives and legislations, the grey wolf (Canis lupus) is recolonising its historic range in Europe. However, wolves have never been extirpated across large areas in Eastern Europe but are often constrained to remote and inaccessible places due to centuries of persecution. This study aimed to identify the potentially suitable wolf habitats in Polesia, a massive cross-border lowland region extending over southern Belarus and northern Ukraine, which are often neglected in large carnivore studies at the continental scale. We hypothesized that anthropogenic rather than environmental factors govern wolf habitat suitability. We used a dataset of 4191 GPS locations obtained from radio-collared wolves (n = 26) and confirmed observations (n = 231) during 2014–2021 and applied maximum entropy method to estimate relative habitat suitability for wolves in Polesia. Artificial light at night (ALAN), proportion of cropland and tree cover were the most important factors affecting wolf habitat suitability. Road densities contributed poorly to predicting habitat suitability for wolves. Our models predicted a quarter of Polesia as suitable habitat and revealed priority areas connecting the important source populations in the Chornobyl Exclusion Zone in the east and the Bialowieza Forest in the west and thus essential for long-term wolf conservation. Our results provide the bases for effective, long-term wolf monitoring and management programs in both Belarus and Ukraine. However, national and transboundary wolf management in Polesia has been extremely challenging since 2022 due to the ongoing war and subsequent habitat degradation in this part of Europe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

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

Nanodomain Formation and Temperature-Dependent Diffusion in Deep Eutectic Solvents Revealed by Single-Molecule Tracking

Deep eutectic solvents (DESs) are typically regarded as homogeneous liquids; however, recent work shows that many exhibit nanoscale structural heterogeneity. Most studies attribute these nanoscale features to short-range chemical interactions. It is still unclear whether a long-range physical mechanism also plays a role. Here, in this study, we examined the nanoscale structure in two hydrophobic DESs, 1:3 tetrabutylammonium bromide: l-menthol (DES-butyl) and 1:3 tetraoctylammonium bromide: l-menthol (DES-octyl). The notation 1:3 represents the molar ratio of the hydrogen bond acceptors to hydrogen bond donors used in the synthesis of the DESs. Single-molecule tracking (SMT) coupled with maximum entropy method (MEM) analysis was used to measure the number of diffusion populations of a dilute concentration of an added fluorescent probe. The presence of more than one population of diffusion coefficients indicates the existence of multiple local environments for the fluorescent probe (i.e., nanoscale structures in the DES). DES-butyl showed a relatively narrow diffusion coefficient distribution centered at 0.55 μm 2 /s, whereas DES-octyl displayed two distinct diffusing populations at 20 °C, with diffusion coefficients of 0.12 μm 2 /s and 0.53 μm 2 /s for the slow and fast populations, respectively. As DES-octyl was heated, the slow-diffusing population steadily diminished and disappeared above ∼30 °C, indicating that the nanodomains present at lower temperatures collapse as the liquid becomes more thermodynamically mixed. This temperature-dependent homogenization is consistent with a physical mechanism of nanostructure formation, for example, liquid–liquid phase separation (LLPS), wherein the structure is not driven solely by specific chemical interactions. The SMT-MEM results suggest that a long-range physical mechanism is the most plausible origin of the measured nanoscale structure in DES-octyl.

Opare-Addo, Jemima [Ames Laboratory (AMES), Ames, ↗
DOE OSTI2025

Investigation of Americium-Containing Phosphates, Silicates, Borates, Molybdates, and Fluorides Synthesized via High-Temperature Flux Crystal Growth

The crystal chemistry of americium-containing extended structures was investigated, and several classes of americium-containing solid-state oxide materials were obtained in single-crystal form via high-temperature flux crystal growth. This enabled the structural characterization of rare examples of ternary, quaternary, and penternary americium-containing silicates K 3 Am (Si 2 O 7 ) and Cs 6 Am 2 Si 21 O 48 , phosphates Na 3 Am (PO 4 ) 2 and K 3 Am (PO 4 ) 2 , borates Ba 3 Am 2 (BO 3 ) 4 and AmBO 3 , borate halides Ca 5 Am(BO 3 ) 4 Cl, molybdates Li 0.5 Am 0.5 MoO 4 , and fluorides CsAm 2 F 7 . Using these crystallographic data, the ionic radii of Am 3+ with coordination numbers of six (0.975 Å), seven (1.052 Å), and nine (1.162 Å) were established. A maximum entropy method (MEM) analysis was performed on the single-crystal X-ray diffraction data that were collected for K 3 Nd(PO 4 ) 2 /K 3 Am(PO 4 ) 2 , K 3 NdSi 2 O 7 /K 3 AmSi 2 O 7 , and NdBO 3 /AmBO 3 , to qualitatively compare the ionicities of the Nd–O and Am–O bonds. In conclusion, Raman spectroscopy data were collected on single crystals of K 3 Am(PO 4 ) 2 and compared to the calculated Raman spectrum of K 3 Am(PO 4 ) 2 obtained from DFT calculations.

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

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

PNNL-CompBio/BoltzmannMFX

BoltzmannMFX is a biological simulation code that solves chemical reaction networks using maximum entropy methods. It uses modules from MFiX-Exa and is based on the AMReX framework for massively parallel block-structured adaptive mesh applications.

Palmer, Bruce↗
DOE OSTI2024

High Magnetic Anisotropy and Magnetocaloric Effects in Single-Crystal Cr 2 Te 3

Here, we report a systematic investigation of anisotropic magnetocaloric effects in single-crystal Cr 2 Te 3 . Single-crystal samples are synthesized by chemical vapor transport and characterized by X-ray and Laue diffraction methods. The maximum magnetic entropy change –ΔS M max is 4.50 J kg –1 K –1 for the easy c-axis (3.36 J kg –1 K –1 for the hard axis along ab-plane), and the relative cooling power (RCP) is 296.7 J kg –1 for the easy c-axis (183.84 J kg –1 for the hard axis along ab-plane) for a magnetic field change of 9 T near the Curie temperature. The magneto-crystalline anisotropy constant K u is estimated to be 580.12 kJ m –3 at 140 K, decreasing to 148.60 kJ m –3 at 168 K. Meanwhile, the maximum of the rotational magnetic entropy change –ΔS M R (T, H) between the c-axis and the ab-plane is about 1.14 J kg –1 K –1 for magnetic-field change of 9 T. The critical exponents are estimated by analyzing magnetocaloric effects, which indicate a 2D-Ising type magnetic system. The accuracy of estimated critical exponents is verified by scaling analysis. The maximum magnetic entropy change –ΔS M max ≈ 5.25 J kg –1 K –1 (along the c-axis) and the corresponding adiabatic temperature change ΔT ad ≈ 3.31 K (along the c-axis) are estimated by analyzing heat capacity measurements with a magnetic field up to 9 T.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Animal movement estimation and network-based epidemic modeling: Illustration for the swine industry in Iowa (US)

Animal movement plays a critical role in disease transmission between farms. However, in the United States, the lack of available animal shipment data, sometimes coupled with a lack of detailed information about farm demographics and characteristics, presents great challenges for epidemic modeling and prediction. In this study, we proposed a new method based on the maximum entropy to generate “synthetic” animal movement networks, considering available statistics about the premises operation type, operation size, and the distance between premises. We illustrated our method for the swine movement networks in Iowa and performed network analyses to gain insights into the swine industry. We then applied the generated networks to a network-based epidemic model to identify potential system vulnerabilities in terms of disease transmission. The model was parameterized for African Swine Fever (ASF) as the US swine industry is quite concerned about this disease. Results show that premises with a central role in the network are more vulnerable to disease outbreaks and play an important role in disease spread. Simulations with outbreaks starting from random farms reveal no significant large outbreaks, indicating the system’s relative robustness against arbitrary disease introductions. However, outbreaks originating from high out-degree farms can lead to large epidemic sizes. This underscores the importance for stakeholders and policymakers to continue improving animal movement records and traceability programs in the US and the value of making that data available to epidemiologists and modelers to better understand risk and inform strategies aimed to cost-effectively prevent and control disease transmission. Our approach could be easily adapted to estimate movement networks in other animal production systems and to inform disease spread models for various infectious diseases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Off-policy deep reinforcement learning with automatic entropy adjustment for adaptive online grid emergency control

Electric overloading conditions and contingencies put modern power systems at risk of voltage collapse and blackouts. Load shedding is crucial to maintain voltage stability for grid emergency control. However, the rule- or model-based schemes rely on accurate dynamic system models and face considerable challenges in adapting to various operating conditions and uncertain event occurrences. Here, to address these issues, this paper proposes a novel deep reinforcement learning (DRL)-based voltage stability control algorithm with automatic entropy adjustment (AEA) for grid emergency control. Various dynamic network components for complex system operations are modeled to construct the DRL environment. An off-policy soft actor-critic architecture is developed to maximize the expected reward and policy entropy simultaneously. The AEA mechanism is proposed to facilitate the policy maximum entropy procedure, and the proposed method can automatically provide effective discrete and continuous actions against various fault scenarios. Our approach accomplishes high sampling efficiency, scalability, and auto-adaptivity of the control policies under high uncertainties. Comparative studies with the existing DRL-based control methods in IEEE benchmarks indicate salient performance improvement of the proposed method for dynamic system emergency control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

SmoQyDEAC.jl: A differential evolution package for the analytic continuation of imaginary time correlation functions

We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Codebase release r1.1 for SmoQyDEAC.jl

We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package’s output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.

Neuhaus, James (ORCID:0000000169048510)↗
DOE OSTI2024

Uncertainty-quantification-enabled inversion of nuclear responses

Nuclear quantum many-body methods rely on integral transform techniques to infer properties of electroweak response functions from ground-state expectation values. Retrieving the energy dependence of these responses is highly nontrivial, especially for quantum Monte Carlo methods, as it requires inverting the Laplace transform, a notoriously ill-posed problem. Here, in this work, we propose an artificial neural network architecture suitable for accurate response function reconstruction with precise estimation of the uncertainty of the inversion. We demonstrate the capabilities of this new architecture benchmarking it against maximum entropy and previously developed neural network methods designed for a similar task, paying particular attention to its robustness noise in the Euclidean

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Modified Maximum Entropy Inverse Reinforcement Learning Approach for Microgrid Energy Scheduling

Increasing popularity of integrating distributed energy resources (DERs) into the power system brings a challenge to optimize the microgrid dispatch policy. The reinforcement learning methods suffer from a long-time problem with the theoretical assumption of the objective/reward function for the microgrid system. Although the traditional inverse reinforcement learning (IRL) approaches can solve this problem to some extent, they encounter a limitation of complex computations for state visitation frequency in the large and continuous state space. To alleviate this limitation, we propose a modified maximum entropy IRL (MMIRL) method to extract the reward function from the expert demonstrations for solving the microgrid energy scheduling problem. The proposed MMIRL algorithm is promising in recovering the reward function and learning the dispatch policy compared to conventional approaches. Case studies are performed in an energy arbitrage problem and a microgrid system with DERs. Results substantiate that the proposed MMIRL approach can learn the dispatch policy with more than 99% efficiency and outperforms other comparative methods.

artificial intelligence, reinforcement learning, m↗
DOE OSTI2022

Quasar Identification Using Multivariate Probability Density Estimated from Nonparametric Conditional Probabilities

Nonparametric estimation for a probability density function that describes multivariate data has typically been addressed by kernel density estimation (KDE). A novel density estimator recently developed by Farmer and Jacobs offers an alternative high-throughput automated approach to univariate nonparametric density estimation based on maximum entropy and order statistics, improving accuracy over univariate KDE. This article presents an extension of the single variable case to multiple variables. The univariate estimator is used to recursively calculate a product array of one-dimensional conditional probabilities. In combination with interpolation methods, a complete joint probability density estimate is generated for multiple variables. Good accuracy and speed performance in synthetic data are demonstrated by a numerical study using known distributions over a range of sample sizes from 100 to 10 6 for two to six variables. Performance in terms of speed and accuracy is compared to KDE. The multivariate density estimate developed here tends to perform better as the number of samples and/or variables increases. As an example application, measurements are analyzed over five filters of photometric data from the Sloan Digital Sky Survey Data Release 17. The multivariate estimation is used to form the basis for a binary classifier that distinguishes quasars from galaxies and stars with up to 94% accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗