SEARCH · Engineering Papers

Results for “PHASE SPACE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

DOE OSTI2017

EXPLICIT SYMPLECTIC-LIKE INTEGRATORS WITH MIDPOINT PERMUTATIONS FOR SPINNING COMPACT BINARIES

We refine the recently developed fourth-order extended phase space explicit symplectic-like methods for inseparable Hamiltonians using Yoshida’s triple product combined with a midpoint permuted map. The midpoint between the original variables and their corresponding extended variables at every integration step is readjusted as the initial values of the original variables and their corresponding extended ones at the next step integration. The triple-product construction is apparently superior to the composition of two triple products in computational efficiency. Above all, the new midpoint permutations are more effective in restraining the equality of the original variables and their corresponding extended ones at each integration step than the existing sequent permutations of momenta and coordinates. As a result, our new construction shares the benefit of implicit symplectic integrators in the conservation of the second post-Newtonian Hamiltonian of spinning compact binaries. Especially for the chaotic case, it can work well, but the existing sequent permuted algorithm cannot. When dissipative effects from the gravitational radiation reaction are included, the new symplectic-like method has a secular drift in the energy error of the dissipative system for the orbits that are regular in the absence of radiation, as an implicit symplectic integrator does. In spite of this, it is superior to the same-order implicit symplectic integrator in accuracy and efficiency. The new method is particularly useful in discussing the long-term evolution of inseparable Hamiltonian problems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

PREDICTING A THIRD PLANET IN THE KEPLER-47 CIRCUMBINARY SYSTEM

We have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system is the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by integrating the entire five bodies (binary + planets) for 10 Myr. We identified several stable regions between the two known planets as well as a region beyond the orbit of Kepler-47c where the orbit of the third planet could be stable. To constrain the orbit of this planet, we used the measured duration of the unexplained transit event (∼4.15 hr) and compared that with the transit duration of the third planet in an ensemble of stable orbits. To remove the degeneracy among the orbits with similar transit durations, we considered the planet to be in a circular orbit and calculated its period analytically. The latter places an upper limit of 424 days on the orbital period of the third planet. Our analysis suggests that if the unexplained transit event detected during the discovery of the Kepler-47 circumbinary system is due to a planetary object, this planet will be in a low eccentricity orbit with a semi-major axis smaller than 1.24 AU. Further constraining of the mass and orbital elements of this planet requires a re-analysis of the entire currently available data, including those obtained post-announcement of the discovery of this system. We present details of our methodology and discuss the implication of the results.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

SU-F-T-79: Monte Carlo Investigation of Optimizing Parameters for Modulated Electron Arc Therapy

Purpose: Electron arc therapy provides excellent dose distributions for treating superficial tumors along curved surfaces. However this modality has not received widespread application due to the lack of needed advancement in electron beam delivery, accurate electron dose calculation and treatment plan optimization. The aim of the current work is to investigate possible parameters that can be optimized for electron arc (eARC) therapy. Methods: The MCBEAM code was used to generate phase space files for 6 and 12MeV electron beam energies from a Varian trilogy machine. An Electron Multi-leaf collimator eMLC of 2cm thickness positioned at 82 cm source collimated distance was used in the study. Dose distributions for electron arcs were calculated inside a cylindrical phantom using the MCSIM code. The Cylindrical phantom was constructed with 0.2cm voxels and a 15cm diameter. Electron arcs were delivered with two different approaches. The first approach was to deliver the arc as segments of very small field widths. In this approach we also tested the impact of the segment size and the arc increment angle. The second approach is to deliver the arc as a sum of large fields each covering the whole target as seen from the beam eye view. Results: In considering 90 % as the prescription isodose line, the first approach showed a region of buildup proceeding before the prescription zone. This build up is minimizing with the second approach neglecting need for bolus. The second approach also showed less x-ray contamination. In both approaches the variation of the segment size changed the size and location of the prescription isodose line. The optimization process for eARC could involve interplay between small and large segments to achieve desired coverage. Conclusion: An advanced modulation of eARCs will allow for tailored dose distribution for superficial curved target as with challenging scalp cases.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

