Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical models”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Precision mass measurement of lightweight self-conjugate nucleus 80 Zr

Protons and neutrons in the atomic nucleus move in shells analogous to the electronic shell structures of atoms. The nuclear shell structure varies due to changes of the nuclear mean field with the number of neutrons N and protons Z, and these variations can be probed by measuring mass differences between nuclei. The N=Z=40 self-conjugate nucleus 80 Zr is of particular interest as its proton and neutron shell structures are expected to be very similar, and its ground state is highly deformed. Here, we provide evidence for the existence of a deformed double shell closure in 80Zr through high precision Penning trap mass measurements of 80-83 Zr. Our mass values show that 80Zr is significantly lighter, and thus more strongly bound than predicted. This can be attributed to the deformed shell closure at N=Z=40 and the large Wigner energy. A statistical Bayesian model mixing analysis employing several global nuclear mass models demonstrates difficulties with reproducing the observed mass anomaly using current theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Eigenvector centrality is a metric of elastomer modulus, heterogeneity, and damage

Here, we present an application of eigenvector centrality to encode the connectivity of polymer networks resolved at the micro- and meso-scopic length scales. This method captures the relative importance of different nodes within the network structure and provides a route toward the development of a statistical mechanics model that correlates connectivity with mechanical response. This scheme may be informed by analytical and semi-analytical models for the network structure, or through direct experimental examination. It may be used to predict the reduction in mechanical performance for heterogeneous materials subjected to specific modes of damage. Here, we develop the method and demonstrate that it leads to the prediction of established trends in elastomers. We also apply the model to the case of a self-healing polymer network reported in the literature, extracting insight about the fraction of bonds broken and re-formed during strain and recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dominant controls of cold-season precipitation variability over the high mountains of Asia

A robust understanding of the sub-seasonal cold season (November–March) precipitation variability over the High Mountains of Asia (HMA) is lacking. Here, we identify dynamic and thermodynamic pathways through which natural modes of climate variability establish their teleconnections over the HMA. First, we identify evaporative sources that contribute to the cold season precipitation over the HMA and surrounding areas. The predominant moisture contribution comes from the mid-latitude regions, including the Mediterranean/Caspian Seas and Mediterranean land. Second, we establish that several tropical and extratropical forcings display a sub-seasonally fluctuating influence on precipitation distribution over the region during the cold season. Many of them varyingly interact, so their impacts cannot be explained independently or at seasonal timescales. Lastly, a single set of evaporative sources is not identifiable as the key determinant in propagating a remote teleconnection because the sources of moisture anomalies depend on the pattern of sub-seasonally varying dynamical forcing in the atmosphere.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Deep Koopman operators for causal discovery

Causal discovery aims to identify cause-effect mechanisms for better scientific understanding, explainable decision-making, and more accurate modeling. Standard statistical frameworks, such as Granger causality, lack the ability to quantify causal relationships in nonlinear dynamics due to the presence of complex feedback mechanisms, timescale mixing, and nonstationarity. Thus, applying these methods to study causal dynamics in real-world systems, such as the Earth, is a major challenge. Addressing this shortcoming, we leverage deep learning and a Koopman operator-theoretic formalism to present a class of causal discovery algorithms. Kausal uses deep Koopman operator methods to approximate nonlinear dynamics in a linearized vector space in which traditional causal inference methods such as Granger causality can be more easily applied. Our idealized experiments demonstrate Kausal’s superior ability in discovering and characterizing causal signals compared to existing deep learning and non-deep learning state-of-the-art approaches. Finally, the successful identification of major El Niño and La Niña events in observations showcases Kausal’s skill to handle real-world applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Design principles for the ultimate gas deliverable capacity material: nonporous to porous deformations without volume change

