Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

RadSim: Particle Transport Wrapper

RadSim is being developed to provide the capability to simulate radiation source emissions, interpolate results from radiation transport tools into a common format to prepare incident flux, and model radiation detector response. This software package includes wrapper API code which is utilized in the development of the second (transport) task of the RadSim project. This software package only contains tools to interface RadSim with radiation transport tools (GEANT4 and MCNP) and does not include any part of the transport source code itself. The users will need prior access to the transport tools (GEANT4 and MCNP), and can use this software package to interface with the transport tool of their choice. This software package takes the geometry file for a radiation transport problem and converts it into an input deck for the transport tools (GEANT4 and MCNP). Additionally, it also contains filters that convert the output of the transport tools to a common flux format.

Cheung, Hoi Sing↗
DOE OSTI2018

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Simulating the SuperCDMS Dark Matter Detector Response and Readout

Dark matter is one of the remaining great mysteries of modern physics. Cosmological and astrophysical evidence hints strongly at a new particle, or set of particles, unknown to the Standard Model, comprising some 80% of all of the matter in the universe. SuperCDMS is one of the world-leading direct-detection experiments attempting to directly measure interactions between dark and regular matter. SuperCDMS uses solid state crystalline detectors in a nearly-zero background environment to attempt to measure the particle properties of dark matter. The detectors themselves, however, are quite complex; and a very detailed understanding of the microscopic physics is helpful in analyzing the very rare events that occur within them. Furthermore, better understanding and modeling of the detectors can aid in the design and optimization of future iterations of the experiment. The SuperCDMS detectors are well-understood in ideal cases. However, an accurate understanding of the detector response to outlier events or under non-ideal conditions is important for a low-statistics experiment. This work describes the design and implementation of a low temperature condensed matter physics simulation library built on top of the popular Geant4 particle tracking framework. The library, named \Geant4 Condensed Matter Physics" or G4CMP, introduces several solid state concepts to the Geant4 framework such as crystal lattices, phonon quasiparticles, non-scalar elective masses, and implements several physics processes relevant to cryogenic temperature crystals. In addition to the physics library, which is intended for general use, this work also describes a full Monte Carlo simulation package for the SuperCDMS iZIP detectors which utilizes G4CMP at its core and also fully simulates the detector readout sensors. The document will begin with the canonical motivations for the cold dark matter hypothesis, followed by a description of the SuperCDMS experiment. It will then describe the microscopic physics inside of the SuperCDMS detectors, as well as details on their implementation in G4CMP. The full detector simulation with corresponding sensor readout will follow, along with plans for future work; including both scientific studies and improvements, as well as computational improvements and added efficiencies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Towards Realistic and High Fidelity Models for Nuclear Reactor Power Synthesis Simulation with Self-Powered Neutron Detectors

As presented in this report, a weighting function–based inferencing method is being applied to synthesize the power distribution in next-generation and university research reactors based on simulated self power neutron detector (SPND) responses. The overall goal is to assess the impacts of sensor uncertainty and true power distribution perturbations on the error in the synthesized power distribution. Regarding sensor uncertainty, the NuScale Small Modular Reactor (SMR) and the Westinghouse AP1000 serve as testbeds for analyzing the impact of varying the sensor uncertainty, as well as varying the number of sensors per sensor string in the reactor core. The reactor models are informed by Monte Carlo N-Particle (MCNP) neutron flux tallies. For the NuScale SMR and Westinghouse AP1000, the SPND response functions (i.e., the response of the SPNDs to individual segments of fuel) were determined homogeneously. Regarding an analysis of power distribution perturbation detection, the Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor was used as a demonstration case with one particular arrangement of SPNDs; the response functions for this reactor model were determined heterogeneously, making this a uniquely high-fidelity demonstration of perturbation detection. Finally, SPND models generated in the Geometry and Tracking 4 (Geant4) code have been generated and tested for comparison with traditionally implemented analytical SPND models, with the intent for Geant4 integration with the full methodological framework. SPND current outputs as a function of distance from some fuel assembly segment in the NuScale SMR are compared with the analytically determined currents. Results from the sensor uncertainty simulations for the NuScale SMR and AP1000 indicate that the average error in the inferred power distribution on the fuel assembly segment level is reasonably low, being slightly less than the random uncertainty applied to all respective SPNDs in both cores. For example, if all SPNDs in the core have a random uncertainty of 5%, then the corresponding fuel assembly segment level error (i.e. difference between the true and inferred local power) is ~2–3%. However, the maximum error in the inferred power distribution on the fuel assembly segment level can be considerably high (>15%) when SPND random uncertainties start to exceed ~3%. In general, the average and maximum errors in the inferred power distribution were slightly higher in the AP1000 as opposed to the NuScale SMR for the sensor string configurations considered herein. Another result determined from analysis of the sensor uncertainty simulations was that increasing the number of SPNDs per string does not clearly reduce inferred power distribution error and can in fact make the error large in some cases; however, this assessment may skewed due to imposed iteration limits. Results from the perturbation detection demonstration using the high-fidelity TAMU TRIGA model indicate that, given the arrangement of 17 SPND strings and 4 SPNDs per string considered herein, there is a clear, provable ability to infer a localized Gaussian-type peak perturbation in the 3D power distribution. Such a perturbation was detected with an average fuel assembly segment level error of 0.19%, and the general visualization of the detected perturbation clearly indicates that the magnitude and shape were appropriately resolved. Finally, the electrical current output generated by the Geant4 modeled SPND indicates significant magnitude differences than the analytically modeled SPND, demonstrating the need for accurate SPND models which account for finite sensor geometry effects to inform the power synthesis work described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

