Search NASASearch

SEARCH · Search NASA

Results for “Interpretability”

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

Interpretations of the ATLAS measurements of Higgs boson production and decay rates and differential cross-sections in pp collisions at $\sqrt{s}$ = 13 TeV

Measurements of the Higgs boson production times decay rates and differential cross-sections have recently been performed by the ATLAS experiment in several decay channels using up to 139 fb -1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded at the Large Hadron Collider. This paper presents multiple interpretations of these Higgs boson measurements. Measurements of production-mode cross-sections, simplified template cross-sections and fiducial differential cross-sections in different decay channels are reparameterised in terms of the impact of Standard Model effective field theory operators, and constraints are reported on the corresponding Wilson coefficients. Production and decay rate measurements are interpreted in UV-complete extensions of the Standard Model, namely the two-Higgs-doublet model (2HDM) near the alignment limit and the Minimal Supersymmetric Standard Model (MSSM) for various MSSM benchmark scenarios. The constraints on the 2HDM parameters (cos(β - α), tanβ) and the MSSM parameters (m A , tanβ) are complementary to those obtained from direct searches for additional Higgs bosons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Interpretation of high-harmonic fast-wave propagation in the scrape-off layer of NSTX-U as a geometrically bounded-waveguide mode

This paper revisits the onset of high-harmonic fast-wave (HHFW) propagation in the scrape-off layer (SOL) plasma of the NSTX/NSTX-U spherical tokamak, motivated by past HHFW heating and current drive experiments and modelling. Previously, the fast-wave propagation in the SOL was correlated with the opening (suppression) of the fast-wave right-hand cutoff layer in front of the antenna. In this work, the SOL propagation is interpreted as a geometric waveguide mode guided by the SOL geometry, whose radial width can be comparable to the wave’s perpendicular wavelength in the HHFW regime. A two-dimensional circular model is first employed using full-wave solvers to characterise the poloidal eigenmode structures and to clarify their relationship to annulus resonance. By progressively adding a tokamak-like magnetic-field configuration, starting from a uniform axial field, the effects of magnetic-field gradients and pitch on the SOL waveguide mode are characterised. The poloidal mode numbers supported by the SOL plasma agree well with analytic estimates, indicating that an anisotropic plasma in a bounded geometry selectively supports and amplifies the resonant poloidal mode number. Additionally, two-dimensional axisymmetric NSTX-U simulations demonstrate that the SOL eigenmode features identified in the circular model persist in the experimentally relevant configurations. A control approach based on lengthening the wave perpendicular wavelength relative to the SOL width, including a higher magnetic (B)-field operation, is discussed. The analysis here shows the key role of a bounded geometry in interpreting HHFW eigenmode coupling and propagation in the NSTX-U SOL plasma.

plasma heating

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE

Learning Physically Interpretable Atmospheric Models From Data With WSINDy

The multiscale and turbulent nature of Earth's atmosphere has historically rendered accurate weather modeling a hard problem. Recently, there has been an explosion of interest surrounding data-driven approaches to weather modeling, which in many cases show improved forecasting accuracy and computational efficiency when compared to traditional methods. However, many of the current data-driven approaches employ highly parameterized neural networks, often resulting in uninterpretable models and limited gains in scientific understanding. In this work, we address the interpretability problem by explicitly discovering partial differential equations governing atmospheric phenomena, identifying symbolic mathematical models with direct physical interpretations. The purpose of this paper is to demonstrate that, in particular, the weak-form sparse identification of nonlinear dynamics (WSINDy) algorithm can learn effective atmospheric models from both simulated and assimilated data. Our approach adapts the standard WSINDy algorithm to work with high-dimensional fluid data of arbitrary spatial dimension.

58 GEOSCIENCES

Evaluating the factors influencing accuracy, interpretability, and reproducibility in the use of machine learning classifiers in biology to enable standardization

The complexity and variability of biological data has promoted the increased use of machine learning methods to understand processes and predict outcomes. These same features complicate reliable, reproducible, interpretable, and responsible use of such methods, resulting in questionable relevance of the derived. outcomes. Here we systematically explore challenges associated with applying machine learning to predict and understand biological processes using a well- characterized in vitro experimental system. We evaluated factors that vary while applying machine learning classifers: (1) type of biochemical signature (transcripts vs. proteins), (2) data curation methods (pre- and post-processing), and (3) choice of machine learning classifier. Using accuracy, generalizability, interpretability, and reproducibility as metrics, we found that the above factors significantly mod- ulate outcomes even within a simple model system. Our results caution against the unregulated use of machine learning methods in the biological sciences, and strongly advocate the need for data standards and validation tool-kits for such studies.

