Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical field theory”

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 397 records · Page 22

Measurement of the muon anomalous precession frequency in runs 4, 5, and 6 of the muon ${g}-2$ Experiment at Fermilab

The Fermilab E989 Muon $g-2$ experiment measures the muon's anomalous magnetic moment to a precision of 127 parts per billion, as reported in June 2025. The value is proportional to the difference between the muon's cyclotron frequency and the spin precession frequency in the presence of a uniform magnetic field, for muons contained within the $g-2$ storage ring. Spin precession frequency is extracted from the time distribution of the muon's decay positrons recorded by 24 electromagnetic calorimeters positioned around the inner circumference of the storage ring. The anomalous precession frequency is one of the primary experimental inputs necessary to estimate the anomalous magnetic moment, the other being the measurement of the magnetic field. This dissertation details the anomalous precession frequency extraction, including reconstruction, time-distribution fitting, and treatment of systematic uncertainties for the final three data-collection runs: Run-4, Run-5, and Run-6. This data represents a fourfold increase in statistics over the previous analysis release, halving the statistical uncertainty. The residual slow term from previous analyses is now well understood and documented in a systematic treatment. As of the writing of this dissertation, the theoretical prediction for the SM estimate of the muon's anomalous magnetic moment is under debate, with two competing prediction methods, so a definitive comparison with theory is not available. The results submitted for experimental release use the kernel-ratio asymmetry method, contributing 115 parts per billion to the statistical uncertainty and 34 parts per billion to the systematic uncertainty. When combined with the previous analyses in earlier data runs, this thereby improves the measurement beyond the experimental goal and sets the world's most precise measurement of the muon's anomalous magnetic moment.

Israel, Scott Nathan [Boston U.]↗

Measurement of the muon anomalous precession frequency in Runs 4, 5, and 6 of the Muon g-2 experiment at Fermilab

The Fermilab E989 Muon g − 2 experiment measures the muon’s anomalous magnetic moment to a precision of 127 parts per billion, as reported in June 2025. The value is proportional to the difference between the muon’s cyclotron frequency and the spin precession frequency in the presence of a uniform magnetic field, for muons contained within the g − 2 storage ring. Spin precession frequency is extracted from the time distribution of the muon’s decay positrons recorded by 24 electromagnetic calorimeters positioned around the inner circumference of the storage ring. The anomalous precession frequency is one of the primary experimental inputs necessary to estimate the anomalous magnetic moment, the other being the measurement of the magnetic field. This dissertation details the anomalous precession frequency extraction, including reconstruction, time-distribution fitting, and treatment of systematic uncertainties for the final three data-collection runs: Run-4, Run-5, and Run-6. This data represents a fourfold increase in statistics over the previous analysis release, halving the statistical uncertainty. The residual slow term from previous analyses is now well understood and documented in a systematic treatment. As of the writing of this dissertation, the theoretical prediction for the SM estimate of the muon’s anomalous magnetic moment is under debate, with two competing prediction methods, so a definitive comparison with theory is not available. The results submitted for experimental release use the kernel-ratio asymmetry method, contributing 115 parts per billion to the statistical uncertainty and 34 parts per billion to the systematic uncertainty. When combined with the previous analyses in earlier data runs, this thereby improves the measurement beyond the experimental goal and sets the world’s most precise measurement of the muon’s anomalous magnetic moment.

Israel, Scott Nathan [Boston U.]↗

Anyonic Membranes and Pontryagin Statistics

