SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Statistical Study of the Nonthermal Continuum Radiation Beaming Angle Measured by the High Frequency Receiver on Van Allen Probes-A

The nonthermal continuum (NTC) radiation beaming angle is computed over the entire Van Allen Probes-A mission when the spacecraft was in the dawn sector. The conditions in the dawn sector are favorable for the wave vector to lie near/in the spacecraft's spin plane allowing a favorable estimate of the beaming angle, and the dawn sector is also advantageous in that previous studies show NTC occurrence to peak in this sector. We found that scatter plots, over the entire mission, of beaming angle versus magnetic latitude form a distinct inverted-V pattern, with the apex at/near the magnetic equator. This pattern was sharpest for frequencies ( f ) ≲ 100 kHz. Using the NTC beaming formula from the linear mode conversion theory (LMCT), we show that such an inverted-V pattern is expected due to the large variation in the plasmapause location over the entire mission. The theoretical derived pattern qualitatively reproduces the observed pattern but not quantitatively. The lack of quantitative agreement is discussed and is attributed to several factors, one factor is off-centered emissions from the radio window. The qualitative agreement strongly supports LMCT as being the dominant mechanism generating NTC for f ≲ 100 kHz. For f ≳ 100 kHz, the inverted-V pattern becomes less distinct, and strong near-equatorial beaming is observed. Finally, after considering contamination of our selections by left-handed polarized AKR, our study suggests that besides LMCT another unidentified NTC generation mechanism becomes important for f ≳ 100 kHz.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗

DOE OSTI2026

MULTI-FIDELITY MODELING AND UNCERTAINTY QUANTIFICATION OF INVERTER BASED RESOURCES IN INTEGRATED T&D SYSTEMS

Uncertainty quantification plays a pivotal role in improving the accuracy and reliability of inverter operation within modern power systems that are increasingly dominated by inverter-based resources (IBRs). IBRs, especially those operating under grid forming (GFM) control, rely heavily on a complex set of control parameters and system measurements to maintain voltage, frequency, and power balance. Traditional deterministic modeling approaches often fail to capture these parameter deviations, potentially resulting in suboptimal control actions, reduced system stability, or even instability under high penetration of IBRs. In this paper, we demonstrate the application of model calibration and uncertainty quantification (UQ) principles to an integrated transmission and distribution (T&D) model involving a GFM converter and provide a framework for prioritizing control improvements, guiding robust design, and informing adaptive strategies that can accommodate real-time variability in system conditions. The proposed approach could be valuable in enhancing the robustness of current and future power systems under increased IBR penetrations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Chemical and Structural Changes in Corn Stover After Ensiling: Influence on Bioconversion

Production of biofuels, bioproducts, and bioenergy requires a well-characterized, stable, and reasonably uniform biomass supply and well-established supply chains for shipping biomass from farm fields to biorefineries, while achieving year-round production targets. Preserving and stabilizing biomass feedstock during storage is a necessity for cost-effective and sustainable biofuel production. Ensiling is a common storage method used to preserve and even improve forage quality; however, the impact of ensiling on biomass physical and chemical properties that influence bioconversion processes has been variable. Our objective in this work was to determine the effects of ensiling on lignocellulosic feedstock physicochemical properties and how that influences bioconversion requirements. We observed statistically significant decreases (p < 0.05) in the content of two major structural carbohydrates (glucan and xylan) of 5 and 8%, respectively, between the ensiled and non-ensiled materials. We were unable to detect differences in sugar yields from structural carbohydrates after pretreatment and enzymatic hydrolysis of the ensiled materials compared to non-ensiled controls. Based on this work, we conclude that ensiling the corn stover did not change the bioconversion requirements compared to the control samples and incurred losses of structural carbohydrates. At the light microscopy level, ensiled corn stover exhibited little structural change or relocation of cell wall components as detected by immunocytochemistry. However, more subtle structural changes were revealed by electron microscopy, as ensiled cell walls exhibit ultrastructural characteristics such as wall delimitation intermediate between non-ensiled and dilute-acid-pretreated cell walls. These findings suggest that alternative methods of conversion, such as deacetylation and mechanical refining, could take advantage of lamellar defects and may be more effective than dilute acid or hot water pretreatment for biomass conversion of ensiled materials.

09 BIOMASS FUELS↗

DOE OSTI2021

Complete polarization of electronic spins in OLEDs