A review on the relativistic effective field theory with parameterized couplings for nuclear matter and neutron stars

Nuclear science has developed many excellent theoretical models for many-body systems in the domain of the baryon-meson strong interaction for the nucleus and nuclear matter at low, medium and high densities. However, a full microscopic understanding of nuclear systems in the extreme density domain of compact stars is still lacking. The aim of this contribution is to shed some light on open questions facing the nuclear many-body problem at the very high density domain. Here we focus our attention on the conceptual issue of naturalness and its role in shaping the baryon-meson phase space dynamics in the description of the equation of state (EoS) of nuclear matter and neutrons stars. In particular, in order to stimulate possible new directions of research, we discuss relevant aspects of a recently developed relativistic effective theory for nuclear matter within Quantum Hadrodynamics (QHD) with genuine many-body forces and derivative natural parametric couplings. Among other topics we discuss in this work the connection of this theory with other known effective QHD models of the literature and its potentiality in describing a new physics for dense matter. The model with parameterized couplings exhausts the whole fundamental baryon octet (n, p, Σ{sup −}, Σ{sup 0}, Σ{sup +}, Λ, Ξ{sup −}, Ξ{sup 0}) and simulates n-order corrections to the minimal Yukawa baryon couplings by considering nonlinear self-couplings of meson fields and meson-meson interaction terms coupled to the baryon fields involving scalar-isoscalar (σ, σ∗), vector-isoscalar (ω, Φ), vector-isovector (ϱ) and scalar-isovector (δ) virtual sectors. Following recent experimental results, we consider in our calculations the extreme case where the Σ{sup −} experiences such a strong repulsion that its influence in the nuclear structure of a neutron star is excluded at all. A few examples of calculations of properties of neutron stars are shown and prospects for the future are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2015

SU-E-T-644: QuAArC: A 3D VMAT QA System Based On Radiochromic Film and Monte Carlo Simulation of Log Files

Purpose: VMAT involves two main sources of uncertainty: one related to the dose calculation accuracy, and the other linked to the continuous delivery of a discrete calculation. The purpose of this work is to present QuAArC, an alternative VMAT QA system to control and potentially reduce these uncertainties. Methods: An automated MC simulation of log files, recorded during VMAT treatment plans delivery, was implemented in order to simulate the actual treatment parameters. The linac head models and the phase-space data of each Control Point (CP) were simulated using the EGSnrc/BEAMnrc MC code, and the corresponding dose calculation was carried out by means of BEAMDOSE, a DOSXYZnrc code modification. A cylindrical phantom was specifically designed to host films rolled up at different radial distances from the isocenter, for a 3D and continuous dosimetric verification. It also allows axial and/or coronal films and point measurements with several types of ion chambers at different locations. Specific software was developed in MATLAB in order to process and evaluate the dosimetric measurements, which incorporates the analysis of dose distributions, profiles, dose difference maps, and 2D/3D gamma index. It is also possible to obtain the experimental DVH reconstructed on the patient CT, by an optimization method to find the individual contribution corresponding to each CP on the film, taking into account the total measured dose, and the corresponding CP dose calculated by MC. Results: The QuAArC system showed high reproducibility of measurements, and consistency with the results obtained with the commercial system implemented in the verification of the evaluated treatment plans. Conclusion: A VMAT QA system based on MC simulation and high resolution dosimetry with film has been developed for treatment verification. It shows to be useful for the study of the real VMAT capabilities, and also for linac commissioning and evaluation of other verification devices.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

Purpose: The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice. Methods: The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called “biophysical” map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process. Results: Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements. Conclusions: A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-F-BRD-02: Application of ARCHERRT-- A GPU-Based Monte Carlo Dose Engine for Radiation Therapy -- to Tomotherapy and Patient-Independent IMRT

