Search NASA⌕ Search

SEARCH · Search NASA

Results for “GEANT4”

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

An Updated Simulation of the Booster Neutrino Beam

For accelerator neutrino experiments, an accurate prediction of the incoming neutrino flux is crucial for reducing uncertainties for all physics measurements. In this exciting period for the Short-Baseline Neutrino program at Fermilab, with far detector (ICARUS) already operating and the near detector (SBND) nearing operation, an updated flux model for the Booster Neutrino Beam (BNB) is presented. The BNB delivers 8 GeV protons to a beryllium target, subsequently producing neutrinos from the decay of the secondary beam of hadrons. A updated Monte-Carlo simulation of beam production in GEANT4 has been created, allowing predictions to be made for detectors with different baselines, offsets and sizes. This new simulation contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, with a view to storing all resulting particles - including neutral mesons - from the proton-Beryllium scatter to allow the study of exotic BSM scenarios. Results are presented, with comparisons to the original flux simulated for the MiniBooNE experiment.

Paton, Josephine↗
DOE OSTI2024

NA22 Quarterly Report

We continue to make progress toward determining the relative $\gamma$-ray intensities of 111 Ag from the secondary experimental configuration. This includes obtaining simulated detection efficiencies from CYLTRAN and Geant4 for $\gamma$-ray and β-$\gamma$ coincidence events, respectively.

42 ENGINEERING↗
DOE OSTI2024

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Simulating the CMS High Granularity Calorimeter with ML

Detector simulation is a key component of physics analysis and related activities in CMS. In the upcoming High Luminosity LHC era, simulation will be required to use a smaller fraction of computing in order to satisfy resource constraints. At the same time, CMS will be upgraded with the new High Granularity Calorimeter (HGCal), which requires significantly more resources to simulate than the existing CMS calorimeters. This computing challenge motivates the use of generative machine learning models as surrogates to replace full physics-based simulation. We study the application of state-of-the-art diffusion models to simulate particle showers in the CMS HGCal. We will discuss methods to overcome the challenges posed by the high-dimensional, irregular geometry of the HGCal. The quality of the showers produced by the diffusion model will be assessed by comparison to the full GEANT4-based simulation. The increase in simulation throughput will be quantified and methods to accelerate the diffusion model inference will also be discussed.

Amram, Oz↗
DOE OSTI2024

Refining Fast Simulation Using Machine Learning

A growing reliance on the fast Monte Carlo (FastSim) will accompany the high luminosity and detector granularity expected in Phase 2. FastSim is roughly 10 times faster than equivalent GEANT4-based full simulation (FullSim). However, reduced accuracy of the FastSim affects some analysis variables and collections. To improve its accuracy, FastSim is refined using regression-based neural networks trained with ML. The status of FastSim refinement is presented. The results show improved agreement with the FullSim output and an improvement in correlations among output observables and external parameters.

Güngördü, Acelya Deniz↗
DOE OSTI2024

Simulating Superconducting Devices in Novel Materials Using the G4CMP Simulation Toolkit

Understanding the kinematics of phonons and charge propagation in superconducting devices is crucial for conducting low-threshold dark matter searches and minimizing correlated errors in superconducting qubits. For nearly a decade, the Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, has been limited to simulating charge and phonon transport in silicon and germanium materials. In this work, we have expanded the capabilities of G4CMP to include novel substrate materials such as sapphire (Al2O3), Gallium Arsenide (GaAs), Lithium Fluoride (LiF), Calcium Tungstate (CaWO4), and Calcium Fluoride (CaF2). We demonstrate the use of this toolkit in generating phonon transport properties for these materials and compare the results with experimentally determined values where available. Additionally, we explore the performance of superconducting devices utilizing these materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Characterization of the Multi-Coincident Bismuth Germanate Oxide Spectrometer at the Nuclear Counting Facility

A multi-detector γ-ray spectroscopy system including 3 high-purity germanium (HPGe) detectors has been developed at the Nuclear Counting Facility at Lawrence Livermore National Laboratory. This system operates not only in conventional singles γ-ray detection mode but also incorporates features like coincidence and anti-coincidence techniques. These advanced detection modes can significantly enhance the system’s sensitivity, making it a powerful tool for isotope identification and quantification, especially for complex samples with high background and overlapping peaks. In this re port, we present a comprehensive performance characterization of the Multi-Coincident Bismuth Oxide Spectrometer (MCBOS) using a variety of sample types, including 60 Co, 48 Sc, and mixed fission product samples. Each sample provides unique tests for the MCBOS capabilities. Additionally, a GEANT4 model of the MCBOS has been generated to aid with the characterization of the system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Finite-Geometry Corrections for Neutron Scattering Cross Section Measurements