Understanding the fundamental limits of gas deliverable capacity in porous materials is of critical importance as it informs whether technical targets (e.g., for on-board vehicular storage) are feasible. High-throughput screening studies of rigid materials, for example, have shown they are not able to achieve the original ARPA-E methane storage targets, yet an interesting question remains: what is the upper limit of deliverable capacity in flexible materials? In this work we develop a statistical adsorption model that specifically probes the limit of deliverable capacity in intrinsically flexible materials. The resulting adsorption thermodynamics indicate that a perfectly designed, intrinsically flexible nanoporous material could achieve higher methane deliverable capacity than the best benchmark systems known to date with little to no total volume change. Density functional theory and grand canonical Monte Carlo simulations identify a known metal–organic framework (MOF) that validates key features of the model. Therefore, this work (1) motivates a continued, extensive effort to rationally design a porous material analogous to the adsorption model and (2) calls for continued discovery of additional high deliverable capacity materials that remain hidden from rigid structure screening studies due to nominal non-porosity.

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

Exploration of lattice Hamiltonians for functional and structural discovery via Gaussian process-based exploration–exploitation

Statistical physics models ranging from simple lattice to complex quantum Hamiltonians are one of the mainstays of modern physics that have allowed both decades of scientific discovery and provided a universal framework to understand a broad range of phenomena from alloying to frustrated and phase separated materials to quantum systems. Traditionally, exploration of the phase diagrams corresponding to multidimensional parameter spaces of Hamiltonians was performed using a combination of basic physical principles, analytical approximations, and extensive numerical modeling. However, exploration of complex multidimensional parameter spaces is subject to the classic dimensionality problem, and the behaviors of interest concentrated on low dimensional manifolds remain undiscovered. Here, we demonstrate that a combination of exploration and exploration–exploitation with Gaussian process modeling and Bayesian optimization allows effective exploration of the parameter space for lattice Hamiltonians and effectively maps the regions at which specific macroscopic functionalities or local structures are maximized. We argue that this approach is general and can be further extended well beyond the lattice Hamiltonians to effectively explore the parameter space of more complex off-lattice and dynamic models.

42 ENGINEERING↗
DOE OSTI2022

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Temporal structure of blobs in NSTX

The time dependence of the blob pulse shape and the waiting time between blobs was found using data from the gas puff imaging diagnostic in NSTX. The database used was of 103 shots from 2010 as described in a previous paper (Zweben et al., Phys. Plasmas 29, 012505 (2022).]. The blob pulse shape was well fit by an exponential rise and fall where the average rise time was τ r = 9.0 ± 2.7 μs and the average fall time was τ f = 16.6 ± 5.8 μs. The waiting times between blob pulses above a threshold of three times the mean had a broad distribution with an average of τ w = 1.2 ± 0.85 ms over the database. The blob intermittency parameter γ b = τ d /τ w , where the blob pulse duration was τ d = τ r + τ f , ranged from γ b ~1% to 5% for shots in this database and increased almost linearly with the blob fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Incorporating Beam Attenuation Calculations into an Integrated Data Analysis Model for Ion Effective Charge

Here, we have created a forward model for charge-exchange impurity density measurements that incorporates neutral beam attenuation measurements self-consistently for the purpose of determining the ion-effective charge Z eff . The model is constructed within an integrated data analysis framework to include a self-consistent calculation of neutral beam attenuation due to multiple impurity species into the measurement of a single impurity density. The model includes measurements of the beam Doppler-shift spectrum and shine-through particle flux to determine the neutral beam particle density which is attenuated by ion collisions. Synthetic data are generated from the diagnostic forward model using statistical and calibration uncertainties. These “noisy” data are used in the analysis to evaluate how accurately Z eff is determined. Methods of experimental design are employed to calculate the information gained from different diagnostic combinations. The analysis shows that while attenuation measurements alone do not provide a unique impurity density measurement in the case of a multispecies inhomogeneous plasma, they do provide an effective measurement of the Z eff profile and place constraints on the impurity density profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

EVAP: A two-photon imaging tool to study conformational changes in endogenous Kv2 channels in live tissues

