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 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2020

2D Ising Model for Adsorption–induced Enantiopurification of Racemates

Mechanisms for the spontaneous transformation of achiral chemical systems into states of enantiomeric purity have important ramifications in modern pharmacology and potential relevance to the origins of homochirality in life on Earth. Such mechanisms for enantiopurification are needed for production of chiral pharmaceuticals and other bioactive compounds. Previously proposed chemical mechanisms leading from achiral systems to near homochirality are initiated by a symmetry–breaking step resulting in a minor excess of one enantiomer via statistical fluctuations in enantiomer concentrations. Subsequent irreversible processes then amplify the majority enantiomer concentration while simultaneously suppressing minority enantiomer production. Herein, equilibrium adsorption of amino acid enantiomer mixtures onto chiral and achiral surfaces reveals amplification of surface enantiomeric excess relative to the gas phase; i. e. enantiopurification of chiral adsorbates by adsorption. Furthermore, this adsorption–induced amplification of enantiomeric excess is shown to be well–describe by the 2D Ising model. More importantly, the 2D–Ising model predicts formation of homochiral monolayers from adsorption of racemic mixtures or prochiral molecules on achiral surfaces; i.e. enantiopurification with no apparent chiral driving force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

A Bayesian Framework for Spectral Reprojection

Abstract Fourier partial sum approximations yield exponential accuracy for smooth and periodic functions, but produce the infamous Gibbs phenomenon for non-periodic ones. Spectral reprojection resolves the Gibbs phenomenon by projecting the Fourier partial sum onto a Gibbs complementary basis, often prescribed as the Gegenbauer polynomials. Noise in the Fourier data and the Runge phenomenon both degrade the quality of the Gegenbauer reconstruction solution, however. Motivated by its theoretical convergence properties, this paper proposes a new Bayesian framework for spectral reprojection, which allows a greater understanding of the impact of noise on the reprojection method from a statistical point of view. We are also able to improve the robustness with respect to the Gegenbauer polynomials parameters. Finally, the framework provides a mechanism to quantify the uncertainty of the solution estimate.

Li, Tongtong (ORCID:0000000276644764)↗

DOE OSTI2025

Critical impact of experimentally-driven strut level anisotropic material models in advanced stress analysis of additively manufactured lattice structures

The rapid acceleration in materials discovery may overshadow the importance of thoroughly understanding the mechanical performance of newly developed materials in demanding environments. The recent interest in combining parametric studies with machine learning techniques to explore how changes in specific processing parameters or model inputs affect the overall behavior of a material system can only be truly beneficial if the governing constitutive relations describing material behavior are accurately established. In this study, we demonstrate the critical impact of accurately representing strut-level anisotropic material behavior in advanced stress analysis of additively manufactured lattice structures (AMLS). We introduce a systematic experimental and modeling approach for developing strut-level anisotropic elastoplastic material models that account for the influence of microstructural features such as porosity, texture, and surface roughness on the development of local anisotropic mechanical properties, which vary with strut orientation relative to the build direction (BD). As a result the presented material model captures and relates the statistics of spatially varying struts’ microstructural features to the local stress distribution. Our findings suggest that incorporating strut-level anisotropic material behavior into unit cell analysis significantly influences the load distribution and evolution of local stresses within the structure. Therefore, accounting for this anisotropy is critical for developing an understanding of unit cell behavior and performance, including subsequent topology/component design optimization based on this analysis.

Sahoo, Subhadip [University of Arizona]↗

DOE OSTI2024

The Impact of Lithium Anode Interface on Capacity Fade in Polymer Electrolyte-Based Solid-State Batteries

This study investigates the Li stripping-plating morphology and failure mechanisms in full cells consisting of a solid polymer electrolyte (SPE) with two commercial Li anodes: Li chip and Li foil. The primary identified failure mechanism of the SPE cell is capacity fade, regardless of the Li manufacturer. While the cathode’s role in capacity fade is evident, the Li anode significantly influences cycling performance, with Li foil cells cycling 50% longer than Li chip cells, a statistical difference. Further, post-mortem scanning electron microscopy and X-ray photoelectron spectroscopy results attribute the Li chip’s faster capacity fade to a loss of contact and continuous growth of the solid electrolyte interphase (SEI). Conversely, Li foil maintains consistent contact with the solid polymer, displaying a thin and stable SEI. Additionally, failure mechanisms between a gel electrolyte in previous work and the dry SPE are compared.