At low temperatures and high magnetic fields, electron and hole spins in an organic light-emitting diode become polarized so that recombination preferentially forms molecular triplet excited-state species. For low device currents, magnetoelectroluminescence perfectly follows Boltzmann activation, implying a virtually complete polarization outcome. As the current increases, the magnetoelectroluminescence effect is reduced because spin polarization is suppressed by the reduction in carrier residence time within the device. Under these conditions, an additional field-dependent process affecting the spin-dependent recombination emerges, possibly related to the build-up of triplet excitons and their interaction with free charge carriers. Suppression of the EL alone does not prove electronic spin polarization. We therefore probe changes in the spin statistics of recombination directly in a dual singlet-triplet emitting material, which shows a concomitant rise in phosphorescence intensity as fluorescence is suppressed. Finite spin-orbit coupling in these materials gives rise to a microscopic distribution in effective g-factors of electrons and holes, Δg, i.e., a distribution in Larmor frequencies. This Δg effect in the pair, which mixes singlet and triplet, further suppresses singlet-exciton formation at high fields in addition to thermal spin polarization of the individual carriers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

AGC-2 Specimen Post Irradiation Data Package Report

This report documents results of the post-irradiation examination material property testing of the creep, control, and piggyback specimens from the irradiation creep capsule Advanced Graphite Creep (AGC)-2 are reported. This is the second of a series of six irradiation test trains planned as part of the AGC experiment to fully characterize the neutron irradiation effects and radiation creep behavior of current nuclear graphite grades. The AGC-2 capsule was irradiated in the Idaho National Laboratory Advanced Test Reactor at a nominal temperature of 600°C and to a peak dose of 5 dpa (displacements per atom). One half of the creep specimens were subjected to mechanical stresses (an applied stress of either 13.8, 17.2, or 20.7 MPa) to induce irradiation creep. All post-irradiation testing and measurement results are reported with the exception of the irradiation mechanical strength testing, which is the last destructive testing stage of the irradiation testing program. Material property tests were conducted on specimens from 15 nuclear graphite grades using a similar loading configuration as the first AGC capsule (AGC-1) to provide easy comparison between the two capsules. However, AGC-2 contained an increased number of specimens (i.e., 487 total specimens irradiated) and replaced specimens of the minor grade 2020 with the newer grade 2114. The data reported include specimen dimensions for both stressed and unstressed specimens to establish the irradiation creep rates, mass and volume data necessary to derive density, elastic constants (Young’s modulus, shear modulus, and Poisson’s ratio) from ultrasonic time of flight velocity measurements, Young’s modulus from the fundamental frequency of vibration, electrical resistivity, and thermal diffusivity and thermal expansion data from 100–500°C. No data outliers were determined after all measurements were completed. A brief statistical analysis was performed on the irradiated data and a limited comparison between pre- and post-irradiation properties is presented. A more complete evaluation of trends in the material property changes, as well as irradiation-induced creep due to irradiation, temperature, and applied load on specimens will be discussed in later AGC-2 post-irradiation examination analysis reports.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Stochastic modeling and statistical calibration with model error and scarce data

This paper introduces a procedure to assess the predictive accuracy of stochastic models subject to model error and sparse data. Model error is introduced as uncertainty on the coefficients of appropriate polynomial chaos expansions (PCE). The error associated with finite sample size allows us to conceive of these coefficients as statistics of the data that we describe as random variables whose influence on output quantities of interest is evaluated through the extended polynomial chaos expansion (EPCE). A Bayesian data assimilation scheme is introduced to update these expansions by considering the resulting nested chaos expansion as a hierarchical probabilistic model. Stochastic models of quantities of interest (QoI) are thus constructed and efficiently evaluated. Here, the Metropolis–Hastings Markov chain Monte Carlo procedure is used to sample the posterior. Two illustrative analytical and numerical problems are used to demonstrate the proposed approach.

Bayesian inference↗

DOE OSTI2020

Photo-Production of eta Mesons

In this work we present the GlueX experiment extensive high statistics measurements of the eta meson photo-production differential cross sections at beam energies between 3 GeV and 11.6 GeV. We compared the low energy (2.9 GeV< Ey < 5.9 GeV) differential cross section (d-sigma/d-Omega) results to previous experimental results from the Continuous Electron Beam Accelerator Facility (CEBAF) Large Acceptance Spectrometer (CLAS) collaboration. These (d-sigma/d-Omega) results are also compared to the theoretical isobar model EtaMAID2018 by Tiator et. al., and the Regge model by Nys et. al. The differential cross section d-sigma/dt is also presented and compared to the previous experimental results, and to the above models, in addition to a model by J. M. Laget. The high energy (6.2 GeV< Ey < 11.6 GeV) cross section results (d-sigma/d-Omega, and d-sigma/dt) are compared to the isobar model EtaMAID2018. These results will provide new constraints on eta photo-production mechanisms and high energy t-channel processes.