Purpose: As a module of ARCHER -- Accelerated Radiation-transport Computations in Heterogeneous EnviRonments, ARCHER{sub RT} is designed for RadioTherapy (RT) dose calculation. This paper describes the application of ARCHERRT on patient-dependent TomoTherapy and patient-independent IMRT. It also conducts a 'fair' comparison of different GPUs and multicore CPU. Methods: The source input used for patient-dependent TomoTherapy is phase space file (PSF) generated from optimized plan. For patient-independent IMRT, the open filed PSF is used for different cases. The intensity modulation is simulated by fluence map. The GEANT4 code is used as benchmark. DVH and gamma index test are employed to evaluate the accuracy of ARCHER{sub RT} code. Some previous studies reported misleading speedups by comparing GPU code with serial CPU code. To perform a fairer comparison, we write multi-thread code with OpenMP to fully exploit computing potential of CPU. The hardware involved in this study are a 6-core Intel E5-2620 CPU and 6 NVIDIA M2090 GPUs, a K20 GPU and a K40 GPU. Results: Dosimetric results from ARCHER{sub RT} and GEANT4 show good agreement. The 2%/2mm gamma test pass rates for different clinical cases are 97.2% to 99.7%. A single M2090 GPU needs 50~79 seconds for the simulation to achieve a statistical error of 1% in the PTV. The K40 card is about 1.7∼1.8 times faster than M2090 card. Using 6 M2090 card, the simulation can be finished in about 10 seconds. For comparison, Intel E5-2620 needs 507∼879 seconds for the same simulation. Conclusion: We successfully applied ARCHER{sub RT} to Tomotherapy and patient-independent IMRT, and conducted a fair comparison between GPU and CPU performance. The ARCHER{sub RT} code is both accurate and efficient and may be used towards clinical applications.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Action-based Dynamical Modeling for the Milky Way Disk: The Influence of Spiral Arms

RoadMapping is a dynamical modeling machinery developed to constrain the Milky Way’s (MW) gravitational potential by simultaneously fitting an axisymmetric parametrized potential and an action-based orbit distribution function (DF) to discrete 6D phase-space measurements of stars in the Galactic disk. In this work, we demonstrate RoadMapping's robustness in the presence of spiral arms by modeling data drawn from an N-body simulation snapshot of a disk-dominated galaxy of MW mass with strong spiral arms (but no bar), exploring survey volumes with radii 500 pc⩽r{sub max}⩽5 kpc. The potential constraints are very robust, even though we use a simple action-based DF, the quasi-isothermal DF. The best-fit RoadMapping model always recovers the correct gravitational forces where most of the stars that entered the analysis are located, even for small volumes. For data from large survey volumes, RoadMapping finds axisymmetric models that average well over the spiral arms. Unsurprisingly, the models are slightly biased by the excess of stars in the spiral arms. Gravitational potential models derived from survey volumes with at least r{sub max}=3 kpc can be reliably extrapolated to larger volumes. However, a large radial survey extent, r{sub max}∼5 kpc, is needed to correctly recover the halo scale length. In general, the recovery and extrapolability of potentials inferred from data sets that were drawn from inter-arm regions appear to be better than those of data sets drawn from spiral arms. Our analysis implies that building axisymmetric models for the Galaxy with upcoming Gaia data will lead to sensible and robust approximations of the MW’s potential.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Wave dispersion in the hybrid-Vlasov model: Verification of Vlasiator

Vlasiator is a new hybrid-Vlasov plasma simulation code aimed at simulating the entire magnetosphere of the Earth. The code treats ions (protons) kinetically through Vlasov's equation in the six-dimensional phase space while electrons are a massless charge-neutralizing fluid [M. Palmroth et al., J. Atmos. Sol.-Terr. Phys. 99, 41 (2013); A. Sandroos et al., Parallel Comput. 39, 306 (2013)]. For first global simulations of the magnetosphere, it is critical to verify and validate the model by established methods. Here, as part of the verification of Vlasiator, we characterize the low-β plasma wave modes described by this model and compare with the solution computed by the Waves in Homogeneous, Anisotropic Multicomponent Plasmas (WHAMP) code [K. Rönnmark, Kiruna Geophysical Institute Reports No. 179, 1982], using dispersion curves and surfaces produced with both programs. The match between the two fundamentally different approaches is excellent in the low-frequency, long wavelength range which is of interest in global magnetospheric simulations. The left-hand and right-hand polarized wave modes as well as the Bernstein modes in the Vlasiator simulations agree well with the WHAMP solutions. Vlasiator allows a direct investigation of the importance of the Hall term by including it in or excluding it from Ohm's law in simulations. This is illustrated showing examples of waves obtained using the ideal Ohm's law and Ohm's law including the Hall term. Our analysis emphasizes the role of the Hall term in Ohm's law in obtaining wave modes departing from ideal magnetohydrodynamics in the hybrid-Vlasov model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Signatures of discrete breathers in coherent state quantum dynamics

