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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ionoacoustic characterization of the proton Bragg peak with submillimeter accuracy

Purpose: Range verification in ion beam therapy relies to date on nuclear imaging techniques which require complex and costly detector systems. A different approach is the detection of thermoacoustic signals that are generated due to localized energy loss of ion beams in tissue (ionoacoustics). Aim of this work was to study experimentally the achievable position resolution of ionoacoustics under idealized conditions using high frequency ultrasonic transducers and a specifically selected probing beam. Methods: A water phantom was irradiated by a pulsed 20 MeV proton beam with varying pulse intensity and length. The acoustic signal of single proton pulses was measured by different PZT-based ultrasound detectors (3.5 and 10 MHz central frequencies). The proton dose distribution in water was calculated by Geant4 and used as input for simulation of the generated acoustic wave by the matlab toolbox k-WAVE. Results: In measurements from this study, a clear signal of the Bragg peak was observed for an energy deposition as low as 10{sup 12} eV. The signal amplitude showed a linear increase with particle number per pulse and thus, dose. Bragg peak position measurements were reproducible within ±30 μm and agreed with Geant4 simulations to better than 100 μm. The ionoacoustic signal pattern allowed for a detailed analysis of the Bragg peak and could be well reproduced by k-WAVE simulations. Conclusions: The authors have studied the ionoacoustic signal of the Bragg peak in experiments using a 20 MeV proton beam with its correspondingly localized energy deposition, demonstrating submillimeter position resolution and providing a deep insight in the correlation between the acoustic signal and Bragg peak shape. These results, together with earlier experiments and new simulations (including the results in this study) at higher energies, suggest ionoacoustics as a technique for range verification in particle therapy at locations, where the tumor can be localized by ultrasound imaging. This acoustic range verification approach could offer the possibility of combining anatomical ultrasound and Bragg peak imaging, but further studies are required for translation of these findings to clinical application.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Dosimetric characterization of two radium sources for retrospective dosimetry studies

Purpose: During the first part of the 20th century, {sup 226}Ra was the most used radionuclide for brachytherapy. Retrospective accurate dosimetry, coupled with patient follow up, is important for advancing knowledge on long-term radiation effects. The purpose of this work was to dosimetrically characterize two {sup 226}Ra sources, commonly used in Sweden during the first half of the 20th century, for retrospective dose–effect studies. Methods: An 8 mg {sup 226}Ra tube and a 10 mg {sup 226}Ra needle, used at Radiumhemmet (Karolinska University Hospital, Stockholm, Sweden), from 1925 to the 1960s, were modeled in two independent Monte Carlo (MC) radiation transport codes: GEANT4 and MCNP5. Absorbed dose and collision kerma around the two sources were obtained, from which the TG-43 parameters were derived for the secular equilibrium state. Furthermore, results from this dosimetric formalism were compared with results from a MC simulation with a superficial mould constituted by five needles inside a glass casing, placed over a water phantom, trying to mimic a typical clinical setup. Calculated absorbed doses using the TG-43 formalism were also compared with previously reported measurements and calculations based on the Sievert integral. Finally, the dose rate at large distances from a {sup 226}Ra point-like-source placed in the center of 1 m radius water sphere was calculated with GEANT4. Results: TG-43 parameters [including g{sub L}(r), F(r, θ), Λ, and s{sub K}] have been uploaded in spreadsheets as additional material, and the fitting parameters of a mathematical curve that provides the dose rate between 10 and 60 cm from the source have been provided. Results from TG-43 formalism are consistent within the treatment volume with those of a MC simulation of a typical clinical scenario. Comparisons with reported measurements made with thermoluminescent dosimeters show differences up to 13% along the transverse axis of the radium needle. It has been estimated that the uncertainty associated to the absorbed dose within the treatment volume is 10%–15%, whereas uncertainty of absorbed dose to distant organs is roughly 20%–25%. Conclusions: The results provided here facilitate retrospective dosimetry studies of {sup 226}Ra using modern treatment planning systems, which may be used to improve knowledge on long term radiation effects. It is surely important for the epidemiologic studies to be aware of the estimated uncertainty provided here before extracting their conclusions.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-T-221: Computational Evaluation of Proton Induced Gain in a Portable Faraday Cup