A primary goal of molecular physiology is to understand how conformational changes of proteins affect the function of cells, tissues, and organisms. Here, we describe an imaging method for measuring the conformational changes of the voltage sensors of endogenous ion channel proteins within live tissue, without genetic modification. We synthesized GxTX-594, a variant of the peptidyl tarantula toxin guangxitoxin-1E, conjugated to a fluorophore optimal for two-photon excitation imaging through light-scattering tissue. We term this tool EVAP (Endogenous Voltage-sensor Activity Probe). GxTX-594 targets the voltage sensors of Kv2 proteins, which form potassium channels and plasma membrane–endoplasmic reticulum junctions. GxTX-594 dynamically labels Kv2 proteins on cell surfaces in response to voltage stimulation. To interpret dynamic changes in fluorescence intensity, we developed a statistical thermodynamic model that relates the conformational changes of Kv2 voltage sensors to degree of labeling. We used two-photon excitation imaging of rat brain slices to image Kv2 proteins in neurons. We found puncta of GxTX-594 on hippocampal CA1 neurons that responded to voltage stimulation and retain a voltage response roughly similar to heterologously expressed Kv2.1 protein. Our findings show that EVAP imaging methods enable the identification of conformational changes of endogenous Kv2 voltage sensors in tissue.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A data-driven approach to real-time vertical position estimation for NSTX-U vertical stability control

In this paper, a database of 77 996 plasma equilibrium reconstructions from 727 discharges during the initial operation of the NSTX-U spherical tokamak is analyzed to develop a statistically robust model of the plasma vertical position for real-time control. A variety of regression models are developed and tested, ranging in complexity from linear models to deep neural networks, and including input signals ranging from the four pairs of flux loops used historically on NSTX-U up to the full set of 389 real-time signals available to the plasma control system. A linear model based on 140 real-time magnetics signals is found to offer excellent accuracy, with a coefficient of determination R 2 = 0.906. The robustness of this model to limited training data, new operating scenarios, and signal errors is tested, and a procedure is demonstrated to tune the model parameters to optimize its robustness. A time-dependent plasma equilibrium solver, TokaMaker, is used to simulate vertical stability control in NSTX-U, demonstrating that it should be possible to iteratively tune the parameters of a linear vertical position model to stabilize both positive and negative triangularity plasmas in future experiments.

magnetic diagnostics↗
DOE OSTI2022

Alternate Analysis of an In Vitro Solubility Study on the Lung Dissolution Rate of 238PuO2 Material Involved in the 2020 Incident at Los Alamos National Laboratory

Here this work presents an alternate analysis of an in vitro solubility study on the lung dissolution rate of 238 PuO 2 material involved in a recent inhalation incident at Los Alamos National Laboratory (LANL). The original dataset used in this work was retrieved from a recently published report. The present work shows an analysis of the same dataset by modeling the dissolution in separate time intervals rather than modeling the cumulative dissolution.

238Pu↗
DOE OSTI2025

Deviations from the Porter-Thomas Distribution due to Nonstatistical 𝛾 Decay below the 150 Nd Neutron Separation Threshold

We introduce a new method for the study of fluctuations of partial transition widths based on nuclear resonance fluorescence experiments with quasimonochromatic linearly polarized photon beams below particle separation thresholds. It is based on the average branching of decays of 𝐽=1 states of an even-even nucleus to the 2$^{+}_{1}$ state in comparison to the ground state. Between 5 and 7 MeV, a constant average branching ratio for 𝛾 decays from 1 − states of 0.490(16) is observed for the nuclide 150 Nd. Assuming 𝜒 2 -distributed partial transition widths, this average branching ratio is related to a degree of freedom of 𝜈 = 1.93⁢(12), rejecting the validity of the Porter-Thomas distribution, requiring 𝜈 = 1. The observed deviation can be explained by nonstatistical effects in the 𝛾-decay behavior with contributions in the range of 9.4(10)% up to 94(10)%.

150 ≤ A ≤ 189↗
DOE OSTI2019

Obstacles to quantum annealing in a planar embedding of XORSAT

