Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

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

TU-CD-304-06: Using FFF Beams Improves Tumor Control in Radiotherapy of Lung Cancers

Purpose: Electron disequilibrium at the lung-tumor interface results in an under-dosage of tumor regions close to its surface. This under-dosage is known to be significant and can compromise tumor control. Previous studies have shown that in FFF beams, disequilibrium effects are less pronounced, which is manifested in an increased skin dose. In this study we investigate the improvement in tumor dose coverage that can be achieved with FFF beams. The significance of this improvement is evaluated by comparing tumor control probabilities of FFF beams and conventional flattened beams. Methods: The dosimetric coverage was investigated in a virtual phantom representing the chest wall, lung tissue and the tumor. A range of tumor sizes was investigated, and two tumor locations – central and adjacent to the chest wall. Calculations were performed with BEAMnrc Monte Carlo code. Parallel-opposed and multiple coplanar 6-MV beams were simulated. The tumor control probabilities were calculated using the logistic model with parameters derived from clinical data for non-small lung cancer patients. Results: FFF beams were not entirely immune to disequilibrium effects. They nevertheless consistently delivered more uniform dose distribution throughout the volume of the tumor, and eliminated up to ∼15% of under-dosage in the most affected by disequilibrium 1-mm thick surface region of the tumor. A voxel-by-voxel comparison of tumor control probabilities between FFF and conventional flattened beams showed an advantage of FFF beams that, depending on the set up, was from a few to ∼9 percent. Conclusion: A modest improvement in tumor control probability on the order of a few percent can be achieved by replacing conventional flattened beams with FFF beams. However, given the large number of lung cancer patients undergoing radiotherapy, these few percent can potentially prevent local tumor recurrence for a significant number of patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Qualification of JEFF3.1.1 library for high conversion reactor calculations using the ERASME/R experiment

With its low CO{sub 2} production, Nuclear Energy appears to be an efficient solution to the global warming due to green-house effect. However, current LWR reactors are poor uranium users and, pending the development of Fast Neutron Reactors, alternative concepts of PWR with higher conversion ratio (HCPWR) are being studied again at CEA, first studies dating from the middle 80's. In these French designs, low moderation ratio has been performed by tightening the lattice pitch, achieving a conversion ratio of 0.8-0.9 with a MOX fuel coming from PWR UOX recycling. Theses HCPWRs are characterized by a harder neutron spectrum and the calculation uncertainties on the fundamental neutronics parameters are increased by a factor 3 regarding a standard PWR lattice, due to the major contribution of the Plutonium isotopes and of the epithermal energy range to the reaction rates. In order to reduce these uncertainties, a 3-year experimental validation program called ERASME has been performed by CEA from 1984 to 1986 in the EOLE reactor. Monte Carlo analysis of the ERASME/R experiments with the Monte Carlo code TRIPOLI4 allowed the qualification of the recommended JEFF.3.1.1 library for major neutronics parameters. K{sub eff} of the MOX under-moderated lattice is over-predicted by 440 {+-} 830 pcm (2{sigma}); the conversion ratio, indicator of the good use of uranium, is also slightly over-predicted: 2 % {+-} 4 % (2{sigma}) and the same for B4C absorber rods worth and soluble boron worth, over-predicted by 2 %, both in the 2 standard deviations range. The radial fission maps of heterogeneities (water-holes, B4C and fertile rods) are well reproduced: maximal (C-E)/E dispersion is 1.3 %, maximal power peak error is 2.7 %. The void reactivity worth is the only parameter poorly calculated with an overprediction of +12.4% {+-} 1.5%. ERASME/R analysis of MOX reactivity, void effect and spectral indexes will contribute to the reevaluation of {sup 241}Am and Plutonium isotopes nuclear data for the next library JEFF3.2. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Assessment of the mechanical performance of the Westinghouse BWR control rod CR 99 at high depletion levels