Purpose: The design of a new Portable Faraday Cup (PFC) used to calibrate proton accelerators was evaluated for energies between 50 and 220 MeV. Monte Carlo simulations performed in Geant4–10.0 were used to evaluate experimental results and reduce the relative detector error for this vacuum-less and low mass system, and invalidate current MCNP releases. Methods: The detector construction consisted of a copper conductor coated with an insulator and grounded with silver. Monte Carlo calculations in Geant4 were used to determine the net charge per proton input (gain) as a function of insulator thickness and beam energy. Kapton was chosen as the insulating material and was designed to capture backscattered electrons. Charge displacement from/into Kapton was assumed to follow a linear proportionality to the origin/terminus depth toward the outer ground layer. Kapton thicknesses ranged from 0 to 200 microns, proton energies were set to match empirical studies ranging from 70 to 250 MeV. Each setup was averaged over 1 million events using the FTFP-BERT 2.0 physics list. Results: With increasing proton energy, the gain of Cu+KA gradually converges to the limit of pure copper, with relative error between 1.52% and 0.72%. The Ag layer created a more diverging behavior, accelerating the flux of negative charge into the device and increasing relative error when compared to pure copper from 1.21% to 1.63%. Conclusion: Gain vs. beam energy signatures were acquired for each device. Further analysis reveals proportionality between insulator thickness and measured gain, albeit an inverse proportionality between beam energy and in-flux of electrons. Increased silver grounding layer thickness also decreases gain, though the relative error expands with beam energy, contrary to the Kapton layer.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-E-T-239: Design and Evaluation of a Nozzle Shielding System in Passively Scattered Proton Therapy

Purpose: To design, implement and evaluate a shielding system that will reduce out-of-field dose experienced by the patient and associated electronic systems in passively scattered proton therapy treatment. Methods: A multi-stage neutron shielding system was retrofitted to the Gantry 1 treatment nozzle at Loma Linda University Medical Center. The system uses multiple borated polyethylene plates staged after the primary beam modifying devices to attenuate and absorb neutrons produced by such devices. This arrangement locates increasing levels of shielding between the sources of secondary particles in the nozzle and the patient. Additionally, the design of this shielding structure allows it to be easily retrofitted to an existing proton nozzle system without impacting design or treatment beam characteristics. The effectiveness of the shielding was evaluated both through experimental measurements and Geant4 Monte Carlo simulations. Results: Measurements were completed with Landauer Luxel+ dosimeters that use optically stimulated luminescence and CR-39 to detect fast neutrons, thermal neutrons, protons, photons and beta particles. Measurements of a 250 MeV proton beam indicated that the shielding system reduced out-of-field dose to the patient by almost half with dose equivalent values at 50 and 40 cm from the field edge decreasing from 0.965 and 1.262 mSv/Gy to 0.596 and 0.777 mSv/Gy respectively. The installation of the multi-stage shielding system also reduced dose equivalent experienced by electronic systems installed in the treatment room by up to 80%. Geant4 simulations were also used to evaluate the neutron fluence at various positions in the treatment room as well as provide information on microdosimetry spectra within the patient and treatment room. Conclusion: The shielding system described above proved to be an effective an inexpensive method of reducing out-of-field doses to the patient and electronic systems and can be easily retrofitted to existing passive scattering nozzles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-AB-204-12: Dosimetry at the Sub-Cellular Scale of Auger-Electron Emitter 99m-Tc in a Mouse Single Thyroid Follicle Model

Purpose: To perform a dosimetry study at the sub-cellular scale of Auger-electron emitter 99m-Tc using a mouse single thyroid cellular model to investigate the contribution of the 99m-Tc Auger-electrons to the absorbed dose and possible link to the thyroid stunning in in vivo experiments in mice, recently reported in literature. Methods: The simulation of S-values for Auger-electron emitting radionuclides was performed using both the recent MCNP6 software and the Geant4-DNA extension of the Geant4 toolkit. The dosimetric calculations were validated through comparison with results from literature, using a simple model of a single cell consisting of two concentric spheres of unit density water and for six Auger-electron emitting radionuclides. Furthermore, the S-values were calculated using a single thyroid follicle model for uniformly distributed 123-I and 125-I radionuclides and compared with published S-values. After validation, the simulation of the S-values was performed for the 99m-Tc radionuclide within the several mouse thyroid follicle cellular compartments, considering the radiative and non-radiative transitions of the 99m-Tc radiation spectrum. Results: The calculated S-values using MCNP6 are in good agreement with the results from literature, validating its use for the 99m-Tc S-values calculations. The most significant absorbed dose corresponds to the case where the radionuclide is uniformly distributed in the follicular cell’s nucleus, with a S-value of 7.8 mGy/disintegration, due mainly to the absorbed Auger-electrons. The results show that, at a sub-cellular scale, the emitted X-rays and gamma particles do not contribute significantly to the absorbed dose. Conclusion: In this work, MCNP6 was validated for dosimetric studies at the sub-cellular scale. It was shown that the contribution of the Auger-electrons to the absorbed dose is important at this scale compared to the emitted photons’ contribution and can’t be neglected. The obtained S-values of Auger-electron emitting 99m-Tc radionuclide will be presented and discussed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Discontinuous finite element space-angle treatment of the first order linear Boltzmann transport equation with magnetic fields: Application to MRI-guided radiotherapy

