Zdenka Kuncic

Theory of type II radio emission from the foreshock of an interplanetary shock

We present an analytical model for type II solar radio bursts and then apply it to an observed type II event. Electron beams are produced in the foreshock of an interplanetary shock via shock drift acceleration. Reflection is treated in the de Hoffman-Teller frame with efficiencies modeled by a losscone that incorporates the effects of the static cross-shock potential ϕ. Stochastic growth theory is used to treat electron beam driven Langmuir wave growth in the type II foreshock. Nonlinear wave-wave interactions are used as the mechanisms for converting Langmuir wave energy into freely propagating radio emission. The electron beams produced in the foreshock have a wide range of speeds and number densities. These electron beams are qualitatively consistent with observations in a type II foreshock as well as earlier theoretical predictions, and observations in Earths foreshock. Significant levels of Langmuir waves and ƒp and 2ƒp emission are predicted. In particular, the predicted volume emissivities are similar to those predicted for type III bursts. The simple model developed for the source environment of the type II event on August 26, 1998, produces fluxes in reasonable agreement with observation.

Kilovoltage Intrafraction Monitoring for Prostate Intensity Modulated Arc Therapy: First Clinical Results

PURPOSE Most linear accelerators purchased today are equipped with a gantry-mounted kilovoltage X-ray imager which is typically used for patient imaging prior to therapy. A novel application of the X-ray system is kilovoltage intrafraction monitoring (KIM), in which the 3-dimensional (3D) tumor position is determined during treatment. In this paper, we report on the first use of KIM in a prospective clinical study of prostate cancer patients undergoing intensity modulated arc therapy (IMAT). METHODS AND MATERIALS Ten prostate cancer patients with implanted fiducial markers undergoing conventionally fractionated IMAT (RapidArc) were enrolled in an ethics-approved study of KIM. KIM involves acquiring kV images as the gantry rotates around the patient during treatment. Post-treatment, markers in these images were segmented to obtain 2D positions. From the 2D positions, a maximum likelihood estimation of a probability density function was used to obtain 3D prostate trajectories. The trajectories were analyzed to determine the motion type and the percentage of time the prostate was displaced ≥ 3, 5, 7, and 10 mm. Independent verification of KIM positional accuracy was performed using kV/MV triangulation. RESULTS KIM was performed for 268 fractions. Various prostate trajectories were observed (ie, continuous target drift, transient excursion, stable target position, persistent excursion, high-frequency excursions, and erratic behavior). For all patients, 3D displacements of ≥ 3, 5, 7, and 10 mm were observed 5.6%, 2.2%, 0.7% and 0.4% of the time, respectively. The average systematic accuracy of KIM was measured at 0.46 mm. CONCLUSIONS KIM for prostate IMAT was successfully implemented clinically for the first time. Key advantages of this method are (1) submillimeter accuracy, (2) widespread applicability, and (3) a low barrier to clinical implementation. A disadvantage is that KIM delivers additional imaging dose to the patient.

Dynamics and Energetics of Turbulent, Magnetized Disk Accretion around Black Holes: A First-Principles Approach to Disk-Corona-Outflow Coupling

We present an analytic description of turbulent, magnetohydrodynamic (MHD) disk accretion around black holes that specifically addresses the relationship between radial and vertical mean field transport of mass, momentum, and energy, thereby complementing and extending numerical simulations. The azimuthal-vertical component of the magnetic stress is fundamental to an understanding of disk-corona-outflow coupling: when it is important for driving the angular momentum transport and mass accretion in the disk, it also has an important influence on the disk-corona-outflow energy budget. The Poynting flux derived from the product of this term with the Keplerian velocity also dominates the Poynting flux into the corona. The ratio of the coronal Alfven velocity to the Keplerian velocity is an important parameter in disk-corona-outflow physics. If this parameter is greater than unity, then energetically significant winds and Poynting flux into the corona occur. However, significant effects could also occur when this parameter is much less than unity. A limiting solution describing the case of angular momentum transport solely by the vertical-azimuthal stress has the property that all of the accretion power is channeled into a wind, some of which would be dissipated in the corona. More realistic solutions in which there is both radial and vertical transport of angular momentum would have different fractions of the accretion power emitted by the disk and corona, respectively. These results have important implications for existing accretion disk theory and for our interpretation of high-energy emission and nuclear outflows from the central engines of active galactic nuclei and galactic black hole candidates.