In this report, we describe the corrections necessary for neutron scattering cross section measurements in experiments which use finite-geometry scattering samples. A new Geant4 based framework was developed to calculate these corrections which has been benchmarked against a well-measured test case. This framework also calculates the outgoing neutron energy dependence of the correction factors, which is crucial when resolving individual final states in the scatterer is not possible.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Measurement of mean excitation energies of neutrino-relevant materials

Modern neutrino experiments require precision reconstruction of events. A crucial component of this reconstruction is the stopping power for charged particles, calculated using the Bethe equation. The main free parameter of the Bethe equation is the mean excitation energy (the "I-value"), which in most cases cannot be calculated, but must be measured for each substance. In many cases, the values are derived from very old experiments with large quoted uncertainties, or worse, small quoted uncertainties and inadequate treatment of systematics. Even if the tabulated values were reliable to the necessary degree, the I-value is affected by the phase of the substance and by chemical bonding, and only rough heuristics have been developed to convert measurements of single elements into I-values for compounds or from one phase to another. Modern neutrino experiments which need to measure an absolute energy scale, while lacking calibration sources, suffer the most from uncertainties in the I-value. DUNE is a primary example. We are performing measurements with the 400MeV Fermilab LINAC beam at the Irradiation Test Facility. We describe our efforts to measure the I-value of liquid argon, as well as those for several other materials used in past, present and potential future neutrino experiments, including water, NOvA scintillator, MINOS steel, zirconium, and molybdenum. These measurements use a set of degraders to scan the beam energy around the Bragg peak for each substance to determine the proton range. Experimental results are compared to results from Geant4 and FLUKA to determine the I-value.