EM Physics

Geant4 provides a comprehensive set of electromagnetic (EM) processes and models for electron/positron, gamma and long-lived charged particles, spanning energies from 100 eV to 100 TeV. Covering diverse energy regions often requires multiple models, which can be constructed using pre-packaged or user-defined EM physics constructors. Geant4 also supports detailed low-energy EM physics through models like Livermore, Penelope, and ICRU73, offering extensive data for elements across a wide energy range (250 eV–100 GeV). Application domains include space, medical, and radiobiology simulations. Additionally, Geant4 provides robust options for simulating complex optical photon production and transportation processes, enhancing its versatility in physics research and applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

D$e^+e^-$ffusion: Capturing the Beam-Beam Physics of $e^+e^-$ Collisions with Diffusion Models

Beam-induced backgrounds at high-luminosity $e^+e^-$ colliders, such as the FCC-ee, are dominated by incoherent pair creation (IPC), and require computationally expensive simulations with dedicated Monte Carlo (MC) event generators. Reliable detector and machine-detector interface studies necessitate event samples that are several orders of magnitude larger than what is practically attainable with existing MC. To alleviate this bottleneck, we present D$e^+e^-$ffusion, a denoising diffusion probabilistic model that operates as a permutation-equivariant, set-valued surrogate for fast IPC simulation. Trained on a small GuineaPig++ sample, D$e^+e^-$ffusion faithfully reproduces the marginal and joint kinematic, angular, and positional distributions of all three IPC production processes. In addition, we assess the fidelity at the detector level by propagating both Geant4 and D$e^+e^-$ffusion events through a Geant4 simulation of the CLD vertex detector and by training a transformer-based two-sample classifier; the classifier achieves an area under the ROC curve of $0.553 \pm 0.016$. The trained model generates events nearly four orders of magnitude faster than Geant4, paving the way for a fast-simulation surrogate for FCC-ee design studies.

Chahine, Antonio [Imperial Coll., London]↗
DOE OSTI2018

Muon Borehole Detector Design for Use in 4-D Density Overburden Monitoring

The increase of atmospheric CO2 concentrations has initiated research into carbon sequestration methods. One possibility is to store CO2 in subsurface porous reservoirs. The monitoring of the injected CO2 plume is vital because escaping CO2 poses health and environmental risks. Typically, seismic reflection methods are used to determine the change in density due to the replacement of brine by CO2 in the reservoir, but this is expensive and not continuous. A potential alternative is to use cosmic muon tomography to measure density changes in the reservoir as a function of time. This paper describes the development of a muon detector that will be capable of being deployed in boreholes. The detector has the required dimensions, an angular resolution of approximately 2 degrees, and is robust. The design is based on alternating layers of scintillating rods, which can provide 4-D tomographic data to detect small changes in density at depths up to approximately 2 km. Geant4, a Monte Carlo simulation code, is being used to develop models to guide the design of the physical configuration. Preliminary testing and measurements have been performed to validate the simulation predictions and optimize the performance of the system. The simulation and preliminary experimental results are presented here.

Borehole Detector, CO¬2 Sequestration, Geant4, Muo↗
DOE OSTI2025

Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions

The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

G EANT 4 atomic relaxation data for transfermium nuclei (Z = 101–104)

