Search NASA⌕ Search

SEARCH · Search NASA

Results for “Montecarlo”

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.

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exploring unsupervised top tagging using Bayesian inference

Recognizing hadronically decaying top-quark jets in a sample of jets, or even its total fraction in the sample, is an important step in many LHC searches for Standard Model and Beyond Standard Model physics as well. Although there exists outstanding top-tagger algorithms, their construction and their expected performance rely on Montecarlo simulations, which may induce potential biases. For these reasons we develop two simple unsupervised top-tagger algorithms based on performing Bayesian inference on a mixture model. In one of them we use as the observed variable a new geometrically-based observable \tilde{A}_{3} A ̃ 3 , and in the other we consider the more traditional \tau_{3}/\tau_{2} τ 3 / τ 2 N N -subjettiness ratio, which yields a better performance. As expected, we find that the unsupervised tagger performance is below existing supervised taggers, reaching expected Area Under Curve AUC \sim 0.80-0.81 ∼ 0.80 − 0.81 and accuracies of about 69% - − 75% in a full range of sample purity. However, these performances are more robust to possible biases in the Montecarlo that their supervised counterparts. Our findings are a step towards exploring and considering simpler and unbiased taggers.

Alvarez, Ezequiel↗
DOE OSTI2023

ECCE sensitivity studies for single hadron transverse single spin asymmetry measurements

Here, we performed feasibility studies for various single transverse spin measurements that are related to the Sivers effect, transversity and the tensor charge, and the Collins fragmentation function. The processes studied include semi-inclusive deep inelastic scattering (SIDIS) where single hadrons (pions and kaons) were detected in addition to the scattered DIS lepton. The data were obtained in pythia6 and geant4 simulated e+p collisions at 18 GeV on 275 GeV, 18 on 100, 10 on 100, and 5 on 41 that use the ECCE detector configuration. Typical DIS kinematics were selected, most notably Q 2 >1 GeV 2 , and cover the x range from 10 -4 to 1. The single spin asymmetries were extracted as a function of x and Q 2 , as well as the semi-inclusive variables z, which corresponds to the momentum fraction the detected hadron carries relative to the struck parton, and P T , which corresponds to the transverse momentum of the detected hadron relative to the virtual photon. They are obtained in azimuthal moments in combinations of the azimuthal angles of the hadron transverse momentum and transverse spin of the nucleon relative to the lepton scattering plane. In order to extract asymmetries, the initially unpolarized MonteCarlo was re-weighted in the true kinematic variables, hadron types and parton flavors based on global fits of fixed target SIDIS experiments and e + e – annihilation data. The expected statistical precision of such measurements is extrapolated to 10 fb –1 and potential systematic uncertainties are approximated given the deviations between true and reconstructed yields. Similar neutron information is obtained by comparing the ECCE e+p pseudo-data with the same from the EIC Yellow Report and scaling the corresponding Yellow Report e+ 3 He pseudo-data uncertainties accordingly. The impact on the knowledge of the Sivers functions, transversity and tensor charges, and the Collins function has then been evaluated in the same phenomenological extractions as in the Yellow Report. Finally, the impact is found to be comparable to that obtained with the parametrized Yellow Report detector and shows that the ECCE detector configuration can fulfill the physics goals on these quantities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

FLUKA simulation validation to measure the beam energy of a cyclotron

FLUKA is a fully integrated particle physics MonteCarlo simulation package which has many applications in high energy experimental physics and engineering. When using a cyclotron for medical radionuclide production, it is important to know that the beam energy is suitable to cause the intended nuclear reaction and minimize unwanted impurities. Here, this work focuses on the comparison of FLUKA as a tool to measure the beam energy of a cyclotron with experimental results using the method of stacked copper foils interspaced by an aluminum degrader. We compared two methods to determine the activity ratio: determining the activity in each foil using the experimentally measured cross-sections from the National Nuclear Data Center (NNDC) cross-section data and calculating the activity in the two Cu foils using the FLUKA program. These results were compared to those reported by Gagnon et al., 2011. The FLUKA results gave energy measurements that were different than the true energy by a percentage in the range of 2.4% to 13%, suggesting that FLUKA is not recommended for pre-calculating the cyclotron beam energy using these low Z materials without experimental confirmation. It also confirms that the experimental cross-section data from the NNDC database is superior for this application giving energy values that ranged from 0.1% to 1.6% from the true energy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Laser calibration system and lost muons correction in the g-2 experiment