25 ENERGY STORAGE↗

DOE OSTI2021

Machine learning to accelerate screening for Marcus reorganization energies

Understanding and predicting the charge transport properties of π-conjugated materials is an important challenge for designing new organic electronic devices, such as solar cells, plastic transistors, light-emitting devices, and chemical sensors. A key component of the hopping mechanism of charge transfer in these materials is the Marcus reorganization energy which serves as an activation barrier to hole or electron transfer. While modern density functional methods have proven to accurately predict trends in intramolecular reorganization energy, such calculations are computationally expensive. In this work, we outline active machine learning methods to predict computed intramolecular reorganization energies of a wide range of polythiophenes and their use toward screening new compounds with low internal reorganization energies. Our models have an overall root mean square error (RMSE) of ±0.113 eV, but a much smaller RMSE of only ±0.036 eV on the new screening set. Since the larger error derives from high-reorganization energy compounds, the new method is highly effective to screen for compounds with potentially efficient charge transport parameters.

14 SOLAR ENERGY↗

DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Asymmetries in cloud microphysical properties ascribed to sea ice leads via water vapour transport in the central Arctic

To investigate the influence of sea ice openings like leads on wintertime Arctic clouds, the air mass transport is exploited as a heat and humidity feeding mechanism which can modify Arctic cloud properties. Cloud microphysical properties in the central Arctic are analysed as a function of sea ice conditions during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020. The Cloudnet classification algorithm is used to characterize the clouds based on remote sensing observations and the atmospheric thermodynamic state from the observatory on board the research vessel (RV) Polarstern. To link the sea ice conditions around the observational site with the cloud observations, the water vapour transport (WVT) being conveyed towards RV Polarstern has been utilized as a mechanism to associate upwind sea ice conditions with the measured cloud properties. This novel methodology is used to classify the observed clouds as coupled or decoupled to the WVT based on the location of the maximum vertical gradient of WVT height relative to the cloud-driven mixing layer. Only a conical sub-sector of sea ice concentration (SIC) and the lead fraction (LF) centred on the RV Polarstern location and extending up to 50 km in radius and with an azimuth angle governed by the time-dependent wind direction measured at the maximum WVT is related to the observed clouds. We found significant asymmetries for cases when the clouds are coupled or decoupled to the WVT and selected by LF regimes. Liquid water path of low-level clouds is found to increase as a function of LF, while the ice water path does so only for deep precipitating systems. Clouds coupled to WVT are found to generally have a lower cloud base and larger thickness than decoupled clouds. Thermodynamically, for coupled cases the cloud-top temperature is warmer and accompanied by a temperature inversion at the cloud top, whereas the decoupled cases are found to be closely compliant with the moist adiabatic temperature lapse rate. The ice water fraction within the cloud layer has been found to present a noticeable asymmetry when comparing coupled versus decoupled cases. This novel approach of coupling sea ice to cloud properties via the WVT mechanism unfolds a new tool to study Arctic surface–atmosphere processes. With this formulation, long-term observations can be analysed to enforce the statistical significance of the asymmetries. Furthermore, our results serve as an opportunity to better understand the dynamic linkage between clouds and sea ice and to evaluate its representation in numerical climate models for the Arctic system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1 -2 kcal · mol -1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal · mol -1 and -6.3 to +3.8 kcal · mol -1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. In conclusion, the results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Overall Performance Losses and Activated Mechanisms in Double Glass and Glass-backsheet Photovoltaic Modules with Monofacial and Bifacial PERC Cells, under Accelerated Exposures

