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,513 records · Page 84

A Study of the Gamma-Ray Burst Fundamental Plane

Long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obey a 3D relation, between the rest-frame time at the end of the plateau, T {sub a} , its corresponding X-ray luminosity, L {sub a} , and the peak luminosity in the prompt emission, L {sub peak}. This 3D relation identifies a GRB fundamental plane whose existence we here confirm. Here we include the most recent GRBs observed by Swift to define a “gold sample” (45 GRBs) and obtain an intrinsic scatter about the plane compatible within 1 σ with the previous result. We compare GRB categories, such as short GRBs with extended emission (SEE), X-ray flashes, GRBs associated with supernovae, a sample of only long-duration GRBs (132), selected from the total sample by excluding GRBs of the previous categories, and the gold sample, composed by GRBs with light curves with good data coverage and relatively flat plateaus. We find that the relation planes for each of these categories are not statistically different from the gold fundamental plane, with the exception of the SSE, which are hence identified as a physically distinct class. The gold fundamental plane has an intrinsic scatter smaller than any plane derived from the other sample categories. Thus, the distance of any particular GRB category from this plane becomes a key parameter. We computed the several category planes with T {sub a} as a dependent parameter obtaining for each category smaller intrinsic scatters (reaching a reduction of 24% for the long GRBs). The fundamental plane is independent from several prompt and afterglow parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum statistics and spacetime topology: Quantum surgery formulas

To formulate the universal constraints of quantum statistics data of generic long-range entangled quantum systems, we introduce the geometric-topology surgery theory on spacetime manifolds where quantum systems reside, cutting and gluing the associated quantum amplitudes, specifically in 2+1 and 3+1 spacetime dimensions. First, we introduce the fusion data for worldline and worldsheet operators capable of creating anyonic excitations of particles and strings, well-defined in gapped states of matter with intrinsic topological orders. Second, we introduce the braiding statistics data of particles and strings, such as the geometric Berry matrices for particle–string Aharonov–Bohm, 3-string, 4-string, or multi-string adiabatic loop braiding process, encoded by submanifold links, in the closed spacetime 3-manifolds and 4-manifolds. Third, we derive new “quantum surgery” formulas and constraints, analogous to Verlinde formula associating fusion and braiding statistics data via spacetime surgery, essential for defining the theory of topological orders, 3d and 4d TQFTs and potentially correlated to bootstrap boundary physics such as gapless modes, extended defects, 2d and 3d conformal field theories or quantum anomalies. This article is meant to be an extended and further detailed elaboration of our previous work Wang, Wen and Yau (2016) [1] and Chapter 6 of Wang (2015) [2]. Our theory applies to general quantum theories and quantum mechanical systems, also applicable to, but not necessarily requiring the quantum field theory description.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Disordered vortex matter out of equilibrium: a Langevin molecular dynamics study

We discuss the use of Langevin molecular dynamics in the investigation of the non-equilibrium properties of disordered vortex matter. Here, our special focus is set on values of system parameters that are realistic for disordered high-T c superconductors such as YBCO. Using a discretised elastic line model, we study different aspects of vortices far from thermal equilibrium. On the one hand, we investigate steady-state properties of driven magnetic flux lines in a disordered environment, namely the current–voltage characteristics, the gyration radius and the pinning time statistics. On the other hand, we study the complex relaxation processes and glassy-like dynamics that emerge in type-II superconductors due to the intricate competition between the long-range vortex–vortex repulsion and flux pinning due to randomly placed point defects. To this end, we consider different types of sudden perturbations: temperature, magnetic field and external current quenches.

74 ATOMIC AND MOLECULAR PHYSICS↗

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗

Quantum Liouville theorem based on Haar measure