59 BASIC BIOLOGICAL SCIENCES

Informed total-error-minimizing priors: Interpretable cosmological parameter constraints despite complex nuisance effects

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, for example, when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of informed total-error-minimizing (ITEM) priors to address this problem. An ITEM prior is a prior distribution on a set of nuisance parameters, such as those describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters (e.g., cosmological parameters). Our method works as follows. For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut, we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we applied our method to the density split statistics measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and that they allow for sharpened yet robust constraints on the parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS

Interpreting scattered neutron spectra on omega cryogenic implosions with measured backgrounds and higher-order scattering

Measurements of the areal density (⁠ρR⁠) of inertially confined implosions are critical to evaluate their performance. On OMEGA, ρR is inferred from measurements of the scattered neutron spectrum via neutron-time-of-flight (nTOF) and magnetic recoil spectrometer detectors. The nTOF measurements, in particular, have the ability to measure a wide range of neutron energies and thus scattering angles, with high precision. However, the neutrons that are backscattered into the detector must transit through the entire dense fuel assembly, and as a result are subject to rescattering effects, which have been heretofore neglected in the interpretation of OMEGA neutron spectra. At the backscatter edge, neglecting rescattering on OMEGA can lead to a ~10 to 20% reduction of the apparent areal density and therefore must be included in the analysis. The low ρR s on OMEGA also mean that backgrounds from non-target physics effects, such as scattering in the target chamber or slow scintillator decays, can significantly alter the measured signal. Here, we discuss how the Monte Carlo neutron spectrum postprocessor IRIS is used to include second-order scattering effects in the interpretation of the scattered neutron spectrum in OMEGA implosion experiments, and how dedicated implosion experiments are used to measure nTOF backgrounds in order to infer backscattered areal densities on OMEGA.

Deuterium

Interpreting AI for fusion: An application to plasma profile analysis for tearing mode stability

Artificial intelligence models have demonstrated strong predictive capabilities for various instabilities in fusion devices such as Tokamaks, including tearing modes (TM), edge localized modes, and disruptive events, but their opaque nature raises concerns about safety and trustworthiness when applied to fusion power plants. Here, we present a physics-based interpretation framework using a TM prediction model as a demonstration that is validated through a dedicated DIII-D TM avoidance experiment. By applying Shapley analysis, we identify how profiles such as rotation, temperature, and density contribute to the model's prediction of TM stability. Our analysis shows that in our experimental scenario, core electron temperature and rotation peaking play the primary role in TM stability, while density changes have smaller effects on stability. We show that off-axis ion temperature stabilizes TMs, suggesting that off-axis neutral beam heating can further stabilize this scenario. This work presents a generalizable ML-based event prediction methodology, from training to physics-driven interpretation, bridging the gap between physics understanding and opaque ML models.