A long-term program assessing the mechanical performance of the Westinghouse BWR control rod CR 99 at high depletion levels has been performed. The scope of the program has mainly been based on the operation of four CR 99 Generation 2 control rods in demanding positions during 6 and 7 cycles in the Leibstadt Nuclear Power Plant (KKL) and on the detailed visual inspections and blade wing thickness measurements that were performed after the rods were discharged. By correlating statistically the blade wing thickness measurements to the appearance of irradiation-assisted stress corrosion cracking (IASCC), the probability of IASCC appearance as function of the blade wing swelling was estimated. In order to correlate the IASCC probability of a CR 99 to its depletion, the {sup 10}B depletion of the studied rods was calculated in detail on a local level with the stochastic Monte Carlo code MCNP in combination with the Westinghouse nodal code system PHOENIX4/POLCA7. Using this information coupled to the blade wing measurement data, a finite element model describing the blade wing swelling of an arbitrary CR 99 design as function of {sup 10}B depletion could then be generated. In the final step, these relationships were used to quantify the probability of IASCC appearance as function of the {sup 10}B depletion of the CR 99 Generations 2 and 3. Applying this detailed mapping of the CR 99 behavior at high depletion levels and using an on-line core monitoring system with explicit {sup 10}B depletion tracking capabilities will enable a reliable prediction of the probability for IASCC appearance, thus enhancing the optimized design and the sound operation of the CR 99 control rod. Another important outcome of the program was that it was clearly shown that no significant amount of boron leakage did occur through any of the detected IASCC cracks, despite the very high depletion levels achieved. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

AN EMBEDDED ACTIVE NUCLEUS IN THE OH MEGAMASER GALAXY IRAS16399–0937

We present a multiwavelength study of the OH megamaser galaxy IRAS16399–0937, based on new Hubble Space Telescope (HST)/Advanced Camera for Surveys F814W and Hα+[N II] images and archive data from HST, Two Micron All Sky Survey, Spitzer, Herschel and the Very Large Array. This system has a double nucleus, whose northern (IRAS16399N) and southern (IRAS16399S) components have a projected separation of ∼6'' (3.4 kpc) and have previously been identified based on optical spectra as a low ionization nuclear emission line region (LINER) and starburst nucleus, respectively. The nuclei are embedded in a tidally distorted common envelope, in which star formation is mostly heavily obscured. The infrared spectrum is dominated by strong polycyclic aromatic hydrocarbon, but deep silicate and molecular absorption features are also present, and are strongest in the IRAS16399N nucleus. The 0.435-500 μm spectral energy distribution was fitted with a model including stellar, interstellar medium and active galactic nucleus (AGN) torus components using our new Markov Chain Monte Carlo code, CLUMPYDREAM. The results indicate that the IRAS16399N contains an AGN (L {sub bol} ∼ 10{sup 44} erg s{sup –1}) deeply embedded in a quasi-spherical distribution of optically thick clumps with a covering fraction ≈1. We suggest that these clumps are the source of the OHM emission in IRAS16399–0937. The high torus covering fraction precludes AGN photoionization as the origin of the LINER spectrum, however, the spectrum is consistent with shocks (v ∼ 100-200 km s{sup –1}). We infer that the ∼10{sup 8} M {sub ☉} black hole in IRAS16399N is accreting at a small fraction (∼1%) of its Eddington rate. The low accretion rate and modest nuclear star formation rates suggest that while the gas-rich major merger forming the IRAS16399–0937 system has triggered widespread star formation, the massive gas inflows expected from merger simulations have not yet fully developed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Feasibility study of noise analysis methods on virtual thermal reactor subcriticality monitoring