In classical mechanics, discrete breathers (DBs) - a spatial time-periodic localization of energy - are predicted in a large variety of nonlinear systems. Motivated by a conceptual bridging of the DB phenomena in classical and quantum mechanical representations, we study their signatures in the dynamics of a quantum equivalent of a classical mechanical point in phase space - a coherent state. In contrast to the classical point that exhibits either delocalized or localized motion, the coherent state shows signatures of both localized and delocalized behavior. The transition from normal to local modes have different characteristics in quantum and classical perspectives. Here, we get an insight into the connection between classical and quantum perspectives by analyzing the decomposition of the coherent state into system's eigenstates, and analyzing the spacial distribution of the wave-function density within these eigenstates. We find that the delocalized and localized eigenvalue components of the coherent state are separated by a mixed region, where both kinds of behavior can be observed. Further analysis leads to the following observations. Considered as a function of coupling, energy eigenstates go through avoided crossings between tunneling and non-tunneling modes. The dominance of tunneling modes in the high nonlinearity region is compromised by the appearance of new types of modes - high order tunneling modes - that are similar to the tunneling modes but have attributes of non-tunneling modes. Certain types of excitations preferentially excite higher order tunneling modes, allowing one to study their properties. Since auto-correlation functions decrease quickly in highly nonlinear systems, short-time dynamics are sufficient for modeling quantum DBs. This work provides a foundation for implementing modern semi-classical methods to model quantum DBs, bridging classical and quantum mechanical signatures of DBs, and understanding spectroscopic experiments that involve a coherent state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Spitzer Microlensing Parallax for OGLE-2017-BLG-0896 Reveals a Counter-rotating Low-mass Brown Dwarf

The kinematics of isolated brown dwarfs in the Galaxy, beyond the solar neighborhood, is virtually unknown. Microlensing has the potential to probe this hidden population, as it can measure both the mass and five of the six phase-space coordinates (all except the radial velocity) even of a dark isolated lens. However, the measurements of both the microlens-parallax and finite-source effects are needed in order to recover the full information. Here, we combine the Spitzer satellite parallax measurement with the ground-based light curve, which exhibits strong finite-source effects, of event OGLE-2017-BLG-0896. We find two degenerate solutions for the lens (due to the known satellite-parallax degeneracy), which are consistent with each other except for their proper motion. The lens is an isolated brown dwarf with a mass of either 18 ± 1 M {sub J} or 20 ± 1 M {sub J}. This is the lowest isolated-object mass measurement to date, only ∼45% more massive than the theoretical deuterium-fusion boundary at solar metallicity, which is the common definition of a free-floating planet. The brown dwarf is located at either 3.9 ± 0.1 kpc or 4.1 ± 0.1 kpc toward the Galactic bulge, but with proper motion in the opposite direction of disk stars, with one solution suggesting it is moving within the Galactic plane. While it is possibly a halo brown dwarf, it might also represent a different, unknown population.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

MULTI-DIMENSIONAL FEATURES OF NEUTRINO TRANSFER IN CORE-COLLAPSE SUPERNOVAE