Anyons, unique to two spatial dimensions, underlie extraordinary phenomena such as the fractional quantum Hall effect, but their generalization to higher dimensions has remained elusive. The topology of Eilenberg-MacLane spaces constrains the loop statistics to be only bosonic or fermionic in any dimension. In this work, we introduce the novel anyonic statistics for membrane excitations in four dimensions. Analogous to the $\mathbb{Z}_N$-particle exhibiting $\mathbb{Z}_{N\times \gcd(2,N)}$ anyonic statistics in two dimensions, we show that the $\mathbb{Z}_N$-membrane possesses $\mathbb{Z}_{N\times \gcd(3,N)}$ anyonic statistics in four dimensions. Given unitary volume operators that create membrane excitations on the boundary, we propose an explicit 56-step unitary sequence that detects the membrane statistics. We further analyze the boundary theory of $(5{+}1)$D 1-form $\mathbb{Z}_N$ symmetry-protected topological phases and demonstrate that their domain walls realize all possible anyonic membrane statistics. We then show that the $\mathbb{Z}_3$ subgroup persists in all higher dimensions. In addition to the standard fermionic $\mathbb{Z}_2$ membrane statistics arising from Stiefel-Whitney classes, membranes also exhibit $\mathbb{Z}_3$ statistics associated with Pontryagin classes. We explicitly verify that the 56-step process detects the nontrivial $\mathbb{Z}_3$ statistics in 5, 6, and 7 spatial dimensions. Furthermore, in 7 and higher dimensions, the statistics of membrane excitations stabilize to $\mathbb{Z}_{2} \times \mathbb{Z}_{3}$, with the $\mathbb{Z}_3$ sector consistently captured by this process.

Abstract algebra↗

The Response of Ribbed Panels to Reverberant Acoustic Fields

A statistical method for estimating the response of ribbed panels to acoustic excitation is discussed. It is shown that the acceleration spectrum of the vibrational field is related to the pressure spectrum by a coupling factor which is a simple function of the radiation and mechanical resistance of the structure. The radiation resistance of a ribbed panel is studied as a function of frequency. The analysis predicts that ribbing increases the radiation resistance of the panel and hence its coupling to the acoustic field. The effect of various panel-rib boundary conditions is also considered. The results of experiments which were conducted to test the theory are reported. The agreement between theory and experiments is shown to be satisfactory.

BOUNDARY↗

Timing of solar cycles by rigid internal rotations

The so-called 11-year cycle of solar activity is really more complex and contains many periods of greatly different lengths. Periods as long as 178 years and as short as 3.1 years are predicted by a theory based on beats between rigidly rotating, inertially oscillating g-modes inside the sun. Most of the beat periods are then confirmed to about 1 percent accuracy in sunspot observations. Since the agreement is of high statistical significance, one can conclude that approximate alignment of major solar oscillation modes contributes to high solar activity. The theory receives further support when tested against an independent class of observations - the large-scale magnetic sector structure. Predicted rotation rates of at least four solar oscillation modes are detected in the sector data with discrepancies all less than 0.3 percent. As a by-product of these successful fits to observation, the mean rotation of the entire solar mass becomes known. Its rotation frequency is 4.49 by 10 to the -7th power Hz, which is a sidereal period of 25.8 days. Magnetic fields have played no role in calculating the length of any of these solar cycles.

Wolff, C. L.↗

An update on the correlation between the cosmic radiation intensity and the geomagnetic AA index

A statistical study between the cosmic ray intensity, as observed by a neutron monitor, and of the geomagnetic aa index, as representative of perturbations in the plasma and interplanetary magnetic field in the heliosphere, has been updated to specifically exclude time periods around the reversal of the solar magnetic field. The results of this study show a strong negative correlation for the period 1960 through 1968 with a correlation coefficient of approximately -0.86. However, there is essentially no correlation between the cosmic ray intensity and the aa index for the period 1972-1979 (i.e. correlation coefficient less than 0.16). These results would appear to support the theory of preferential particle propagation into the heliosphere vis the ecliptic during the period 1960-1968 and via the solar polar regions during 1972-1979.

Shea, M. A.↗

Backbone Stitching in Bottlebrush Copolymer Mesodomains and the Impact of Side Chain Crystallization

