Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical physics”

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,585 records · Page 88

Model comparisons for two-temperature plasma equations of state

When a plasma is generated in the laboratory, energy is often deposited preferentially into either the electrons or the ions, giving rise to a quasiequilibrium state in which the two species, electrons and ions, are well described by two effective temperatures, T e and T i . Accurate hydrodynamic modeling of such a two-temperature plasma requires an equation of state that captures the relevant many-body physics without assuming a strict local thermodynamic equilibrium. Several models have been proposed within the literature, which extend conventional statistical approaches, each employing a different combination of assumptions for modifying the equilibrium equations. Here, in this work, we compare the predictions for several models, presenting derivations of the internal energy and pressure for each microscopic model within a unified framework so that the assumptions of each model may be more easily compared to one another. We find that for sufficiently weak coupling, all models agree with one another. However, as the coupling strength is increased, the disagreement between the models becomes more pronounced. Moreover, the relative sizes of the corrections predicted by each model depend on which species has the higher temperature, T e >T i vs T e < T i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancing Interpretability in Generative Modeling: Statistically Disentangled Latent Spaces Guided by Generative Factors in Scientific Datasets

This study addresses the challenge of statistically extracting generative factors from complex, high-dimensional datasets in unsupervised or semi-supervised settings. We investigate encoder-decoder-based generative models for nonlinear dimensionality reduction, focusing on disentangling low-dimensional latent variables corresponding to independent physical factors. Introducing Aux-VAE, a novel architecture within the classical Variational Autoencoder framework, we achieve disentanglement with minimal modifications to the standard VAE loss function by leveraging prior statistical knowledge through auxiliary variables. These variables guide the shaping of the latent space by aligning latent factors with learned auxiliary variables. We validate the efficacy of Aux-VAE through comparative assessments on multiple datasets, including astronomical simulations.

97 MATHEMATICS AND COMPUTING↗

Most theoretical approaches used in nuclear astrophysics to model the nucleosynthesis of heavy elements incorporate the so-called statistical model in order to describe the excitation and decay properties of atomic nuclei. One of the basic assumptions of this model is the validity of the Brink–Axel hypothesis and the related concept of so-called photon strength functions to describe γ-ray transition probabilities. We present a novel experimental approach that allows for the first time to experimentally determine the photon strength function simultaneously in two independent ways by a unique combination of quasi-monochromatic photon beams and a newly implemented γ–γ coincidence setup. This technique does not assume a priori the validity of the Brink–Axel hypothesis and sets a benchmark in terms of the detection sensitivity for measuring decay properties of photo-excited states below the neutron separation energy. The data for the spherical off-shell nucleus 128 Te were obtained for γ-ray beam-energy settings between 3 MeV and 9 MeV in steps of 130 keV for the lower beam energies and in steps of up to 280 keV for the highest beam settings. We present a quantitative analysis on the consistency of the derived photon strength function with the Brink–Axel hypothesis. The data clearly demonstrate a discrepancy of up to a factor of two between the photon strength functions extracted from the photoabsorption and photon emission process, respectively. In addition, we observe that the photon strength functions are not independent of the excitation energy, as usually assumed. Thus, we conclude, that the Brink–Axel hypothesis is not strictly fulfilled in the excitation-energy region below the neutron separation threshold (S n = 8.78 MeV) for the studied case of 128 Te.

Not Available

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Persona of Transition Metal Ions in Solids: A Statistical Learning on Local Structures of Transition Metal Oxides

Abstract The local structure of a transition metal (TM) ion is a function of cation elements and valence states. More than that, in this work, by employing a trove of first‐principles data of TM oxides, the local structures of TM cations are statistically analyzed to extract detailed information about cation site preference, bond length, site structural distortion, and cation magnetization. It is found that cation radius alone poorly describes the local structure of a transition metal oxide, while the statistics of coordination number as well as the TMO bond length distribution, especially that of the 3d TMs, can provide comprehensive knowledge for understanding the behavior of TM elements. Based on these statistics, the interplay of site distortion due to the Jahn–Teller effect, cation site similarity, and a new set of ionic radii are all obtained to chart the “persona” of transition metal ions in solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced Method Optimization for Sampling and Analysis Instrumentation