The Muon g-2 experiment at Fermilab has the main goal to measure the muon anomalous magnetic moment $a_{µ}$ = ($g$ − 2)/2 to a precision of 0.14 parts per million (ppm), which means 4 times improvement in precision with respect to the final result from BNL: $a_{µ}$ (expt. BNL) = 11659208.0(6.3) × 10$^{−10}$ (0.54 ppm) In this report I summarize the work done at the $g$-2 experiment during my summer internship at Fermilab. In the first two weeks my first task has been to replace the NIM logic used in the Laser Calibration System with a new FPGA. Then, I was involved in the Lost Muons analysis studying both real data and MonteCarlo simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Search for Long-Lived Particles with Muon Detector Shower Signature in the CMS Run-3 data

Many beyond standard model theories predict the existence of long-lived particles (LLPs). These LLPs can have sizable lifetimes and decay several meters from their production vertex. In this poster/talk, we present the analysis strategy for searching for LLPs using the Compact Muon Solenoid (CMS) Experiment. The proton-proton collision data used in this analysis were collected from 2022 to 2024 at a center-of-mass energy of 13.6 TeV, corresponding to an integrated luminosity of 170 fb^-1. The LLP decays are reconstructed as a high-multiplicity cluster of detector hits in the cathode strip chambers (CSC) of the muon system endcap. This signature is referred to as the Muon Detector Showers (MDS). This search requires events to contain at least one MDS cluster. The analysis focuses on LLP hadronic decays and LLP masses up to a few tens of GeV. We present the signal properties in MonteCarlo simulation, event selection, background modeling, and evaluation of the expected sensitivity. The results are interpreted under the Twin Higgs model benchmark.

Agyemang-Duah, Andrews [Grambling State U.]↗
NASA NTRS2003

Infrared Spectroscopy of Star Formation in Galactic and Extragalactic Regions

In this program we proposed to perform a series of spectroscopic studies, including data analysis and modeling, of star formation regions using an ensemble of archival space-based data from the Infrared Space Observatory's Long Wavelength Spectrometer and Short Wavelength Spectrometer, and to take advantage of other spectroscopic databases including the first results from SIRTF. Our emphasis has been on star formation in external, bright IR galaxies, but other areas of research have included young, low or high mass pre-main sequence stars in star formation regions, and the galactic center. The OH lines in the far infrared were proposed as one key focus of this inquiry, because the Principal Investigator (H. Smith) had a full set of OH IR lines from IS0 observations. It was planned that during the proposed 2-1/2 year timeframe of the proposal other data (including perhaps from SIRTF) would become available, and we intended to be responsive to these and other such spectroscopic data sets. The program has the following goals: 1) Refine the data analysis of IS0 observations to obtain deeper and better SNR results on selected sources. The IS0 data itself underwent pipeline 10 reductions in early 2001, and the more 'hands-on data reduction packages' have been released. The IS0 Fabry-Perot database is particularly sensitive to noise and can have slight calibration errors, and improvements are anticipated. We plan to build on these deep analysis tools and contribute to their development. Model the atomic and molecular line shapes, in particular the OH lines, using revised montecarlo techniques developed by the Submillimeter Wave Astronomy Satellite (SWAS) team at the Center for Astrophysics. 2) 3) Use newly acquired space-based SIRTF or SOFIA spectroscopic data as they become available, and contribute to these observing programs as appropriate. 4) Attend scientific meetings and workshops. 5) E&PO activities, especially as related to infrared astrophysics and/or spectroscopy.

Smith, Howard A.↗
NASA NTRS2004

The Far Infrared Lines of OH as Molecular Cloud Diagnostics

