Search NASA⌕ Search

SEARCH · Search NASA

Results for “Maximum Entropy Model”

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 OSTI2017

An in situ USAXS–SAXS–WAXS study of precipitate size distribution evolution in a model Ni-based alloy

Intermetallic γ' precipitates typically strengthen nickel-based superalloys. The shape, size and spatial distribution of strengthening precipitates critically influence alloy strength, while their temporal evolution characteristics determine the high-temperature alloy stability. Combined ultra-small-, small- and wide-angle X-ray scattering (USAXS–SAXS–WAXS) analysis can be used to evaluate the temporal evolution of an alloy's precipitate size distribution (PSD) and phase structure during in situ heat treatment. Analysis of PSDs from USAXS–SAXS data employs either least-squares fitting of a preordained PSD model or a maximum entropy (MaxEnt) approach, the latter avoiding a priori definition of a functional form of the PSD. However, strong low- q scattering from grain boundaries and/or structure factor effects inhibit MaxEnt analysis of typical alloys. Lastly, this work describes the extension of Bayesian–MaxEnt analysis methods to data exhibiting structure factor effects and low- q power law slopes and demonstrates their use in an in situ study of precipitate size evolution during heat treatment of a model Ni–Al–Si alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2015

CORRELATION ANALYSIS BETWEEN TIBET AS-γ TeV COSMIC RAY AND WMAP NINE-YEAR DATA

The WMAP team subtracted template-based foreground models to produce foreground-reduced maps, and masked point sources and uncertain sky regions directly; however, whether foreground residuals exist in the WMAP foreground-reduced maps is still an open question. Here, we use Pearson correlation coefficient analysis with AS-γ TeV cosmic ray (CR) data to probe possible foreground residuals in the WMAP nine-year data. The correlation results between the CR and foreground-contained maps (WMAP foreground-unreduced maps, WMAP template-based, and Maximum Entropy Method foreground models) suggest that: (1) CRs can trace foregrounds in the WMAP data; (2) at least some TeV CRs originate from the Milky Way; (3) foregrounds may be related to the existence of CR anisotropy (loss-cone and tail-in structures); (4) there exist differences among different types of foregrounds in the decl. range of <15°. Then, we generate 10,000 mock cosmic microwave background (CMB) sky maps to describe the cosmic variance, which is used to measure the effect of the fluctuations of all possible CMB maps to the correlations between CR and CMB maps. Finally, we do correlation analysis between the CR and WMAP foreground-reduced maps, and find that: (1) there are significant anticorrelations; and (2) the WMAP foreground-reduced maps are credible. However, the significant anticorrelations may be accidental, and the higher signal-to-noise ratio Planck SMICA map cannot reject the hypothesis of accidental correlations. We therefore can only conclude that the foreground residuals exist with ∼95% probability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Complex Structural Disorder in a Polar Orthorhombic Perovskite Observed through the Maximum Entropy Method/Rietveld Technique

Ambient pressure stable perovskite oxides with all Bi3+ on the A-site are rare, with only four examples known. Due to the lone pair on Bi 3+ , these materials are seen as the best alternative to Pb-based piezoelectrics, which are used widely in society. The industry standard piezoelectric, Pb (Zr 1 – x Ti x )O 3 , relies on the [001] polarization of PbTiO 3 , but there are currently no ambient pressure stable Bi-based perovskites with this polarization vector, preventing the creation of an analogous system. We present the full structural analysis of the orthorhombic phase of (1 – x)Bi (Ti 3/8 Fe 2/8 Mg 3/8 )O 3 – x CaTiO 3 , which crystallizes in Pna2 1 symmetry with [001] polarization. This symmetry is rare and has only been reported twice for perovskites at ambient conditions. Analysis of maximum entropy method (MEM) models using synchrotron radiation powder X-ray diffraction reveals a disordered A-site configuration, and the MEM/Rietveld technique generates a structural model of this extreme disorder. Combined Rietveld analysis of X-ray and neutron diffraction data yields an accurate description of the local A-site configuration, which we use to understand our dielectric, ferroelectric, and piezoelectric measurements. These results give insight into how to stabilize this unique symmetry and inspire new design principles for Bi-based piezoelectrics.

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

