Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,171 records · Page 65

Positivity from Cosmological Correlators

Abstract Effective field theories in flat space and in anti-de Sitter space are constrained by causality and unitarity, often in the form of positivity bounds. Similar bounds have been harder to demonstrate in cosmological backgrounds, where the roles of unitarity and causality are more obscure. Fortunately, the expansion of the universe ensures that late-time cosmological correlators are effectively classical and the role of unitarity is played by classical statistical inequalities. For multi-field inflation, the resulting positivity constraints have long been known in terms of the Suyama-Yamaguchi inequality. In this paper, we demonstrate that similar statistical bounds imply nontrivial constraints for massive fields in the early universe. We show that any real anomalous dimensions for principal series fields in de Sitter space must be positive. We also derive a limit on the amplitude of particular oscillatory signals from inflation, including those arising in cosmological collider physics. Finally, we demonstrate that these constraints manifest themselves directly in the two-point statistics of matter and galaxies that will be measured in upcoming surveys.

Physics↗

Where Lies the Quantumness Behind Detecting Electromagnetic Waves for Frequencies from Infrared and Up?

The synthesis of Newtonian concept of corpuscles during emission with Huygens’ concept of secondary-wavelets during propagation implies that all EM radiations from quantized atoms and molecules are released as discrete amount of energies. However, they propagate out as time-finite Maxwellian light pulses. Huygens also underscored that his secondary wavelets keep propagating as independent pulses in the absence of any interacting medium, or until intercepted by an interface with a medium or a detector. Then we use the Superposition Principle and the coherence theory to derive Einstein’s photoelectric equation by summing innumerable random time-finite pulses. This process driven approach should yield the characteristic statistical variations of photoelectron current pulses, as generated by photodetectors for different kinds of light sources. Lamb & Scully originally proposed this semiclassical approach without assuming that light actually consists of time finite pulses. The quantumness remains confined within the excitation and de-excitations processes in the material particles.

Electromagnetic Waves↗

Optimized setup for two-dimensional convection experiments in thin liquid films

We present a novel experimental setup to investigate two-dimensional thermal convection in a freestanding thin liquid film. Such films can be produced in a controlled way on the scale of 5–1000 nm. Our primary goal is to investigate convection patterns and the statistics of reversals in Rayleigh-Bénard convection with varying aspect ratio. Additionally, questions regarding the physics of liquid films under controlled conditions can be investigated, like surface forces, or stability under varying thermodynamical parameters. The film is suspended in a frame which can be adjusted in height and width to span an aspect ratio range of Γ = 0.16–10. The top and bottom frame elements can be set to specific temperature within T = 15 °C to 55 °C. A thickness to area ratio of approximately 10{sup 8} enables only two-dimensional fluid motion in the time scales relevant for turbulent motion. The chemical composition of the film is well-defined and optimized for film stability and reproducibility and in combination with carefully controlled ambient parameters allows the comparison to existing experimental and numerical data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Statistical model of exotic rotational correlations in emergent space-time

A statistical model is formulated to compute exotic rotational correlations that arise as inertial frames and causal structure emerge on large scales from entangled Planck scale quantum systems. Noncommutative quantum dynamics are represented by random transverse displacements that respect causal symmetry. Entanglement is represented by covariance of these displacements in Planck scale intervals defined by future null cones of events on an observer's world line. Light that propagates in a nonradial direction inherits a projected component of the exotic rotational correlation that accumulates as a random walk in phase. A calculation of the projection and accumulation leads to exact predictions for statistical properties of exotic Planck scale correlations in an interferometer of any configuration. The cross-covariance for two nearly co-located interferometers is shown to depart only slightly from the autocovariance. Specific examples are computed for configurations that approximate realistic experiments, and show that the model can be rigorously tested.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗

Fractional excitations in foliated fracton phases