Liouville theorem (L theorem) reveals robust incompressibility of the distribution function in phase space, given arbitrary potentials. However, its quantum generalization, Wigner flow, is compressible, i.e., L theorem is only conditionally true (e.g., for perfect Harmonic potential). Here, we develop quantum L theorem (rigorous incompressibility) for arbitrary potentials (interacting or not) in Hamiltonians. Haar measure, instead of symplectic measure dp$\bigwedge$dq used in Wigner’s scheme, plays a central role. The argument is based on general measure theory, independent of specific spaces or coordinates. Comparison of classical and quantum is made: for instance, here we address why Haar measure and metric preservation do not work in the classical case. Applications of the theorems in statistics, topological phase transition, ergodic theory, etc., are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Release of MCNP6.2 and WHISPER-1.1 - Guidance for NCS users

In 2017, a new version of the MCNP6 Monte Carlo code package was released: MCNP6.2. This is the first version of MCNP to include the WHISPER package to support sensitivity-uncertainty based nuclear criticality safety validation. The coding changes to MCNP6.2 physics are relatively insignificant. Corrections to the S(α,β) thermal scattering numerics are generally negligible relative to problem statistics (or, in rare cases, prevent aborts). Similarly, the changes to adjoint-weighting for computing kinetics parameters may result in small differences, generally negligible compared to problem statistics. The change to the MCNP6.2 geometry treatment to correctly handle coincident surfaces in problems with universe/fill features will produce different round-off in the geometry tracking. This will produce differences in results relative to previous versions, but those differences should be small relative to problem statistics, and are not a concern. An immediate benefit of using MCNP6.2 (rather than MCNP6.1) is that the new version is typically 1.5-2 times faster.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Duality and topology

Mappings between models may be obtained by unitary transformations with preservation of the spectra but in general a change in the states. Non-canonical transformations in general also change the statistics of the operators involved. In these cases one may expect a change of topological properties as a consequence of the mapping. Here we consider some dualities resulting from mappings, by systematically using a Majorana fermion representation of spin and fermionic problems. We focus on the change of topological invariants that results from unitary transformations taking as examples the mapping between a spin system and a topological superconductor, and between different fermionic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified “local” Brink-Axel hypothesis for Gamow-Teller transitions for pf-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semiconverged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing results to previous calculations of astrophysical rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

FLC Nomination: Stephen R. Lee

STEM Mentorship Award: Given to an individual who may or may not be a direct employee of a government laboratory or facility in the Mid Continent Region, but has demonstrated exceptional character, mentorship, and creativity in the areas of science and math or by promoting the development of young people in grades K-12 or higher in technical career fields. Contributions must be unique and show measurable benefit to the laboratory, school system/students and must have occurred within the past three calendar years. Stephen R. Lee is the division leader of the Computer, Computational, and Statistical Sciences Division (CCS) at Los Alamos National Laboratory (LANL or the Laboratory). He holds a Bachelor of Science in Physics from New Mexico State University and a Master of Science in Physics from Kansas State University and began his 27-year career at the Laboratory in 1986. As division leader of CCS, Stephen is responsible for the management and oversight of CCS and all of its activities, as well as approximately 200 technical staff engaged in computer and computational, and statistical sciences. His experience and interests include computational algorithm development for largescale parallel computing (peta-scale and beyond), predictive simulations for nuclear weapons, simulation capabilities for non-nuclear weapons applications such as bioscience, homeland defense, nanotechnology, advanced computational platform research and development.

97 MATHEMATICS AND COMPUTING↗

Impact of model improvements on 80 m wind speeds during the second Wind Forecast Improvement Project (WFIP2)

