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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Quantum statistical theory of semiconductor junctions in thermal equilibrium

Free carrier and electric field distributions of one-dimensional semiconductor junctions are evaluated using a quantum mechanical phase-space distribution and its corresponding Boltzmann equation. Attention is given to quantum and exchange corrections in cases of high doping concentrations when carrier densities become degenerate. Quantitative differences between degenerate and classical junction characteristics, e.g., maximum electric field and built-in voltage and carrier concentration within the transition region, are evaluated numerically.

Von Roos, O.↗
NASA NTRS1990

Wave propagation through random media: A local method of small perturbations based on the Helmholtz equation

Propagation of sound through the turbulent atmosphere is a statistical problem. The randomness of the refractive index field causes sound pressure fluctuations. Although no general theory to predict sound pressure statistics from given refractive index statistics exists, there are several approximate solutions to the problem. The most common approximation is the parabolic equation method. Results obtained by this method are restricted to small refractive index fluctuations and to small wave lengths. While the first condition is generally met in the atmosphere, it is desirable to overcome the second. A generalization of the parabolic equation method with respect to the small wave length restriction is presented.

Grosse, Ralf↗
NASA NTRS1999

Characterization of Sound Radiation by Unresolved Scales of Motion in Computational Aeroacoustics

Evaluation of the sound sources in a high Reynolds number turbulent flow requires time-accurate resolution of an extremely large number of scales of motion. Direct numerical simulations will therefore remain infeasible for the forseeable future: although current large eddy simulation methods can resolve the largest scales of motion accurately the, they must leave some scales of motion unresolved. A priori studies show that acoustic power can be underestimated significantly if the contribution of these unresolved scales is simply neglected. In this paper, the problem of evaluating the sound radiation properties of the unresolved, subgrid-scale motions is approached in the spirit of the simplest subgrid stress models: the unresolved velocity field is treated as isotropic turbulence with statistical descriptors, evaluated from the resolved field. The theory of isotropic turbulence is applied to derive formulas for the total power and the power spectral density of the sound radiated by a filtered velocity field. These quantities are compared with the corresponding quantities for the unfiltered field for a range of filter widths and Reynolds numbers.

Rubinstein, Robert↗
NASA NTRS2000

Scale-Invariant Forms of Conservation Equations in Reactive Fields and a Modified Hydro-Thermo-Diffusive Theory of Laminar Flames

A scale-invariant model of statistical mechanics is applied to present invariant forms of mass, energy, linear, and angular momentum conservation equations in reactive fields. The resulting conservation equations at molecular-dynamic scale are solved by the method of large activation energy asymptotics to describe the hydro-thermo-diffusive structure of laminar premixed flames. The predicted temperature and velocity profiles are in agreement with the observations. Also, with realistic physico-chemical properties and chemical-kinetic parameters for a single-step overall combustion of stoichiometric methane-air premixed flame, the laminar flame propagation velocity of 42.1 cm/s is calculated in agreement with the experimental value.

Sohrab, Siavash H.↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Diagnosis of 3D magnetic field and mode composition in MHD turbulence with Y -parameter

ABSTRACT Magnetic fields are crucial in numerous astrophysical processes within the interstellar medium (ISM). However, the detailed determination of magnetic field geometry is notoriously challenging. Based on the modern magnetohydrodynamic (MHD) turbulence theory, we introduce a novel statistical technique, the ‘Y-parameter’, to decipher the magnetic field inclination in the ISM and identify dominant turbulence modes. The Y-parameter, calculated as the ratio of anisotropies of different Stokes parameter combinations, displays contrasting trends with the mean-field inclination angle in Alfvénic and compressible turbulence modes. A Y-parameter value around 1.5 ± 0.5 provide a statistical boundary to determine the dominant MHD turbulence modes. We have discovered specific correlations between the Y-parameter value and the inclination angle that unveil the dominant turbulence mode. This methodology, when applied to future radio polarization surveys such as LOFAR and SKA, promises to significantly enhance our knowledge of 3D magnetic field in the ISM and improve our understanding of interstellar turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

Interplanetary Alfvenic fluctuations - A stochastic model

The concept of minimum variance is investigated for nonplanar interplanetary Alfvenic fluctuations in which the field direction varies randomly. The theory of the random wandering of a vector of constant length is developed as a representation of the magnetic field, and it is found that the minimum variance tends to coincide with the mean field directions over the correlation time of the fluctuations. The Fokker-Planck limit of the theory is then developed and used to analyze the statistic distribution of field directions with and without a reflecting barrier. Results suggest that the tendency of the Alfvenic fluctuations to have a direction of minimum variance statistically aligned with the mean magnetic field may be purely a consequence of the randomness of the fluctuation and not imply that the fluctuations are necessarily plane waves. Extensive statistical studies of the observed directional variations of the interplanetary magnetic field are necessary to test this hypothesis.

Barnes, A.↗
NASA NTRS2011

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗
DOE OSTI2025

Wormholes with ends of the world

We study classical wormhole solutions in 3D gravity with end-of-the-world (EOW) branes, conical defects, kinks, and punctures. These solutions compute statistical averages of an ensemble of boundary conformal field theories (BCFTs) related to universal asymptotics of OPE data extracted from the 2D conformal bootstrap. Conical defects connect BCFT bulk operators; branes join BCFT boundary intervals with identical boundary conditions; kinks (1D defects along branes) link BCFT boundary operators; and punctures (0D defects) are endpoints where conical defects terminate on branes. We provide evidence for a correspondence between the gravity theory and the ensemble. In particular, the agreement of the g-function dependence results from an underlying topological aspect of the on-shell EOW brane action, from which a BCFT analog of the Schlenker-Witten theorem also follows.