Kamel, Mahmoud↗

DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗

DOE OSTI2024

Electrical Responses of Modified Mineral Surfaces as Observed With Spectral Induced Polarization and Atomic Force Microscopy

Abstract Atomic force microscopy (AFM) and spectral induced polarization (SIP) are widely used to investigate the electrical properties of mineral surfaces at vastly different scales of measurement. We compare AFM and SIP measurements made on two different materials (glass beads and silica gel) subjected to etching, deposition of iron oxide particles, and inclusion of calcite grains. We found that the treatments produced qualitatively consistent behaviors in the AFM and SIP data. Direct AFM measurements of surface charge density for silica and calcite surfaces were quantitatively compared to values estimated from the SIP results using a grain polarization model. No statistically significant difference (at a 95% confidence level) was found between the surface charge density of silica estimated by AFM (2.3 ± 6.6 mC/m 2 for glass beads and 1.6 ± 0.1 mC/m 2 for silica gel) versus SIP (5.4 ± 4.4 mC/m 2 for glass beads and 1.6 ± 0.5 mC/m 2 for silica gel). The surface charge density for calcite determined by AFM (43.5 ± 12.9 mC/m 2 ) was approximately 19 times higher than that found for silica. While the charge density of calcite surfaces determined by SIP was also generally higher than that found for silica, different treatments produced significantly different values between 4.7 and 258 mC/m 2 (with a maximum 95% CI of ±8.7 mC/m 2 ). Several possible explanations exist for the range of the observed SIP measurements, including aging of the calcite surfaces. Overall, this study suggests the potential for the complementary use of AFM and SIP measurements to constrain future investigations of polarization mechanisms in porous media.

Geochemistry & Geophysics↗

DOE OSTI2020

Comparative differential cuproproteomes of Rhodobacter capsulatus reveal novel copper homeostasis related proteins

Abstract Copper (Cu) is an essential, but toxic, micronutrient for living organisms and cells have developed sophisticated response mechanisms towards both the lack and the excess of Cu in their environments. In this study, we achieved a global view of Cu-responsive changes in the prokaryotic model organism Rhodobacter capsulatus using label-free quantitative differential proteomics. Semi-aerobically grown cells under heterotrophic conditions in minimal medium (∼0.3 μM Cu) were compared with cells supplemented with either 5 μM Cu or with 5 mM of the Cu-chelator bathocuproine sulfonate. Mass spectrometry based bottom-up proteomics of unfractionated cell lysates identified 2430 of the 3632 putative proteins encoded by the genome, producing a robust proteome dataset for R. capsulatus. Use of biological and technical replicates for each growth condition yielded high reproducibility and reliable quantification for 1926 of the identified proteins. Comparison of cells grown under Cu-excess or Cu-depleted conditions to those grown under minimal Cu-sufficient conditions revealed that 75 proteins exhibited statistically significant (p < 0.05) abundance changes, ranging from 2- to 300-fold. A subset of the highly Cu-responsive proteins was orthogonally probed using molecular genetics, validating that several of them were indeed involved in cellular Cu homeostasis.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Short-term warming increased soil heterotrophic respiration due to enhanced active microbial community

Soil microbes, particularly active microbes, play a crucial role in conserving soil carbon under climate change, especially in forest ecosystems, constituting over 50% of global soil organic carbon. Nevertheless, the response mechanisms of active microbial community to climate warming and their subsequent impacts on soil heterotrophic respiration (Rh) remain insufficiently understood. To resolve this mechanistic uncertainty, we implemented a 3-year soil translocation experiment to investigate soil warming effects on active microbial community and Rh. We used 16S rRNA gene amplicon sequencing, Metatranscriptomics sequencing, coupled with correlation analyses to explore the linkages between warming-induced shifts in Rh and the active microbial community. Our results demonstrated that warming of soil induced a 55% enhancement in Rh. Concurrently, the Shannon index and Richness of active microbial community increased by 20% and 117%, respectively. Warming significantly altered active microbial community composition, inducing a compositional shift characterized by 0.5-to 2-fold increases in the abundances of Proteobacteria, Chloroflexi, and Actinobacteria. Our study also revealed 92–100% increases in the abundance of C-degradation functional genes encoding starch, cellulose, and lignin decomposition pathways within active microbial community under warming. Statistical analyses identified significant positive correlations between Rh and the abundance of Proteobacteria and Actinobacteria, as well as expression levels of functional genes associated with lignin and cellulose decomposition pathways. Furthermore, our results suggested that short-term warming increased Rh through altering diversity, species composition, and C degradation functional genes of active microbial community providing insights into the influence of microbial communities on soil C-climate feedbacks under climate warming.