Theoretically predicted properties of type II radio emission from an interplanetary foreshock

[1] We investigate the causes of variations in interplanetary type II radio bursts using an analytic model which predicts the emission generated by electron beams in the foreshock regions upstream of an interplanetary shock wave. Trends in source region characteristics and remotely observed radio fluxes are calculated as a function of a number of solar wind and shock parameters. Calculations are performed for a single three-dimensional ripple on the global shock surface. Radio-loud ripples are predicted to have larger shock speeds relative to the solar wind speed, higher levels of nonthermal electrons, larger radii of curvature, and be moving through higher density regions than radio-quiet ripples. These predictions are qualitatively consistent with available observations. The predicted emission depends most sensitively on the speed of the shock relative to the solar wind. Strong correlations are found between the intensity of fundamental emission and the level of nonthermal electrons present in the tail of the incident solar wind electron distribution. Harmonic emission is found to be most sensitive to variations in the electron temperature T e of the incident solar wind. These results indicate that the bursty nature of typical type II observations can be accounted for by a shock propagating through an inhomogeneous solar wind.

The water equivalence of solid phantoms for low energy photon beams.

PURPOSE To compare and evaluate the dosimetric water equivalence of several commonly used solid phantoms for low energy photon beams. METHODS A total of ten different solid phantom materials was used in the study. The PENELOPE Monte Carlo code was used to calculate depth doses and beam profiles in all the phantom materials as well as the dose to a small water voxel at the surface of the solid phantom. These doses were compared to the corresponding doses calculated in a water phantom. The primary photon beams used ranged in energy from 50 to 280 kVp. RESULTS A number of phantom materials had excellent agreement in dose compared to water for all the x-ray beam energies studied. RMI457 Solid Water, Virtual Water, PAGAT, A150, and Plastic Water DT all had depth doses that agreed with those in water to within 2%. For these same phantom materials, the dose changes in the water voxel at the surface of the solid phantom were within 2%, except for A150, which agreed to within 2.7%. By comparison, the largest differences in depth doses occurred for Plastic Water (-21.7%) and polystyrene (17.6%) for the 50 kVp energy photon beam and 8 cm diameter field size. Plastic Water gave the largest difference in the normalized beam profiles with differences of up to 3.5% as compared to water. Surface dose changes, due to the presence of the solid phantom acting as the backscatter material, were found to be up to 9.1% for polystyrene with significant differences also found for Plastic Water, PMMA, and RW3 phantoms. CONCLUSIONS The following solid phantoms can be considered water equivalent and are recommended for relative dosimetry of low energy photon beams: A150, PAGAT, Plastic Water DT, RMI457 Solid Water, and Virtual Water. However, the following solid phantoms give significant differences, compared to water, in depth doses, profiles, and/or in surface doses due to backscatter changes: Plastic Water, PMMA, polystyrene, PRESAGE, and RW3.

Advances in kilovoltage x-ray beam dosimetry

This topical review provides an up-to-date overview of the theoretical and practical aspects of therapeutic kilovoltage x-ray beam dosimetry. Kilovoltage x-ray beams have the property that the maximum dose occurs very close to the surface and thus, they are predominantly used in the treatment of skin cancers but also have applications for the treatment of other cancers. In addition, kilovoltage x-ray beams are used in intra operative units, within animal irradiators and in on-board imagers on linear accelerators and kilovoltage dosimetry is important in these applications as well. This review covers both reference and relative dosimetry of kilovoltage x-ray beams and provides recommendations for clinical measurements based on the literature to date. In particular, practical aspects for the selection of dosimeter and phantom material are reviewed to provide suitable advice for medical physicists. An overview is also presented of dosimeters other than ionization chambers which can be used for both relative and in vivo dosimetry. Finally, issues related to the treatment planning and the use of Monte Carlo codes for solving radiation transport problems in kilovoltage x-ray beams are presented.