This paper presents the analysis results of Rossi-alpha, cross-correlation, Feynman-alpha, and Feynman difference methods applied to the subcriticality monitoring of nuclear reactors. A thermal spectrum Godiva model has been designed for the analysis of the four methods. This Godiva geometry consists of a spherical core containing the isotopes of H-l, U-235 and U-238, and the H{sub 2}O reflector outside the core. A Monte Carlo code, McCARD, is used in real time mode to generate virtual detector signals to analyze the feasibility of the four methods. The analysis results indicate that the four methods can be used with high accuracy for the continuous monitoring of subcriticality. In addition to that, in order to analyze the impact of the random noise contamination on the accuracy of the noise analysis, the McCARD-generated signals are contaminated with arbitrary noise. It is noticed that, even when the detector signals are contaminated, the four methods can predict the subcriticality with reasonable accuracy. Nonetheless, in order to reduce the adverse impact of the random noise, eight detector signals, rather than a single signal, are generated from the core, one signal from each equally divided eighth part of the core. The preliminary analysis with multiple virtual detector signals indicates that the approach of using many detectors is promising to improve the accuracy of criticality prediction and further study will be performed in this regard. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

SU-C-BRB-05: Investigation of Conformal Arc Therapy Utilizing Newly Designed Cobalt 60 Machine

Purpose: A new modernized design for cobalt 60 (Co-60) machines is being developed with a ring type gantry. In this study we investigate the beneficial outcome of the new design for conformal arc therapy for various clinical sites. The new modality was evaluated based on isodose distributions and dose volume histograms as compared to 6MV photon beams from conventional linear accelerators. Methods: Computed tomographic images of seven different SBRT patients were selected from our patient database. All of these cases were previously planned on the Eclipse treatment planning system. New plans for these patients were generated with a modified conformal arc technique using both 6MV and Co-60 beams. The conformal arc was created by the delivery of treatment fields conformal to the target cross-section at every 5 or 10 degrees. The field shape was modified or turned off when it initially passed through a critical structure. Monte Carlo codes, MCBEAM and MCPLAN, were used for the machine head simulation and phantom/patient dose calculation, respectively. In the new Co-60 machine design, the source-to-isocenter distance was 60cm and the treatment head included the Co-60 source, primary collimator, jaws and MLC. Results: For all cases investigated, conformal arc plans utilizing Co-60 beams achieved similar conformity (mean conformity index=1.19) comparing to 6MV photon beams. Isodose distributions were tailored similarly around the PTV; both Co-60 and 6MV plans met our clinical acceptance criteria for the target coverage, and the maximum and minimum target doses. The DVH for the Co-60 plans showed slightly lower doses to the critical structures although the differences were small in most cases. Conclusion: There were no significant dosimetric differences between 6MV and Co-60 plans. Our results confirmed that this new Co-60 design could be a cost-effective machine for advanced radiotherapy due to its low cost, low maintenance and high up time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-533: LET Dependence Correction of Radiochromic Films for Application in Low Energy Proton Irradiation

Purpose: Many cell irradiation experiments with low-energy laser-driven ions rely on radiochromic films (RCF), because of their dose-rate independent response and superior spatial resolution. RCF dosimetry in low-energy ion beams requires a correction of the LET dependent film response. The relative efficiency (RE), the ratio of photon to proton dose that yields the same film darkening, is a measure for the film’s LET dependence. A direct way of RE determination is RCF irradiation with low-energy mono-energetic protons and hence, well-defined LET. However, RE is usually determined using high energy proton depth dose measurements where RE corrections require knowledge of the average LET in each depth, which can be either track (tLET) or dose (dLET) averaged. The appropriate LET concept has to be applied to allow a proper film response correction. Methods: Radiochromic EBT2 and EBT3 films were irradiated in clinical photon and proton beams. For each depth of the 200 MeV proton depth dose curve, tLET and dLET were calculated by special user routines from the Monte Carlo code FLUKA. Additional irradiations with mono-energetic low energy protons (4–20 MeV) serve as reference for the RE determination. Results: The difference of dLET and tLET increases with depth, with the dLET being almost twice as large as the tLET for the maximum depth. The comparison with mono-energetic measurements shows a good agreement of the RE for the dLET concept, while a considerably steeper drop in RE is observed when applying the tLET. Conclusion: RCF can be used as reference dosimeter for biomedical experiments with low-energy proton beams if appropriate LET corrections are applied. When using depth dose measurements from clinical proton accelerators for these corrections, the concept of dLET has to be applied. Acknowledgement: This work was funded by the DFG Cluster of Excellence ‘Munich-Centre for Advanced Photonics’ (MAP). This work was funded by the DFG Cluster of Excellence Munich-Centre for Advanced Photonics (MAP)

61 RADIATION PROTECTION AND DOSIMETRY↗
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 OSTI2015

Monte Carlo simulations of patient dose perturbations in rotational-type radiotherapy due to a transverse magnetic field: A tomotherapy investigation

Purpose: Several groups are exploring the integration of magnetic resonance (MR) image guidance with radiotherapy to reduce tumor position uncertainty during photon radiotherapy. The therapeutic gain from reducing tumor position uncertainty using intrafraction MR imaging during radiotherapy could be partially offset if the negative effects of magnetic field-induced dose perturbations are not appreciated or accounted for. The authors hypothesize that a more rotationally symmetric modality such as helical tomotherapy will permit a systematic mediation of these dose perturbations. This investigation offers a unique look at the dose perturbations due to homogeneous transverse magnetic field during the delivery of Tomotherapy{sup ®} Treatment System plans under varying degrees of rotational beamlet symmetry. Methods: The authors accurately reproduced treatment plan beamlet and patient configurations using the Monte Carlo code GEANT4. This code has a thoroughly benchmarked electromagnetic particle transport physics package well-suited for the radiotherapy energy regime. The three approved clinical treatment plans for this study were for a prostate, head and neck, and lung treatment. The dose heterogeneity index metric was used to quantify the effect of the dose perturbations to the target volumes. Results: The authors demonstrate the ability to reproduce the clinical dose–volume histograms (DVH) to within 4% dose agreement at each DVH point for the target volumes and most planning structures, and therefore, are able to confidently examine the effects of transverse magnetic fields on the plans. The authors investigated field strengths of 0.35, 0.7, 1, 1.5, and 3 T. Changes to the dose heterogeneity index of 0.1% were seen in the prostate and head and neck case, reflecting negligible dose perturbations to the target volumes, a change from 5.5% to 20.1% was observed with the lung case. Conclusions: This study demonstrated that the effect of external magnetic fields can be mitigated by exploiting a more rotationally symmetric treatment modality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Black holes in young stellar clusters

We present theoretical models for stellar black hole (BH) properties in young, massive star clusters. Using a Monte Carlo code for stellar dynamics, we model realistic star clusters with N ≅ 5 × 10{sup 5} stars and significant binary fractions (up to 50%) with self-consistent treatments of stellar dynamics and stellar evolution. We compute the formation rates and characteristic properties of single and binary BHs for various representative ages, cluster parameters, and metallicities. Because of dynamical interactions and supernova (SN) kicks, more single BHs end up retained in clusters compared to BHs in binaries. We also find that the ejection of BHs from a cluster is a strong function of initial density. In low-density clusters (where dynamical effects are negligible), it is mainly SN kicks that eject BHs from the cluster, whereas in high-density clusters (initial central density ρ {sub c}(0) ∼ 10{sup 5} M {sub ☉} pc{sup –3} in our models) the BH ejection rate is enhanced significantly by dynamics. Dynamical interactions of binary systems in dense clusters also modify the orbital period and eccentricity distributions while increasing the probability of a BH having a more massive companion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Sci—Thur AM: YIS - 04: Gold Nanoparticle Enhanced Arc Radiotherapy: A Monte Carlo Feasibility Study

Introduction: The use of gold nanoparticles (GNPs) in radiotherapy has shown promise for therapeutic enhancement. In this study, we explore the feasibility of enhancing radiotherapy with GNPs in an arc-therapy context. We use Monte Carlo simulations to quantify the macroscopic dose-enhancement ratio (DER) and tumour to normal tissue ratio (TNTR) as functions of photon energy over various tumour and body geometries. Methods: GNP-enhanced arc radiotherapy (GEART) was simulated using the PENELOPE Monte Carlo code and penEasy main program. We simulated 360° arc-therapy with monoenergetic photon energies 50 – 1000 keV and several clinical spectra used to treat a spherical tumour containing uniformly distributed GNPs in a cylindrical tissue phantom. Various geometries were used to simulate different tumour sizes and depths. Voxel dose was used to calculate DERs and TNTRs. Inhomogeneity effects were examined through skull dose in brain tumour treatment simulations. Results: Below 100 keV, DERs greater than 2.0 were observed. Compared to 6 MV, tumour dose at low energies was more conformai, with lower normal tissue dose and higher TNTRs. Both the DER and TNTR increased with increasing cylinder radius and decreasing tumour radius. The inclusion of bone showed excellent tumour conformality at low energies, though with an increase in skull dose (40% of tumour dose with 100 keV compared to 25% with 6 MV). Conclusions: Even in the presence of inhomogeneities, our results show promise for the treatment of deep-seated tumours with low-energy GEART, with greater tumour dose conformality and lower normal tissue dose than 6 MV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-90: Accuracy of Calibration of Lithium-6 and -7 Enriched LiF TLDs for Neutron Measurements in High Energy Radiotherapy

Purpose: To determine the potential error involved in the interpretation of neutron measurements from medical linear accelerators (linacs) using TLD-600H and TLD-700H if standard AmBe and {sup 252}Cf neutron sources are used for calibration without proper inclusion of neutron energy spectrum information. Methods: The Kerma due to neutrons can be calculated from the energy released by various nuclear interactions (elastic and inelastic scatter, (n,α), (n,p), (n,d), (n,t), (n,2n), etc.). The response of each TLD can be considered the sum of the neutron and gamma components; each proportional to the Kerma. Using the difference between the measured TLD responses and the ratio of the calculated Kerma for each material, the neutron component of the response can be calculated. The Monte Carlo code MCNP6 has been used to calculate the neutron energy spectra resulting from photonuclear interactions in a Varian 21EX linac. TLDs have been exposed to the mixed (γ-n) field produced by a linac and AmBe and {sup 252}Cf standard neutron sources. Results: For dosimetry of neutrons from AmBe or {sup 252}Cf sources, assuming TLD-700H insensitivity to neutrons will Result in 10% or 20% overestimation of neutron doses respectively.For dosimetry of neutrons produced in a Varian 21EX, applying a calibration factor derived from a standard AmBe or {sup 252}Cf source will Result in an overestimation of neutron fluence, by as much as a factor of 47.The assumption of TLD-700H insensitivity to neutrons produced by linacs leads to a negligible error due to the extremely high Kerma ratio (600H/700H) of 3000 for the assumed neutron spectrum. Conclusion: Lithium-enriched TLDs calibrated with AmBe and/or {sup 252}Cf neutron sources are not accurate for use under the neutron energy spectrum produced by a medical linear accelerator.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-146: Comparing Monte Carlo Simulations with Commissioning Beam Data for Mevion S250 Proton Therapy System

Purpose: To evaluate the Monte Carlo simulated beam data with the measured commissioning data for the Mevion S250 proton therapy system. Method: The Mevion S250 proton therapy system utilizes a passive double scattering technique with a unique gantry mounted superconducting accelerator and offers effective proton therapy in a compact design concept. The field shaping system (FSS) includes first scattering foil, range modulator wheel (RMW), second scattering foil and post absorber and offers two field sizes and a total of 24 treatment options from proton range of 5 cm to 32 cm. The treatment nozzle was modeled in detail using TOPAS (TOolkit for PArticle Simulation) Monte Carlo code. The timing feathers of the moving modulator wheels were also implemented to generate the Spread Out Bragg Peak (SOBP). The simulation results including pristine Bragg Peak, SOBP and dose profiles were compared with the data measured during beam commissioning. Results: The comparison between the measured data and the simulation data show excellent agreement. For pristine proton Bragg Peaks, the simulated proton range (depth of distal 90%) values agreed well with the measured range values within 1 mm accuracy. The differences of the distal falloffs (depth from distal 80% to 20%) were also found to be less than 1 mm between the simulations and measurements. For the SOBP, the widths of modulation (depth of proximal 95% to distal 90%) were also found to agree with the measurement within 1 mm. The flatness of the simulated and measured lateral profiles was found to be 0.6 % and 1.1 %, respectively. Conclusion: The agreement between simulations and measurements demonstrate that TOPAS could be used as a viable platform to proton therapy applications. The matched simulation results offer a great tool and open opportunity for variety of applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MAPS OF EVOLVING CLOUD STRUCTURES IN LUHMAN 16AB FROM HST TIME-RESOLVED SPECTROSCOPY

WISE J104915.57-531906.1 is the nearest brown dwarf binary to our solar system, consisting of two brown dwarfs in the L / T transition: Luhman 16A and B. In this paper, we present the first map of Luhman 16A, and maps of Luhman 16B for two epochs. Our maps were created by applying Aeolus , a Markov-Chain Monte Carlo code that maps the top-of-the-atmosphere (TOA) structure of brown dwarf and other ultracool atmospheres, to light curves of Luhman 16A and B using the Hubble Space Telescope ’s G141 and G102 grisms. Aeolus retrieved three or four spots in the TOA of Luhman 16A and B, with a surface coverage of 19%–32% (depending on an assumed rotational period of 5 hr or 8 hr) or 21%–38.5% (depending on the observational epoch), respectively. The brightness temperature of the spots of the best-fit models was ∼200 K hotter than the background TOA. We compared our Luhman 16B map with the only previously published map. Interestingly, our map contained a large TOA spot that was cooler (Δ T ∼ 51 K) than the background, which lay at low latitudes, in agreement with the previous Luhman 16B map. Finally, we report the detection of a feature reappearing in Luhman 16B light curves that are separated by tens of hundreds of rotations from each other. We speculate that this feature is related to TOA structures of Luhman 16B.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-I-37: Eye Lens Dose Reduction From CT Scan Using Organ Based Tube Current Modulation

Purpose: To investigate the eye lens dose reduction by CT scan with organ based tube current modulation (OBTCM) using GPU Monte Carlo code ARCHER-CT. Methods: 36 X-ray sources and bowtie filters were placed around the patient head with the projection angle interval of 10° for one rotation of CT scan, each projection was simulated respectively. The voxel eye models with high resolution(0.1mm*0.1mm*0.1mm) were used in the simulation and different tube voltage including 80kVp, 100kVp, 120kVp and 140kVp were taken into consideration. Results: The radiation doses to the eye lens increased with the tube voltage raised from 80kVp to 140kVp, and the dose results from 0° (AP) direction are much higher than those from 180° (PA) direction for all the 4 different tube voltage investigated. This 360° projection dose characteristic enables organ based TCM, which can reduce the eye lens dose by more than 55%. Conclusion: As the eye lens belongs to superficial tissues, its radiation dose to external exposure like CT is direction sensitive, and this characteristic feature makes organ based TCM to be an effective way to reduce the eye lens dose, so more clinical use of this technique were recommended. National Nature Science Foundation of China(No.11475047)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Shielded radiography with a laser-driven MeV-energy X-ray source

