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,423 records · Page 79

Statistical Tests for Convergence in Monte Carlo Criticality Calculations

This report describes methods for using statistical tests to automatically determine convergence of Monte Carlo (MC) iterations in criticality calculations. The methods described herein are intended to replace the traditional approach used for the past 60 years: make a trial run; determine convergence based on plots of k-effective and Shannon entropy vs cycle; adjust the input parameters for controlling the iterations; make a final run to obtain results and statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

COG Beta-Effective Benchmarks

COG 11.3 calculations with the “new method” using ENDF/B-VIII.0 data are shown to produce results with |C/E – 1| < 3σ in all cases, whereas the other methods produce a few minor outliers. Furthermore, the “new method” using ENDF/B-VIII.0 data is observed to produce “over” and “under” results consistent with 50/50 expectations whereas the other methods appear to be slightly over or under predicting. JEFF-3.3 results are slightly larger overall than ENDF/B-VIII.0 results. In all but a few cases, it is also observed that $F_d/F_t ≥ (vF)_d/(vF)_t ≥ 1 – (k_{prompt}/k_{eff})$ and all estimators produce equivalent results for the most thermal systems. In those other cases, this relation holds within the given statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Progress Towards an Investigation of the MiniBooNE Low Energy Excess Using Neutral-Current Delta-Like Single Photons in MicroBooNE with Wire-Cell 3D Reconstruction Algorithms

This note presents progress towards a Neutral-Current (NC) Delta-like single photon search in MicroBooNE using Wire-Cell 3D reconstruction and pattern recognition algorithms. This analysis will help address the question of the MiniBooNE Low Energy Excess (LEE) [1]. The selections are similar to the previously published Pandora-based single photon selections [2] in many ways, but here we use independent reconstruction and selection tools and examine a larger phase space (including charged pions and multiple protons). In the future, we will significantly reduce systematic uncertainties via the use of constraining sideband observations in order to statistically test specific models of the MiniBooNE LEE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Progress towards an investigation of the MiniBooNE Low Energy Excess using neutral-current Delta-like single photons in MicroBooNE with Wire-Cell 3D reconstruction algorithms

This note presents progress towards a Neutral-Current (NC) Delta-like single photon search in MicroBooNE using Wire-Cell 3D reconstruction and pattern recognition algorithms. This analysis will help address the question of the MiniBooNE Low Energy Excess (LEE). The selections are similar to the previously published Pandora-based single photon selections in many ways, but here we use independent reconstruction and selection tools and examine a larger phase space (including charged pions and multiple protons). In the future, we will significantly reduce systematic uncertainties via the use of constraining sideband observations in order to statistically test specific models of the MiniBooNE LEE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probing stochastic inter-galactic magnetic fields using blazar-induced gamma ray halo morphology

Inter-galactic magnetic fields can imprint their structure on the morphology of blazar-induced gamma ray halos. We show that the halo morphology arises through the interplay of the source's jet and a two-dimensional surface dictated by the magnetic field. Through extensive numerical simulations, we generate mock halos created by stochastic magnetic fields with and without helicity, and study the dependence of the halo features on the properties of the magnetic field. We propose a sharper version of the Q-statistics and demonstrate its sensitivity to the magnetic field strength, the coherence scale, and the handedness of the helicity. We also identify and explain a new feature of the Q-statistics that can further enhance its power.

79 ASTRONOMY AND ASTROPHYSICS↗

CAN PLANETESIMALS FORM BY COLLISIONAL FUSION?

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk, we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin. Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities, we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution timescales.

79 ASTRONOMY AND ASTROPHYSICS↗

Diverse properties of interstellar medium embedding gamma-ray bursts at the epoch of reionization

Analysis is performed on ultra-high-resolution large-scale cosmological radiation-hydrodynamic simulations to quantify, for the first time, the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high-density (≥10{sup 4} cm{sup –3}) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density (∼10{sup –2.5} cm{sup –3}) high-temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature, and metallicity of intervening absorbers, vary strongly from sight line to sight line. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density, and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium (ISM) of the host galaxy based on a very small number (usually one) of sight lines would be precarious. Utilizing high-z GRBs to probe the ISM and intergalactic medium should be undertaken properly, taking into consideration the physical diversities of the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling baryonic physics in future weak lensing surveys

Modifications of the matter power spectrum due to baryonic physics are one of the major theoretical uncertainties in cosmological weak lensing measurements. Developing robust mitigation schemes for this source of systematic uncertainty increases the robustness of cosmological constraints, and may increase their precision if they enable the use of information from smaller scales. Here we explore the performance of two mitigation schemes for baryonic effects in weak lensing cosmic shear: the principal component analysis (PCA) method and the halo-model approach in hmcode. We construct mock tomographic shear power spectra from four hydrodynamical simulations, and run simulated likelihood analyses with cosmolike assuming LSST-like survey statistics. With an angular scale cut of ℓmax < 2000, both methods successfully remove the biases in cosmological parameters due to the various baryonic physics scenarios, with the PCA method causing less degradation in the parameter constraints than hmcode. For a more aggressive ℓmax = 5000, the PCA method performs well for all but one baryonic physics scenario, requiring additional training simulations to account for the extreme baryonic physics scenario of Illustris; hmcode exhibits tensions in the 2D posterior distributions of cosmological parameters due to lack of freedom in describing the power spectrum for $k \gt 10\ h^{-1}\, \mathrm{Mpc}$. We investigate variants of the PCA method and improve the bias mitigation through PCA by accounting for the noise properties in the data via Cholesky decomposition of the covariance matrix. Our improved PCA method allows us to retain more statistical constraining power while effectively mitigating baryonic uncertainties even for a broad range of baryonic physics scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Report of the theory panel