We introduce a planar embedding of the k-regular k-XORSAT problem, in which solutions are encoded in the ground state of a classical statistical mechanics model of reversible logic gates arranged on a square grid and acting on bits that represent the Boolean variables of the problem. The special feature of this embedding is that the resulting model lacks a finite-temperature phase transition, thus bypassing the first-order thermodynamic transition known to occur in the random graph representation of XORSAT. In spite of this attractive feature, the thermal relaxation into the ground state displays remarkably slow glassy behavior. The question addressed in this paper is whether this planar embedding can afford an efficient path to solution of k-regular k-XORSAT via quantum adiabatic annealing. We first show that our model bypasses an avoided level crossing and consequent exponentially small gap in the limit of small transverse fields. We then present quantum Monte Carlo results for our embedding of the k-regular k-XORSAT that strongly support a picture in which second-order and first-order transitions develop at a finite transverse field for k = 2 and k = 3, respectively. Furthermore, this translates into power-law and exponential dependences in the scaling of energy gaps with system size, corresponding to times-to-solution which are, respectively, polynomial and exponential in the number of variables. We conclude that neither classical nor quantum annealing can efficiently solve our reformulation of XORSAT, even though the original problem can be solved in polynomial time by Gaussian elimination.

2-dimensional systems↗
DOE OSTI2020

Mean-field entanglement transitions in random tree tensor networks

Entanglement phase transitions in quantum chaotic systems subject to projective measurements and in random tensor networks have emerged as a new class of critical points separating phases with different entanglement scaling. We propose a mean-field theory of such transitions by studying the entanglement properties of random tree tensor networks. As a function of bond dimension, we find a phase transition separating area-law from logarithmic scaling of the entanglement entropy. Using a mapping onto a replica statistical mechanics model defined on a Cayley tree and the cavity method, we analyze the scaling properties of such transitions. Our approach provides a tractable, mean-field-like example of an entanglement transition. Furthermore, we verify our predictions numerically by computing directly the entanglement of random tree tensor network states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Role of baryon number conservation in measurements of fluctuations

Here, I discuss the role and impact of net-baryon number conservation in measurements of net-proton fluctuations in heavy-ion collisions. I show that the magnitude of the fluctuations is entirely determined by the strength of two particle correlations. At LHC and top RHIC energy, this implies the fluctuations are proportional to the integral of the balance function (BF), B pp¯ of protons and antiprotons, while in the context of the RHIC beam energy scan (BES), one must also account for correlations of stopped protons. The integral of B pp¯ measured in a 4π detector depends on the relative cross sections of processes yielding pp¯ and those balancing the proton baryon number via the production of other antibaryons. The accepted integral of B pp¯ further depends on the shape and width of the BF relative to the width of the acceptance. The magnitude of the measured second-order cumulant of net-proton fluctuations thus has much less to do with QCD susceptibilities than with the creation/transport of baryons and antibaryons in heavy-ion collisions, and most particularly the impact of radial flow on the width of the BF. I thus advocate that net-proton fluctuations should be studied by means of differential BF measurements rather than with integral correlators. I also derive an expression of net-baryon fluctuations in terms of integrals of balance functions of identified baryon pairs and argue that measurements of such balance functions would enable a better understanding of the collision system expansion dynamics, the hadronization chemistry, and an experimental assessment of the strength of net-baryon fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Isoscaling and nuclear reaction dynamics

Isoscaling parameters α and β have been explored as a function of breakup angle in binary excited projectile-like fragment decays produced in collisions of 70 Zn+ 70 Zn and 64 Zn+ 64 Zn at 35 MeV per nucleon. In this analysis, focus was placed on isoscaling the second heaviest fragment with 4 ≤ Z L ≤ 8 emitted from the excited projectile-like fragment in events that contained a heavy fragment with Z H ≥ 12. The breakup orientation θ prox was defined as the angle between the heavy and light fragments’ center-of-mass velocity and the fragment pairs’ relative velocity. Breakups between 0° < θ prox ≤ 80° have been shown to be dominated by dynamical contributions, while break-up angles of θ prox > 100° are predominantly statistical. Historically, isoscaling has often been understood and applied in a statistical context, assuming that the fragments are produced after statistical equilibrium is achieved. Finally, studying isoscaling parameters as a function of θ prox reveals the sensitivity of α and β to the mechanism of fragment production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