Statistical mechanics of collisionless orbits. IV. Distribution of angular momentum

It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L {sup 2}, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E, L {sup 2}). Instead, we require that when integrating N(E, L {sup 2}) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E, L {sup 2}) we show how the distribution of particles in L {sup 2} is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, β(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and β(r), must have linear or convex distributions in L {sup 2}, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step toward deriving the full distribution function for spherically symmetric dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

MENT-Flow: maximum-entropy phase space tomography using normalizing flows

Generative models can be trained to reproduce low-dimensional projections of high-dimensional phase space distributions. Normalizing flows are generative models that parameterize invertible transformations, allowing exact probability density evaluation and sampling. Consequently, flows are unbiased entropy estimators and could be used to solve the high-dimensional maximum-entropy tomography (MENT) problem. In this work, we evaluate a flow-based MENT solver (MENT-Flow) against exact maximum-entropy solutions and Minerbo's iterative MENT algorithm in two dimensions.

Hoover, Austin↗
DOE OSTI2014

Electron density distribution and disordered crystal structure of 8H-SiAlON, Si{sub 3−x}Al{sub 1+x}O{sub x}N{sub 5−x} (x∼2.2)

The 8H-SiAlON crystal with general formula Si{sub 3−x}Al{sub 1+x}O{sub x}N{sub 5−x} was characterized using laboratory X-ray powder diffraction (CuKα{sub 1}), transmission electron microscopy and energy dispersive X-ray spectroscopy. The [Si: Al] molar ratios were determined to be [0.21(1): 0.79(1)], corresponding to x=2.2(1). The Si{sub 0.8(1)}Al{sub 3.2(1)}O{sub 2.2(1)}N{sub 2.8(1)} compound is hexagonal with space group P6{sub 3}/mmc (Z=2). The unit-cell dimensions are a=0.298877(8) nm, c=2.30872(5) nm and V=0.178602(8) nm{sup 3}. The initial structural model was successfully derived by the charge-flipping method and further refined by the Rietveld method. The final structural model showed the positional disordering of two of the three (Si,Al) sites. The maximum-entropy method-based pattern fitting (MPF) method was used to confirm the validity of the split-atom model, in which conventional structure bias caused by assuming intensity partitioning was minimized. The reliability indices calculated from the MPF were R{sub wp}=4.74%, S (=R{sub wp}/R{sub e})=1.30, R{sub p}=3.40%, R{sub B}=1.04% and R{sub F}=0.74%. The disordered crystal structure was successfully described by overlapping four types of domains with ordered atom arrangements. The distribution of atomic positions in each of the domains can be achieved in the space group P6{sub 3}mc. Two of the four types of domains are related by a pseudo-symmetry inversion, and the two remaining domains also have each other the inversion pseudo-symmetry. - Graphical abstract: Bird's eye view of electron densities up to 49.3% (0.071 nm{sup −3}) of the maximum on the plane parallel to (110) with the corresponding atomic arrangements of Si{sub 0.8}Al{sub 3.2}O{sub 2.2}N{sub 2.8}. - Highlights: • Crystal structure of Si{sub 0.8}Al{sub 3.2}O{sub 2.2}N{sub 2.8} is determined by laboratory X-ray powder diffraction. • The atom arrangements are represented by the split-atom model. • The maximum-entropy method-based pattern fitting method is used to confirm the validity of the model. • The disordered structure is described by overlapping four types of domains with ordered atom arrangements.

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

Electron density distribution and disordered crystal structure of 15R-SiAlON, SiAl{sub 4}O{sub 2}N{sub 4}

The crystal structure of SiAl{sub 4}O{sub 2}N{sub 4} was characterized by laboratory X-ray powder diffraction (CuKα{sub 1}). The title compound is trigonal with space group R3-bar m. The hexagonal unit-cell dimensions (Z=3) are a=0.301332(3) nm, c=4.18616(4) nm and V=0.3291825(5) nm{sup 3}. The initial structural model was successfully derived by the charge-flipping method and further refined by the Rietveld method. The final structural model showed the positional disordering of one of the three (Si,Al) sites. The maximum-entropy method-based pattern fitting (MPF) method was used to confirm the validity of the split-atom model, in which conventional structure bias caused by assuming intensity partitioning was minimized. The reliability indices calculated from the MPF were R{sub wp}=5.05%, S (=R{sub wp}/R{sub e})=1.21, R{sub p}=3.77%, R{sub B}=1.29% and R{sub F}=1.01%. The disordered crystal structure was successfully described by overlapping three types of domains with ordered atom arrangements. The distribution of atomic positions in one of the three types of domains can be achieved in the space group R3-bar m. The atom arrangements in the other two types of domains are noncentrosymmetrical with the space group R3m. These two structural configurations are related by the pseudo-symmetry inversion. -- Graphical abstract: A bird's eye view of electron densities up to 75.3% (0.133 nm{sup −3}) of the maximum on the plane parallel to (110) with the corresponding atomic arrangements of SiAl{sub 4}O{sub 2}N{sub 4}. Highlights: • Crystal structure of SiAl{sub 4}O{sub 2}N{sub 4} is determined by laboratory X-ray powder diffraction. • The atom arrangements are represented by the split-atom model. • The maximum-entropy method-based pattern fitting method is used to confirm the validity of the model. • The disordered structure is described by overlapping three types of domains with ordered atom arrangements.

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

ENTROPY PRODUCTION IN COLLISIONLESS SYSTEMS. III. RESULTS FROM SIMULATIONS

The equilibria formed by the self-gravitating, collisionless collapse of simple initial conditions have been investigated for decades. We present the results of our attempts to describe the equilibria formed in N-body simulations using thermodynamically motivated models. Previous work has suggested that it is possible to define distribution functions for such systems that describe maximum entropy states. These distribution functions are used to create radial density and velocity distributions for comparison to those from simulations. A wide variety of N-body code conditions are used to reduce the chance that results are biased by numerical issues. We find that a subset of initial conditions studied lead to equilibria that can be accurately described by these models, and that direct calculation of the entropy shows maximum values being achieved.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

General multi-group macroscopic modeling for thermo-chemical non-equilibrium gas mixtures

This paper opens a new door to macroscopic modeling for thermal and chemical non-equilibrium. In a game-changing approach, we discard conventional theories and practices stemming from the separation of internal energy modes and the Landau-Teller relaxation equation. Instead, we solve the fundamental microscopic equations in their moment forms but seek only optimum representations for the microscopic state distribution function that provides converged and time accurate solutions for certain macroscopic quantities at all times. The modeling makes no ad hoc assumptions or simplifications at the microscopic level and includes all possible collisional and radiative processes; it therefore retains all non-equilibrium fluid physics. We formulate the thermal and chemical non-equilibrium macroscopic equations and rate coefficients in a coupled and unified fashion for gases undergoing completely general transitions. All collisional partners can have internal structures and can change their internal energy states after transitions. The model is based on the reconstruction of the state distribution function. The internal energy space is subdivided into multiple groups in order to better describe non-equilibrium state distributions. The logarithm of the distribution function in each group is expressed as a power series in internal energy based on the maximum entropy principle. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients succinctly to any order. The model’s accuracy depends only on the assumed expression of the state distribution function and the number of groups used and can be self-checked for accuracy and convergence. We show that the macroscopic internal energy transfer, similar to mass and momentum transfers, occurs through nonlinear collisional processes and is not a simple relaxation process described by, e.g., the Landau-Teller equation. Unlike the classical vibrational energy relaxation model, which can only be applied to molecules, the new model is applicable to atoms, molecules, ions, and their mixtures. Numerical examples and model validations are carried out with two gas mixtures using the maximum entropy linear model: one mixture consists of nitrogen molecules undergoing internal excitation and dissociation and the other consists of nitrogen atoms undergoing internal excitation and ionization. Results show that the original hundreds to thousands of microscopic equations can be reduced to two macroscopic equations with almost perfect agreement for the total number density and total internal energy using only one or two groups. We also obtain good prediction of the microscopic state populations using 5-10 groups in the macroscopic equations.

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

Modelling bivariate extreme precipitation distribution for data-scarce regions using Gumbel-Hougaard copula with maximum entropy estimation

A new method of parameter estimation in data scarce regions is valuable for bivariate hydrological extreme frequency analysis. Here, this paper proposes a new method of parameter estimation (maximum entropy estimation, MEE) for both Gumbel and Gumbel–Hougaard copula in situations when insufficient data are available. MEE requires only the lower and upper bounds of two hydrological variables. To test our new method, two experiments to model the joint distribution of the maximum daily precipitation at two pairs of stations on the tributaries of Heihe and Jinghe River, respectively, were performed and compared with the method of moments, correlation index estimation, and maximum likelihood estimation, which require a large amount of data. Both experiments show that for the Ye Niugou and Qilian stations, the performance of MEE is nearly identical to those of the conventional methods. For the Xifeng and Huanxian stations, MEE can capture information indicating that the maximum daily precipitation at the Xifeng and Huanxian stations has an upper tail dependence, whereas the results generated by correlation index estimation and maximum likelihood estimation are unreasonable. Moreover, MEE is proved tobe generally reliable and robust by many simulations under three different situations. TheGumbel–Hougaard copula with MEE can also be applied to the bivariate frequency analysis ofother extreme events in data-scarce regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

MULTI-PRED: A Software Module for Predictive Modeling of Coupled Multi-Physics Systems: User's Manual

This User’s Manual describes the code module MULTI-PRED, written in FORTRAN which implements the methodology for “predictive modeling of coupled multi-physics systems (PM-CMPS)” formulated by Cacuci (2014). This methodology fully takes into account the coupling terms between the systems but requires only the computational resources that would be needed to perform predictive modeling on each system separately. The PM-CMPS methodology uses the maximum entropy principle to construct an optimal approximation of the unknown a priori distribution based on a priori known mean values and uncertainties characterizing the experimental and computational parameters and results of interest responses, called for the multi-physics models under consideration. This “maximum entropy” a priori distribution is combined, using Bayes’ theorem, with the “likelihood” provided by the multi-physics simulation models to obtain a formal posterior distribution. Subsequently, the posterior distribution thus obtained is evaluated using the saddle-point method to obtain analytical expressions for the optimally predicted values for the multi-physics models parameters and responses along with corresponding reduced uncertainties. Noteworthy, the predictive modeling methodology for the coupled systems is constructed such that the systems can be considered sequentially rather than simultaneously, while preserving exactly the same results as if the systems were treated simultaneously. Consequently, very large coupled systems, which could perhaps exceed available computational resources if treated simultaneously, can be treated with the PM-CMPS methodology presented in this work sequentially and without any loss of generality or information, requiring just the resources that would be needed if the systems were treated sequentially. Three illustrative demonstration problems are also provided. The first problem presents the application of the PM-CMPS methodology to a simple particle diffusion problem which admits a closed-form analytical solution which facilitates a rapid understanding of this methodology and its predicted results. The second demonstration problem presents the application of the PM-CMPS methodology to the problem of inverse prediction, from detector responses in the presence of counting uncertainties, of the thickness of a homogeneous slab of material containing uniformly distributed gamma-emitting sources, for optically thin and thick slabs. This problem highlights the essential role played by the relative uncertainties (or, conversely, accuracies) of measured and computed responses. The third demonstration problem presents the application of the PM-CMPS methodology to the F-area cooling towers at the Savannah River National Lab. This problem demonstrates that the PM-CMPS methodology reduces the predicted response uncertainties not only at locations where measurements are available, but also at locations where measurements are not available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Upscaling ab-initio chemistry models to non-equilibrium flow simulations

Thermochemical nonequilibrium is generated by hypersonic cruise vehicles, high speed combustion and detonations. Nonequilibrium interactions affect reaction rates and energy transfer that are critical for many engineering applications. Ab-initio potential energy surfaces can be used in quasi-classical trajectory (QCT) calculations to calculate state-specific dissociation rates and cross sections, as well as state-to-state energy exchange rates and cross sections. These rates and cross sections are needed for simulating such nonequilibrium flows. However, there are two main challenges in extending QCT calculations to modeling all chemical processes found in typical engineering flowfields. First, ab-initio surfaces for many collision partners are not available and when available their fidelity is not quantified. Second, state-to-state cross section datasets are typically on the order of 1000 4, , and are too large to be efficiently implemented in engineering simulations. This work focuses on O 2 +O interactions because of prevalence of oxygen in many engineering application; however, all of the models developed here may be applied to other diatom-atom interactions and can be extended to other diatom-diatom combinations. It is shown in this work that at temperatures between 2,500 and 20,000K, state-specific cross sections are not highly sensitive to potentials and a Morse additive pairwise (MAP) potential is proposed for modeling collisions for which ab-initio potentials are not available. Next, this work suggests an approach to model state-to-state cross sections based on maximum entropy considerations. The calculated cross sections are fit to a new model with just 8 parameters, which can be easily scaled up to flow simulations.

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

Model-Free Approach for Profiling of Polydisperse Soft Matter Using Small Angle Scattering

A strategy for determining the size polydispersity of systems from their small angle coherent scattering is outlined. Here, using the method of moment expansion, we show that the various central moments representing the average particle size, variance of particle size, and skewness of size distribution function (SDF) for polydisperse systems can be extracted from spectral analysis without bias. When the degree of polydispersity is moderate, SDF can be further reconstructed based on the maximum entropy principle. Numerical benchmarking of a model study over a wide range of size nonuniformity demonstrates the validity of this analytical approach for quantifying the size distribution of general soft matter systems in a model-free manner. Furthermore, the efficacy of this method was validated by successfully applying it to the fitting of small-angle neutron scattering data obtained from L64 Pluronic micelles using various form factor models. The numerical and experimental verification underscores the reliability and versatility of this method in accurately characterizing the size distribution of complex soft matter systems.

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 OSTI2022

An ecological niche model to predict the geographic distribution of Haemagogus janthinomys, Dyar, 1921 a yellow fever and Mayaro virus vector, in South America

Yellow fever virus (YFV) has a long history of impacting human health in South America. Mayaro virus (MAYV) is an emerging arbovirus of public health concern in the Neotropics and its full impact is yet unknown. Both YFV and MAYV are primarily maintained via a sylvatic transmission cycle but can be opportunistically transmitted to humans by the bites of infected forest dwelling Haemagogus janthinomys Dyar, 1921. To better understand the potential risk of YFV and MAYV transmission to humans, a more detailed understanding of this vector species’ distribution is critical. This study compiled a comprehensive database of 177 unique Hg . janthinomys collection sites retrieved from the published literature, digitized museum specimens and publicly accessible mosquito surveillance data. Covariate analysis was performed to optimize a selection of environmental (topographic and bioclimatic) variables associated with predicting habitat suitability, and species distributions modelled across South America using a maximum entropy (MaxEnt) approach. Our results indicate that suitable habitat for Hg . janthinomys can be found across forested regions of South America including the Atlantic forests and interior Amazon.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Fast and Efficient Stochastic Optimization for Analytic Continuation

In this analytic continuation of imaginary-time quantum Monte Carlo data to extract real-frequency spectra remains a key problem in connecting theory with experiment. Here we present a fast and efficient stochastic optimization method (FESOM) as a more accessible variant of the stochastic optimization method introduced by Mishchenko et al. [Phys. Rev. B 62, 6317 (2000)], and we benchmark the resulting spectra with those obtained by the standard maximum entropy method for three representative test cases, including data taken from studies of the two-dimensional Hubbard model. Genearally, we find that our FESOM approach yields spectra similar to the maximum entropy results. In particular, while the maximum entropy method yields superior results when the quality of the data is strong, we find that FESOM is able to resolve fine structure with more detail when the quality of the data is poor. In addition, because of its stochastic nature, the method provides detailed information on the frequency-dependent uncertainty of the resulting spectra, while the maximum entropy method does so only for the spectral weight integrated over a finite frequency region. Therefore, we believe that this variant of the stochastic optimization approach provides a viable alternative to the routinely used maximum entropy method, especially for data of poor quality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