Fractional excitations in fracton models exhibit novel features not present in conventional topological phases: their mobility is constrained, there are an infinitude of types, and they bear an exotic sense of ‘braiding’. Hence, they require a new framework for proper characterization. Based on our definition of foliated fracton phases in which equivalence between models includes the possibility of adding layers of gapped 2D states, we propose to characterize fractional excitations in these phases up to the addition of quasiparticles with 2D mobility. That is, two quasiparticles differing by a set of quasiparticles that move along 2D planes are considered to be equivalent; likewise, ‘braiding’ statistics are measured in a way that is insensitive to the attachment of 2D quasiparticles. The fractional excitation types and statistics defined in this way provide a universal characterization of the underlying foliated fracton order which can subsequently be used to establish phase relations. We demonstrate as an example the equivalence between the X-cube model and the semionic X-cube model both in terms of fractional excitations and through an exact mapping.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Risk-Adaptive Experimental Design for High-Consequence Systems: LDRD Final Report

Constructing accurate statistical models of critical system responses typically requires an enormous amount of data from physical experiments or numerical simulations. Unfortunately, data generation is often expensive and time consuming. To streamline the data generation process, optimal experimental design determines the 'best' allocation of experiments with respect to a criterion that measures the ability to estimate some important aspect of an assumed statistical model. While optimal design has a vast literature, few researchers have developed design paradigms targeting tail statistics, such as quantiles. In this project, we tailored and extended traditional design paradigms to target distribution tails. Our approach included (i) the development of new optimality criteria to shape the distribution of prediction variances, (ii) the development of novel risk-adapted surrogate models that provably overestimate certain statistics including the probability of exceeding a threshold, and (iii) the asymptotic analysis of regression approaches that target tail statistics such as superquantile regression. To accompany our theoretical contributions, we released implementations of our methods for surrogate modeling and design of experiments in two complementary open source software packages, the ROL/OED Toolkit and PyApprox.

97 MATHEMATICS AND COMPUTING↗

Gamma-ray emission in dissipative pulsar magnetospheres: from theory to Fermi observations

We compute the patterns of γ-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed γ-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model γ-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and γ-ray emission on the γ-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the γ-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars↗

Data-driven probability concentration and sampling on manifold

A new methodology is proposed for generating realizations of a random vector with values in a finite-dimensional Euclidean space that are statistically consistent with a dataset of observations of this vector. The probability distribution of this random vector, while a priori not known, is presumed to be concentrated on an unknown subset of the Euclidean space. A random matrix is introduced whose columns are independent copies of the random vector and for which the number of columns is the number of data points in the dataset. The approach is based on the use of (i) the multidimensional kernel-density estimation method for estimating the probability distribution of the random matrix, (ii) a MCMC method for generating realizations for the random matrix, (iii) the diffusion-maps approach for discovering and characterizing the geometry and the structure of the dataset, and (iv) a reduced-order representation of the random matrix, which is constructed using the diffusion-maps vectors associated with the first eigenvalues of the transition matrix relative to the given dataset. The convergence aspects of the proposed methodology are analyzed and a numerical validation is explored through three applications of increasing complexity. The proposed method is found to be robust to noise levels and data complexity as well as to the intrinsic dimension of data and the size of experimental datasets. Both the methodology and the underlying mathematical framework presented in this paper contribute new capabilities and perspectives at the interface of uncertainty quantification, statistical data analysis, stochastic modeling and associated statistical inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fast Quantum Algorithm for Predicting Descriptive Statistics of Stochastic Processes

Stochastic processes are used as a modeling tool in several sub-fields of physics, biology, and finance. Analytic understanding of the long term behavior of such processes is only tractable for very simple types of stochastic processes such as Markovian processes. However, in real world applications more complex stochastic processes often arise. In physics, the complicating factor might be nonlinearities; in biology it might be memory effects; and in finance is might be the non-random intentional behavior of participants in a market. In the absence of analytic insight, one is forced to understand these more complex stochastic processes via numerical simulation techniques. In this paper we present a quantum algorithm for performing such simulations. In particular, we show how a quantum algorithm can predict arbitrary descriptive statistics (moments) of N-step stochastic processes in just O(square root of N) time. That is, the quantum complexity is the square root of the classical complexity for performing such simulations. This is a significant speedup in comparison to the current state of the art.