Commercial PV modules have various packaging choices nowadays, which influence their long-term reliability. This study compared the degradation behaviors of sixteen module variants from two brands with varying encapsulant materials (EVA or POE), encapsulant types, module architectures (GB or DG), and cell types (monofacial or bifacial) using null hypothesis testing to determine statistical significant findings. The modules were exposed for 2,520 hours under two accelerated exposures: modified damp heat (mDH) and modified damp heat with full-spectrum light (mDH+FSL). For both brands, two DG module variants with UV-Cutoff rear encapsulant are found to have significantly lower average power loss than the module variants of EVA+GB with opaque rear encapsulant after each accelerated exposure. Metallization interconnect corrosion is identified as the primary degradation mechanism. Furthermore, unsupervised hierarchical clustering finds that the degradation behaviors of modules from one brand with a more strict manufacturing quality control depends on module architectures only.

14 SOLAR ENERGY↗

DOE OSTI2023

Coherent Random Lasing in Subwavelength Quasi‐2D Perovskites

Abstract Quasi‐2D lead halide perovskites have garnered increasing interest as lasing gain media. Relatively simple fabrication, high refractive index, and unique quantum well structure encourage their use in traditional cavity lasers and cavity‐free systems called random lasers (RLs). Despite tremendous advances reported thus far, coherent random lasing in quasi‐2D perovskite subwavelength films has not been reported. Consequently, coherent optical feedback mechanisms in quasi‐2D perovskite systems are still unexplored. Here, this work reports the observation of coherent random lasing in subwavelength quasi‐2D perovskite films. Statistical analysis of spectral measurements reveals Lévy‐like intensity fluctuations, replica symmetry breaking confirms random lasing, and the coherent modes are studied with spectral and spatial correlation techniques. The observed coherent lasing modes are found to be extended states that arise from the random crystal grain structure during fabrication and span the entire pump volume. These modes out‐compete diffusive lasing due to their coherence.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

DOE OSTI2025

Measurement of the branching fraction of Λ + 𝑐 →𝑝⁢𝐾 0 𝑆⁢ 𝜋 0 at Belle

We report a precise measurement of the ratio of branching fractions ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 )/ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ) using 980 fb −1 of 𝑒 + ⁢𝑒 − data from the Belle experiment. We obtain a value of ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 )/ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ) =0.339 ±0.002 ±0.009, where the first and second uncertainties are statistical and systematic, respectively. This Belle result is consistent with the previous measurement from the CLEO experiment but has a fivefold improvement in precision. By combining our result with the world average ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ), we obtain the absolute branching fraction ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 ) =(2.12 ± 0.01 ±0.05 ±0.10)%, where the uncertainties are statistical, systematic, and the uncertainty in the absolute branching fraction scale ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ), respectively. This measurement can shed light on hadronic decay mechanisms in charmed baryon decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Measurement and memory in the periodically driven complex Ginzburg-Landau equation

In the present work we illustrate that classical but nonlinear systems may possess features reminiscent of quantum ones, such as memory, upon suitable external perturbation. As our prototypical example, we use the two-dimensional complex Ginzburg-Landau equation in its vortex glass regime. Here, we impose an external drive as a perturbation mimicking a quantum measurement protocol, with a given “measurement rate” (the rate of repetition of the drive) and “mixing rate” (characterized by the intensity of the drive). Using a variety of measures, we find that the system may or may not retain its coherence, statistically retrieving its original glass state, depending on the strength and periodicity of the perturbing field. The corresponding parametric regimes and the associated energy cascade mechanisms involving the dynamics of vortex waveforms and domain boundaries are discussed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Experimentally Inferred Fusion Yield Dependencies of OMEGA Inertial Confinement Fusion Implosions

Statistical modeling of experimental and simulation databases has enabled the development of an accurate predictive capability for deuterium-tritium layered cryogenic implosions at the OMEGA laser. Here, a physics-based statistical mapping framework is described and used to uncover the dependencies of the fusion yield. This model is used to identify and quantify the degradation mechanisms of the fusion yield in direct-drive implosions on OMEGA. The yield is found to be reduced by the ratio of laser beam to target radius, the asymmetry in inferred ion temperatures from the $l$ = 1 mode, the time span over which tritium fuel has decayed, and parameters related to the implosion hydrodynamic stability. When adjusted for tritium decay and $l$ = 1 mode, the highest yield in OMEGA cryogenic implosions is predicted to exceed 2 × 10 14 fusion reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Exact quantum scars in the chiral nonlinear Luttinger liquid