AdS-CFT Correspondence↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Precision Computations in Strongly Coupled Conformal Field Theories (Final Technical Report)

Conformal Field Theories (CFTs) are quantum field theories that are invariant under the conformal symmetry group (which includes translations and rotations, but also local rescalings of spacetime). They are building blocks of general quantum field theories, and appear in many areas of physics, including statistical physics, condensed matter physics, particle physics, and quantum gravity. Because of their extra symmetries, the mathematical structure of CFTs is tightly constrained, and this leads to the idea of the ``conformal bootstrap," which is to use these mathematical structures to constrain, and in some cases determine, CFT observables. A new numerical implementation of the conformal bootstrap idea appeared in 2008 with the work of Rattazzi, Rychkov, Tonni, and Vichi. Their observation was that certain bootstrap constraints (conformal symmetry and unitarity) could be combined to yield a convex optimization problem that constraints CFT data. By solving this convex optimization problem on a computer, one could obtain bounds on observables like critical exponents and operator product expansion (OPE) coefficients. Over the course of this award, the PI has improved numerical bootstrap techniques by optimizing known algorithms and finding new ones for performing the required convex optimization computations. The PI has applied these techniques to compute high-precision observables in several important strongly-coupled systems. The PI has also explored both analytical and numerical bootstrap methods for constraining the space of low energy effective field theories of quantum gravity, and developed new analytical techniques for CFT and QFT more broadly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1986

The statistics of peaks of Gaussian random fields

A set of new mathematical results on the theory of Gaussian random fields is presented, and the application of such calculations in cosmology to treat questions of structure formation from small-amplitude initial density fluctuations is addressed. The point process equation is discussed, giving the general formula for the average number density of peaks. The problem of the proper conditional probability constraints appropriate to maxima are examined using a one-dimensional illustration. The average density of maxima of a general three-dimensional Gaussian field is calculated as a function of heights of the maxima, and the average density of 'upcrossing' points on density contour surfaces is computed. The number density of peaks subject to the constraint that the large-scale density field be fixed is determined and used to discuss the segregation of high peaks from the underlying mass distribution. The machinery to calculate n-point peak-peak correlation functions is determined, as are the shapes of the profiles about maxima.

Bardeen, J. M.↗
DOE OSTI2022

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Loops in 4+1d topological phases

2+1d topological phases are well characterized by the fusion rules and braiding/exchange statistics of fractional point excitations. In 4+1d, some topological phases contain only fractional loop excitations. What kind of loop statistics exist? We study the 4+1d gauge theory with 2-form \mathbb{Z}_2 ℤ 2 gauge field (the loop-only toric code) and find that while braiding statistics between two different types of loops can be nontrivial, the self “exchange” statistics are all trivial. In particular, we show that the electric, magnetic, and dyonic loop excitations in the 4+1d toric code are not distinguished by their self-statistics. They tunnel into each other across 3+1d invertible domain walls which in turn give explicit unitary circuits that map the loop excitations into each other. The SL(2, \mathbb{Z}_2 ℤ 2 ) symmetry that permutes the loops, however, cannot be consistently gauged and we discuss the associated obstruction in the process. Moreover, we discuss a gapless boundary condition dubbed the “fractional Maxwell theory” and show how it can be Higgsed into gapped boundary conditions. We also discuss the generalization of these results from the \mathbb{Z}_2 ℤ 2 gauge group to \mathbb{Z}_N ℤ N .

Chen, Xie↗
DOE OSTI2023

Beyond Generalized Eigenvalues

Two analysis techniques, the generalized eigenvalue method (GEM) or Prony's method (PM), are commonly used to analyze statistical estimates of correlation functions produced in lattice quantum field theory calculations. GEM takes full advantage of the matrix structure of correlation functions but only considers individual pairs of time separations when much more data exists. PM can be applied to many time separations and many individual matrix elements simultaneously but does not fully exploit the matrix structure of the correlation function. We combine both these methods into a single framework based on matrix polynomials which we call block Prony method (BPM).

Fleming, George T.↗
DOE OSTI2023

Chemical controls on the propagation and healing of subcritical fractures

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

36 MATERIALS SCIENCE↗
NASA NTRS1966

On a logical basis for division of responsibilities in statistical practice

The purpose of this paper is to explain principles for division of responsibilities between the statistician and the people that he works with, and reasons why this division of responsibilities is important -- that is, how it improves the performance of both statistician and expert in subject-matter. The aim is to find and illustrate principles of practice by which statisticians may make effective use of their knowledge of theory. The specialist in statistical methods may find himself applying the same basic theory in a dozen different fields in a week, rotating through the same projects the next week. Or, he may work day after day primarily in a single substantive field. Either way, he requires rules of practice. A statement of statistical reliability should present any information that might help the reader to form his own opinion concerning the validity of conclusions likely to be drawn from the results. The aim of a statistical report is to protect the client from seeing merely what he would like to see; to protect him from losses that could come from misuse of results. A further aim is to forestall unwarranted claims of accuracy that the client's public might otherwise accept.

Deming, W. Edwards↗