We synthesized bottlebrush statistical copolymers (BSCPs) having poly­(ethylene oxide) (PEO) and poly­(dimethylsiloxane) (PDMS) side chains attached to a polynorbornene backbone. Small-angle X-ray scattering analysis showed that for densely grafted BSCPs, the scattering length density gradually transitions between the PEO and PDMS domains. For loosely grafted BSCPs, the polymer backbone formed a distinct mesodomain, with a lower electron and mass density than both the PEO and PDMS domains. The bottlebrush backbone essentially “stitches” the PEO and PDMS side chains, looping back and forth from the PEO to PDMS domains with the backbone segments oriented normal to the domain interfaces. Self-consistent field theory (SCFT) calculations validated the stitching of the backbone driven by the microphase separation of PEO and PDMS, along with a strong segmental order of the side chains in the melt. The reduced birefringence upon PEO crystallization suggests the disruption of the strong segmental order by the crystallization. Both the static intrinsic and the form birefringences of the BSCPs decreased upon PEO crystallization. Solid-state NMR confirmed the rigidity of PEO crystallites and the bottlebrush backbone. Self-assembly of BSCPs containing polyhedral oligomeric silsesquioxane (POSS) pendent groups was also evaluated by X-ray scattering, showing the formation of lamellar microdomains that inhibited POSS crystallization.

Hu, Mingqiu↗

Short-term climate variability and atmospheric teleconnections from satellite-observed outgoing longwave radiation. I Simultaneous relationships. II - Lagged correlations

Attention is given to the low-frequency variability of outgoing longwave radiation (OLR) fluctuations, their possible correlations over different parts of the globe, and their relationships with teleconnections obtained from other meteorological parameters, for example, geopotential and temperature fields. Simultaneous relationships with respect to the Southern Oscillation (Namais, 1978; Barnett, 1981) signal and the reference OLR fluctuation over the equatorial central Pacific are investigated. Emphasis is placed on the relative importance of the Southern Oscillation (SO) signal over preferred regions. Using lag cross-correlation statistics, possible lagged relationships between the tropics and midlatitudes and their relationships with the SO are then investigated. Only features that are consistent with present knowledge of the dynamics of the system are emphasized. Certain features which may not meet rigorous statistical significance tests but yet are either expected a priori from independent observations or are predicted from dynamical theories are also explored.

Lau, K.-M.↗

Discrete Dirac equation on a finite half-integer lattice

The formulation of the Dirac equation on a discrete lattice with half-integer spacing and periodic boundary conditions is investigated analytically. The importance of lattice formulations for problems in field theory and quantum mechanics is explained; the concept of half-integer Fourier representation is introduced; the discrete Dirac equation for the two-dimensional case is derived; dispersion relations for the four-dimensional case are developed; and the spinor formulation for the Dirac fields on the half-integer lattice and the discrete time variable for the four-dimensional time-dependent Dirac equation are obtained. It is argued that the half-integer lattice, because it takes the Dirac Lagrangian into account, is more than a mere relabeling of the integer lattice and may have fundamental physical meaning (e.g., for the statistics of fermions). It is noted that the present formulation does not lead to species doubling, except in the continuum limit.

Smalley, L. L.↗

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗

Propagation model for the Land Mobile Satellite channel in urban environments

This paper presents the major characteristics of a simulation package capable of performing a complete narrow and wideband analysis of the mobile satellite communication channel in urban environments for any given orbital configuration. The wavelength-to-average urban geometrical dimension ratio has required the use of the Geometrical Theory of Diffraction (GTD). For the RF frequency range, the model has been designed to be (1 up to 60 GHz) extended to include effects of non-perfect conductivity and surface roughness. Taking advantage of the inherent capabilities of such a high frequency method, we are able to provide a complete description of the electromagnetic field at the mobile terminal. Using the information made available at the ray-tracer and GTD solver outputs, the Land Mobile Satellite (LMS) urban model can also give a detailed description of the communication channel in terms of power delay profiles, Doppler spectra, channel scattering functions, and so forth. Statistical data, e.g. cumulative distribution functions, level crossing rates or distributions of fades are also provided. The user can access the simulation tool through a Design-CAD user-friendly interface by means of which she can effectively design her own urban layout and run consequently all the envisaged routines. The software is optimized in its execution time so that numerous runs can be achieved in a considerably short time.