Purpose: The advent of magnetic resonance imaging (MRI) guided radiotherapy systems demands the incorporation of the magnetic field into dose calculation algorithms of treatment planning systems. This is due to the fact that the Lorentz force of the magnetic field perturbs the path of the relativistic electrons, hence altering the dose deposited by them. Building on the previous work, the authors have developed a discontinuous finite element space-angle treatment of the linear Boltzmann transport equation to accurately account for the effects of magnetic fields on radiotherapy doses. Methods: The authors present a detailed description of their new formalism and compare its accuracy to GEANT4 Monte Carlo calculations for magnetic fields parallel and perpendicular to the radiation beam at field strengths of 0.5 and 3 T for an inhomogeneous 3D slab geometry phantom comprising water, bone, and air or lung. The accuracy of the authors’ new formalism was determined using a gamma analysis with a 2%/2 mm criterion. Results: Greater than 98.9% of all points analyzed passed the 2%/2 mm gamma criterion for the field strengths and orientations tested. The authors have benchmarked their new formalism against Monte Carlo in a challenging radiation transport problem with a high density material (bone) directly adjacent to a very low density material (dry air at STP) where the effects of the magnetic field dominate collisions. Conclusions: A discontinuous finite element space-angle approach has been proven to be an accurate method for solving the linear Boltzmann transport equation with magnetic fields for cases relevant to MRI guided radiotherapy. The authors have validated the accuracy of this novel technique against GEANT4, even in cases of strong magnetic field strengths and low density air.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Backscatter dose effects for high atomic number materials being irradiated in the presence of a magnetic field: A Monte Carlo study for the MRI linac

Purpose: To quantify and explain the backscatter dose effects for clinically relevant high atomic number materials being irradiated in the presence of a 1.5 T transverse magnetic field. Methods: Interface effects were investigated using Monte Carlo simulation techniques. We used GPUMCD (v5.1) and GEANT4 (v10.1) for this purpose. GPUMCD is a commercial software written for the Elekta AB, MRI linac. Dose was scored using GPUMCD in cubic voxels of side 1 and 0.5 mm, in two different virtual phantoms of dimensions 20 × 20 × 20 cm and 5 × 5 × 13.3 cm, respectively. A photon beam was generated from a point 143.5 cm away from the isocenter with energy distribution sampled from a histogram representing the true Elekta, MRI linac photon spectrum. A slab of variable thickness and position containing either bone, aluminum, titanium, stainless steel, or one of the two different dental filling materials was inserted as an inhomogeneity in the 20 × 20 × 20 cm phantom. The 5 × 5 × 13.3 cm phantom was used as a clinical test case in order to explain the dose perturbation effects for a head and neck cancer patient. The back scatter dose factor (BSDF) was defined as the ratio of the doses at a given depth with and without the presence of the inhomogeneity. Backscattered electron fluence was calculated at the inhomogeneity interface using GEANT4. A 1.5 T magnetic field was applied perpendicular to the direction of the beam in both phantoms, identical to the geometry in the Elekta MRI linac. Results: With the application of a 1.5 T magnetic field, all the BSDF’s were reduced by 12%–47%, compared to the no magnetic field case. The corresponding backscattered electron fluence at the interface was also reduced by 45%–64%. The reduction in the BSDF at the interface, due to the application of the magnetic field, is manifested in a different manner for each material. In the case of bone, the dose drops at the interface contrary to the expected increase when no magnetic field is applied. In the case of aluminum, the dose at the interface is the same with and without the presence of the aluminum. For all of the other materials the dose increases at the interface. Conclusions: The reduction in dose at the interface, in the presence of the magnetic field, is directly related to the reduction in backscattered electron fluence. This reduction occurs due to two different reasons. First, the electron spectrum hitting the interface is changed when the magnetic field is turned on, which results in changes in the electron scattering probability. Second, some electrons that have curved trajectories due to the presence of the magnetic field are absorbed by the higher density side of the interface and no longer contribute to the backscattered electron fluence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-217: A Comprehensive Monte-Carlo Study of Out-Of-Field Secondary Neutron Spectra in a Scanned-Beam Proton Therapy Treatment Room

