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 289 records · Page 16

A new look at decomposition of turbulence forcing field and the structural response

Measured cross-spectrum of a turbulence field usually shows some decay in the statistical correlation in addition to convection at a characteristic velocity. It is shown that a decaying turbulence can be decomposed into frozen-pattern components thus permitting a simpler way to calculate the structural response. This procedure also provides a relationship whereby the measured input spectra can be incorporated. The theory is applied to an infinite beam which is backed on one side by a fluid filled cavity and is exposed on the other side by the turbulence excitation. The effect of the free stream velocity is also taken into consideration.

Lin, Y. K.↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Digital image profilers for detecting faint sources which have bright companions, phase 2

A breadboard image profiling system developed for the first phase of this project demonstrated the potential for detecting extremely faint optical sources in the presence of light companions. Experimental data derived from laboratory testing of the device supports the theory that image profilers of this type may approach the theoretical limit imposed by photon statistics. The objective of Phase 2 of this program is the development of a ground-based multichannel image profiling system capable of detecting faint stellar objects slightly displaced from brighter stars. We have finalized the multichannel image profiling system and attempted three field tests.

Morris, Elena↗

Information-Theoretic Assessment of Sample Imaging Systems

By rigorously extending modern communication theory to the assessment of sampled imaging systems, we develop the formulations that are required to optimize the performance of these systems within the critical constraints of image gathering, data transmission, and image display. The goal of this optimization is to produce images with the best possible visual quality for the wide range of statistical properties of the radiance field of natural scenes that one normally encounters. Extensive computational results are presented to assess the performance of sampled imaging systems in terms of information rate, theoretical minimum data rate, and fidelity. Comparisons of this assessment with perceptual and measurable performance demonstrate that (1) the information rate that a sampled imaging system conveys from the captured radiance field to the observer is closely correlated with the fidelity, sharpness and clarity with which the observed images can be restored and (2) the associated theoretical minimum data rate is closely correlated with the lowest data rate with which the acquired signal can be encoded for efficient transmission.

Huck, Friedrich O.↗

Stochastic Growth Theory of Spatially-Averaged Distributions of Langmuir Fields in Earth's Foreshock

Langmuir-like waves in the foreshock of Earth are characteristically bursty and irregular, and are the subject of a number of recent studies. Averaged over the foreshock, it is observed that the probability distribution is power-law P(bar)(log E) in the wave field E with the bar denoting this averaging over position, In this paper it is shown that stochastic growth theory (SGT) can explain a power-law spatially-averaged distributions P(bar)(log E), when the observed power-law variations of the mean and standard deviation of log E with position are combined with the log normal statistics predicted by SGT at each location.

Boshuizen, Christopher R.↗

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning BPS spectra and the gap conjecture

We explore statistical properties of Bogomol’nyi-Prasad-Sommerfield q-series for strongly coupled supersymmetric theories that correspond to a particular family of three-manifolds. We discover that gaps between exponents in the -series are statistically more significant at the beginning of the -series compared to gaps that appear in higher powers of. Our observations are obtained by calculating saliencies of -series features used as input data for principal component analysis, which is a standard example of an explainable machine learning technique that allows for a direct calculation and a better analysis of feature saliencies.

97 MATHEMATICS AND COMPUTING↗

Observation of the modification of quantum statistics of plasmonic systems

Abstract For almost two decades, researchers have observed the preservation of the quantum statistical properties of bosons in a large variety of plasmonic systems. In addition, the possibility of preserving nonclassical correlations in light-matter interactions mediated by scattering among photons and plasmons stimulated the idea of the conservation of quantum statistics in plasmonic systems. It has also been assumed that similar dynamics underlie the conservation of the quantum fluctuations that define the nature of light sources. So far, plasmonic experiments have been performed in nanoscale systems in which complex multiparticle interactions are restrained. Here, we demonstrate that the quantum statistics of multiparticle systems are not always preserved in plasmonic platforms and report the observation of their modification. Moreover, we show that optical near fields provide additional scattering paths that can induce complex multiparticle interactions. Remarkably, the resulting multiparticle dynamics can, in turn, lead to the modification of the excitation mode of plasmonic systems. These observations are validated through the quantum theory of optical coherence for single- and multi-mode plasmonic systems. Our findings unveil the possibility of using multiparticle scattering to perform exquisite control of quantum plasmonic systems.

36 MATERIALS SCIENCE↗

Production of neutron-rich isotopes for 𝑍 ≥ 98 in the 238 U + 248 Cm reaction