Active microbial community↗

DOE OSTI2020

Dynamic microbial assembly processes correspond to soil fertility in sustainable paddy agroecosystems

1. Microorganisms underpin the sustainability of agroecosystems, but we have a poor understanding of how the assembly processes that govern microbial community structure interact through time to maintain soil fertility. 2. Ancient rice paddies have a history of cultivation spanning hundreds or thousands of years and provide an ideal opportunity to explore relationships between community assembly processes and long-term agroecosystem sustainability. We used statistical null models to evaluate temporal patterns of bacterial community assembly processes across rice cultivation years and their associations with increases in soil fertility in paddy agroecosystems. 3. We found that stochasticity was a persistent force structuring paddy soil microbial communities, putatively due to frequent flooding that facilities dispersal. However, determinism gradually increased over centuries of agriculture and corresponded to increases in soil fertility. Accordingly, multifunctionality, calculated using enzymes catalyzing key soil functions in carbon, nitrogen, phosphorous and sulfur cycling significantly increased along rice cultivation years. We therefore propose that balanced assembly processes preserve a diverse array of traits that can dominate under various circumstances (stochasticity) while allowing for traits relevant to soil fertility to be selected for by soil fertility (determinism), thus supporting multifunctionality that underlies sustainable agroecosystems. 4. In total, our results highlight that knowledge on temporal microbial assembly processes, community composition, and ecological function in context of soil fertility is key to understanding ecological mechanisms maintaining agroecosystems sustainability.

Liu, Wenjing↗

DOE OSTI2025

Widening of Wind Stress Anomalies Amplifies ENSO in a Warming Climate

Abstract Climate change simulations generally indicate the strengthening of El Niño–Southern Oscillation (ENSO) sea surface temperature (SST) variability through the twenty-first century, yet a robust physical mechanism explaining this change across different models is still lacking, and the projections for ENSO amplitude exhibit a large spread. Most commonly, changes in the background state of the tropical Pacific are invoked to explain these changes of ENSO. Here, we show that changes in the structure of wind stress anomalies associated with ENSO are potentially as important as these background state changes. Specifically, changes in the magnitude, meridional width, and zonal structure of wind stress anomalies can explain approximately 53% of the intermodel variance in the projected change of ENSO magnitude through the twenty-first century as well as 43% in ENSO periodicity changes. Among these changes in the wind structure, the meridional widening of wind anomalies plays the most important role. To demonstrate that these changes are indeed critical, we develop a hybrid model of ENSO based on the Community Earth System Model, version 2, which incorporates a dynamical ocean coupled to a simplified statistical atmosphere within the tropical Pacific. In the absence of external forcing and corresponding mean-state changes, the imposed changes in wind stress anomalies in this hybrid model result in an increase of ENSO amplitudes of nearly 10% along the equator. Our results are also theoretically supported by a recharge-oscillator model that incorporates the meridional wind structure. Thus, changes in the structure of wind stress anomalies, together with changes in the mean state, likely play a critical role in the projected strengthening of ENSO.

Stuivenvolt-Allen, Jacob [Yale University, New Hav↗

DOE OSTI2022

Why it is Unfortunate that Linear Machine Learning “Works” so well in Electromechanical Switching of Ferroelectric Thin Films

Machine learning (ML) is relied on for materials spectroscopy. It is challenging to make ML models fail because statistical correlations can mimic the physics without causality. Here, using a benchmark band-excitation piezoresponse force microscopy polarization spectroscopy (BEPS) dataset the pitfalls of the so-called “better”, “faster”, and “less-biased” ML of electromechanical switching are demonstrated and overcome. Using a toy and real experimental dataset, it is demonstrated how linear nontemporal ML methods result in physically reasonable embedding (eigenvalues) while producing nonsensical eigenvectors and generated spectra, promoting misleading interpretations. A new method of unsupervised multimodal hyperspectral analysis of BEPS is demonstrated using long-short-term memory (LSTM) β-variational autoencoders (β-VAEs) . By including LSTM neurons, the ordinal nature of ferroelectric switching is considered. Further, to improve the interpretability of the latent space, a variational Kullback–Leibler-divergency regularization is imposed . Finally, regularization scheduling of β as a disentanglement metric is leveraged to reduce user bias. Combining these experiment-inspired modifications enables the automated detection of ferroelectric switching mechanisms, including a complex two-step, three-state one. Ultimately, this work provides a robust ML method for the rapid discovery of electromechanical switching mechanisms in ferroelectrics and is applicable to other multimodal hyperspectral materials spectroscopies.

36 MATERIALS SCIENCE↗