Abstract. During the second Wind Forecast Improvement Project (WFIP2; October 2015–March 2017, held in the Columbia River Gorge and Basin area of eastern Washington and Oregon states), several improvements to the parameterizations used in the High Resolution Rapid Refresh (HRRR – 3 km horizontal grid spacing) and the High Resolution Rapid Refresh Nest (HRRRNEST – 750 m horizontal grid spacing) numerical weather prediction (NWP) models were tested during four 6-week reforecast periods (one for each season). For these tests the models were run in control (CNT) and experimental (EXP) configurations, with the EXP configuration including all the improved parameterizations. The impacts of the experimental parameterizations on the forecast of 80 m wind speeds (wind turbine hub height) from the HRRR and HRRRNEST models are assessed, using observations collected by 19 sodars and three profiling lidars for comparison. Improvements due to the experimental physics (EXP vs. CNT runs) and those due to finer horizontal grid spacing (HRRRNEST vs. HRRR) and the combination of the two are compared, using standard bulk statistics such as mean absolute error (MAE) and mean bias error (bias). On average, the HRRR 80 m wind speed MAE is reduced by 3 %–4 % due to the experimental physics. The impact of the finer horizontal grid spacing in the CNT runs also shows a positive improvement of 5 % on MAE, which is particularly large at nighttime and during the morning transition. Lastly, the combined impact of the experimental physics and finer horizontal grid spacing produces larger improvements in the 80 m wind speed MAE, up to 7 %–8 %. The improvements are evaluated as a function of the model's initialization time, forecast horizon, time of the day, season of the year, site elevation, and meteorological phenomena. Causes of model weaknesses are identified. Finally, bias correction methods are applied to the 80 m wind speed model outputs to measure their impact on the improvements due to the removal of the systematic component of the errors.

17 WIND ENERGY↗

Proton Isovector Helicity Distribution on the Lattice at Physical Pion Mass

We present a state-of-the-art calculation of the isovector quark helicity Bjorken-x distribution in the proton using lattice-QCD ensembles at the physical pion mass. We compute quasi-distributions at proton momenta P z ϵ {2.2; 2.6; 3.0} GeV on the lattice, and match them systematically to the physical parton distribution using large-momentum effective theory (LaMET).We reach an unprecedented precision through high statistics in simulations, large-momentum proton matrix elements, and control of excited-state contamination. The resulting distribution with combined statistical and systematic errors is in agreement with the latest phenomenological analysis of the spin-dependent experimental data; in particular, Δ$\bar{u}$(x) > Δ$\bar{d}$(x).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗

Lattice QCD Calculation of the Subtraction Function in Forward Compton Amplitude

The subtraction function plays a pivotal role in calculations involving the forward Compton amplitude, which is crucial for predicting the Lamb shift in muonic atoms, as well as the proton-neutron mass difference. In this Letter, we present a lattice QCD calculation of the subtraction function using two domain wall fermion gauge ensembles near the physical pion mass. We utilize a recently proposed subtraction point, demonstrating its advantage in mitigating statistical and systematic uncertainties by eliminating the need for ground-state subtraction. Our results reveal significant contributions from N⁢π intermediate states to the subtraction function. Incorporating these contributions, we compute the proton, neutron, and nucleon isovector subtraction functions at photon momentum transfer Q 2 ϵ[0,2] GeV 2 . For the proton subtraction function, we compare our lattice results with chiral perturbation theory prediction at low Q 2 and with the results from the perturbative operator-product expansion at high Q 2 . Finally, using these subtraction functions as input, we determine their contribution to two-photon exchange effects in the Lamb shift and isovector nucleon electromagnetic self-energy.

74 ATOMIC AND MOLECULAR PHYSICS↗

Urban forms, physical activity and body mass index: a cross-city examination using ISS Earth Observation photographs

Johnson Space Center has archived thousands of astronauts acquired Earth images. Some spectacular images have been widely used in news media and in k-12 class room, but their potential utilizations in health promotion and disease prevention have relatively untapped. The project uses daytime ISS photographs to define city forms and links them to city or metropolitan level health data in a multicity context. Road connectivity, landuse mix and Shannon's information indices were used in the classification of photographs. In contrast to previous remote-sensing studies, which tend to focus on a single city or a portion of a city, this project utilized photographs of 39 U.S. cities. And in contrast to previous health-promotion studies on the built environment, which tend to rely on survey respondents' responses to evaluate road connectivity or mixed land use for a single study site, the project examined the built environments of multiple cities based on ISS photos. It was found that road connectivity and landuse mix were not statistically significant by themselves, but the composite measure of the Shannon index was significantly associated with physical activity, but not BMI. Consequently, leisure-time physical activity seems to be positively associated with the urban complexity scale. It was also concluded that unless they are planned or designed in advance, photographs taken by astronauts generally are not appropriate for a study of a single-site built environment nor are they appropriate for a study of infectious diseases at a local scale. To link urban built environment with city-wide health indicators, both the traditional nadir view and oblique views should be emphasized in future astronauts' earth observation photographs.

