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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pion valence quark distribution at physical pion mass of N f = 2 + 1 + 1 lattice QCD

Abstract We present a state-of-the-art calculation of the unpolarized pion valence-quark distribution in the framework of large-momentum effective theory (LaMET) with improved handling of systematic errors as well as two-loop perturbative matching. We use lattice ensembles generated by the MILC collaboration at lattice spacinga≈ 0.09 fm, lattice volume 64 3 × 96,N f = 2 + 1 + 1 flavors of highly-improved staggered quarks and a physical pion mass. The LaMET matrix elements are calculated with pions boosted to momentumP z ≈ 1.72 GeV with high-statistics ofO(10 6 ) measurements. We study the pion PDF in both hybrid-ratio and hybrid-regularization-independent momentum subtraction (hybrid-RI/MOM) schemes and also compare the systematic errors with and without the addition of leading-renormalon resummation (LRR) and renormalization-group resummation (RGR) in both the renormalization and lightcone matching. The final lightcone PDF results are presented in the modified minimal-subtraction scheme at renormalization scaleμ= 2.0 GeV. We show that thex-dependent PDFs are compatible between the hybrid-ratio and hybrid-RI/MOM renormalization with the same improvements. We also show that systematics are greatly reduced by the simultaneous inclusion of RGR and LRR and that these methods are necessary if improved precision is to be reached with higher-order terms in renormalization and matching.

Physics↗
DOE OSTI2022

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

On the origin of red spirals: does assembly bias play a role?

The formation of the red spirals is a puzzling issue in the standard picture of galaxy formation and evolution. Most studies attribute the colour of the red spirals to different environmental effects. We analyze a volume limited sample from the SDSS to study the roles of small-scale and large-scale environments on the colour of spiral galaxies. We compare the star formation rate, stellar age and stellar mass distributions of the red and blue spirals and find statistically significant differences between them at 99.9% confidence level. The red spirals inhabit significantly denser regions than the blue spirals, explaining some of the observed differences in their physical properties. However, the differences persist in all types of environments, indicating that the local density alone is not sufficient to explain the origin of the red spirals. Using an information theoretic framework, we find a small but non-zero mutual information between the colour of spiral galaxies and their large-scale environment that are statistically significant (99.9% confidence level) throughout the entire length scale probed. Such correlations between the colour and the large-scale environment of spiral galaxies may result from the assembly bias. Thus both the local environment and the assembly bias may play essential roles in forming the red spirals. The spiral galaxies may have different assembly history across all types of environments. We propose a picture where the differences in the assembly history may produce spiral galaxies with different cold gas content. Such a difference would make some spirals more susceptible to quenching. Finally, in all environments, the spirals with high cold gas content could delay the quenching and maintain a blue colour, whereas the spirals with low cold gas fractions would be easily quenched and become red.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Potential Constraints to Neutrino-Nucleus Interactions Based on Electron Scattering Data

A thorough understanding of neutrino-nucleus interactions physics is crucial to achieving precision goals in broader neutrino physics programs. The complexity of nuclei comprising the detectors and limited understanding of their weak response constitutes one of the biggest systematic uncertainties in neutrino experiments - both at intermediate energies affecting the short- and long-baseline neutrino programs as well as at lower energies affecting coherent scattering neutrino programs. While electron and neutrino interactions are different at the primary vertex, many underlying relevant physical processes in the nucleus are the same in both cases, and electron scattering data collected with precisely controlled kinematics, large statistics and high precision allows one to constrain nuclear properties and specific interaction processes. To this end, electron-nucleus scattering experiments provide vital complementary information to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. In fact, for many decades, the study of electron scattering off a nucleus has been used as a tool to probe the properties of that nucleus and its electromagnetic response. While previously existing electron scattering data provide important information, new and proposed measurements are tied closely to what is required for the neutrino program in terms of expanding kinematic reach, the addition of relevant nuclei and information on the final states hadronic system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Range-separated stochastic resolution of identity: Formulation and application to second-order Green’s function theory

We develop a range-separated stochastic resolution of identity (RS-SRI) approach for the four-index electron repulsion integrals, where the larger terms (above a predefined threshold) are treated using a deterministic RI and the remaining terms are treated using a SRI. The approach is implemented within a second-order Green’s function formalism with an improved O(N 3 ) scaling with the size of the basis set, N. Moreover, the RS approach greatly reduces the statistical error compared to the full stochastic version [T. Y. Takeshita et al., J. Chem. Phys. 151, 044114 (2019)], resulting in computational speedups of ground and excited state energies of nearly two orders of magnitude, as demonstrated for hydrogen dimer chains and water clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Indirectly probing the structural change in ion-irradiated Zr-Based metallic glasses from small scale mechanical tests