We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six-dimensional phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolution to stationary states of neutrino distribution functions using our Boltzmann solver. Taking advantage of the multi-angle and multi-energy feature realized by the S {sub n} method in our code, we reveal the genuine characteristics of spatially three-dimensional neutrino transfer, such as nonradial fluxes and nondiagonal Eddington tensors. In addition, we assess the ray-by-ray approximation, turning off the lateral-transport terms in our code. We demonstrate that the ray-by-ray approximation tends to propagate fluctuations in thermodynamical states around the neutrino sphere along each radial ray and overestimate the variations between the neutrino distributions on different radial rays. We find that the difference in the densities and fluxes of neutrinos between the ray-by-ray approximation and the full Boltzmann transport becomes ∼20%, which is also the case for the local heating rate, whereas the volume-integrated heating rate in the Boltzmann transport is found to be only slightly larger (∼2%) than the counterpart in the ray-by-ray approximation due to cancellation among different rays. These results suggest that we should carefully assess the possible influences of various approximations in the neutrino transfer employed in current simulations of supernova dynamics. Detailed information on the angle and energy moments of neutrino distribution functions will be profitable for the future development of numerical methods in neutrino-radiation hydrodynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

SU-E-T-578: MCEBRT, A Monte Carlo Code for External Beam Treatment Plan Verifications

Purpose: Present a new Monte Carlo code (MCEBRT) for patient-specific dose calculations in external beam radiotherapy. The code MLC model is benchmarked and real patient plans are re-calculated using MCEBRT and compared with commercial TPS. Methods: MCEBRT is based on the GEPTS system (Med. Phys. 29 (2002) 835–846). Phase space data generated for Varian linac photon beams (6 – 15 MV) are used as source term. MCEBRT uses a realistic MLC model (tongue and groove, rounded ends). Patient CT and DICOM RT files are used to generate a 3D patient phantom and simulate the treatment configuration (gantry, collimator and couch angles; jaw positions; MLC sequences; MUs). MCEBRT dose distributions and DVHs are compared with those from TPS in absolute way (Gy). Results: Calculations based on the developed MLC model closely matches transmission measurements (pin-point ionization chamber at selected positions and film for lateral dose profile). See Fig.1. Dose calculations for two clinical cases (whole brain irradiation with opposed beams and lung case with eight fields) are carried out and outcomes are compared with the Eclipse AAA algorithm. Good agreement is observed for the brain case (Figs 2-3) except at the surface where MCEBRT dose can be higher by 20%. This is due to better modeling of electron contamination by MCEBRT. For the lung case an overall good agreement (91% gamma index passing rate with 3%/3mm DTA criterion) is observed (Fig.4) but dose in lung can be over-estimated by up to 10% by AAA (Fig.5). CTV and PTV DVHs from TPS and MCEBRT are nevertheless close (Fig.6). Conclusion: A new Monte Carlo code is developed for plan verification. Contrary to phantombased QA measurements, MCEBRT simulate the exact patient geometry and tissue composition. MCEBRT can be used as extra verification layer for plans where surface dose and tissue heterogeneity are an issue.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-E-J-09: A Monte Carlo Analysis of the Relationship Between Cherenkov Light Emission and Dose for Electrons, Protons, and X-Ray Photons

Purpose: A number of recent studies have proposed that light emitted by the Cherenkov effect may be used for a number of radiation therapy dosimetry applications. Here we investigate the fundamental nature and accuracy of the technique for the first time by using a theoretical and Monte Carlo based analysis. Methods: Using the GEANT4 architecture for medically-oriented simulations (GAMOS) and BEAMnrc for phase space file generation, the light yield, material variability, field size and energy dependence, and overall agreement between the Cherenkov light emission and dose deposition for electron, proton, and flattened, unflattened, and parallel opposed x-ray photon beams was explored. Results: Due to the exponential attenuation of x-ray photons, Cherenkov light emission and dose deposition were identical for monoenergetic pencil beams. However, polyenergetic beams exhibited errors with depth due to beam hardening, with the error being inversely related to beam energy. For finite field sizes, the error with depth was inversely proportional to field size, and lateral errors in the umbra were greater for larger field sizes. For opposed beams, the technique was most accurate due to an averaging out of beam hardening in a single beam. The technique was found to be not suitable for measuring electron beams, except for relative dosimetry of a plane at a single depth. Due to a lack of light emission, the technique was found to be unsuitable for proton beams. Conclusions: The results from this exploratory study suggest that optical dosimetry by the Cherenkov effect may be most applicable to near monoenergetic x-ray photon beams (e.g. Co-60), dynamic IMRT and VMAT plans, as well as narrow beams used for SRT and SRS. For electron beams, the technique would be best suited for superficial dosimetry, and for protons the technique is not applicable due to a lack of light emission. NIH R01CA109558 and R21EB017559.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