Williams Colin P.↗

Closing the Gap Between Modeling and Experiments in the Self-assembly of Biomolecules at Interfaces and in Solution

Molecular self-assembly is a powerful tool in materials design, wherein non-covalent interactions like electrostatic, hydrophobic, hydrogen bonding, and van der Waals can be exploited to produce supramolecular nanostructures that are functional and highly tunable. Biomolecules are attractive building blocks, as they are biocompatible, biodegradable and adopt a wide array of higher order structures. Moreover, naturally occurring protein systems display a manifold of structures and interactions that can be replicated in synthetic biomolecules. In this perspective, we highlight advances in multiscale simulation techniques across broad spatiotemporal scales that can aid in characterizing self-assembly of hybrid and hierarchical bionanomaterial systems, with an emphasis on physics-based simulation approaches currently employed to study biomolecules at mineral interfaces. The power of these approaches is highlighted across a few recent areas where molecular simulations have advanced our understanding of self-assembly spanning peptides to protein self-assembly. Looking forward, we discuss how in the near future emerging methods in statistical and machine learning will advance this research field in all areas from expanding the capabilities of physics-based simulation methods to enabling new analyses of high throughput experiments. These advances will pave the way for understanding the molecular recognition patterns in systems that are dictated by self-assembly - biomineralizing peptides, hierarchical peptoids, and large protein assemblies, and will aid in the development of a new synthesis science for achieving precise molecular control in materials design

Sampath, Janani↗

Two-Dimensional Triplet Magnetoexcitons and a Magnetofermionic Condensate in the GaAs/AlGaAs Heterostructures

A fundamentally new collective state, namely, the magnetofermionic condensate, is discovered during photoexcitation of a sufficiently dense gas of long-lived triplet cyclotron magnetoexcitons in a twodimensional Hall insulator with a high electron mobility, a filling factor of ν = 2, and temperatures of T < 1 K. The condensed phase coherently interacts with an external electromagnetic field, exhibits superradiant properties in the recombination of correlated condensate electrons with heavy holes in the valence band, and spreads nondissipatively in the layer of a two-dimensional electron gas to macroscopical large distances, transferring an integer spin. The observed effects are explained in terms of a coherent condensate in a nonequilibrium system of two-dimensional fermions with a fully quantized energy spectrum, in which a degenerate ensemble of long-lived triplet magnetoexcitons obeying the Bose statistics is present.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nuclear EMC effect in non-extensive statistical model

In the present work, we attempt to describe the nuclear EMC effect by using the proton structure functions obtained from the non-extensive statistical quark model. We record that such model has three fundamental variables, the temperature T, the radius, and the Tsallis parameter q. By combining different small changes, a good agreement with the experimental data may be obtained. Another interesting point of the model is to allow phenomenological interpretation, for instance, with q constant and changing the radius and the temperature or changing the radius and q and keeping the temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A regression-based approach to estimating retrofit savings using the Building Performance Database

Retrofitting building systems is known to provide cost-effective energy savings. This article addresses how the Building Performance Database is used to help identify potential savings. Currently, prioritizing retrofits and computing their expected energy savings and cost/benefits can be a complicated, costly, and an uncertain effort. Prioritizing retrofits for a portfolio of buildings can be even more difficult if the owner must determine different investment strategies for each of the buildings. Meanwhile, we are seeing greater availability of data on building energy use, characteristics, and equipment. These data provide opportunities for the development of algorithms that link building characteristics and retrofits empirically. In this paper we explore the potential of using such data for predicting the expected energy savings from equipment retrofits for a large number of buildings. We show that building data with statistical algorithms can provide savings estimates when detailed energy audits and physics-based simulations are not cost- or time-feasible. We develop a multivariate linear regression model with numerical predictors (e.g., operating hours, occupant density) and categorical indicator variables (e.g., climate zone, heating system type) to predict energy use intensity. The model quantifies the contribution of building characteristics and systems to energy use, and we use it to infer the expected savings when modifying particular equipment. We verify the model using residual analysis and cross-validation. We demonstrate the retrofit analysis by providing a probabilistic estimate of energy savings for several hypothetical building retrofits. We discuss the ways understanding the risk associated with retrofit investments can inform decision making. The contributions of this work are the development of a statistical model for estimating energy savings, its application to a large empirical building dataset, and a discussion of its use in informing building retrofit decisions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Status and Prospects of SBND