Ion irradiation was applied to tailor the structural heterogeneities in Zr-based metallic glasses at room temperature. Experimental methods of X-ray diffraction, nanoindentation, and micropillar compression were conducted to examine the irradiation effects on their structural and mechanical property changes. It is found that the irradiated materials retained amorphous structure after room-temperature Ni ion irradiation. The reduction of elastic modulus and hardness measured by nanoindentation indicated the irradiation-induced mechanical degradation. A unified statistic model was employed to quantitatively predict the density and strength of irradiation defects, although their specific structural and physical nature is not explicitly included in this model. The transition of non-intersecting shear bands to multiple intersecting shear bands was observed on compression tests of irradiated micropillars with the increase of irradiation dose. Furthermore, the analysis of the displacement excursion of micropillar compression tests indicated a different deformation mode from the unirradiated state. These results from nanoindentation pop-in and micro-pillar compression tests suggested that irradiation eventually leads to a new state with different types and characteristics of structural heterogeneities from violent displacement cascades and non-equilibrium energy deposition/dissipation processes, which also proves ion irradiation as an effective method to tune the structure and mechanical properties of metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Full event particle-level unfolding with variable-length latent variational diffusion

The measurements performed by particle physics experiments must account for the imperfect response of the detectors used to observe the interactions. One approach, unfolding, statistically adjusts the experimental data for detector effects. Recently, generative machine learning models have shown promise for performing unbinned unfolding in a high number of dimensions. However, all current generative approaches are limited to unfolding a fixed set of observables, making them unable to perform full-event unfolding in the variable dimensional environment of collider data. A novel modification to the variational latent diffusion model (VLD) approach to generative unfolding is presented, which allows for unfolding of high- and variable-dimensional feature spaces. The performance of this method is evaluated in the context of semi-leptonic t\bar{t} t t ‾ production at the Large Hadron Collider.

Shmakov, Alexander↗
DOE OSTI2020

End-Use Load Profiles for the U.S. Building Stock

The End-Use Load Profiles for the U.S. Building Stock project uses New ResStock and ComStock models to statistically represent the energy use of U.S. buildings. The project's hybrid approach combines best-available ground truth data, such as submetering studies and statistical disaggregation of whole-building interval meter data, with the reach, cost-effectiveness, and granularity of physics-based and data-driven building stock modeling to deliver a nationally-comprehensive data set at a fraction of the historical cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Simulation of Particle-Material Interactions

This paper gives an overview of the particle transport theory essentials, the basics of particle-material interaction simulation, physical quantities needed to simulate particle transport and interactions in materials, Monte Carlo simulation flow, response of additive detectors, statistical weights and other techniques to minimize statistical errors. Effects in materials under irradiation, materials response related to component lifetime and performance are considered with a focus on high-energy and high-power accelerator applications. Implementation of simulation of particle-material interactions in the modern Monte Carlo codes along with the code s main features and results of recent benchmarking are described.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Analysis of the multicomponent relativistic Boltzmann equation for electron scattering in big bang nucleosynthesis

Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multicomponent relativistic Boltzmann equation. However, the full multicomponent relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce nonthermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a nonthermal injection. Further, we show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Equations of State in Computational Physics: What, Why, and How

Summary: 1. Equations of state furnish thermodynamic properties of a material. 2. Usually involve P, p (or V), and T, or P, r, and E 3. Necessary ingredient to solve compressible flow problems. 4. Varieties: i. Analytical expressions (ideal gas, Grüneisen) – fast, convenient for verification and test problems ii. Graphs and charts – good for hand calculations iii. Tables and databases – generally most accurate, require interpolation i. Steam tables ii. NIST database iii. SESAME EOS Library 5. Constructed by using semi-empirical models from Statistical Mechanics and fitting to available data

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data Selection Improvement For MicroBooNE

Data selection is an extremely important part of data analysis for any experiment. Finding a physics result is often the result of sifting through a massive amount of data, keeping data that we believe to be signal, and throwing out data we do not. This process is called data selection. Creating a selection algorithm is an intensive process that must balance keeping enough data to have statistics and maximizing the signal purity of that data. We also need to choose the right reconstruction method, a tool to take raw data from the detector and convert it into physics results. In this study, we used three different reconstruction tools, Pandora, WireCell, and LANTERN, for the MicroBooNE experiment in conjunction to improve the selection algorithm for analysis. For the case of this study, we look into the charged current N proton 0 pions (CCNp0$\pi$) interaction channel. This is the dominant channel for the Short Baseline Neutrino (SBN) program and is expected to be a large contributor to the Deep Underground Neutrino Experiment (DUNE). We first investigated each of the three tools to find out more about their strengths and weaknesses as reconstructions, and compared them to the truth information directly from the MicroBooNE simulation pipeline. We then put together a direct comparison of the three methods to find which method or combination of methods would return the best result for us. While the study is ongoing, we have learned a lot about data selection for the experiment and the differences between the reconstruction tools.

Dillon, Brayden [Fermilab]↗
DOE OSTI2023

Methodology for physics-informed generation of synthetic neutron time-of-flight measurement data

Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

New physics in single resonant top quarks

Abstract Searches for new physics in the top quark sector are of great theoretical interest, yet some powerful avenues for discovery remain unexplored. We characterize the expected statistical power of the LHC dataset to constrain the single production of heavy top partnersTdecaying to a top quark and a photon or a top quark and a gluon. We describe an effective interaction which could generate such production, though the limits apply to a range of theoretical models. We find sensitivity to cross sections in the 10 2 − 10 5 fb range, forTmasses between 300 and 1000 GeV, depending on decay mode.

Physics↗
DOE OSTI2025

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