The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

The effect of organ motion of radiotherapy plans

Organ motion as a result of respiration is an important field of research for medical physics. Knowledge of magnitude and direction of this motion is necessary to allow for more accurate radiotherapy treatment planning. This will result in higher doses to the tumour whilst sparing healthy tissue. This project involved human trials, where the radiation therapy patient's kidneys were CT scanned under three different conditions; whilst free breathing (FB), breath-hold at normal tidal inspiration (BHIN), and breath-hold at normal tidal expiration (BHEX). The magnitude of motion was measured by recording the outline of the kidney from a Beam's Eye View (BEV). The centre of mass of this 2D shape was calculated for each set using "ImageJ" software and the magnitude of movement determined from the change in the centroid's coordinates between the BHIN and BHEX scans. The movement ranged from, for the left and right kidneys, 4-46mm and 2-44mm in the superior/inferior (axial) plane, 1-21mm and 2- 16mm in the anterior/posterior (coronal) plane, and 0-6mm and 0-8mm in the lateral/medial (sagittal) plane. From exhale to inhale, the kidneys tended to move inferiorly, anteriorly and laterally. A standard radiotherapy plan, designed to treat the para-aortics with opposed lateral fields was performed on the free breathing (planning) CT set. The field size and arrangement was set up using the same parameters for each subject. The prescription was to deliver 45 Gray in 25 fractions. This field arrangement and prescription was then copied over to the breath hold CT sets, and the dosimetric differences were compared using Dose Volume Histograms (DVH). The point of comparison for the three sets was recorded as the percentage volume of kidney receiving less than or equal to 10 Gray. The QUASAR respiratory motion phantom was used with the range of motion determined from the human study. The phantom was imaged, planned and treated with a linear accelerator with dose determined by film. The effect of the motion was measured by the change in the penumbra of the film and compared to the penumbra from the treatment planning system.

The complete mitochondrial genome of the yellow coaster, Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini): sequence, gene organization and a unique tRNA translocation event

In this paper, the complete mitochondrial genome of Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini) is reported; a circular molecule of 15,245 bp in size. For A. issoria, genes are arranged in the same order and orientation as the complete sequenced mitochondrial genomes of the other lepidopteran species, except for the presence of an extra copy of tRNAIle(AUR)b in the control region. All protein-coding genes of A. issoria mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that COI gene uses TTG as its initial codon and terminates in a single T residue. All tRNA genes possess the typical clover leaf secondary structure except for tRNASer(AGN), which has a simple loop with the absence of the DHU stem. The sequence, organization and other features including nucleotide composition and codon usage of this mitochondrial genome were also reported and compared with those of other sequenced lepidopterans mitochondrial genomes. There are some short microsatellite-like repeat regions (e.g., (TA)9, polyA and polyT) scattered in the control region, however, the conspicuous macro-repeats units commonly found in other insect species are absent.

The complete mitochondrial genome of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae)

We determined the nucleotide sequence of the mitochondrial genome (mtgenome) of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae). The entire closed circular molecule is 15,368 bp and contains 37 genes with the typical gene complement and order for lepidopteran mtgenomes. All tRNAs except tRNASer(AGN) can be folded into the typical cloverleaf secondary structures. The protein-coding genes (PCGs) have typical mitochondrial start codons, with the exception of COI, which uses the unusual CGA one as is found in all other Lepidoptera sequenced to date. In addition, six of 13 PCGs harbor the incomplete termination codons, a single T. The A+T-rich region contains some conserved structures that are similar to those found in other lepidopteran mtgenomes, including a structure combining the motif 'ATAGA', a 19-bp poly(T) stretch and three microsatellite (AT)n elements which are part of larger 122+ bp macrorepeats. This is the first report of macrorepeats in a lepidopteran mtgenome.

A comparative analysis of Mitochondrial genomes in Coleoptera (Arthropoda: Insecta) and genome descriptions of six new beetles

Coleoptera is the most diverse group of insects with over 360,000 described species divided into four suborders: Adephaga, Archostemata, Myxophaga, and Polyphaga. In this study, we present six new complete mitochondrial genome (mtgenome) descriptions, including a representative of each suborder, and analyze the evolution of mtgenomes from a comparative framework using all available coleopteran mtgenomes. We propose a modification of atypical cox1 start codons based on sequence alignment to better reflect the conservation observed across species as well as findings of TTG start codons in other genes. We also analyze tRNA-Ser(AGN) anticodons, usually GCU in arthropods, and report a conserved UCU anticodon as a possible synapomorphy across Polyphaga. We further analyze the secondary structure of tRNA-Ser(AGN) and present a consensus structure and an updated covariance model that allows tRNAscan-SE (via the COVE software package) to locate and fold these atypical tRNAs with much greater consistency. We also report secondary structure predictions for both rRNA genes based on conserved stems. All six species of beetle have the same gene order as the ancestral insect. We report noncoding DNA regions, including a small gap region of about 20 bp between tRNA-Ser(UCN) and nad1 that is present in all six genomes, and present results of a base composition analysis.