This work presents a generalized approach for analytical method optimization that branches the gap between techniques historically employed and accurate modern optimization techniques suitable for various applications. The novelty of the described strategy is the utilization of multivariate, multiobjective optimization with Karush-Kuhn-Tucker conditions to bound the optimization space to solutions within the physical limitations of instrumentation. Briefly, the basic steps outlined in this paper are to (1) determine the objective(s) that should be maximized or minimized based on the goals of the analytical application, (2) conduct a screening experiment, (3) perform ANOVA to determine the parameters which have a statistically significant effect on the objective, (4) conduct an experiment (e.g., Box-Behnken design) to collect data for fitting the objective equation, and (5) determine the physical constraints of the parameters and solve the Lagrangian to determine the optimal method parameters. A broad approach to optimization target selection allows for robust method tuning to develop improved data sets amenable for chemometrics and machine learning algorithm development. Gas chromatography-mass spectrometry was selected as a use case due to its broad use across scientific fields and time-consuming method development involving numerous parameters. In conclusion, this strategy can reduce the cost of research, improve data quality, and enable the rapid development of new analytical technique.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗

Subtropical Gyre Variability as Seen from Satellites

A satellite multi-sensor approach is used to analyse the biological response of open ocean regions of the subtropical gyres to changes in physical forcing. Thirteen years (1998-2010) of SeaWiFS chlorophyll a (Chl-a), combined with concurrent satellite records of sea-surface temperature (SST) and sea level height, were analysed to investigate the seasonal and interannual variability of Chl-a concentration within these immense so-called ocean deserts. The seasonal variability of Chl-a within the gyres is driven mostly by the warming/cooling of surface waters. Summer warming promotes shallower mixed layers and lower Chl-a due to a reduction of vertical mixing and consequently a decrease in nutrient supply. The opposite happens during the winter cooling period. Therefore, long-term trends in SST have the potential to cause an impact on the interannual variability of Chl-a. Our analyses show that, during the 13 whole years of SeaWiFS data record, the North Pacific, Indian Ocean, and North Atlantic gyres experienced a decrease in Chl-a of 9%, 12%, and 11%, respectively, with corresponding SST increases of 0.27 C, 0.42 C, and 0.32 C. The South Pacific and South Atlantic gyres also showed warming trends but with weak positive trends in Chl-a that are not statistically significant. We hypothesize that the warming of surface waters in these two gyres are counterbalanced by other interacting physical and biological driving mechanisms, as indicated in previous studies.

Signorini, Sergio R.↗

GenASiS Basics: Object-oriented utilitarian functionality for large-scale physics simulations (Version 3)

GenASiS Basics provides Fortran 2003 classes furnishing extensible object-oriented utilitarian functionality for large-scale physics simulations on distributed memory supercomputers. This functionality includes physical units and constants; display to the screen or standard output device; message passing; I/O to disk; and runtime parameter management and usage statistics. This revision – Version 3 of Basics – includes a significant name change, some minor additions to functionality, and a major addition to functionality: infrastructure facilitating the offloading of computational kernels to devices such as GPUs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

q-coherent states associated with the noncommutative complex plane C{sub q{sup 2}} for the Biedenharn–Macfarlane q-oscillator

The q-coherent states associated with the noncommutative complex plane C{sub q{sup 2}} for the Biedenharn–Macfarlane q-oscillator are constructed by using a q-generalization of the Glauber’s unitary displacement operator. The q-deformed counterpart of the eigenvalue equation of undeformed harmonic oscillator coherent states is obtained. Some properties of the q-coherent states are considered. It is shown that super-Poissonian statistics and saturation of the uncertainty relation occur simultaneously with the nonclassical behaviors including antibunching effects with nonnegative values less than unity and the negative values, for the second-order intensity q-correlation function in the case of q>1. Furthermore, the nonclassical sub-Poissonian statistics is demonstrated simultaneously with bunching and antibunching effects such that in some regions of the coherence parameter the uncertainty inequality is satisfied, and violated elsewhere. We show that another nonclassical effect, the so-called strong squeezing effect, is exhibited by both quadrature variances in the latter issue.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calibration strategy of the PROSPECT-II detector with external and intrinsic sources