Sforza, M.↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

The microwave spectrum of solar millisecond spikes

The microwave radiation from solar flares sometimes shows short and intensive spikes which are superimposed on the burst continuum. In order to determine the upper frequency limit of their occurrence and the circular polarization, a statistical analysis was performed on digital microwave observations from 3.2 to 92.5 GHz. Additionally, fine structures were investigated with a fast 32-channel spectrometer at 3.47 GHz. It was found that about 10 percent of the bursts show fine structures at 3.2 and 5.2 GHz, whereas none occurred above 8.4 GHz. Most of the observed spikes were very short and their bandwidth varied from below 0.5 MHz to more than 200 MHz. Simultaneous observations at two further frequencies showed no coincident spikes at the second and third harmonic. The observations can be explained by the theory of electron cyclotron masering if the observed bandwidths are determined by magnetic field inhomogeneities or if the rise times are independent of the source diameters. The latter would imply source sizes between 50 and 100 km.

Staehli, M.↗

Forecast of an exceptionally large even-numbered solar cycle

Using the 'dynamo theory' method to predict solar activity, an accurate prediction was made for solar cycle 21 by Schatten et al. (1978). Using the same dynamo technique for solar cycle 22, a value for the smoothed sunspot number of 170 + or - 25 is obtained. This large sunspot number is expected to peak in 1990 + or - 1 year. The F(10.7) radio flux is expected to reach a smoothed value of 210 + or - 25 flux units. Since this value is larger than values obtained with prediction schemes based upon 'statistical' and 'periodicity' methods, it provides a useful test for the current methodology, based upon the strength of the sun's polar field near solar minimum. The predicted degree of solar activity is expected to enhance the density and temperature of the earth's thermosphere to values somewhat larger than those found in solar cycle 21. This will have an impact on the orbital lifetime of low altitude satellites.

Schatten, Kenneth H.↗

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗

Microstructural Influence on Deformation and Fatigue Life of Composites Using the Generalized Method of Cells

A fully coupled deformation and damage approach to modeling the response of composite materials and composite laminates is presented. It is based on the semi-­‐analytical generalized method of cells (GMC) micromechanics model as well as its higher fidelity counterpart, HFGMC, both of which provide closed-form constitutive equations for composite materials as well as the micro scale stress and strain fields in the composite phases. The provided constitutive equations allow GMC and HFGMC to function within a higher scale structural analysis (e.g., finite element analysis or lamination theory) to represent a composite material point, while the availability of the micro fields allow the incorporation of lower scale sub­‐models to represent local phenomena in the fiber and matrix. Further, GMC's formulation performs averaging when applying certain governing equations such that some degree of microscale field accuracy is surrendered in favor of extreme computational efficiency, rendering the method quite attractive as the centerpiece in a integrated computational material engineering (ICME) structural analysis; whereas HFGMC retains this microscale field accuracy, but at the price of significantly slower computational speed. Herein, the sensitivity of deformation and the fatigue life of graphite/epoxy PMC composites, with both ordered and disordered microstructures, has been investigated using this coupled deformation and damage micromechanics based approach. The local effects of fiber breakage and fatigue damage are included as sub‐models that operate on the microscale for the individual composite phases. For analysis of laminates, classical lamination theory is employed as the global or structural scale model, while GMC/HFGMC is embedded to operate on the microscale to simulate the behavior of the composite material within each laminate layer. A key outcome of this study is the statistical influence of microstructure and micromechanics idealization (GMC or HFGMC) on the overall accuracy of unidirectional and laminated composite deformation and fatigue response.

Microscale Models↗

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