Lin, Ge↗

Ground testing at MSFC

The large-scale field collection of acoustic data and their analysis are described in this review. The entire sound-field program has the purpose of solving complex problems of excessive noise propagation associated with the static test firing of large space-booster powerplants. The site of data collection has been Mississippi Test Facility (MTF), where far-field sound-propagation and correlated meteorological measurements were initiated in 1962. The statistical evaluation of sound-measurement data shows that the polar distribution of sound in the springtime (at MTF) is determined mainly by the temperature field of the atmosphere. Greater deformation of sound-intensity polar distribution occurs in the autumn as a result of the influence of the wind field. The data on sound attenuation as a function of distance from its source revealed attenuation in excess of the expected magnitude. The excess was mainly between 1 and 10 kilometers from the source, with maximum attenuation (17 to 18 decibels) near 6 kilometers. Excess sound attenuation was compared with calculated sound dispersion. For all sound frequencies in all seasons the deviation between reference and empirical attenuation values formed a bell-shaped distribution on a semilogarithmic plot. The physical causes for the sound- attenuation behavior are considered to be atmospheric diffraction and scattering effects. Statistical analysis of sod-intensity extremes indicated amplification of the sound signals by focusing which could be expected along selected azimuths. Empirical amplification values showed the maximum occurring in winter (50 decibels at 22 kilometers). For the analysis of correlation between empirical sound data and local atmospheric sound- propagation conditions, #e concept of analytical classification of sound-velocity profiles has been developed. Examination of over a quarter-million individual sound-velocity profiles indicates that 33 profile types for heights up to 3 kilometers constitute the empirical possibilities of sound-propagation conditions in southeastern United States. Since the focusing qualities of each sound- velocity profile can be determined as a class characteristic, the overall probabilities of focusing occurrence can he obtained. The severity of focusing, as based upon lateral spread of focusing areas, has been determined. The data for MTF indicate focusing is most severe in winter, with the direction toward east being most affected. Sound propagation is a fluctuating phenomenon because atmospheric conditions vary continuously. A measuring system for analyzing the nonstationary sound-propagation phenomena has been established. The analysis of the acoustic fluctuation data is sup- posed to provide "dynamic” perturbation functions of sound, which are to be superimposed on the related characteristic sound-velocity-profile type. The ultimate result of the program will be the establishment of reliable contingency tables showing the relationship between atmospheric parameters and far-field sound intensities. These data will be the basis for more accurate sound-propagation forecasting for static firing tests.

STATIC TESTING↗

The Ginibre evolution in the large-N limit

We analyse statistics of the real eigenvalues of gl(N, R)-valued Brownian motion (the Ginibre evolution) in the limit of large N. In particular, we calculate the limiting two-time correlation function of spin variables associated with real eigenvalues of the Ginibre evolution. We also show how the formalism of spin variables can be used to compute the fixed time correlation functions of real eigenvalues discovered originally by Forrester and Nagao [“Eigenvalue statistics of the real Ginibre ensemble,” Phys. Rev. Lett. 99(5), 050603 (2007)] and Borodin and Sinclair [“The Ginibre ensemble of real random matrices and its scaling limits,” Commun. Math. Phys. 291(1), 177–224 (2009)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Gaussian wave packet phase-space representation of quantum canonical statistics

We present a mapping of quantum canonical statistical averages onto a phase-space average over thawed Gaussian wave-packet (GWP) parameters, which is exact for harmonic systems at all temperatures. The mapping invokes an effective potential surface, experienced by the wave packets, and a temperature-dependent phase-space integrand, to correctly transition from the GWP average at low temperature to classical statistics at high temperature. Numerical tests on weakly and strongly anharmonic model systems demonstrate that thermal averages of the system energy and geometric properties are accurate to within 1% of the exact quantum values at all temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