This paper presents an energy calibration scheme for an upgraded reactor antineutrino detector for the Precision Reactor Oscillation and Spectrum Experiment (PROSPECT). The PROSPECT collaboration is preparing an upgraded detector, PROSPECT-II (P-II), to advance capabilities for the investigation of fundamental neutrino physics, fission processes and associated reactor neutrino flux, and nuclear security applications. P-II will expand the statistical power of the original PROSPECT (P-I) dataset by at least an order of magnitude. The new design builds upon previous P-I design and focuses on improving the detector robustness and long-term stability to enable multi-year operation at one or more sites. The new design optimizes the fiducial volume by elimination of dead space previously occupied by internal calibration channels, which in turn necessitates the external deployment. In this paper, we describe a calibration strategy for P-II. Here, the expected performance of externally deployed calibration sources is evaluated using P-I data and a well-benchmarked simulation package by varying detector segmentation configurations in the analysis. The proposed external calibration scheme delivers a compatible energy scale model and achieves comparable performance with the inclusion of an additional AmBe neutron source, in comparison to the previous internal arrangement. Most importantly, the estimated uncertainty contribution from the external energy scale calibration model meets the precision requirements of the P-II experiment.

47 OTHER INSTRUMENTATION↗

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neutron and Gamma Correlations using CGM in MCNP 6.2.0

The general-purpose Monte Carlo radiation-transport code MCNP6{sup TM} came from the integration of MCNP5{sup TM} and MCNPX{sup TM}. The most recent release of MCNP, version 6.2.0, is scheduled for release in early 2017, which includes a variety of bug fixes and new features to the code. MCNP traditionally models nuclear reactions using Monte Carlo sampling techniques on measured and modeled cross-sectional data contained within ACE (A Compact Evaluated Nuclear Data File) libraries. In addition to reaction cross-sections, ACE libraries contain information relating to the production of secondary particles, such as secondary neutrons and photons. For example, an incident neutron on a nucleus can produce a nuclear reaction in the form (n, M{sub n}n'), where n is the incident neutron and M{sub n} is the number of secondary neutrons n' produced. Using this notation, different reactions can be described based on their neutron multiplicity, such as neutron capture (M{sub n} = 0), elastic and inelastic scattering (M{sub n} = 1), and reactions with multiple secondary neutrons (M{sub n} > 1). Accounting for secondary γ-rays emitted from the residual nuclei, this reaction can further be generalized as n, Mnn'M to describe a reaction producing M{sub n} neutrons and M -rays. These values M{sub n} and M are referred to as the neutron and multiplicities, respectively. ACE libraries contain the cross-sectional data for each reaction, and their corresponding statistically averaged multiplicities for neutrons and photons, M-bar{sub n} and M-bar{sub γ}. However, there are several limitations to reaction sampling using ACE libraries. First, the data contained within ACE libraries is limited by the accuracy of the model and/or experiment, which out of necessity can greatly simplify the true physics behind a reaction. For example, the value M{sub γ} in a (n, M{sub n}n'M{sub γ}γ) reaction is determined by the statistically averaged multiplicity M-bar{sub γ} instead of a distribution. Consequently, MCNP, unable to produce a continuous distribution, will utilizes a binary sampling around M-bar{sub γ}. While this sampling method provides a statistically average M{sub γ} over a large simulation, it is an inaccurate representation of an individual reaction and its corresponding secondary particles. Secondary neutrons are affected in a similar manner, such that different kinds of secondary neutrons can be sampled from the same interaction (eg., an inelastic neutron with an [n, 2n] neutron). Lastly, the random sampling technique for secondary particles remains completely independent of the sampled reaction, eliminating the ability to correlate secondary particles for a given reaction. To supplement these data libraries, additional codes and models describing certain nuclear reactions are employed. These codes can fill in the information missing from ACE, provide a better model to the true physics distribution, and can provide true correlated secondary particles, which has many useful applications in research and industry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Legendre submanifolds in contact manifolds as attractors and geometric nonequilibrium thermodynamics