Starburst in the ultraluminous galaxy Arp 220: Constraints from observations of radio recombination lines and continuum

We present observations of radio recombination lines (RRL) from the starburst galaxy Arp 220 at 8.1 GHz (H92 alpha) and 1.4 GHz (H167 alpha and H165 alpha) and at 84 GHz (H42 alpha), 96 GHz (H40 alpha) and 207 GHz (H31 alpha) using the Very Large Array and the IRAM 30 m telescope, respectively. RRLs were detected at all the frequencies except 1.4 GHz, where a sensitive upper limit was obtained. We also present continuum flux measurements at these frequencies as well as at 327 MHz made with the VLA. The continuum spectrum, which has a spectral index alpha similar to -0.6 (S-nu proportional to nu(alpha)) between 5 and 10 GHz, shows a break near 1.5 GHz, a prominent turnover below 500 MHz, and a flatter spectral index above 50 GHz. We show that a model with three components of ionized gas with different densities and area covering factors can consistently explain both RRL and continuum data. The total mass of ionized gas in the three components is 3.2 x 10(7) M., requiring 3 x 10(5) O5 stars with a total Lyman continuum production rate N-Lyc similar to 1.3 x 10(55) photons s(-1). The ratio of the expected to observed Br alpha and Br gamma fluxes implies a dust extinction A(V) similar to 45 mag. The derived Lyman continuum photon production rate implies a continuous star formation rate (SFR) averaged over the lifetime of OB stars of similar to 240 M yr(-1). The Lyman continuum photon Production rate of similar to 3% associated with the high-density H II regions implies a similar SFR at recent epochs (t < 10(5) yr). An alternative model of high-density gas, which cannot be excluded on the basis of the available data, predicts 10 times higher SFR at recent epochs. If confirmed, this model implies that star formation in Arp 220 consists of multiple starbursts of very high SFR (few times 10(3) M. yr(-1)) and short duration (similar to 10(5) yr). The similarity of IR excess, L-IR/L-Ly alpha similar to 24, in Arp 220 to values observed in starburst galaxies shows that most of the high luminosity of Arp 220 is due to the ongoing starburst rather than to a hidden active galactic nucleus (AGN). A comparison of the IR excesses in Arp 220, the Galaxy, and M33 indicates that the starburst in Arp 220 has an initial mass function that is similar to that in normal galaxies and has a duration longer than 107 yr. If there was no infall of gas during this period, then the star formation efficiency (SFE) in Arp 220 is similar to 50%. The high SFR and SFE in Arp 220 is consistent with their known dependences on mass and density of gas in star-forming regions of normal galaxies.

Effect of a QSO on the interstellar medium of its host galaxy

We consider the hydrodynamic evolution of gas in the interstellar medium of the host galaxy of a quasar due to Compton heating by the QSO radiation. We show that a Lagrangean formulation of the problem is necessary. It is found that the "hydrodynamic time scale" becomes important compared to the Compton heating time scale. We also relax the "single fluid" approximation by considering the existence of clouds and taking into account the mass loss from stars. The results predict star burst activity, and thus we explain the blue colors of the active galaxies.

The modulation of radiation in an electron-positron plasma

The modulational instability of a large-amplitude, linearly polarized electromagnetic wave propagating in an electron-positron plasma is considered, including the combined effect of relativistic mass variation of the plasma particles, harmonic generation, and the non-resonant, finite-frequency electrostatic density perturbations, all caused by the large-amplitude radiation field. The radiation from many strong sources, such as AGN and pulsars, has been observed to vary over a host of time-scales. It is possible that the extremely rapid variations in the non-thermal continuum of AGN, as well as in the non-thermal radio radiation from pulsars, can be accounted for by the modulational instabilities to which radiation may be subjected during its propagation out of the emission region.

Very large array detection of radio recombination lines from the radio nucleus of NGC 253: Ionization by a weak active galactic nucleus, an obscured super star cluster, or a compact supernova remnant?