The ultimate goal of this research is to develop an understanding which is sufficiently comprehensive to allow realistic predictions of the behavior of the physical systems. Theory has a central role to play in the quest for this understanding. The level of theoretical description is dependent on three constraints: (1) the available computer hardware may limit both the number and the size of physical processes the model system can describe; (2) the fact that some natural systems may only be described in a statistical manner; and (3) the fact that some natural systems may be observable only through remote sensing which is intrinsically limited by spatial resolution and line of sight integration. From this the report discusses present accomplishments and future goals of theoretical space physics. Finally, the development and use of new supercomputer is examined.

Ashourabdalla, Maha↗

Latest results of searches for new physics in the D0 experiment

Results on searches for particles and phenomena beyond the Standard Model (new-physics effects) in the D0 experiment at the Tevatron accelerator (FNAL, USA) are examined. These results were obtained on the basis of a statistical sample corresponding to an integrated luminosity of 5 to 7 fb{sup -1} and were published in 2010 and 2011.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hadronic contributions to (g - 2) µ

The Muon g-2 Experiment at Fermilab, which recently started running, plans to reduce the uncertain- ties on the already very precisely measured anomalous magnetic moment of the muon by a factor of four. The goal of this effort is to probe the observed difference of more than three standard deviations between Standard-Model theory and experiment, one of the few persistent hints for physics beyond the Standard Model. The Fermilab experiment collected data from their first run last year with statistics comparable to BNL E821. They expect to release the measurement result in 2019. On the theoretical side, because the muon g - 2 arises from quantum- mechanical loop contributions in the Standard Model, it is sensitive to virtual effects of new particles, and places important constraints on Standard-Model extensions. To leverage the anticipated reduction in experimental errors, and determine unambiguously whether or not new-physics effects contribute to this quantity, the theoretical errors must be made more reliable and reduced to a commensurate precision. The Muon g-2 Theory Initiative was created to facilitate this development. The dominant sources of uncertainty in the Standard-Model prediction of the muon g -2 are from the hadronic contributions. The hadronic vacuum polarization (HVP) provides the leading correction followed by hadronic light-by-light (HLbL) scattering. There are a number of complementary theoretical efforts underway to better understand and quantify these contributions, including dispersive and data driven methods, lattice QCD, and effective field theories. Given the precision goals and the phenomenological importance, it is important to have more than one independent method for each of the two hadronic corrections, each with fully quantified uncertainties. Fostering the development of such methods is a prime goal of the initiative, as this will enable critical cross checks, and, upon combination, may yield gains in precision, to maximize the impact of E989. An important aspect of the Muon g-2 Theory Initiative’s activities are providing platforms that facilitate interactions between the different groups, as well as between the theoretical and experimental g - 2 communities. To this end, several workshops were organized in 2017 and 2018. The first meeting, held at Fermilab (June 3–6, 2017, St. Charles, IL, USA), served to kick-off the Initiative’s activities. Two meetings in early 2018 were focused respectively on the HVP and HLbL corrections. The HVP meeting was held at KEK (February 12–14, 2018, Tsukuba, Japan) and the HLbL meeting at the University of Connecticut (March 12–14, 2018, Storrs, CT, USA). The most recent workshop, which served as the second plenary meeting of the Initiative, was held at the University of Mainz (June 18–22, 2018, Mainz, Germany). An important outcome of these meetings are concrete plans for a first white paper, which is currently being written. We aim to post the white paper just prior to the release of the first E989 measurement, to present a clean theoretical prediction. The first white paper is focused on assessing and improving the reliability of the SM prediction. The purpose of the INT workshop in September 2019 is to start the next stage of the Initiative, focusing on the development of strategies to improve the theory uncertainties beyond the current level towards the E989 precision goal. We aim to accelerate theoretical developments on the hadronic contributions to the muon g - 2 so that the Standard-Model theory error can be brought to the needed precision, again in advance of the next release from E989. Hence the workshop will provide crucial theory support for a US experiment with broad impact. The Muon g-2 Theory Initiative relies on input from representatives of all the different communities that are engaged in this effort. It is therefore important that all these areas are properly represented. The funds from this grant will be used to enable more people, especially early-career scientists, to participate and make essential contributions to the workshop discussions. Since the workshop’s main goal is to kick-off the next stage of the theory initiative’s activities, support for this workshop from the DOE will help the theory community provide crucial theoretical support to the Fermilab Muon g-2 Experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stacked reverberation mapping of high-redshift quasars in DESI. I. Feasibility analysis