Investigation of radiological properties and water equivalency of PRESAGE ® dosimeters

PURPOSE PRESAGE is a dosimeter made of polyurethane, which is suitable for 3D dosimetry in modern radiation treatment techniques. Since an ideal dosimeter is radiologically water equivalent, the authors investigated water equivalency and the radiological properties of three different PRESAGE formulations that differ primarily in their elemental compositions. Two of the formulations are new and have lower halogen content than the original formulation. METHODS The radiological water equivalence was assessed by comparing the densities, interaction probabilities, and radiation dosimetry properties of the three different PRESAGE formulations to the corresponding values for water. The relative depth doses were calculated using Monte Carlo methods for 50, 100, 200, and 350 kVp and 6 MV x-ray beams. RESULTS The mass densities of the three PRESAGE formulations varied from 5.3% higher than that of water to as much as 10% higher than that of water for the original formulation. The probability of photoelectric absorption in the three different PRESAGE formulations varied from 2.2 times greater than that of water for the new formulations to 3.5 times greater than that of water for the original formulation. The mass attenuation coefficient for the three formulations is 12%-50% higher than the value for water. These differences occur over an energy range (10-100 keV) in which the photoelectric effect is the dominant interaction. The collision mass stopping powers of the relatively lower halogen-containing PRESAGE formulations also exhibit marginally better water equivalency than the original higher halogen-containing PRESAGE formulation. Furthermore, the depth dose curves for the lower halogen-containing PRESAGE formulations are slightly closer to that of water for a 6 MV beam. In the kilovoltage energy range, the depth dose curves for the lower halogen-containing PRESAGE formulations are in better agreement with water than the original PRESAGE formulation. CONCLUSIONS Based on the results of this study, the new PRESAGE formulations with lower halogen content are more radiologically water equivalent overall than the original formulation. This indicates that the new PRESAGE formulations are better suited to clinical applications and are more accurate dosimeters and phantoms than the original PRESAGE formulation. While correction factors are still needed to convert the dose measured by the dosimeter to an absorbed dose in water in the kilovoltage energy range, these correction factors are considerably smaller for the new PRESAGE formulations compared to the original PRESAGE and the existing polymer gel dosimeters.

X-ray polarization in relativistic jets

We investigate the polarization properties of Comptonized X-rays from relativistic jets in active galactic nuclei (AGN) using Monte Carlo simulations. We consider three scenarios commonly proposed for the observed X-ray emission in AGN: Compton scattering of blackbody photons emitted from an accretion disc; scattering of cosmic microwave background (CMB) photons and self-Comptonization of intrinsically polarized synchrotron photons emitted by jet electrons. Our simulations show that for Comptonization of disc and CMB photons, the degree of polarization of the scattered photons increases with the viewing inclination angle with respect to the jet axis. In both cases, the maximum linear polarization is ≈20 per cent. In the case of synchrotron self-Comptonization (SSC), we find that the resulting X-ray polarization depends strongly on the seed synchrotron photon injection site, with typical fractional polarizations P ≈ 10–20 per cent when synchrotron emission is localized near the jet base, while P ≈ 20–70 per cent for the case of uniform emission throughout the jet. These results indicate that X-ray polarimetry may be capable of providing unique clues to identify the location of particle acceleration sites in relativistic jets. In particular, if synchrotron photons are emitted quasiuniformly throughout a jet, then the observed degree of X-ray polarization may be sufficiently different for each of the competing X-ray emission mechanisms (synchrotron, SSC or external Comptonization) to determine which is the dominant process. However, X-ray polarimetry alone is unlikely to be able to distinguish between disc and CMB Comptonization.