Purpose: To simulate secondary neutron radiation-fields produced at different positions during phantom irradiation inside a scanning proton therapy gantry treatment room. Further, to identify origin, energy distribution and angular emission as function of proton beam energy. Methods: GEANT4 and FLUKA Monte-Carlo codes were used to model the relevant parts of the treatment room in a gantry-equipped pencil beam scanning proton therapy facility including walls, floor, metallic gantry-components, patient table and the homogeneous PMMA target. The proton beams were modeled based on experimental beam ranges in water and spot shapes in air. Neutron energy spectra were simulated at 0°, 45°, 90° and 135° relative to the beam axis at 2m distance from isocenter, as well as 11×11 cm2 fields for 75MeV, 140MeV, 200MeV and for 118MeV with 5cm PMMA range-shifter. The total neutron energy distribution was recorded for these four positions and proton energies. Additionally, the room-components generating secondary neutrons in the room and their contributions to the total spectrum were identified and quantified. Results: FLUKA and GEANT4 simulated neutron spectra showed good general agreement in the whole energy range of 10{sup −}9 to 10{sup 2} MeV. Comparison of measured spectra with the simulated contributions of the various room components helped to limit the complexity of the room model, by identifying the dominant contributions to the secondary neutron spectrum. The iron of the bending magnet and counterweight were identified as sources of secondary evaporation-neutrons, which were lacking in simplified room models. Conclusion: Thorough Monte-Carlo simulations have been performed to complement Bonner-sphere spectrometry measurements of secondary neutrons in a clinical proton therapy treatment room. Such calculations helped disentangling the origin of secondary neutrons and their dominant contributions to measured spectra, besides providing a useful validation of widely used Monte-Carlo packages in comparison to experimental data. Cluster of Excellence of the German Research Foundation (DFG) “Munich-Centre for Advanced Photonics (MAP)”.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-658: Out-Of-Field Dose Comparison for TrueBeam Low Energy Beams for Extended Distances: Measurement Vs Monte Carlo Simulation

Purpose: Patient dose far from the treatment field is comprised of scatter from within the patient, and treatment head leakage. We quantify the treatment head leakage for TrueBeam linear accelerator for 6X and 6X-FFF beams by comparing measurements to Monte Carlo simulations for a variety of jaw sizes and collimator rotations. This work is conceptually similar to that of Kry et al. (Medical Physics 2006; 33: 4405), who considered a Clinac linear accelerator. Methods: Measurements were made using an EXRADIN A101 ion chamber positioned in the patient plane, at distances up to 100 cm from isocenter. Simulations were done using VirtuaLinac, the GEANT4-based Monte Carlo model of the TrueBeam treatment head, and an in-house (U. Virginia) GEANT4-based code. In-house code modelled an ion chamber with build-up, based on a CT scan of the chamber. VirtuaLinac included a detailed model of the treatment head shielding, and was run on the Amazon Web Services cloud to generate spherical phase space files surrounding the treatment head. These phase space files were imported into the in-house code. Results: Initial comparisons between measurements and simulation showed an excess of dose in the in-plane direction, away from the gantry, in the simulations. This was traced to an incomplete model of the shielding—specifically, the component holding the primary collimator was smaller in the model than in the TrueBeam. Modifications were made to VirtuaLinac to more closely match the engineering drawings. In the in-plane direction, the lowest out of field dose was away from gantry (negative abscissa values) at around 60 cm from isocenter, for fields smaller than 10×10 cm2. Out of field dose decreased with decreasing jaw size. Flattening-filter free beam produced out-of-field doses as low as 65% of those with flattened beam. Conclusion: Doses determined from simulation and measurement were in close agreement. Funding support from the Jefferson Trust Foundation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-IeP2-13: Toward Heavy Ion Computed Tomography with Carbon Ions: A Monte Carlo Study