Advanced theoretical methods can accurately calculate various atomic observables and predict electronic structure. Still, systematic computations of the radiative and non-radiative transition probabilities and energies are missing for the actinides and all the transfermium elements. However, these compilations are needed for comprehensive Monte-Carlo simulations (such as GEANT4) of the radioactive decay of transfermium nuclei. These simulations can forma basis for data analysis of experiments, especially with complex detection setups. Investigation of the transfermium nuclei is crucial for understanding the nature of the nuclear force. In this study, simulations based on data from the Jena Atomic Calculator (JAC) and the data from the Evaluated Atomic Data Library (EADL) present in GEANT4 were found compatible for the three elements Ba(Z = 56), U(Z = 92), and Fm(Z = 100), thus, validating the JAC calculations. For Z> 100, we also found sound agreement between simulations that used data generated with JAC and experimental results involving No(Z = 102) and Rf(Z = 104) isotopes. In conclusion, these results demonstrate that JAC can produce reliable atomic data sets for transfermium elements, which will assist in analyzing nuclear-decay-spectroscopy experiments.

GEANT4↗
DOE OSTI2021

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Characterization of a gamma-ray tracking array: A comparison of GRETINA and Gammasphere using a 60 Co source

In this paper, we provide a formalism for the characterization of the tracking arrays with emphasis on the proper corrections required to extract their photopeak efficiencies and peak-to-total ratios. The methods are first applied to Gammasphere, a well characterized 4pi array based on the principle of Compton suppression, and subsequently to GRETINA. The tracking efficiencies are then discussed and some guidelines as to what clustering angle to use in the tracking algorithm are presented. It was possible, using GEANT4 simulations, to scale the measured efficiencies up to the expected values for the full 4pi implementation of GRETA.

efficiency measurements↗
DOE OSTI2025

Neutron capture measurements at UMass lowell research reactor

A new facility was designed around the thermal column beam port of the 1 MW Research Reactor at the University of Massachusetts Lowell. Thermal neutrons are collimated to a 1-inch diameter beam and incident on samples to induce the radiative neutron capture. New measurements of capture gamma rays are planned for Mn, Cu, Ni, Cr, and Gd samples in the next few years. The experiments will be carried out in close collaboration with the ENSDF (Evaluated Nuclear Structure Data File) evaluation group at Brookhaven National Laboratory. The gamma rays are measured using an array of high-resolution HPGe detectors. Some of the HPGe detectors use active Compton shields that were designed from scintillation detectors and improve significantly the signal-to-background in the measured gamma-ray spectra with HPGe. The experimental spectra will be validated by the GEANT4 simulations of the array and theoretical models of the emission of gamma rays from the compound nucleus. In conclusion, the new experimental results on the capture gamma ray intensities will be incorporated in future ENSDF evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

59 Co(p,X) spallation reaction cross sections for 250 MeV to 2 GeV protons

Cobalt is an advantageous target for probing the physics of nuclear spallation, because it is a naturally mono-isotopic element ( 59 Co), and its per-nucleon binding energy (BE/A = 8.768 MeV) is near the maximum value for all nuclei. We measured nuclear spallation cross sections for the 59 Co(p,X) reaction at five kinetic energies ranging from 250 MeV to 2 GeV. Cross sections for the production of 58 Co, 57 Co, 57 Mn, 56 Co, 56 Mn, 56 Cr, 55 Fe, 53 Fe, 54 Mn, 52 Mn, 51 Cr, 49 Cr, 48 V, 47 Sc, 46 Sc, 44 Sc, and 44 Scm are reported. Where comparable data exist in the EXFOR reaction database, we find that our measured cross sections generally agree. In many cases, we provide data for reactions or energies not currently reported in EXFOR. Our cross sections also provide evidence for the presence of α-clusters within the 59 Co nucleus, a surprising result given the asymmetry in Z (27) and N (32) for this nucleus. Finally, we use our measurements to evaluate the accuracy of spallation cross section simulations from GEANT4-based radiation transport toolkit, performed with the INCLXX-, Bertini-, and Binary-ion-cascade (BIC) based physics lists. This benchmarking activity revealed that the simulations overestimated the cross sections by a factor of ∼2–4 on average, and that the INCL-XX physics list provides the most reliable results. This evaluation informs the selection of the GEANT4 physics lists used for the analysis of data from NASA’s Psyche mission, which will measure γ rays and neutrons resulting from spallation reactions occurring on the surface of an asteroid whose surface is thought to be rich in iron-nickel metal.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Current and Future Performance of the CMS Simulation