While the chiral linear Luttinger liquid is integrable via bosonization, its nonlinear counterpart does not admit for an analytic solution. In this work, we find a subextensive number of exact eigenstates for a large family of density-density interaction terms. These states are embedded in a continuum of strongly correlated excited states. The real-space entanglement entropy of some exact states scales logarithmically with system size while that of others has volume-law scaling. We introduce momentum-space entanglement as an unambiguous differentiator between these exact states and the remaining excited states. With regard to momentum space, the exact states behave as bona fide quantum many-body scars: they exhibit identically zero momentum-space entanglement, while typical eigenstates behave thermally. We corroborate this finding by a level statistics analysis. Here we detail the general formalism for systematically finding all interaction terms and associated exact states, and present a number of infinite exact state sequences extending to arbitrarily high energies. Unlike many previous examples of quantum many-body scars, the exact states uncovered here do not lie at equidistant energies and do not follow from a special operator algebra. Instead, they are uniquely enabled by the interplay of Fermi statistics and chirality.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Self-organized criticality in multi-pulse gamma-ray bursts

The variability in multi-pulse gamma-ray bursts (GRBs) may help to reveal the mechanism of underlying processes from the central engine. To investigate whether the self-organized criticality (SOC) phenomena exist in the prompt phase of GRBs, we statistically study the properties of GRBs with more than 3 pulses in each burst by fitting the distributions of several observed physical variables with a Markov Chain Monte Carlo approach, including the isotropic energy E iso , the duration time T, and the peak count rate P of each pulse. Our sample consists of 454 pulses in 93 GRBs observed by the CGRO/BATSE satellite. The best-fitting values and uncertainties for these power-law indices of the differential frequency distributions are: α$^{d}_{E}$ = 1.54±0.09, α$^{d}_{T}$ = 1.82±$^{+0.14}_{0.15}$ and α$^{d}_{P}$ = 2.09$^{+0.18}_{–0.19}$, while the power-law indices in the cumulative frequency distributions are: α$^{c}_{E}$ = 1.44$^{+0.08}_{–0.10}$, α$^{c}_{T}$ = 1.75$^{+0.141}_{–0.13}$ and α$^{c}_{P}$ = 1.99$^{+0.16}_{–0.19}$. We find that these distributions are roughly consistent with the physical framework of a Fractal-Diffusive, Self-Organized Criticality (FD-SOC) system with the spatial dimension S = 3 and the classical diffusion β=1. The results of this study support that the jet responsible for the GRBs should be magnetically dominated and magnetic instabilities (e.g., kink model, or tearing-model instability) lead the GRB emission region into the SOC state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

En echelon faults reactivated by wastewater disposal near Musreau Lake, Alberta

We use machine-learning and cross-correlation techniques to enhance earthquake detectability by two magnitude units for the earthquake sequence near Musreau Lake, Alberta, which is induced by wastewater disposal. This deep catalogue reveals a series of en echelon ~N–S oriented strike-slip faults that are favourably oriented for reactivation. These faults require only ~0.6 MPa overpressure for triggering to occur. Earthquake activity occurs in bursts, or episodes; episodes restricted to the largest fault tend to have earthquakes starting near the southern end (distant from injectors) and progressing northwards (towards the injectors). While most events are concentrated along these ~N–S oriented faults, we also delineate smaller faults. Together, these findings suggest pore pressure as the triggering mechanism, where a time-dependent increase in pore pressure likely caused these faults to progressively reawaken. Analysis of the ‘next record-breaking event’, a statistical model that forecasts the sequencing of earthquake magnitudes, suggests that the next largest event would be M L ~4.3. The seismically illuminated length of the largest fault indicates potential magnitudes as large as M w 5.3.

58 GEOSCIENCES↗