The broad-line region of quasars has long been probed by reverberation mapping techniques that measure time lags between continuum and broad emission-line variations. Stacked reverberation mapping has been proposed as a less observationally expensive alternative to traditional methods. This ensemble approach also reduces biases from small-number statistics. The Dark Energy Spectroscopic Instrument (DESI) is conducting the most extensive spectroscopic survey of quasars to date. We create mock light curves emulating expected DESI quasar observations at redshifts $1.48\lt z\lt 5.2$ and luminosities $44.68 \le \log \lambda L_{1350 \mathring{\rm A}{}} / \mathrm{erg\, s^{-1}} \le 45.99$ to test stacked reverberation mapping feasibility using sparse spectroscopic data paired with well-sampled photometric data. The pipeline, using the lag estimation code JAVELIN (Just Another Vehicle for Estimating Lags In Nuclei), successfully recovers the simulated C IV lags within 1σ of the true values using spectroscopic light curves composed of only a few spectral epochs (2–10) with irregular cadences. We investigate how observational factors, including C IV flux error magnitude, number of stacked quasars, and spectral epoch count, affect performance. This work motivates a pathway for future stacked reverberation mapping projects with large-scale spectroscopic surveys of quasars having $\ge 2$ spectroscopic observations. Our results suggest an economical alternative for constraining and extending the radius–luminosity relation to higher redshifts and luminosities. Subsequently, this relation can be employed more reliably in single-epoch black hole mass measurements and quasar cosmology in these distant regimes.

quasars: general, quasars: supermassive black hole↗

Quantum state-resolved differential cross sections for complex-forming chemical reactions: Asymmetry is the rule, symmetry the exception

We argue that statistical theories are generally unable to accurately predict state-resolved differential cross sections for triatomic bimolecular reactions studied in beam experiments, even in the idealized limit where the dynamics are fully chaotic. The basic reason is that quenching of interferences between partial waves is less efficient than intuitively expected, especially around the poles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bubble Departure Diameter and Bubble Release Frequency Measurement from TAMU Subcooled Flow Boiling Experiment

The bubble departure diameter and bubble release frequency were obtained through the analysis of TAMU subcooled flow boiling experimental data. The numerous images of bubbles at departure were analyzed for each experimental condition to achieve the reliable statistics of the measured bubble parameters. The results are provided in this report with simple discussion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Use of dMLC for implementation of dynamic respiratory-gated radiation therapy

Purpose: To simulate and evaluate the use of dynamic multileaf collimators (dMLC) in respiratory gating to compensate for baseline drift.Methods: Tumor motion tracking data from 30 lung tumors over 322 treatment fractions was analyzed with the finite state model. A dynamic respiratory gating window was established in real-time by determining the average positions during the previous two end-of-expiration breathing phases and centering the dMLC aperture on a weighted average of these positions. A simulated dMLC with physical motion constraints was used in dynamic gating treatment simulations. Fluence maps were created to provide a statistical description of radiation delivery for each fraction. Duty cycle was also calculated for each fraction.Results: The average duty cycle was 2.3% greater under dynamic gating conditions. Dynamic gating also showed higher fluences and less tumor obstruction. Additionally, dynamic gating required fewer beam toggles and each delivery period was longer on average than with static gating.Conclusions: The use of dynamic gating showed better performance than static gating and the physical constraints of a dMLC were shown to not be an impediment to dynamic gating.

61 RADIATION PROTECTION AND DOSIMETRY↗

Algorithms for Non-Negative Matrix Factorization on Noisy Data With Negative Values

Non-negative matrix factorization (NMF) is a dimensionality reduction technique that has shown promise for analyzing noisy data, especially astronomical data. For these datasets, the observed data may contain negative values due to noise even when the true underlying physical signal is strictly positive. Prior NMF work has not treated negative data in a statistically consistent manner, which becomes problematic for low signal-to-noise data with many negative values. In this paper we present two algorithms, Shift-NMF and Nearly-NMF, that can handle both the noisiness of the input data and also any introduced negativity. Both of these algorithms use the negative data space without clipping or masking and recover non-negative signals without any introduced positive offset that occurs when clipping or masking negative data. We demonstrate this numerically on both simple and more realistic examples, and prove that both algorithms have monotonically decreasing update rules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Comprehensive Monte Carlo Framework for Jet-Quenching

This article presents the motivation for developing a comprehensive modeling framework in which different models and parameter inputs can be compared and evaluated for a large range of jet-quenching observables measured in relativistic heavy-ion collisions at RHIC and the LHC. Here, the concept of a framework is discussed within the context of recent efforts by the JET Collaboration, the authors of JEWEL, and the JETSCAPE collaborations. The framework ingredients for each of these approaches is presented with a sample of important results from each. The role of advanced statistical tools in comparing models to data is also discussed, along with the need for a more detailed accounting of correlated errors in experimental results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