Multinucleon transfer (MNT) reactions in actinide systems are a promising method to synthesize transuranium neutron-rich elements. Appropriate theoretical approaches are needed to understand the mechanisms behind MNT. We employ a microscopic approach to calculate neutron-rich isotope production in the reaction 238 U + 248 Cm system. Here, the stochastic mean-field (SMF) approach is used to calculate the primary cross sections in MNT reactions based on the quasifission and inverse quasifission processes, and a statistical de-excitation model with GEMINI ++ code to calculate the secondary fragment cross sections. The calculated cross sections using SMF and GEMINI ++ explain available experimental results for the 238 U + 248 Cm system at 𝐸 c.m. = 898.7 MeV energy. This shows the effectiveness and applicability of the quantal diffusion approach, based on the SMF theory, in heavy-ion collisions. Production of transuranium neutron-rich elements with a proton number up to 𝑍=101 is obtained with sizable cross sections. Theoretical results calculated for the 𝑍=102–105 region, for which there are no experimental data, show that the cross-section values would be lower than the microbarn level. SMF theory does not contain any adjustable parameters other than the standard parameters of the energy density functional used in the TDHF theory and is an important approach for the microscopic understanding of reaction mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collective dynamics of polarized spin-half fermions in relativistic heavy-ion collisions

Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Steady induction effects in geomagnetism. Part 1A: Steady motional induction of geomagnetic chaos

Geomagnetic effects of magnetic induction by hypothetically steady fluid motion and steady magnetic flux diffusion near the top of Earth's core are investigated using electromagnetic theory, simple magnetic earth models, and numerical experiments with geomagnetic field models. The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation indicated by broad-scale models of the observed geomagnetic field is examined and solved. In Part 1, the steady surficial core flow estimation problem is solved in the context of the source-free mantle/frozen-flux core model. In the first paper (IA), the theory underlying such estimates is reviewed and some consequences of various kinematic and dynamic flow hypotheses are derived. For a frozen-flux core, fluid downwelling is required to change the mean square normal magnetic flux density averaged over the core-mantle boundary. For surficially geostrophic flow, downwelling implies poleward flow. The solution of the forward steady motional induction problem at the surface of a frozen-flux core is derived and found to be a fine, easily visualized example of deterministic chaos. Geomagnetic effects of statistically steady core surface flow may well dominate secular variation over several decades. Indeed, effects of persistent, if not steady, surficially geostrophic core flow are described which may help explain certain features of the present broad-scale geomagnetic field and perhaps paleomagnetic secular variation.

Voorhies, Coerte V.↗

Structure of turbulence at high shear rate

The structure of homogeneous turbulence subject to high shear rate has been investigated by using three-dimensional, time-dependent numerical simulations of the Navier-Stokes equations. This study indicates that high shear rate alone is sufficient for generation of the streaky structures, and that the presence of a solid boundary is not necessary. Evolution of the statistical correlations is examined to determine the effect of high shear rate on the development of anisotropy in turbulence. It is shown that the streamwise fluctuating motions are enhanced so profoundly that a highly anisotropic turbulence state with a 'one-component' velocity field and 'two-component' vorticity field develops asymptotically as total shear increases. Because of high-shear rate, rapid distortion theory predicts remarkably well the anisotropic behavior of the structural quantities.

Lee, Moon Joo↗

Noniterative finite amplitude methods for giant resonances and the application to the neutron radiative capture cross sections

We calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances with noniterative finite amplitude methods and demonstrate how the fully microscopic density functional theory predicts the giant resonances without any phenomenological parameters. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the result of E1 and M1 transitions and find that the capture cross sections for deformed nuclei are enhanced due to the contribution from the low energy M1 scissors mode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗

Measurement of the muon anomalous precession frequency at the Muon g − 2 Experiment at Fermilab

The anomalous magnetic moment of the muon, $a_\mu = \frac{g-2}{2}$ is the fractional deviation of the muon $g$-factor from the value of 2. It emerges as the cumulative effect of the virtual particles participating in the muon interaction with a magnetic field via quantum loop corrections. Its value encodes all the possible interactions between the virtual particles and, for this reason, represents an important test of the Standard Model (SM). In particular, any deviation from the SM theoretical evaluation could be due to new physics contributions. The new Muon $g-2$ (E989) Experiment at Fermilab is currently operating to repeat and improve the previous E821 experiment at Brookhaven National Laboratory (BNL), aiming to reduce the experimental error by a factor of 4 to the final accuracy of 140 parts per billion (ppb). On April 7th, 2021, the E989 collaboration published the first result based on the first year of data taking (Run-1), measuring $a_\mu = 0.001~165~920~40(54)$ with a precision of 460 ppb. The measured value is consistent with the BNL measurement and strengthens the long-standing tension with the data-driven SM prediction to a combined discrepancy of 4.2$\sigma$. On the theory side, however, new efforts involving lattice-QCD techniques are starting to question the current consensus on the theoretical prediction, demanding new improvements on both the experimental and theoretical sides. The E989 collaboration is now finalizing the analysis of Run-2 and Run-3 data and a new publication is expected in the first half of 2023 with a combined statistical uncertainty of 200 ppb. The anomalous magnetic moment $a_\mu$ is measured as the ratio between the muon spin anomalous precession frequency, $\omega_a$, and the average magnetic field experienced by the muons as they circulate in the storage ring. This thesis presents a precession frequency analysis of the Run-1 data and an evaluation of the related systematic uncertainties. A new positron reconstruction developed for the analysis of the subsequent data-taking periods, aiming to reduce some of the major systematic uncertainties of the $\omega_a$ measurement, is presented. The author's involvement in the production of the Run-2/5 data and in the precise calibration of the detectors is discussed. Finally, the complete Run-1 $a_\mu$ result is presented.

43 PARTICLE ACCELERATORS↗

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗

Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point

Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