Future IR missions should give some priority to high resolution spectroscopic observations of the set of far-IR transitions of OH. There are 15 far-IR lines arising between the lowest eight rotational levels of OH, and ISO detected nine of them. Furthermore, ISO found the OH lines, sometimes in emission and sometimes in absorption, in a wide variety of galactic and extragalactic objects ranging from AGB stars to molecular clouds to active galactic nuclei and ultra-luminous IR galaxies. The ISO/LWS Fabry-Perot resolved the 119 m doublet line in a few of the strong sources. This set of OH lines provides a uniquely important diagnostic for many reasons: the lines span a wide wavelength range (28.9 m to 163.2 m); the transitions have fast radiative rates; the abundance of the species is relatively high; the IR continuum plays an important role as a pump; the contribution from shocks is relatively minor; and, not least, the powerful centimeter-wave radiation from OH allows comparison with radio and VLBI datasets. The problem is that the large number of sensitive free parameters, and the large optical depths of the strongest lines, make modeling the full set a difficult job. The SWAS montecarlo radiative transfer code has been used to analyze the ISO/LWS spectra of a number of objects with good success, including in both the lines and the FIR continuum; the DUSTY radiative transfer code was used to insure a self-consistent continuum. Other far IR lines including those from H2O, CO, and [OI] are also in the code. The OH lines all show features which future FIR spectrometers should be able to resolve, and which will enable further refinements in the details of each cloud's structure. Some examples are given, including the case of S140, for which independent SWAS data found evidence for bulk flows.

Smith, Howard A.↗
NASA NTRSDate not supplied

Assimilation of microwave satellite observations over the rainbands of tropical cyclones

The potential of microwave observations in improving the weather forecasts is limited by the accuracy of all-sky radiative transfer calculations. We introduce a novel Bayesian MonteCarlo technique to improve the assimilation of microwave observations over the rainbands of tropical cyclones. The BMCI technique eliminates the need for a forward model in the data assimilation system. The technique includes three steps, (i) generating a comprehensive dataset using in-situ cloud measurements and atmospheric profiles, (ii) generating synthetic ATMS observations from the training dataset, and (iii) using real observations to estimate the geophysical variables over the rainbands of tropical cyclones. The retrieved profiles of temperature, relative humidity, and cloud water content as well assurface information such as SST were then assimilated into the model. The results show that assimilating the BMCI retrievals can influence the dynamical features of the cyclone, including a stronger warm core, a symmetric eye, and vertically aligned wind columns.

Monte Carlo↗
DOE OSTI2024

Software engineering to sustain a high-performance computing scientific application: QMCPACK

We provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum MonteCarlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion ofcontinuous integration (CI) targeting CPU, using GitHub Actions runners, and graphics processing units (GPU) in pre-exascalesystems, using self-hosted hardware; (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Dockercontainers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, sustainable maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the sustainability of community HPC codes and scientific discovery at scale.

Godoy, William↗
DOE OSTI2021

Realization of a high luminosity muon beam at the Jefferson Lab

The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.

Fulci, Antonino↗
DOE OSTI2024

ADRIANO2 Calorimeter performance from 2022 Prototypes

A novel high-granularity dual-readout calorimetric technique was developed as part of the T1604 collaboration research program. The building block of ADRIANO2 Calorimeter consists of a pair of optically isolated, small sized tiles made of scintillating plastic and lead glass. The two components of the dual-readout energy compensation technique. Furthermore,Cˇerenkovlightfromtheleadglassareexploitedtoforhighresolutiontimingmeasurements,whilehigh granularity from scintillating plastic can be used to probe the spatial component of the particle shower. This setup works for excellent energy resolution and particle identification for REDTOP as it is crucial for a calorimeter to detect and identify the decay products of eta/eta-prime mesons. Measurements were collected on several ADRIANO2 prototypes between February to December 2022 to evaluate the detector performance. The key metrics extracted from the analyses are presented in this paper are the efficiencies and light-yi eld for various tile configurations. The measured values will be used as input parameters for improved REDTOP Montecarlo simulation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Serpent - Bison - THM Preliminary Multiphysics Modeling of a Nuclear Thermal Propulsion System Fuel Assembly

This work demonstrates the Monte Carlo neutronic and Thermo-Hydraulic coupling scheme using the Serpent code and MOOSE application Bison and Thermo Hydraulic Module. The coupling scheme is then applied to the reference BWX Technologies Nuclear Ther- mal Propulsion system at he fuel assembly level where it’s used to perform an analysis of the isothermal material coefficients and potential material reactivity worth. A method is developed to isolate which feedback effects should be considered for proceeding with reduced order deterministic neutronic modeling where branch off analysis must be con- ducted. The convergence behavior of the coupling scheme is demonstrated where it fol- lows the standard Picard iteration approach. Verification studies for the method of deduc- ing relevant feedback effects is also demonstrated.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