Farre-Kaga, Hiro J. [Princeton Univ., NJ (United S

Combined measurements and interpretations of Higgs boson production and decay in proton–proton collisions at $\sqrt{s}$ = 13 TeV

Combined measurements of Higgs boson production and decay rates are reported, representing the most comprehensive study performed by the CMS Collaboration to date. The included analyses use proton–proton collision data recorded by the CMS experiment at $\sqrt{s}$ = 13 TeV from 2016 to 2018, corresponding to an integrated luminosity of 138 fb –1 . The statistical combination is based on analyses that measure the following decay channels: H → γγ, H → ZZ, H → WW, H → ππ, H → bb, H → μμ, and H → Zγ. Information in the events from each decay channel is used to target multiple Higgs boson production processes. Searches for invisible Higgs boson decays are also considered, as well as an analysis that measures off-shell Higgs boson production in the H → ZZ → 4$\ell$ ($\ell$ = e, μ) decay channel. The best fit inclusive signal yield is measured to be 1.014$^{+0.055}_{–0.053}$ times the standard model (SM) expectation, for a Higgs boson mass of 125.38 GeV. Measurements in kinematic regions defined by the simplified template cross section framework are also provided, as well as interpretations in the coupling modifier and SM effective field theory frameworks. The coupling modifier interpretation is further used to place constraints on various two-Higgs-doublet models. The results show good compatibility with the SM predictions for the majority of the measured parameters.

CMS

Massive νs through the CNN lens: interpreting the field-level neutrino mass information in weak lensing

Modern cosmological surveys probe the Universe deep into the nonlinear regime, where massive neutrinos suppress cosmic structure. Traditional cosmological analyses, which use the 2-point correlation function to extract information, are no longer optimal in the nonlinear regime, and there is thus much interest in extracting beyond-2-point information to improve constraints on neutrino mass. Quantifying and interpreting the beyond-2-point information is thus a pressing task. We study the field-level information in weak lensing convergence maps using convolution neural networks. We find that the network performance increases as higher source redshifts and smaller scales are considered — investigating up to a source redshift of 2.5 and ℓ max ≃ 10 4 — verifying that massive neutrinos leave a distinct effect on weak lensing. However, the performance of the network significantly drops after scaling out the 2-point information from the maps, implying that most of the field-level information can be found in the 2-point correlation function alone. We quantify these findings in terms of the likelihood ratio and also use Integrated Gradient saliency maps to interpret which parts of the map the network is learning the most from. We find that, in the absence of noise, the network extracts a similar amount of information from the most overdense and underdense regions. However, upon adding noise, the information in underdense regions is distorted as noise disproportionately washes out void-like structures.

Golshan, Malika [University of California, Berkele

R ( 3780 ) Resonance Interpreted as the 1 D 1 3 -Wave Dominant State of Charmonium from Precise Measurements of the Cross Section of e + e − → Hadrons

We report the precise measurements of the cross section of e + e − → hadrons at center-of-mass energies from 3.645 to 3.871 GeV. We thereby perform the most precise study of the cross sections and find a complex system composed of three resonances of R ( 3760 ) , R ( 3780 ) , and R ( 3810 ) . For the first time, we measure the R ( 3810 ) electronic width to be ( 19.4 ± 7.4 ± 12.1 ) eV . For the R ( 3760 ) resonance, we measure the mass to be ( 3751.9 ± 3.8 ± 2.8 ) MeV / c 2 , the total width to be ( 32.8 ± 5.8 ± 8.7 ) MeV , and the electronic width to be ( 184 ± 75 ± 86 ) eV . For the R ( 3780 ) resonance, we measure its mass to be ( 3778.7 ± 0.5 ± 0.3 ) MeV / c 2 , total width to be ( 20.3 ± 0.8 ± 1.7 ) MeV , and electronic width to be ( 265 ± 67 ± 83 ) eV . Forty-seven years ago, the ψ ( 3770 ) resonance was discovered, and was subsequently interpreted as the 1 D 1 3 -wave dominant state of charmonium. However, our analysis of the total-hadron cross sections indicates that the ψ ( 3770 ) is not a single state, but a complex system composed of the R ( 3760 ) , R ( 3780 ) , and R ( 3810 ) resonances. Among these, we interpret the R ( 3780 ) is a resonance dominated by the 1 D 1 3 charmonium state. Published by the American Physical Society 2024

Ablikim, M.

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie

On interpretation of fluctuations of conserved charges at high T

Abstract Fluctuations of conserved charges calculated on the lattice which can be measured experimentally, are well reproduced by a hadron resonanse gas model at temperatures below$$T_{ch} \sim 155$$ T ch ∼ 155 MeV and radically deviate from the hadron resonance gas predictions above the chiral restoration crossover. This behaviour is typically interpreted as an indication of deconfinement in the quark-gluon plasma regime. We present an argument that this interpretation may be too simple. The argument is based on the scaling of quantities with the number of colors: demonstration of deconfinement and QGP requires observable that is sensitive to$$\sim N_c^2$$ ∼ N c 2 gluons while the conserved charges are sensitive only to quarks and above$$T_{ch}$$ T ch scale as$$N_c^1$$ N c 1 . The latter scaling is consistent with the existence of an intermediate regime characterized by restored chiral symmetry and by approximate chiral spin symmetry which is a symmetry of confining interaction. In this regime the energy density, pressure and entropy density scale as$$N_c^1$$ N c 1 . In the large$$N_c$$ N c limit this regime might become a distinct phase separated from the hadron gas and from QGP by phase transitions. A natural observable that associates with deconfinement and is directly sensitive to deconfined$$N_c^2-1$$ N c 2 - 1 gluons is the Polyakov loop; in the$$N_c=3$$ N c = 3 world it remains very close to 0 at temperatures well above chiral crossover, reaches the value$$\sim 0.5$$ ∼ 0.5 around$$\sim 3T_{ch}$$ ∼ 3 T ch and the value close to 1 at temperatures$$\sim 1$$ ∼ 1 GeV.

Physics

Measurement of high-mass $t\bar{t}\ell ^{+}\ell ^{-}$ production and lepton flavour universality-inspired effective field theory interpretations at $\sqrt{s}=13$ $\text {T}\text {e}\hspace{-1.00006pt}\text {V}$ with the ATLAS detector

Measurements of $t\bar{t}\ell ^{+}\ell ^{-}$ production in the region of high dilepton invariant mass with effective field theory (EFT) interpretations are presented. They are performed using final states with three isolated leptons (electrons or muons) and are based on $\sqrt{s} = 13$ TeV proton–proton collision data with an integrated luminosity of $140\,\textrm{fb}^{-1}$, recorded from 2015 to 2018 with the ATLAS detector at the Large Hadron Collider. Measurements of the $t\bar{t}\ell ^{+}\ell ^{-}$ signal strength and cross-section upper-limits are performed inclusively in lepton flavour and separately for electrons and muons. The study also aims to probe anomalous four-fermion interactions including to test for possible lepton flavor universality violation. No significant deviations from the Standard Model predictions are observed and the measurements are interpreted through the EFT formalism to provide new constraints on the relevant operators.

Aad, G. [CNRS/IN2P3] (ORCID:0000000266654934)

Poster Abstract: Leveraging Large Language Models to Reveal Interpretable Cooling Behaviors from Smart Thermostat Data

Frequent heatwaves and hot summers increasingly challenge occupant comfort, health, and energy grid stability. Addressing these challenges requires a detailed understanding of household cooling behaviors, such as thermostat adjustments and adaptive responses to extreme conditions. Traditional analyses often rely on aggregated numerical metrics that overlook subtle but important household-specific variations. In this study, we introduce a generalizable methodology that integrates large language models (LLMs) with vision capabilities to enable scalable and detailed analysis of residential thermostat data. Using Ecobee's Donate Your Data (DYD) dataset—which provides five-minute records of indoor temperatures, thermostat setpoints, and HVAC runtimes—we focus on two U.S. cities with contrasting summer climates : Austin (TX) and Phoenix (AZ). Because raw time-series data are not well suited for direct LLM analysis, we transform them into visual representations, such as daily indoor temperature trajectories and weekly runtime histograms, to better capture behavioral variations. Leveraging LLMs' visual interpretation, we extract descriptive behavioral features, including temperature preferences, time-of-day cooling orientation, anticipatory versus reactive heatwave responses, and behavioral consistency. These semantic features support unsupervised clustering to identify distinct occupant archetypes at scale, revealing differences—such as morning-centric anticipatory coolers versus households that shift toward warmer setpoints during heatwaves—that can inform demand response, resilience planning, and health-aware interventions. By converting raw numerical data into interpretable behavioral patterns, this methodology enables scalable and practical analysis of occupant behavior, supporting actionable insights for comfort, resilience, and energy management.

Nihar, Kopal

Physical Interpretation of Early Battery Life Prediction Models

Early battery life prediction models are most useful for R&D if they help us understand the early changes in battery electrochemical response that correspond with long-term degradation and failure. Linear regression models such as Fused lasso and Partial Least Squares can fit coefficients directly to high-dimensional electrochemical data like capacity-voltage and ΔV–state-of-charge, i.e., Q(V) and ΔV(SOC) curves, learning coefficients that can be physically interpreted. We leverage the ISU-ILCC battery aging data set to learn high-dimensional coefficients for early battery life prediction from traditional slow-rate capacity check data, demonstrating learning on Q(V), d Q· d V −1 , and ΔV(SOC) curves. A thorough study on the dependence of coefficient values on train/test size and data preprocessing methods is made, demonstrating the reliability of high-dimensional regression approaches unless very small amounts of data are used for model training. For this data set, coefficients from Q(V) and d Q· d V −1 models highlight changes in electrode stoichiometry due to lithium loss, while ΔV(SOC) coefficients highlight changes in positive electrode diffusivity due to particle cracking as well as electrode stoichiometry shifts. By directly interpreting the coefficients of a regression model, we make physical insights into battery degradation mechanisms without requiring the assumptions of traditional battery data analysis methods.

25 ENERGY STORAGE

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES