Results for “OPACITY” — Page 41

DOE OSTI2016

A stable 1D multigroup high-order low-order method

The high-order low-order (HOLO) method is a recently developed moment-based acceleration scheme for solving time-dependent thermal radiative transfer problems, and has been shown to exhibit orders of magnitude speedups over traditional time-stepping schemes. However, a linear stability analysis by Haut et al. (2015 Haut, T. S., Lowrie, R. B., Park, H., Rauenzahn, R. M., Wollaber, A. B. (2015). A linear stability analysis of the multigroup High-Order Low-Order (HOLO) method. In Proceedings of the Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method; Nashville, TN, April 19–23, 2015. American Nuclear Society.) revealed that the current formulation of the multigroup HOLO method was unstable in certain parameter regions. Since then, we have replaced the intensity-weighted opacity in the first angular moment equation of the low-order (LO) system with the Rosseland opacity. Furthermore, this results in a modified HOLO method (HOLO-R) that is significantly more stable.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2013

RESOLVED NEAR-INFRARED SPECTROSCOPY OF WISE J104915.57-531906.1AB: A FLUX-REVERSAL BINARY AT THE L DWARF/T DWARF TRANSITION

We report resolved near-infrared spectroscopy and photometry of the recently identified brown dwarf binary WISE J104915.57-531906.1AB, located 2.02 {+-} 0.15 pc from the Sun. Low-resolution spectral data from Magellan/FIRE and IRTF/SpeX reveal strong H{sub 2}O and CO absorption features in the spectra of both components, while the secondary also exhibits weak CH{sub 4} absorption at 1.6 {mu}m and 2.2 {mu}m. Spectral indices and comparison to low-resolution spectral standards indicate component types of L7.5 and T0.5 {+-} 1, the former consistent with the optical classification of the primary. Both sources also have unusually red spectral energy distributions for their spectral types, which we attribute to enhanced condensate opacity (thick clouds). Relative photometry reveals a flux reversal between the J and K bands, with the T dwarf component being brighter in the 0.95-1.3 {mu}m region ({Delta}J = -0.31 {+-} 0.05). As with other L/T transition binaries, this reversal likely reflects the depletion of condensate opacity in the T dwarf, with the contrast enhanced by the thick clouds present in the photosphere of the L dwarf primary. The 1 {mu}m flux from the T dwarf most likely emerges from gaps in its cloud layer, as suggested by the significant optical variability detected from this source by Gillon et al. Component mass measurements of the WISE J1049-5319AB system through astrometric and component radial velocity monitoring may resolve the current debate as to whether the loss of photospheric condensate clouds at the L dwarf/T dwarf boundary is a slow or rapid process, a conceivable endeavor given its proximity, brightness, small separation (3.1 {+-} 0.3 AU), and reasonable orbital period (20-30 yr)

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

The impact of dust evolution and photoevaporation on disk dispersal

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron-sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet (FUV) photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains, leaving the larger particles—which decouple, from the gas—behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV, and X-rays, dust evolution, and radial drift using a one-dimensional (1D) multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, ∼3×10{sup −4} M {sub ⊙} of mass in solids typically remain, comparable to the solids inventory of our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Herschel dust emission as a probe of starless cores mass: MCLD 123.5+24.9 of the polaris flare

We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, and hence masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-infrared grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) ≫ 1 × 10{sup 22} cm{sup −2}), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that {sup 13}CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Numerical simulations of turbulent molecular clouds regulated by reprocessed radiation feedback from nascent super star clusters

Radiation feedback from young star clusters embedded in giant molecular clouds (GMCs) is believed to be important to the control of star formation. For the most massive and dense clouds, including those in which super star clusters (SSCs) are born, pressure from reprocessed radiation exerted on dust grains may disperse a significant portion of the cloud mass back into the interstellar medium. Using our radiation hydrodynamics code, Hyperion, we conduct a series of numerical simulations to test this idea. Our models follow the evolution of self-gravitating, strongly turbulent clouds in which collapsing regions are replaced by radiating sink particles representing stellar clusters. We evaluate the dependence of the star formation efficiency (SFE) on the size and mass of the cloud and κ, the opacity of the gas to infrared (IR) radiation. We find that the single most important parameter determining the evolutionary outcome is κ, with κ≳15 cm{sup 2} g{sup −1} needed to disrupt clouds. For κ=20−40 cm{sup 2} g{sup −1}, the resulting SFE = 50%−70% is similar to empirical estimates for some SSC-forming clouds. The opacities required for GMC disruption likely apply only in dust-enriched environments. We find that the subgrid model approach of boosting the direct radiation force L/c by a “trapping factor” equal to a cloud’s mean IR optical depth can overestimate the true radiation force by factors of ∼4−5. We conclude that feedback from reprocessed IR radiation alone is unlikely to significantly reduce star formation within GMCs unless their dust abundances or cluster light-to-mass ratios are enhanced.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Detection of significant cm to sub-mm band radio and γ-ray correlated variability in Fermi bright blazars

The exact location of the γ-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA programme) and γ-ray (0.1–300 GeV) ~3.5 yr light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. These results reveal: (i) the first highly significant multiband radio and γ-ray correlations (radio lagging γ rays) when averaging over the whole sample, (ii) average time delays (source frame: 76 ± 23 to 7 ± 9 d), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence τr, γ(ν) ∝ ν -1 , in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk γ-ray production region typically located within/upstream of the 3 mm core region (τ3mm, γ = 12 ± 8 d), (iv) mean distances between the region of γ-ray peak emission and the radio ‘τ = 1 photosphere’ decreasing from 9.8 ± 3.0 pc (11 cm) to 0.9 ± 1.1 pc (2 mm) and 1.4 ± 0.8 pc (0.8 mm), (v) 3 mm/γ-ray correlations in nine individual sources at a significance level where one is expected by chance (probability: 4 × 10 -6 ), (vi) opacity and ‘time lag core shift’ estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk γ-ray production region from the supermassive black hole (SMBH) of ~0.8–1.6 pc, i.e. at the outer edge of the broad-line region (BLR) or beyond. A 3 mm τ = 1 surface at ~2–3 pc from the jet base (i.e. well outside the ‘canonical BLR’) finally suggests that BLR material extends to several parsec distances from the SMBH.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

What physics determines the peak of the IMF? Insights from the structure of cores in radiation-magnetohydrodynamic simulations

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical time-scale, the opacity limit for fragmentation, at ~0.01M⊙. However, the observed peak of the initial mass function (IMF) lies a factor of 20-30 higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. Here, we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, ~0.01M⊙, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why ~0.01M⊙ objects are rare: unless an outside agent intervenes (e.g. a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of ~0.2M⊙, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at ~0.2M⊙.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Quantum simulations of nuclei and nuclear pasta with the multiresolution adaptive numerical environment for scientific simulations

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation, and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle three-dimensional (3D) Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of 16 O, 208 Pb, and 238 U nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24 fm with an average density of about ρ = 0.05fm -3 , proton fraction of Yp = 0.3, and temperature of T = 0 MeV. Overall, our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. We compare calculations with and without spin-orbit forces. We find that while subtle differences are present, the pasta phase remains in the waffle geometry. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Effects of Landau quantization on neutrino emission and absorption

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as 10 15 G and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as 10 16 G. In the presence of a strong magnetic field, Leinson and Perez showed and Baiko and Yakovlev confirmed that the energy levels of electrons and protons are quantized and the direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields 𝐵 ≳ 10 16 G, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. Here, we present a fully relativistic calculation of the direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wave functions present and a new semianalytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at subsaturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

astrophysical electromagnetic fields↗

DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

RadVolViz: An Information Display-Inspired Transfer Function Editor for Multivariate Volume Visualization

In volume visualization transfer functions are widely used for mapping voxel properties to color and opacity. Typically, volume density data are scalars which require simple 1D transfer functions to achieve this mapping. If the volume densities are vectors of three channels, one can straightforwardly map each channel to either red, green or blue, which requires a trivial extension of the 1D transfer function editor. Here, we devise a new method that applies to volume data with more than three channels. These types of data often arise in scientific scanning applications, where the data are separated into spectral bands or chemical elements. Our method expands on prior work in which a multivariate information display, RadViz, was fused with a radial color map, in order to visualize multi-band 2D images. In this work, we extend this joint interface to blended volume rendering. The information display allows users to recognize the presence and value distribution of the multivariate voxels and the joint volume rendering display visualizes their spatial distribution. We design a set of operators and lenses that allow users to interactively control the mapping of the multivariate voxels to opacity and color. This enables users to isolate or emphasize volumetric structures with desired multivariate properties. Furthermore, it turns out that our method also enables more insightful displays even for RGB data. We demonstrate our method with three datasets obtained from spectral electron microscopy, high energy X-ray scanning, and atmospheric science.

36 MATERIALS SCIENCE↗

DOE OSTI2013

A prototype fan-beam optical CT scanner for 3D dosimetry

Purpose: The objective of this work is to introduce a prototype fan-beam optical computed tomography scanner for three-dimensional (3D) radiation dosimetry. Methods: Two techniques of fan-beam creation were evaluated: a helium-neon laser (HeNe, {lambda} = 543 nm) with line-generating lens, and a laser diode module (LDM, {lambda} = 635 nm) with line-creating head module. Two physical collimator designs were assessed: a single-slot collimator and a multihole collimator. Optimal collimator depth was determined by observing the signal of a single photodiode with varying collimator depths. A method of extending the dynamic range of the system is presented. Two sample types were used for evaluations: nondosimetric absorbent solutions and irradiated polymer gel dosimeters, each housed in 1 liter cylindrical plastic flasks. Imaging protocol investigations were performed to address ring artefacts and image noise. Two image artefact removal techniques were performed in sinogram space. Collimator efficacy was evaluated by imaging highly opaque samples of scatter-based and absorption-based solutions. A noise-based flask registration technique was developed. Two protocols for gel manufacture were examined. Results: The LDM proved advantageous over the HeNe laser due to its reduced noise. Also, the LDM uses a wavelength more suitable for the PRESAGE{sup TM} dosimeter. Collimator depth of 1.5 cm was found to be an optimal balance between scatter rejection, signal strength, and manufacture ease. The multihole collimator is capable of maintaining accurate scatter-rejection to high levels of opacity with scatter-based solutions (T < 0.015%). Imaging protocol investigations support the need for preirradiation and postirradiation scanning to reduce reflection-based ring artefacts and to accommodate flask imperfections and gel inhomogeneities. Artefact removal techniques in sinogram space eliminate streaking artefacts and reduce ring artefacts of up to {approx}40% in magnitude. The flask registration technique was shown to achieve submillimetre and subdegree placement accuracy. Dosimetry protocol investigations emphasize the need to allow gel dosimeters to cool gradually and to be scanned while at room temperature. Preliminary tests show that considerable noise reduction can be achieved with sinogram filtering and by binning image pixels into more clinically relevant grid sizes. Conclusions: This paper describes a new optical CT scanner for 3D radiation dosimetry. Tests demonstrate that it is capable of imaging both absorption-based and scatter-based samples of high opacities. Imaging protocol and gel dosimeter manufacture techniques have been adapted to produce optimal reconstruction results. These optimal results will require suitable filtering and binning techniques for noise reduction purposes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Sci—Fri PM: Dosimetry—01: Radiation-induced refraction artefacts in the optical CT readout of polymer gel dosimeters

Polymer gel dosimeters (PGDs) are a desirable tool for the verification of advanced radiotherapy treatments. Fully 3D, deformable, and tissue-equivalent, the PGD polymerizes wherever it absorbs dose. To measure the dose absorbed by a PGD, optical computed tomography (CT) can be used to evaluate, in full 3D, the opacity distribution that coincides with polymerization. In addition to an increase in opacity with dose, an increase in refractive index (RI) is also known to occur in irradiated polymer gels. The increase in RI is slight and was previously assumed insignificant. This work reveals the effects that radiation-induced RI changes can have on the optical CT readout of PGDs. A fan-beam optical CT scanner was used to image a cylindrical PGD irradiated by a pair of 3×3 cm{sup 2}, 6 MV photon beams in an orthogonal arrangement. Investigative scans were performed to evaluate refraction errors occurring: i) within the plane, and ii) out of the plane of the fan-beam. In-plane refraction was shown to cause distinct streaking artefacts along dose gradients (i.e. RI gradients) due to higher intensity rays being refracted into more opaque regions. Out-of-plane refraction was shown to produce severe, widespread artefacts due to rays missing the detector array. An iterative Savitzky-Golay filtering technique was developed to reduce both types of artefacts by specifically targeting structured errors in sinogram space. Results introduce a new category of imaging artefacts to be aware of when using optical CT for PGD readout.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

TU-E-201-03: Eye Lens Dosimetry in Radiotherapy Using Contact Lens-Shaped Applicator

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

TU-E-201-00: Eye Lens Dosimetry for Patients and Staff

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

TU-E-201-02: Eye Lens Dosimetry From CT Perfusion Studies

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-J-129: Verification of Geometric and Dosimetric Accuracy of Respiratory Management Systems Using Homemade Phantom

Purpose: Different placements of Infrared Cameras (IRC) in CT and treatment rooms can cause gating window level (GWL) variations leading to differences between GWL used for planning and treatments. Although, Varian Clinac DHX-OBI sytem and CT are equipped with the same kind of IRC, Truebeam STx (TB) has a different type of IRC known as banana type. In this study; geometric and dosimetric accuracy of respiratory management system (RPM) for different machines were investigated with a special homemade phantom. Methods: Special phantom was placed on the respiratory simulator machine and a CT data set was obtained at the end of the expirium phase (EOE). Conformal and IMRT plans were generated on the EOE CT image series for both DHX-OBI and TB LINACs while a VMAT plan was generated only for TB.The acquired respiratory graphs in the CT were directly sent to DHX-OBI system, and they were converted with software before sending to TB. EBT3 films were placed inside the phantom and were irradiated using RPM system with two machines for different plans. Planar dose distributions were compared with gamma analysis (GA) method (3mm, %3) to evaluate planned-measured dose differences. In addition, radio-opac marker was placed in the center of the phantom to evaluate the geometric accuracy of treatment field with gated flouroscopy (GF). Results: There were no shifts detected between planning and treeatment GWL for both DHX-OBI and TB. Difference on the GF image between digital graticule and radio-opac marker was <1mm for TB and 1mm for DHX-OBI. Although, GA agreement was 97% for conformal and IMRT techniques in TB, it was 96% for VMAT technique. While GA agreement was 98% for conformal technique in DHX-OBI, IMRT was 95%.ConclusionThis study showed that RPM can be used accurately in spite of different IRC placements or different types of ICR used.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Branson: A Mini-App for Studying Parallel IMC, Version 1.0

This code solves the gray thermal radiative transfer (TRT) equations in parallel using simple opacities and Cartesian meshes. Although Branson solves the TRT equations it is not designed to model radiation transport: Branson contains simple physics and does not have a multigroup treatment, nor can it use physical material data. The opacities have are simple polynomials in temperature there is a limited ability to specify complex geometries and sources. Branson was designed only to capture the computational demands of production IMC codes, especially in large parallel runs. It was also intended to foster collaboration with vendors, universities and other DOE partners. Branson is similar in character to the neutron transport proxy-app Quicksilver from LLNL, which was recently open-sourced.

Long, Alex↗