Rapidly evolving and luminous transients from Pan-STARRS1

In the past decade, several rapidly evolving transients have been discovered whose timescales and luminosities are not easily explained by traditional supernovae (SNe) models. The sample size of these objects has remained small due, at least in part, to the challenges of detecting short timescale transients with traditional survey cadences. Here we present the results from a search within the Pan-STARRS1 Medium Deep Survey (PS1-MDS) for rapidly evolving and luminous transients. We identify 10 new transients with a time above half-maximum (t {sub 1/2}) of less than 12 days and –16.5 > M > –20 mag. This increases the number of known events in this region of SN phase space by roughly a factor of three. The median redshift of the PS1-MDS sample is z = 0.275 and they all exploded in star-forming galaxies. In general, the transients possess faster rise than decline timescale and blue colors at maximum light (g {sub P1} – r {sub P1} ≲ –0.2). Best-fit blackbodies reveal photospheric temperatures/radii that expand/cool with time and explosion spectra taken near maximum light are dominated by a blue continuum, consistent with a hot, optically thick, ejecta. We find it difficult to reconcile the short timescale, high peak luminosity (L > 10{sup 43} erg s{sup –1}), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of {sup 56}Ni. Rather, we find that many are consistent with either (1) cooling envelope emission from the explosion of a star with a low-mass extended envelope that ejected very little (<0.03 M {sub ☉}) radioactive material, or (2) a shock breakout within a dense, optically thick, wind surrounding the progenitor star. After calculating the detection efficiency for objects with rapid timescales in the PS1-MDS we find a volumetric rate of 4800-8000 events yr{sup –1} Gpc{sup –3} (4%-7% of the core-collapse SN rate at z = 0.2).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

A Grid Intensity Map-Based Source Model for Monte Carlo Modeling of IMRT Field Irradiation - Paper 5

An accurate source model is important for TPS (Treatment Planning System) dose calculation and verification, especially for IMRT (Intensity Modulated Radiotherapy) irregular and inhomogeneous fields. Monte Carlo simulation is the most accurate way to calculate dose distributions in the heterogeneous media like human bodies. A full radiotherapy modeling generally includes Linac modeling and phantom dose calculation. Full Linac modeling is very time-consuming and requires highly detailed knowledge of the accelerator composition and structure, however, which is very difficult to obtain for commercial interests. Often all relevant parameters are not known with sufficient accuracy, and thus a time-consuming tuning of e.g. the electron beam spot size and energy must be performed. A multiple source model (MSM) regards particles coming from different parts of an accelerator as if they are coming from different virtual sources. MSM uses a set of virtual sources representing the treatment beam (e.g. primary photons, head scatter, electron contamination) by a set of analytic functions or histograms with parameters derived from dosimetric measurements or a full MC simulation. Tillikainen et al built a three sub-source model by matching with measurement data using an automatic process to optimize the free model parameters for a Varian 2100 C/D accelerator. Lin et al have built a convenient five-photon-sub-source model for Varian 2300C 6 MV-X beam based on Monte Carlo simulated Linac for ARTS (Accurate Radiotherapy System). The particle generation of MSM is very fast compared to full Linac simulation, and its sampling efficiency is not limited by the size of the phase space file (PSF). However, generally MSM is generated based on regular field PSF and needs an irregular component to form the beam shape meanwhile apt to bringing error Linac modeling and MSM producing from PSF analysis actually both are multi-freedom optimization process by adjusting several parameters (e.g. electron energy and angular divergence). Linac modeling adjusts its parameters from the beginning of electron hitting target through all components to fit the simulation to the measurement. MSM tunes its PSF-derived parameters by optimizing to minimize the deviations between the simulation and the measurement. Since the full Linac modeling and the MSM model both need to fine tune their parameters to match the simulation to the measurement, and sometimes maybe offsets the real accelerator for its multi-parameter optimization process. So, why not fine tune directly all source model parameters just behind all Linac components to fit the simulation to the measurement? A grid intensity-based surface source model was suggested here, which assumed the particle emitting from the bottom surface of Linac last beam trimmer MLC (Multi Leaf Collimator). The particle emitting position is sampled on the grid intensity map optimized by TPS, and the particle weight is adjusted by its original grid intensity. This source model not only takes advantage of the Monte Carlo method in the heterogeneous media, but also avoids effectively the transportation speed shortcoming of Monte Carlo simulation through Linac components. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2014