Strait, Matthew L. [Fermilab] (ORCID:0000000157088↗
DOE OSTI2026

Induced Magnetic Fields in LBNF Horns and DUNE-PRISM Off-Axis Neutrinos

Deep Underground Neutrino Experiment (DUNE) is a next-generation, long-baseline neutrino oscillation experiment that will utilize an intense neutrino beam from Fermilab to measure neutrino oscillation parameters with unprecedented precision. The DUNE-PRISM near detector concept employs an off-axis measurement strategy to mitigate neutrino-nucleus interaction uncertainties; however, this approach relies critically on the accurate characterization of the neutrino flux. As statistical uncertainties are suppressed in the high-intensity Long Baseline Neutrino Facility (LBNF) beam, beamline focusing uncertainties, specifically those arising from focusing horn geometry, become an important systematic uncertainty source. This work presents Geant4 LBNF beam simulations quantifying one of the most important systematic uncertainties for off-axis fluxes: the impact of manufacturing tolerances in the LBNF horns. The focusing horns consist of coaxial inner and outer conductors. Ideally, the region inside the inner conductor is field-free, while the region between the inner and outer conductors serves as the focusing region, where the magnetic field follows a nominal $1/r$ dependence. However, inner conductor deformations such as eccentricity and ellipticity introduce asymmetry, inducing unintended magnetic fields inside the field-free region. Crucially, while the on-axis neutrino flux remains largely unaffected by these induced fields, the resulting flux fractional shifts become pronounced over a particular range of off-axis positions where the DUNE-PRISM program will perform measurements. Consequently, these inner conductor deformations can impact the precision of flux predictions, which could potentially degrade the sensitivity to oscillation parameters. Possible mitigation strategies for these effects will also be presented.

Bae, Yuri [U. Minnesota, Minneapolis (main)] (ORCI↗
DOE OSTI2026

Bayesian Exploration and Surrogate Emulation of Nonlinear Beam-Response Geometry in the LBNF Beamline

Next-generation long-baseline neutrino experiments aim to achieve multi-MW proton beam power while reducing accelerator-induced systematic uncertainties. At Fermilab, the LBNF beamline is designed for 1.2 MW operation with PIP-II and is upgradeable to 2.4 MW. DUNE will probe the three-flavor neutrino paradigm and search for CP violation, requiring precise neutrino-flux normalization and improved control of accelerator-related uncertainties. Within the LBNF beamline, the System for On-Axis Neutrino Detection (SAND) will constrain flux uncertainties using precision near-detector measurements, while the Muon Monitor System (MuMS) will provide beamline diagnostics sensitive to the proton beam, target, and horn configuration. However, the pion phase space relevant for DUNE depends simultaneously on many correlated parameters, including beam centroid, beam width, horn current and alignment, target position, optics shifts, and radiation-induced changes. Consequently, MuMS observables exhibit nonlinear and coupled responses that are difficult to characterize using traditional one-parameter scans. To address this challenge, we are developing a Bayesian Exploration framework coupled to physics-informed surrogate emulators trained on Geant4 beamline simulations. Gaussian-process emulators provide both fast predictions and uncertainty estimates, enabling adaptive selection of new simulation points in beam-parameter space. As an initial demonstration, we construct surrogate emulators for MuMS response observables using a verified simulation campaign spanning proton-beam steering conditions. The emulators reproduce the simulated dependence of MuMS centroid and gradient observables while providing predictive uncertainties, and serve as the foundation for future multidimensional exploration including beam width, horn current, and additional beamline parameters. This work establishes a framework for uncertainty-aware beam monitoring, adaptive simulation campaigns, and rapid beam-response inference for future DUNE operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗
DOE OSTI2022

Research and Development Studies for Reactor Neutrino Experiments in Turkey (RNET)

The program of the Reactor Neutrino Experiments of Turkey includes a small portable Water-based liquid scintillator detector to detect neutrinos from the Akkuyu nuclear power plant, planned to begin operating in 2023. The small near-field detector will weigh about 2-3 tons and will be placed less than 100 meters from the reactor cores. The Reactor Neutrino Experiments of Turkey program also includes a medium-size 30-ton Water-based liquid scintillator detector, which will be placed 1-2 km away from the reactor cores and will be used as a far detector. Both detectors and their responses to neutrino interactions were simulated using a GEANT4-based RAT-PAC simulation package. Here, we will share the technical and physical details of both detectors, and discuss the ongoing R&D effort for neutrino studies in Turkey.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Status of experimental knowledge on the unbound nucleus 13 Be

The structure of the unbound nucleus 13 Be is important for understanding the Borromean, two-neutron halo nucleus 14 Be. The experimental studies conducted over the last four decades are reviewed in the context of the beryllium chain of isotopes and some significant theoretical studies. The focus of this paper is the comparison of new data from a 12 Be(d,p) reaction in inverse kinematics, which was analyzed using Geant4 simulations and a Bayesian fitting procedure, with previous measurements. Two possible scenarios to explain the strength below 1 MeV above the neutron separation energy were proposed in that study: a single p-wave resonance or a mixture of an s-wave virtual state with a weaker p- or d-wave resonance. Comparisons of recent invariant mass and the (d,p) experiments show good agreement between the transfer measurement and the two most recent high-energy nucleon removal measurements.

12Be↗
DOE OSTI2022

The (Un)reasonable Effectiveness of Neural Network in Cherenkov Calorimetry

We report a greater than factor of two improvement in the hadronic energy resolution of a simulated Cherenkov calorimeter by estimating the energy with machine learning over traditional techniques. The prompt signal formation and energy threshold properties of Cherenkov radiation provide identifiable features that machine learning techniques can exploit to produce a superior model for energy reconstruction. We simulated a quartz-fiber calorimeter via the GEANT4 framework to study the reconstruction techniques in single events. We compared the machine learning-based reconstruction performance to the traditional simple sum of signal and dual-readout techniques that use both Cherenkov and scintillation signals. We describe why this game-changing approach to Cherenkov hadron calorimetry excels and our plans for a dedicated beam test to validate these findings with a fast, radiation-hard hadron calorimeter prototype.

Akchurin, Nural↗
DOE OSTI2022

The SiD Digital ECal Based on Monolithic Active Pixel Sensors

The SiD detector concept capitalizes on high granularity in its tracker and calorimeter to achieve the momentum resolution and particle flow calorimetry physics goals in a compact design. The collaboration has had a long interest in the potential for improved granularity in both the tracker and ECal with an application of monolithic active pixel sensors (MAPS) and a study of MAPS in the SiD ECal was described in the ILC TDR. Work is progressing on the MAPS application in an upgraded SiD design with a prototyping design effort for a common SiD tracker/ECal design based on stitched reticules to achieve 10 × 10 cm2 sensors with 25 × 100 micron2 pixels. Application of large area MAPS in these systems would limit delicate and expensive bump-bonding, provide possibilities for better timing, and should be significantly cheaper than the TDR concept due to being a more conventional CMOS foundry process. The small pixels significantly improve shower separation. Recent simulation studies confirm previous performance projections, indicating electromagnetic energy resolution based on digital hit cluster counting provides better performance than the SiD TDR analog design based on 13 mm2 pixels. Furthermore, the two shower separation is excellent down to the millimeter scale. Geant4 simulation results demonstrate these expectations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Simple Model of the Energy Threshold for Snowball Chambers

Cloud and bubble chambers have historically been used for particle detection, capitalizing on supersaturation and superheating, respectively. Here, we present new results from a prototype snowball chamber, in which an incoming particle triggers the crystallization of a purified, supercooled liquid. We demonstrate, for the first time, simulation agreement with our first results from 5 years ago: the higher temperature of the freezing of water and significantly shorter time spent supercooled compared to the control in the presence of a Cf-252 fission neutron source. This is accomplished by combining Geant4 modeling of neutron interactions with the Seitz nucleation model used in superheated bubble chambers, including those seeking dark matter. We explore the possible implications of using this new technology for GeV-scale WIMP searches, especially in terms of spin-dependent proton coupling, and report the first supercooling of WbLS (water-based liquid scintillator).

47 OTHER INSTRUMENTATION↗