Purpose: In the present Monte Carlo study, we investigated the use of Carbon ions for computed tomography (CT) with a relatively low imaging dose. This will enable us to avoid any conversion of X-ray CT numbers to the relative stopping power (or relative electron density) values and the associated uncertainties in Carbon dose calculation. Methods: In the first stage, we studied the propagation of Carbon nuclei through a water phantom using the Geant4 specially to understand their lateral displacement inside the phantom. In the second stage, we used our GPU-based Monte Carlo code, which has been cross validated against Geant4, to create the 2D map of the water equivalent path length (WEPL) inside a human head size phantom by acquiring 240 projections each 1.5° apart. Subsequently the 3D relative electron density map of the phantom was reconstructed from the 2D WEPL map using the Algebraic Reconstruction Technique (ART) coupled with total variation (TV) minimization Results: A high quality image of the relative electron density inside the phantom was reconstructed by ARTTV. The mean relative error between the reconstructed image for low contrast object (PMMA) was about 1.74%. The delivered dose per scan at the center of the phantom was about 0.1 Gy. Conclusion: We have been able to obtain a 3D map of the electron density using a human head size phantom while keeping the delivered dose to relatively low value. Using the GPU capabilities of our simulation engine, we believe that a real time CT with Carbon ions could be a reality in future.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-BRB-04: A New Optimization Method for Pre-Treatment Patient-Specific Stopping-Power by Combining Proton Radiography and X-Ray CT

Purpose: The relative stopping power (RSP) uncertainty is the largest contributor to the range uncertainty in proton therapy. The purpose of this work is to develop a robust and systematic method that yields accurate patient specific RSPs by combining pre-treatment X-ray CT and daily proton radiography. Methods: The method is formulated as a penalized least squares optimization (PLSO) problem min(|Ax-B|). The matrix A represents the cumulative path-length crossed in each material computed by calculating proton trajectories through the X-ray CT. The material RSPs are denoted by x and B is the pRad, expressed as water equivalent thickness. The equation is solved using a convex-conic optimizer. Geant4 simulations were made to assess the feasibility of the method. RSP extracted from the Geant4 materials were used as a reference and the clinical HU-RSP curve as a comparison. The PLSO was first tested on a Gammex RMI-467 phantom. Then, anthropomorphic phantoms of the head, pelvis and lung were studied and resulting RSPs were evaluated. A pencil beam was generated in each phantom to evaluate the proton range accuracy achievable by using the optimized RSPs. Finally, experimental data of a pediatric head phantom (CIRS) were acquired using a recently completed experimental pCT scanner. Results: Numerical simulations showed precise RSP (<0.75%) for Gammex materials except low-density lung (LN-300) (1.2%). Accurate RSP have been obtained for the head (µ=−0.10%, 1.5σ=1.12%), lung (µ=−0.33, 1.5σ=1.02%) and pelvis anthropomorphic phantoms (µ=0.12, 1.5σ=0,99%). The range precision has been improved with an average R80 difference to the reference (µ±1.5σ) of −0.20±0.35%, −0.47±0.92% and −0.06±0.17% in the head, lung and pelvis phantoms respectively, compared to the 3.5% clinical margin. Experimental HU-RSP curve have been produced on the CIRS pediatric head. Conclusion: The proposed PLSO with prior knowledge X-ray CT shows promising potential (R80 σ<1.0% over all sites) to decrease the range uncertainty.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-CAMPUS-TeP3-05: Evaluation of the Microscopic Dose Enhancement in the Nanoparticle-Enhanced Auger Therapy