HDRMC, an accelerated Monte Carlo dose calculator for high dose rate brachytherapy with CT-compatible applicators

Purpose: To present a new accelerated Monte Carlo code for CT-based dose calculations in high dose rate (HDR) brachytherapy. The new code (HDRMC) accounts for both tissue and nontissue heterogeneities (applicator and contrast medium). Methods: HDRMC uses a fast ray-tracing technique and detailed physics algorithms to transport photons through a 3D mesh of voxels representing the patient anatomy with applicator and contrast medium included. A precalculated phase space file for the{sup 192}Ir source is used as source term. HDRM is calibrated to calculated absolute dose for real plans. A postprocessing technique is used to include the exact density and composition of nontissue heterogeneities in the 3D phantom. Dwell positions and angular orientations of the source are reconstructed using data from the treatment planning system (TPS). Structure contours are also imported from the TPS to recalculate dose-volume histograms. Results: HDRMC was first benchmarked against the MCNP5 code for a single source in homogenous water and for a loaded gynecologic applicator in water. The accuracy of the voxel-based applicator model used in HDRMC was also verified by comparing 3D dose distributions and dose-volume parameters obtained using 1-mm{sup 3} versus 2-mm{sup 3} phantom resolutions. HDRMC can calculate the 3D dose distribution for a typical HDR cervix case with 2-mm resolution in 5 min on a single CPU. Examples of heterogeneity effects for two clinical cases (cervix and esophagus) were demonstrated using HDRMC. The neglect of tissue heterogeneity for the esophageal case leads to the overestimate of CTV D90, CTV D100, and spinal cord maximum dose by 3.2%, 3.9%, and 3.6%, respectively. Conclusions: A fast Monte Carlo code for CT-based dose calculations which does not require a prebuilt applicator model is developed for those HDR brachytherapy treatments that use CT-compatible applicators. Tissue and nontissue heterogeneities should be taken into account in modern HDR brachytherapy planning.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

MO-FG-CAMPUS-TeP3-03: Calculation of Proton Pencil Beam Properties with Full Beamline Model in TOPAS

Purpose: Introducing Monte Carlo based dose calculation algorithms into proton therapy planning systems (TPS) leads to improved accuracy. However accurate modelling of the proton pencil beam impinging the patient is necessary. Current approaches rely on measurement-driven reconstruction of phase-space and spectrum properties, typically constrained to analytical model functions. In this study a detailed Monte Carlo model of the complete cyclotron-based delivery system was created with the aim of providing more representative beam properties at treatment position. Methods: A model of the Varian Probeam proton system from the cyclotron exit to isocenter was constructed in the TOPAS Monte Carlo framework. The beam evolution through apertures and magnetic elements was validated using Transport/Turtle calculations and additionally against measurements from the Probeam™ system at Scripps Proton Therapy Center (SPTC) in San Diego, CA. A voxelized water phantom at isocenter allowed for comparison of the dose-depth curve from the Probeam model with that of a corresponding Gaussian beam over the entire energy range (70–240 MeV). Measurements of relative beam fluence cross-profiles and depth-dose curves at and around isocenter were also compared to the MC results. Results: The simulated TOPAS beam envelope was found to agree with both the Transport/Turtle and measurements to within 5% for most of the beamline. The MC predicted energy spectrum at isocenter was found to deviate increasingly from Gaussian at energies below 160 MeV. The corresponding effects on the depth dose curve agreed well with measurements. Conclusion: Given the flexibility of TOPAS and available details of the delivery system, an accurate characterization of a proton pencil beam at isocenter is possible. Incorporation of the MC derived properties of the proton pencil beam can eliminate analytical approximations and ultimately increase treatment plan accuracy and quality. Both authors are employees of Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