The CMS full simulation using Geant4 has delivered billions of simulated events for analysis during Runs 1 and 2 of the LHC. However, the HL-LHC dataset will be an order of magnitude larger, with a similar increase in occupancy per event. In addition, the upgraded CMS detector will be considerably more complex, with an extended silicon tracker and a high granularity calorimeter in the endcap region. Increases in conventional computing resources are subject to both technological and budgetary limitations, so novel approaches are needed to improve software efficiency and to take advantage of new architectures and heterogeneous resources. Several projects are in development to address these needs, including the vectorized geometry library Vec-Geom and the GeantV transport engine, which uses track-level parallelization. The current computing performance of the CMS simulation will be presented as a baseline, along with an overview of the various optimizations already available for Geant4. Finally, the progress and outlook for integrating VecGeom and GeantV in the CMS software framework will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

R&D Adoption and Progress in Full Simulation of the CMS experiment

In this work we report on evolution of usage of Geant4 within CMSSW and adaptation of the newest Geant4 11.2.1, which is expected to be used for CMS simulation production in 2025. Physics validation results and results on CPU performance are reported. For the Phase-2 simulation several R&D are carried out. A significant update for CMS geometry description is performed. Different aspects of geometry description and physics simulation for the new detectors will be discussed. Progress on R&D efforts for the Phase-2 simulation will be presented, which includes reports on experience of application of G4HepEm external library.

Srimanobhas, Norraphat [Chulalongkorn U.]↗
DOE OSTI2013

A compact stilbene crystal neutron spectrometer for EAST D-D plasma neutron diagnostics

A new compact stilbene crystal neutron spectrometer has been investigated and applied in the neutron emission spectroscopy on the EAST tokamak. A new components analysis method is presented to study the anisotropic light output in the stilbene crystal detector. A Geant4 code was developed to simulate the neutron responses in the spectrometer. Based on both the optimal light output function and the fitted pulse height resolution function, a reliable neutron response matrix was obtained by Geant4 simulations and validated by 2.5 MeV and 14 MeV neutron measurements at a 4.5 MV Van de Graaff accelerator. The spectrometer was used to diagnose the ion temperature in plasma discharges with lower hybrid wave injection and ion cyclotron resonance heating on the EAST tokamak.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Review of Monte Carlo simulations for backgrounds from radioactivity

For all experiments dealing with the rare event searches (neutrino, dark matter, neutrino-less double-beta decay), the reduction of the radioactive background is one of the most important and difficult tasks. There are basically two types of background, electron recoils and nuclear recoils. The electron recoil background is mostly from the gamma rays through the radioactive decay. The nuclear recoil background is from neutrons from spontaneous fission, (α, n) reactions and muoninduced interactions (spallations, photo-nuclear and hadronic interaction). The external gammas and neutrons from the muons and laboratory environment, can be reduced by operating the detector at deep underground laboratories and by placing active or passive shield materials around the detector. The radioactivity of the detector materials also contributes to the background; in order to reduce it a careful screening campaign is mandatory to select highly radio-pure materials. In this review I present the status of current Monte Carlo simulations aimed to estimate and reproduce the background induced by gamma and neutron radioactivity of the materials and the shield of rare event search experiment. For the electromagnetic background a good level of agreement between the data and the MC simulation has been reached by the XENON100 and EDELWEISS experiments, using the GEANT4 toolkit. For the neutron background, a comparison between the yield of neutrons from spontaneous fission and (α, n) obtained with two dedicated softwares, SOURCES-4A and the one developed by Mei-Zhang-Hime, show a good overall agreement, with total yields within a factor 2 difference. The energy spectra from SOURCES-4A are in general smoother, while those from MZH presents sharp peaks. The neutron propagation through various materials has been studied with two MC codes, GEANT4 and MCNPX, showing a reasonably good agreement, inside 50% discrepancy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

A simulation of X-ray shielding for a superconducting electron cyclotron resonance ion source

It is generally assumed that large amounts of x-rays are emitted from the ion source of an Electron Cyclotron Resonance (ECR) instrument. The total amount of x-rays should be strictly limited to avoid the extra heat load to the cryostat of the superconducting ECR ion source, since they are partly absorbed by the cold mass into the cryostat. A simulation of x-ray shielding was carried out to determine the effective thickness of the x-ray shield needed via the use of Geant4. X-ray spectra of the 10 GHz Nanogan ECR ion source were measured as a function of the thickness variation in the x-ray shield. The experimental results were compared with Geant4 results to verify the effectiveness of the x-ray shield. Based on the validity in the case of the 10 GHz ECR ion source, the x-ray shielding results are presented by assuming the spectral temperature of the 28 GHz ECR ion source.

43 PARTICLE ACCELERATORS↗