A system for EPID-based real-time treatment delivery verification during dynamic IMRT treatment

PURPOSE To design and develop a real-time electronic portal imaging device (EPID)-based delivery verification system for dynamic intensity modulated radiation therapy (IMRT) which enables detection of gross treatment delivery errors before delivery of substantial radiation to the patient. METHODS The system utilizes a comprehensive physics-based model to generate a series of predicted transit EPID image frames as a reference dataset and compares these to measured EPID frames acquired during treatment. The two datasets are using MLC aperture comparison and cumulative signal checking techniques. The system operation in real-time was simulated offline using previously acquired images for 19 IMRT patient deliveries with both frame-by-frame comparison and cumulative frame comparison. Simulated error case studies were used to demonstrate the system sensitivity and performance. RESULTS The accuracy of the synchronization method was shown to agree within two control points which corresponds to approximately ∼1% of the total MU to be delivered for dynamic IMRT. The system achieved mean real-time gamma results for frame-by-frame analysis of 86.6% and 89.0% for 3%, 3 mm and 4%, 4 mm criteria, respectively, and 97.9% and 98.6% for cumulative gamma analysis. The system can detect a 10% MU error using 3%, 3 mm criteria within approximately 10 s. The EPID-based real-time delivery verification system successfully detected simulated gross errors introduced into patient plan deliveries in near real-time (within 0.1 s). CONCLUSIONS A real-time radiation delivery verification system for dynamic IMRT has been demonstrated that is designed to prevent major mistreatments in modern radiation therapy.

An investigation of backscatter factors for kilovoltage x-rays: a comparison between Monte Carlo simulations and Gafchromic EBT film measurements

Backscatter factors are important parameters in the determination of dose for kilovoltage x-ray beams. However, backscatter factors are difficult to measure experimentally, and tabulated values are based largely on Monte Carlo calculations. In this study we have determined new backscatter factors by both experimental and Monte Carlo methods, and compared them with existing backscatter factors published in the AAPM TG-61 protocol. The purpose of this study is twofold: (1) to evaluate the overall effectiveness of using Gafchromic EBT film for backscatter factor measurements and (2) to determine whether existing Monte Carlo-calculated backscatter factors need to be updated. We measured backscatter factors using Gafchromic EBT film for three field sizes (2, 4 and 6 cm diameter cones) and three kilovoltage beam qualities, including 280 kVp for which similar measurements have not previously been reported. We also present new Monte Carlo-calculated backscatter factors obtained using the EGSnrc/BEAMnrc code system to simulate the Pantak kilovoltage x-ray unit used in our measurements. The results were compared with backscatter factors tabulated in the AAPM TG-61 protocol for kilovoltage x-ray dosimetry. The largest difference between our measured and calculated backscatter factors and the AAPM TG-61 values was found to be 2.5%. This agreement is remarkably good, considering that the AAPM TG-61 values consist of a combination of experimental and Monte Carlo calculations obtained over 20 years ago using different measurement techniques, as well as older Monte Carlo code and cross-section data. Furthermore, our Monte Carlo-calculated backscatter factors agree within 1% with the AAPM TG-61 values for all beam qualities and field sizes. Our Gafchromic film measurements had slightly larger differences with the AAPM TG-61 backscatter factors, up to approximately 2% for the 6 cm diameter cone at a beam quality of 50 kVp. The largest difference in backscatter factors, of 2.5%, was found between Monte Carlo-calculated and Gafchromic film-measured data for the 100 kVp x-ray beam with the 4 cm diameter cone. The differences in backscatter factors between the three data sets (measurements, calculations and published values) are all within the uncertainties from our Gafchromic film measurements and Monte Carlo calculations. Our results demonstrate the suitability of using Gafchromic EBT film to measure equipment-specific backscatter factors for kilovoltage x-ray beams over the entire energy range and also confirm that backscatter factors published in kilovoltage dosimetry protocols still remain valid.