We have imaged the H92alpha and H75alpha radio recombination line (RRL) emissions from the starburst galaxy NGC 253 with a resolution of similar to4 pc. The peak of the RRL emission at both frequencies coincides with the unresolved radio nucleus. Both lines observed toward the nucleus are extremely wide, with FWHMs of similar to200 km s(-1). Modeling the RRL and radio continuum data for the radio nucleus shows that the lines arise in gas whose density is similar to10(4) cm(-3) and mass is a few thousand M., which requires an ionizing flux of (6-20) x 10(51) photons s(-1). We consider a supernova remnant (SNR) expanding in a dense medium, a star cluster, and also an active galactic nucleus (AGN) as potential ionizing sources. Based on dynamical arguments, we rule out an SNR as a viable ionizing source. A star cluster model is considered, and the dynamics of the ionized gas in a stellar-wind driven structure are investigated. Such a model is only consistent with the properties of the ionized gas for a cluster younger than similar to10(5) yr. The existence of such a young cluster at the nucleus seems improbable. The third model assumes the ionizing source to be an AGN at the nucleus. In this model, it is shown that the observed X-ray flux is too weak to account for the required ionizing photon flux. However, the ionization requirement can be explained if the accretion disk is assumed to have a big blue bump in its spectrum. Hence, we favor an AGN at the nucleus as the source responsible for ionizing the observed RRLs. A hybrid model consisting of an inner advection-dominated accretion flow disk and an outer thin disk is suggested, which could explain the radio, UV, and X-ray luminosities of the nucleus.

Isotropic heating of galaxy cluster cores via rapidly reorienting active galactic nucleus jets

Active galactic nucleus (AGN) jets carry more than sufficient energy to stave off catastrophic cooling of the intracluster medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, the ICM must be heated in a near-isotropic fashion and narrow bipolar jets with P-jet = 10(44-45) erg s(-1), typical of radio AGNs at cluster centers, are inefficient in heating the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster cores, the supermassive black holes (SMBHs) will, in addition to accreting gas via radiatively inefficient flows, experience short stochastic episodes of enhanced accretion via thin disks. In general, the orientation of these accretion disks will be misaligned with the spin axis of the black holes (BHs) and the ensuing torques will cause the BH's spin axis (and therefore the jet axis) to slew and rapidly change direction. This model not only explains recent observations showing successive generations of jet-lobes-bubbles in individual cool-core clusters that are offset from each other in the angular direction with respect to the cluster center, but also shows that AGN jets can heat the cluster core nearly isotropically on the gas cooling timescale. Our model does require that the SMBHs at the centers of cool-core clusters be spinning relatively slowly. Torques from individual misaligned disks are ineffective at tilting rapidly spinning BHs by more than a few degrees. Additionally, since SMBHs that host thin accretion disks will manifest as quasars, we predict that roughly 1-2 rich clusters within z < 0.5 should have quasars at their centers.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

Lurking in ULIRGs: Molecular Gas in local Merging Galaxies

Mergers and interacting galaxies are pivotal to the evolution of galaxies in the universe. They are the sites of prodigious star formation and key to understanding the starburst processes: the physical and chemical properties and the dynamics of the molecular gas. ULIRGs or Ultraluminous Infrared Galaxies are a result of many of these mergers. They host extreme starbursts, AGNs, and mergers. They are the perfect laboratory to probe the connection between starbursts, black hole accretion and mergers and to further our understanding of star formation and merging.

NGC 6240 and Arp 220 can be considered the founding members of this very active class of objects. They are in different stages of merging and hence are excellent case studies to further our understanding about the merging process. We have imaged the dense star-forming regions of these galaxies at sub-arcsec resolution with CARMA C and B Configurations (2" and 0.5 - 0.8"). Multi-band imaging allows excitation analysis of HCN, HCO⁺, HNC, and CS along with CO transitions to constrain the properties of the gas. Our dataset is unique in that we have observed these lines at similar resolutions and high sensitivity which can be used to derive line ratios of faint high excitation lines.

Arp 220 has not had confirmed X-ray AGN detections for either nuclei. However, our observations indicate HCN/HNC ratios consistent with the chemistry of X-ray Dominated Regions (XDRs) -- a likely symptom of AGN. We calculated the molecular Hydrogen densities using each of the molecular species and conclude that assuming abundances of HNC and HCO⁺ similar to those in galactic sources are incorrect in the case of ULIRGs. The physical conditions in the dense molecular gas in ULIRGs alter these abundances. The derived H2 volume densities are ~ 5 x 10⁴ cm^-3 in both Arp 220 nuclei and ~ 10⁴ cm^-3 in NGC 6240.

The Kiloparsec-scale structure of 3C 286