Purpose: The aim of this study is to apply Monte Carlo simulations to investigate the nanoparticle dose enhancement for Auger therapy. Methods: Two nanoparticle fabrications were considered: nanoshell and nanosphere. In the first step, a single nanoparticle was irradiated with Auger emitters. The electrons were scored in a phase space at the outer surface of the nanoparticle with Geant4-Penelope. In the second step, the previously recorded phase space was used as a source and placed at the center of a cell-size water phantom. The nanoscale dose was evaluated in water around the nanoparticle with Geant4-DNA. The dose enhancement factor (DEF) is defined as the ratio of doses with and without nanoparticles. The nanoparticles were replaced by corresponding water nanoparticle with the same size and volume source which represents typical situation of Auger emitters without nanoparticle. Various sizes/materials of nanoparticles and isotopes were considered. Results: Nanoshell - Microscopic dose was increased up to 130% at 20 – 100 nm distances from the surface of Au-{sup 125}I nanoshell. However, dose at less than 20 nm distance was reduced due to absorbed low energy electrons in gold nanoshell. The amounts and regions of the dose enhancement were dependent on nanoshell size, materials, and isotopes. Nanosphere - The increased amounts of electrons up to 300% and reduced average energy with nanosphere were observed compared with water nanoparticle. We observed localized dose enhancement (up to a factor 3.6) in the immediate vicinity (< 50 nm) of Au-{sup 125} I nanosphere. The dose enhancement patterns vary according to nanosphere sizes and isotopes. Conclusion: We conclude that Auger therapy with nanoparticles can lead to change of electron energy spectrum and dose enhancements at certain range. The dose enhancement patterns vary according to nanoparticle sizes, materials, and isotopes. This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP: Ministry of Science, ICT and Future Planning) (No. NRF-2013M2B2B1075776)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-H-BRA-04: Biological Geometries for the Monte Carlo Simulation Toolkit TOPASNBio

Purpose: New advances in radiation therapy are most likely to come from the complex interface of physics, chemistry and biology. Computational simulations offer a powerful tool for quantitatively investigating radiation interactions with biological tissue and can thus help bridge the gap between physics and biology. The aim of TOPAS-nBio is to provide a comprehensive tool to generate advanced radiobiology simulations. Methods: TOPAS wraps and extends the Geant4 Monte Carlo (MC) simulation toolkit. TOPAS-nBio is an extension to TOPAS which utilizes the physics processes in Geant4-DNA to model biological damage from very low energy secondary electrons. Specialized cell, organelle and molecular geometries were designed for the toolkit. Results: TOPAS-nBio gives the user the capability of simulating biological geometries, ranging from the micron-scale (e.g. cells and organelles) to complex nano-scale geometries (e.g. DNA and proteins). The user interacts with TOPAS-nBio through easy-to-use input parameter files. For example, in a simple cell simulation the user can specify the cell type and size as well as the type, number and size of included organelles. For more detailed nuclear simulations, the user can specify chromosome territories containing chromatin fiber loops, the later comprised of nucleosomes on a double helix. The chromatin fibers can be arranged in simple rigid geometries or within factual globules, mimicking realistic chromosome territories. TOPAS-nBio also provides users with the capability of reading protein data bank 3D structural files to simulate radiation damage to proteins or nucleic acids e.g. histones or RNA. TOPAS-nBio has been validated by comparing results to other track structure simulation software and published experimental measurements. Conclusion: TOPAS-nBio provides users with a comprehensive MC simulation tool for radiobiological simulations, giving users without advanced programming skills the ability to design and run complex simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Neutron production by cosmic-ray muons in various materials

The results obtained by studying the background of neutrons produced by cosmic-raymuons in underground experimental facilities intended for rare-event searches and in surrounding rock are presented. The types of this rock may include granite, sedimentary rock, gypsum, and rock salt. Neutron production and transfer were simulated using the Geant4 and SHIELD transport codes. These codes were tuned via a comparison of the results of calculations with experimental data—in particular, with data of the Artemovsk research station of the Institute for Nuclear Research (INR, Moscow, Russia)—as well as via an intercomparison of results of calculations with the Geant4 and SHIELD codes. It turns out that the atomic-number dependence of the production and yield of neutrons has an irregular character and does not allow a description in terms of a universal function of the atomic number. The parameters of this dependence are different for two groups of nuclei—nuclei consisting of alpha particles and all of the remaining nuclei. Moreover, there are manifest exceptions from a power-law dependence—for example, argon. This may entail important consequences both for the existing underground experimental facilities and for those under construction. Investigation of cosmic-ray-induced neutron production in various materials is of paramount importance for the interpretation of experiments conducted at large depths under the Earth’s surface.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Neutron flux measurement using activated radioactive isotopes at the Baksan underground scintillation telescope

Preliminary results of a neutron background measurement at the Baksan underground scintillation telescope (BUST) are presented. The external planes of the BUST are fully covered with standard scintillation detectors shielding the internal planes and suppressing thus background events due to cosmogenic and local radioactivity. The shielded internal planes were used as target for the neutron flux registration. The experimental method is based on the delayed coincidences between signals from any of the BUST counters. It is assumed that the first signal is due to inelastic interaction of a neutron with the organic scintillator, while the second signal comes from the decay of an unstable radioactive isotope formed when the fast neutron interacts with the {sup 12}C nuclei. Using the Monte-Carlo method (GEANT4) we also simulated propagation of neutrons through a layer of scintillator. The experimentally found muon induced neutron flux is j =1.3{sup +0.7}{sub -0.3} ×10{sup -10}cm{sup -2}s{sup -1} for neutron energies E ≥ 22MeV, which is in a qualitative agreement with similar measurements of other underground laboratories as well as with predictions of the GEANT4.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