The Short-Baseline Near Detector (SBND) is one of the neutrino detectors positioned along the the Booster Neutrino Beam (BNB) at Fermilab. It is the near detector in the Short-Baseline Neutrino (SBN) Program, as well as having a wide-ranging independent physics program. The detector began taking neutrino data in the summer of 2024 and has already collected over 1 million neutrino events. Thanks to its unique combination of measurement resolution and statistics, SBND will soon carry out a rich program of neutrino interaction measurements and novel searches for physics beyond the Standard Model (BSM). As the SBN near detector, it will enable the full potential of the SBN sterile neutrino program by performing a precise characterisation of the unoscillated event rate and constraining BNB flux and neutrino-argon cross-section systematic uncertainties. In this talk, the physics reach, current status, and future prospects of SBND are discussed.

Paton, Josephine [Fermilab]↗

An efficient optimization of well placement and control for a geothermal prospect under geological uncertainty

This work applies an efficient optimization technique based on a multivariate adaptive regression spline (MARS) technique to determine the optimal design and engineering of a potential geothermal production operation at a prospect near Superstition Mountain in Southern California, USA. The faster MARS-based statistical model is used as a surrogate for higher-fidelity physical models within the intensive optimization process. Its use allows for the exploration of the impacts of specific engineering design parameters in the context of geologic uncertainty as a means to both understand and maximize profitability of the production operation. The MARS model is initially developed from a training dataset generated by a finite set of computationally complex hydrothermal models applied to the prospect. Its application reveals that the optimal engineering design variables can differ considerably assuming different choices of hydrothermal flow properties, which, in turn, indicates the importance of reducing the uncertainty of key geologic properties. The major uncertainty sources in the natural-system are identified and ranked first by an efficient MARS-enabled total order sensitivity quantification, which is then used to assist evaluating the effect of geological uncertainties on optimized results. At the Southern California prospect, this parameter sensitivity analysis suggests that groundwater circulation through high permeable structures, rather than heat conduction through impermeable granite, is the primary heat transfer method during geothermal extraction. Reservoir histories simulated using optimal parameters with different constraints are analyzed and compared to investigate the longevity and maximum profit of the geothermal resources. The comparison shows that the longevity and profit are very likely to be overestimated by optimizations without appropriate constraints on natural conditions. In addition to geothermal energy production, this optimization approach can also be used to manage other geologic resource operations, such as hydrocarbon production or CO 2 sequestration, under uncertain reservoir conditions.

15 GEOTHERMAL ENERGY↗

Inadequacy of current approaches for characterizing membrane transport properties at high salinities

Cost optimal design of osmotic membrane processes requires an accurate estimate of membrane transport parameters across their full operational range. However, standard approaches for estimating these parameters rely on empirical methods, the accuracy of which remains unquantified as a function of temperature, salinity, and measurement error. Herein, we present a systematic accuracy analysis of previously developed methods for estimation of membrane transport properties in reverse osmosis, high-pressure reverse osmosis, forward osmosis, pressure retarded reverse osmosis, and osmotically assisted reverse osmosis. We use a Monte Carlo approach to sample the full range of feasible membrane water permeabilities, salt permeabilities, structural parameters, and operating conditions for these processes. These material and process parameters are then incorporated into a physical transport model for each process. Our analysis shows that the statistical uncertainty of current empirical methods for estimating membrane parameters increases by 5 times from low-salinity to high-salinity conditions. Here, the result of this work demonstrates that empirical methods are inadequate for precisely estimating membrane transport properties at high salinity and highlight a critical need for the development of statistically validated higher accuracy methods.

42 ENGINEERING↗