We report the results of experimental and numerical-simulation studies of shielded radiography using narrowband MeV-energy X-rays from a compact all-laser-driven inverse-Compton-scattering X-ray light source. This recently developed X-ray light source is based on a laser-wakefield accelerator with ultra-high-field gradient (GeV/cm). We demonstrate experimentally high-quality radiographic imaging (image contrast of 0.4 and signal-to-noise ratio of 2:1) of a target composed of 8-mm thick depleted uranium shielded by 80-mm thick steel, using a 6-MeV X-ray beam with a spread of 45% (FWHM) and 10 7 photons in a single shot. The corresponding dose of the X-ray pulse measured in front of the target is ~100 nGy/pulse. Simulations performed using the Monte-Carlo code MCNPX accurately reproduce the experimental results. These simulations also demonstrate that the narrow bandwidth of the Compton X-ray source operating at 6 and 9 MeV leads to a reduction of deposited dose as compared to broadband bremsstrahlung sources with the same end-point energy. The X-ray beam’s inherently low-divergence angle (~mrad) is advantageous and effective for interrogation at standoff distance. Furthermore, these results demonstrate significant benefits of all-laser driven Compton X-rays for shielded radiography.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

SU-F-T-656: Monte Carlo Study On Air Activation Around a Medical Electron Linac

Purpose: In high energy photon therapy, several radiation protection issues result from photonuclear reactions. The activation of air - directly by photonuclear reactions as well as indirectly by capture of photoneutrons generated inside the linac head - is a major point of concern for the medical staff. The purpose of this study was to estimate the annual effective dose to medical workers due to activated air around a medical high energy electron linac by means of Monte Carlo simulations. Methods: The treatment head of a Varian Clinac in 18 MV-X mode as well as the surrounding concrete bunker were modeled and the radiation transport was simulated using the Monte Carlo code FLUKA, starting from the primary electron striking the bremsstrahlung target. The activation yields in air from photo-disintegration of O-16 and N-14 nuclei as well as from neutron capture on Ar-40 nuclei were obtained from the simulations. The activation build-up, radioactive decay and air ventilation were studied using a mathematical model. The annual effective dose to workers was estimated by using published isotope specific conversion factors. Results: The oxygen and nitrogen activation yields were in contrast to the argon activation yield found to be field size dependent. The impact of the treatment room ventilation on the different air activation products was investigated and quantified. An estimate with very conservative assumptions gave an annual effective dose to workers of < 1 mSv/a. Conclusion: From the results of this study it can be concluded that the contribution of air activation to the radiation exposure to medical workers should be negligible in modern photon therapy, especially when it is compared to the dose due to prompt neutrons and the activation of heavy solid materials such as the jaws and the collimators inside the linac head.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A scintillator-based approach to monitor secondary neutron production during proton therapy

Here, the primary objective of this work is to measure the secondary neutron field produced by an uncollimated proton pencil beam impinging on different tissue-equivalent phantom materials using organic scintillation detectors. Additionally, the Monte Carlo code mcnpx-PoliMi was used to simulate the detector response for comparison to the measured data. Comparison of the measured and simulated data will validate this approach for monitoring secondary neutron dose during proton therapy. Proton beams of 155- and 200-MeV were used to irradiate a variety of phantom materials and secondary particles were detected using organic liquid scintillators. These detectors are sensitive to fast neutrons and gamma rays: pulse shape discrimination was used to classify each detected pulse as either a neutron or a gamma ray. The mcnpx-PoliMi code was used to simulate the secondary neutron field produced during proton irradiation of the same tissue-equivalent phantom materials. As a result, an experiment was performed at the Loma Linda University Medical Center proton therapy research beam line and corresponding models were created using the mcnpx-PoliMi code. The authors’ analysis showed agreement between the simulations and the measurements. The simulated detector response can be used to validate the simulations of neutron and gamma doses on a particular beam line with or without a phantom. In conclusion, the authors have demonstrated a method of monitoring the neutron component of the secondary radiation field produced by therapeutic protons. The method relies on direct detection of secondary neutrons and gamma rays using organic scintillation detectors. These detectors are sensitive over the full range of biologically relevant neutron energies above 0.5 MeV and allow effective discrimination between neutron and photon dose. Because the detector system is portable, the described system could be used in the future to evaluate secondary neutron and gamma doses on various clinical beam lines for commissioning and prospective data collection in pediatric patients treated with proton therapy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