It has been proposed that equilibrium thermodynamics is described on Legendre submanifolds in contact geometry. It is shown in this paper that Legendre submanifolds embedded in a contact manifold can be expressed as attractors in phase space for a certain class of contact Hamiltonian vector fields. By giving a physical interpretation that points outside the Legendre submanifold can represent nonequilibrium states of thermodynamic variables, in addition to that points of a given Legendre submanifold can represent equilibrium states of the variables, this class of contact Hamiltonian vector fields is physically interpreted as a class of relaxation processes, in which thermodynamic variables achieve an equilibrium state from a nonequilibrium state through a time evolution, a typical nonequilibrium phenomenon. Geometric properties of such vector fields on contact manifolds are characterized after introducing a metric tensor field on a contact manifold. It is also shown that a contact manifold and a strictly convex function induce a lower dimensional dually flat space used in information geometry where a geometrization of equilibrium statistical mechanics is constructed. Legendre duality on contact manifolds is explicitly stated throughout.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistics of base polytopes in F-theory

We propose a new statistical ensemble of toric bases for elliptic Calabi-Yaus used in F-theory models, by focusing on only the convex hull of the base, i.e., the base polytope. This physically motivated coarse-graining greatly simplifies the combinatorial complexity of the part of the 4d F-theory landscape with toric bases. We develop a Monte Carlo approach that randomly samples the base polytopes within fixed boxes, with proper statistical weights. We first apply the algorithm to the set of 2d base polytopes, generating an enlarged set of toric 2d bases that include certain types of codimension-two (4,6) points, and we validate our approach against exact numbers. We then explore the set of 3d base polytopes which fit in a set of “maximal” 3d boxes, and estimate the total number of inequivalent 3d base polytopes to be 10 85 –10 90 . We provide statistical data such as the distribution of non-Higgsable gauge groups on these bases. Amusingly, a similar method can also be applied to generate reflexive polytopes in various dimensions. In both the reflexive and base polytope cases, the number of relevant polytopes obeys a Gaussian distribution as a function of the number of vertices, which can be understood in terms of other results on random polytopes in the math literature.

Differential and algebraic geometry↗

How Does Pressure Fluctuate in Equilibrium?

We study fluctuations of pressure in equilibrium for classical particle systems. In equilibrium statistical mechanics, pressure for a microscopic state is defined by the derivative of a thermodynamic function or, more mechanically, through the momentum current. We show that although the two expectation values converge to the same equilibrium value in the thermodynamic limit, the variance of the mechanical pressure is in general greater than that of the pressure defined through the thermodynamic relation. We also present a condition for experimentally detecting the difference between them in an idealized measurement of momentum transfer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Gaussian CL s method for searches of new physics

Here we describe a method based on the CL s approach to present results in searches of new physics, under the condition that the relevant parameter space is continuous. Our method relies on a class of test statistics developed for non-nested hypotheses testing problems, denoted by ΔT, which has a Gaussian approximation to its parent distribution when the sample size is large. This leads to a simple procedure of forming exclusion sets for the parameters of interest, which we call the Gaussian CL s method. Our work provides a self-contained mathematical proof for the Gaussian CL s method, that explicitly outlines the required conditions. These conditions are milder than that required by the Wilks' theorem to set confidence intervals (CIs). We illustrate the Gaussian CL s method in an example of searching for a sterile neutrino, where the CL s approach was rarely used before. We also compare data analysis results produced by the Gaussian CL s method and various CI methods to showcase their differences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Estimating Vegetation Parameters from Interferometric and Polarimetric Radar Using Physical Scattering Models

Radar data from vegetated land surfaces depend on many structural and compositional parameters describing the terrain. Because early, noninterferometric radar systems usually constituted an insufficient observation set from which to estimate parameters of the terrain, statistical regression techniques were used which incorporated some level of apriori knowledge or field measurements. With the advent of radar interferometry and polarimetric interferometry, potentially at multiple baselines, the observation set is now approaching that required to quantitatively estimate the parameters describing a vegetated land surface. Quantitative estimation entails formulating a physical scattering model relating the radar observations to the vegetation and surface parameters on which they depend. This paper describes the physics of candidate scattering models, and shows how the models determine the estimable parameter set. It also indicates the measurement accuracy of parameters such as vegetation height, height-to-base-of-live-crown, and surface topography with multibaseline polarimetric interferometry.

Treuhaft, Robert N.↗