We present radio images of the compact steep spectrum (CSS) quasar 3C 286 acquired with the Very Large Array (VLA) at 8.4 and 22.5 GHz. The source exhibits a two-sided core-jet structure with a bright central component and two extended components one to the east (P.A. 100degrees) and another to the southwest (P.A. -116degrees). From the compact core, an extension runs towards the southwest component up to similar to 0.7 arcsecond. The emission between the primary central component and the southwest component exhibits a knotty structure. A gradual change of the jet position angles from -135degrees to -120degrees in the inner southwest jet suggests a local bend. The position angle changes of the major eastern components E2 and E1 suggest that the eastern jet likely follows a curved trace. The bends in the jet trace may be associated with a relativistic precession or some interaction between the jet and the ambient matter. A mean spectral index of alpha(8.4)(22.5) similar to -0.76 (S-nu proportional to nu(alpha)) is estimated for the core component. Steep spectra are also obtained for the extended southwest component (2.6", P.A. -116degrees) and eastern component (0.8", P.A. 100degrees), with alpha(8.4)(22.5) similar to -0.88 and alpha(8.4)(22.5) similar to -1.79, respectively. The radio morphologies and spectral index distributions suggest that the core seen in our images is likely to be the beamed inner jet while the real nucleus is dimmed by it beaming away from us.

The radio counter-jet of the QSO 3C 48

We present multi- frequency radio observational results of the quasar 3C 48. The observations were carried out with the Very Large Array ( VLA) at five frequencies, 0.33, 1.5, 4.8, 8.4, and 22.5 GHz, and with the Multi- Element Radio Linked Interferometer Network ( MERLIN) at the two frequencies of 1.6 and 5 GHz. The source shows a one- sided jet to the north within 1", which then extends to the northeast and becomes diffuse. Two bright components ( N2 and N3), containing most of the flux density, are present in the northern jet. The spectral index of the two components is alpha(N2) similar to -0.99 +/- 0.12 and alpha(N3) similar to - 0.84 +/- 0.23 ( S proportional to nu(alpha)). Our images show the presence of an extended structure surrounding component N2, suggestive of strong interaction between the jet and the interstellar medium ( ISM) of the host galaxy. A steep- spectrum component, labelled S, located 0.25 " southwest to the flat- spectrum component which could be the core of 3C 48, is detected at a significance of > 15 sigma. Both the location and the steepness of the spectrum of component S suggest the presence of a counter- jet in 3C 48.

The development of new methods for high resolution radio astronomy imaging

Very Long Baseline Interferometry (VLBI) polarisation observations of the relativistic jets from Active Galactic Nuclei (AGN) allow the magnetic field environment around the jet to be probed. In particular, multi-wavelength observations of AGN jets allow the creation of Faraday rotation measure maps which can be used to gain an insight into the magnetic field component of the jet along the line of sight. Recent polarisation and Faraday rotation measure maps of many AGN show possible evidence for the presence of helical magnetic fields. The detection of such evidence is highly dependent both on the resolution of the images and the quality of the error analysis and statistics used in the detection. This thesis focuses on the development of new methods for high resolution radio astronomy imaging in both of these areas. An implementation of the Maximum Entropy Method (MEM) suitable for multi-wavelength VLBI polarisation observations is presented and the advantage in resolution it possesses over the CLEAN algorithm is discussed and demonstrated using Monte Carlo simulations. This new polarisation MEM code has been applied to multi-wavelength imaging of the Active Galactic Nuclei 0716+714, Mrk 501 and 1633+382, in each case providing improved polarisation imaging compared to the case of deconvolution using the standard CLEAN algorithm. The first MEM-based fractional polarisation and Faraday-rotation VLBI images are presented, using these sources as examples. Recent detections of gradients in Faraday rotation measure are presented, including an observation of a reversal in the direction of a gradient further along a jet. Simulated observations confirming the observability of such a phenomenon are conducted, and possible explanations for a reversal in the direction of the Faraday rotation measure gradient are discussed. These results were originally published in Mahmud et al. (2013). Finally, a new error model for the CLEAN algorithm is developed which takes into account correlation between neighbouring pixels. Comparison of error maps calculated using this new model and Monte Carlo maps show striking similarities when the sources considered are well resolved, indicating that the method is correctly reproducing at least some component of the overall uncertainty in the images. The calculation of many useful quantities using this model is demonstrated and the advantages it poses over traditional single pixel calculations is illustrated. The limitations of the model as revealed by Monte Carlo simulations are also discussed; unfortunately, the error model does not work well when applied to compact regions of emission.