G4VG

G4VG reads the in-memory Geant4 geometry (the volume hierarchy and physical shape descriptions defined by a user application) and produces an equivalent VecGeom representation. This provides a direct bridge between Geant4 detector descriptions and VecGeom-based applications, eliminating the need for manual export pipelines or specialized GDML workflows. The code is adapted from the Celeritas project and developed as an independent utility for broader reuse.

Johnson, SethR [Oak Ridge National Laboratory] (00↗
DOE OSTI2018

Measurement of the Negative Pion and Positive Kaon Total Hadronic Cross Sections on Argon at the LArIAT Experiment

The Liquid Argon Time Projection Chamber (LArTPC) represents one of the most advanced experimental technologies for physics at the Intensity Frontier due to its full 3D-imaging, excellent particle identication and precise calorimetric energy reconstruction. By deploying a LArTPC in a dedicated calibration test beam line at Fermilab, LArIAT (Liquid Argon In A Testbeam) aims to experimentally calibrate this technology in a controlled environment and to provide physics results key to the neutrino oscillation physics and proton decay searches of the Short Baseline Neutrino (SBND, MicroBooNE, ICARUS) and Long Baseline Neutrino programs (DUNE). LArIAT's physics program entails a vast set of topics with a particular focus on the study of nuclear eects such as pion and kaon characteristic interaction modes. This thesis presents two measurements performed for the rst time: the measure- ment of the negative pion total hadronic cross section on argon in the 100-1050 MeV kinetic energy range and the measurement of the positive kaon total hadronic cross section on argon in the 100-650 MeV kinetic energy range. The analyses devised for these measurements use both the core elements of LArIAT: beamline and LArTPC. The rst step in each analysis is the development of an event selection based on beam- line and TPC information geared towards the identication of the hadron of interest. We then proceed to match the beamline candidate to a suitable TPC track. Finally, we apply the \thin slice method" technique and measure the cross section, correcting for background contamination and detector eects. The thin slice technique is a new method to measure hadron-argon cross sections possible only due to the combination of the tracking and calorimetry capability of the LArTPC technology. Albeit never on argon, the hadronic cross section of pions has been extensively measured before on lighter and heavier dierent elements in thin target experiments, leading to solid predictions for measurements on argon. Through the use of a dierent technique, our measurement of the negative pion total hadronic cross section on argon is in general agreement with the predictions by the Geant4 Bertini Cascade model which are based on data from thin target experiments. On the contrary, cross section measurements for kaons are extremely scarce, thus more dicult to model. Our measurement of the positive kaon total hadronic cross section on argon is mostly in tension with the Geant4 prediction over the explored range of energies and provides new experimental data to properly compare existing interaction models. This thesis also reports two ancillary detector physics measurements necessary for the cross section analyses: the measurements of the LArIAT electric eld and calorimetry constants. We developed a technique to measure the LArIAT electric eld using cathode-anode piercing tracks with cosmic data. We applied a new technique for the measurement of the calorimetry calibration constants based on the particles' momentum measurement. The negative pion and the positive kaon total hadronic cross measurements are the rst physics results of the LArIAT experiment and will be the basis for the future LArIAT measurements of pion and kaon cross sections in the exclusive channels. The outcome of these measurements will ultimately enable to quantify and reduce the systematic associated with the hadronic interaction models in neutrino-argon interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

AI Denoising to Accelerate Detector Simulation

Detector simulation is critical to experimental HEP; however this simulation (commonly done through toolkits such as Geant4) is computationally intensive. Performance can be improved somewhat through technical optimization, but more is needed. Using machine learning (ML) to accelerate simulation is a promising field, however efforts to use generative adversarial networks (GANs) or optimized autoencoders have faced issues. Using convolutional neural networks (CNNs) for denoising has been successful in non-HEP applications such as image processing. This poster investigates the efficacy of using CNNs to denoise Geant4 simulations. This could increase the accuracy of simulations performed under settings designed to increase computational efficiency.

Franklin, Lena↗