Effect of CT electron density identification on Monte Carlo simulations of radiotherapeutic treatments

Introduction: The accurate identification of tissue electron densities is of great importance for Monte Carlo (MC) dose calculations. When converting patient CT data into a voxelised format suitable for MC simulations, however, it is common to simplify the assignment of electron densities so that the complex tissues existing in the human body are categorized into a few basic types. This study examines the effects that the assignment of tissue types and the calculation of densities can have on the results of MC simulations, for the particular case of a Siemen’s Sensation 4 CT scanner located in a radiotherapy centre where QA measurements are routinely made using 11 tissue types (plus air). Methods: DOSXYZnrc phantoms are generated from CT data, using the CTCREATE user code, with the relationship between Hounsfield units (HU) and density determined via linear interpolation between a series of specified points on the ‘CT-density ramp’ (see Figure 1(a)). Tissue types are assigned according to HU ranges. Each voxel in the DOSXYZnrc phantom therefore has an electron density (electrons/cm3) defined by the product of the mass density (from the HU conversion) and the intrinsic electron density (electrons /gram) (from the material assignment), in that voxel. In this study, we consider the problems of density conversion and material identification separately: the CT-density ramp is simplified by decreasing the number of points which define it from 12 down to 8, 3 and 2; and the material-type-assignment is varied by defining the materials which comprise our test phantom (a Supertech head) as two tissues and bone, two plastics and bone, water only and (as an extreme case) lead only. The effect of these parameters on radiological thickness maps derived from simulated portal images is investigated. Results & Discussion: Increasing the degree of simplification of the CT-density ramp results in an increasing effect on the resulting radiological thickness calculated for the Supertech head phantom. For instance, defining the CT-density ramp using 8 points, instead of 12, results in a maximum radiological thickness change of 0.2 cm, whereas defining the CT-density ramp using only 2 points results in a maximum radiological thickness change of 11.2 cm. Changing the definition of the materials comprising the phantom between water and plastic and tissue results in millimetre-scale changes to the resulting radiological thickness. When the entire phantom is defined as lead, this alteration changes the calculated radiological thickness by a maximum of 9.7 cm. Evidently, the simplification of the CT-density ramp has a greater effect on the resulting radiological thickness map than does the alteration of the assignment of tissue types. Conclusions: It is possible to alter the definitions of the tissue types comprising the phantom (or patient) without substantially altering the results of simulated portal images. However, these images are very sensitive to the accurate identification of the HU-density relationship. When converting data from a patient’s CT into a MC simulation phantom, therefore, all possible care should be taken to accurately reproduce the conversion between HU and mass density, for the specific CT scanner used. Acknowledgements: This work is funded by the NHMRC, through a project grant, and supported by the Queensland University of Technology (QUT) and the Royal Brisbane and Women's Hospital (RBWH), Brisbane, Australia. The authors are grateful to the staff of the RBWH, especially Darren Cassidy, for assistance in obtaining the phantom CT data used in this study. The authors also wish to thank Cathy Hargrave, of QUT, for assistance in formatting the CT data, using the Pinnacle TPS. Computational resources and services used in this work were provided by the HPC and Research Support Group, QUT, Brisbane, Australia.

Pixel value to electron density conversion for a Megavoltage Cone Beam CT System

Introduction: The motivation for developing megavoltage (and kilovoltage) cone beam CT (MV CBCT) capabilities in the radiotherapy treatment room was primarily based on the need to improve patient set-up accuracy. There has recently been an interest in using the cone beam CT data for treatment planning. Accurate treatment planning, however, requires knowledge of the electron density of the tissues receiving radiation in order to calculate dose distributions. This is obtained from CT, utilising a conversion between CT number and electron density of various tissues. The use of MV CBCT has particular advantages compared to treatment planning with kilovoltage CT in the presence of high atomic number materials and requires the conversion of pixel values from the image sets to electron density. Therefore, a study was undertaken to characterise the pixel value to electron density relationship for the Siemens MV CBCT system, MVision, and determine the effect, if any, of differing the number of monitor units used for acquisition. If a significant difference with number of monitor units was seen then pixel value to ED conversions may be required for each of the clinical settings. The calibration of the MV CT images for electron density offers the possibility for a daily recalculation of the dose distribution and the introduction of new adaptive radiotherapy treatment strategies. Methods: A Gammex Electron Density CT Phantom was imaged with the MVCB CT system. The pixel value for each of the sixteen inserts, which ranged from 0.292 to 1.707 relative electron density to the background solid water, was determined by taking the mean value from within a region of interest centred on the insert, over 5 slices within the centre of the phantom. These results were averaged and plotted against the relative electron densities of each insert with a linear least squares fit was preformed. This procedure was performed for images acquired with 5, 8, 15 and 60 monitor units. Results: The linear relationship between MVCT pixel value and ED was demonstrated for all monitor unit settings and over a range of electron densities. The number of monitor units utilised was found to have no significant impact on this relationship. Discussion: It was found that the number of MU utilised does not significantly alter the pixel value obtained for different ED materials. However, to ensure the most accurate and reproducible MV to ED calibration, one MU setting should be chosen and used routinely. To ensure accuracy for the clinical situation this MU setting should correspond to that which is used clinically. If more than one MU setting is used clinically then an average of the CT values acquired with different numbers of MU could be utilized without loss in accuracy. Conclusions: No significant differences have been shown between the pixel value to ED conversion for the Siemens MV CT cone beam unit with change in monitor units. Thus as single conversion curve could be utilised for MV CT treatment planning. To fully utilise MV CT imaging for radiotherapy treatment planning further work will be undertaken to ensure all corrections have been made and dose calculations verified. These dose calculations may be either for treatment planning purposes or for reconstructing the delivered dose distribution from transit dosimetry measurements made using electronic portal imaging devices. This will potentially allow the cumulative dose distribution to be determined through the patient’s multi-fraction treatment and adaptive treatment strategies developed to optimize the tumour response.

Calibration of the electron volt spectrometer, a deep inelastic neutron scattering spectrometer at the ISIS pulsed neutron spallation source

The electron Volt Spectrometer (eVS) is an inverse geometry filter difference spectrometer that has been optimised to measure the single atom properties of condensed matter systems using a technique known as Neutron Compton Scattering (NCS) or Deep Inelastic Neutron Scattering (DINS). The spectrometer utilises the high flux of epithermal neutrons that are produced by the ISIS neutron spallation source enabling the direct measurement of atomic momentum distributions and ground state kinetic energies. In this paper the procedure that is used to calibrate the spectrometer is described. This includes details of the method used to determine detector positions and neutron flight path lengths as well as the determination of the instrument resolution. Examples of measurements on 3 different samples are shown, ZrH2, 4He and Sn which show the self-consistency of the calibration procedure.

The effect of hydroxyl groups and surface area of hematite derived from annealing goethite for phosphate removal

Synthetic goethite and thermally treated goethite at different temperatures were used to remove phosphate from sewage. The effect of annealing temperature on phosphate removal over time was investigated. X-ray diffraction(XRD), transmission electron microscopy (TEM), N2 adsorption and desorption (BET), and infrared emission spectrum (FT-IES) were utilized to characterize the phase, morphology, specific surface area, pore distribution, and the surface groups of samples. The results show that annealed products of goethite at temperatures over 250 °C are hematite with the similar morphology as the original goethite with different hydroxyl groups and surface area. Increasing temperature causes the decrease in hydroxyl groups, consequential increase in surface area at first and then experiences a decrease (14.8–110.4–12.6 m2/g) and the subsequent formation of nanoscale pores. The variation rate of hydroxyl groups and surface area based on FT-IES and BET, respectively, are used to evaluate the effect of annealing temperature on phosphate removal. By using all of the characterization techniques, it is concluded that the changes of phosphate removal basically result from the total variation rate between hydroxyl groups and surface area.

Thermal treatment of natural goethite : thermal transformation and physical properties

XRD (X-ray diffraction), XRF (X-ray fluorescence), TG (thermogravimetry), FT-IES (Fourier transform infrared emission spectroscopy), FESEM (field emission scanning electron microscope), TEM (transmission electron microscope) and nitrogen–adsorption–desorption analysis were used to characterize the composition and thermal evolution of the structure of natural goethite. The in situ FT-IES demonstrated the start temperature (250 °C) of the transformation of natural goethite to hematite and the thermodynamic stability of protohematite between 250 and 600 °C. The heated products showed a topotactic relationship to the original mineral based on SEM analysis. Finally, the nitrogen–adsorption–desorption isotherm provided the variation of surface area and pore size distribution as a function of temperature. The surface area displayed a remarkable increase up to 350 °C, and then decreased above this temperature. The significant increase in surface area was attributed to the formation of regularly arranged slit-shaped micropores running parallel to elongated direction of hematite microcrystal. The main pore size varied from 0.99 nm to 3.5 nm when heating temperature increases from 300 to 400 °C. The hematite derived from heating goethite possesses high surface area and favors the possible application of hematite as an adsorbent as well as catalyst carrier.

In situ observation of structural changes in polycrystalline silver catalysts by environmental scanning electron microscopy

Morphology changes induced in polycrystalline silver catalysts as a result of heating in either oxygen, water or oxygen-methanol atmospheres have been investigated by environmental scanning electron microscopy (ESEM), FT-Raman spectroscopy and temperature programmed desorption (TPD). The silver catalyst of interest consisted of two distinct particle types, one of which contained a significant concentration of sub-surface hydroxy species (in addition to surface adsorbed atomic oxygen). Heating the sample to 663 K resulted in the production of 'pin-holes' in the silver structure as a consequence of near-surface explosions caused by sub-surface hydroxy recombination. Furthermore, 'pin-holes' were predominantly found in the vicinity of surface defects, such as platelets and edge structures. Reaction between methanol and oxygen also resulted in the formation of 'pin-holes' in the silver surface, which were inherently associated with the catalytic process. A reaction mechanism is suggested that involves the interaction of methanol with sub-surface oxygen species to form sub-surface hydroxy groups. The sub-surface hydroxy species subsequently erupt through the silver surface to again produce 'pin-holes'.

A combined environmental scanning electron microscopy and Raman microscopy study of methanol oxidation on silver(I) oxide

The techniques of environmental scanning electron microscopy (ESEM) and Raman microscopy have been used to respectively elucidate the morphological changes and nature of the adsorbed species on silver(I) oxide powder, during methanol oxidation conditions. Heating Ag2O in either water vapour or oxygen resulted firstly in the decomposition of silver(I) oxide to polycrystalline silver at 578 K followed by sintering of the particles at higher temperature. Raman spectroscopy revealed the presence of subsurface oxygen and hydroxyl species in addition to surface hydroxyl groups after interaction with water vapour. Similar species were identified following exposure to oxygen in an ambient atmosphere. This behaviour indicated that the polycrystalline silver formed from Ag2O decomposition was substantially more reactive than silver produced by electrochemical methods. The interaction of water at elevated temperatures subsequent to heating silver(I) oxide in oxygen resulted in a significantly enhanced concentration of subsurface hydroxyl species. The reaction of methanol with Ag2O at high temperatures was interesting in that an inhibition in silver grain growth was noted. Substantial structural modification of the silver(I) oxide material was induced by catalytic etching in a methanol/air mixture. In particular, "pin-hole" formation was observed to occur at temperatures in excess of 773 K, and it was also recorded that these "pin- holes" coalesced to form large-scale defects under typical industrial reaction conditions. Raman spectroscopy revealed that the working surface consisted mainly of subsurface oxygen and surface Ag=O species. The relative lack of sub-surface hydroxyl species suggested that it was the desorption of such moieties which was the cause of the "pin-hole" formation.

In situ imaging of catalytic etching on silver during methanol oxidation conditions by environmental scanning electron microscopy

Polycrystalline silver is used to catalytically oxidise methanol to formaldehyde. This paper reports the results of extensive investigations involving the use of environmental scanning electron microscopy (ESEM) to monitor structural changes in silver during simulated industrial reaction conditions. The interaction of oxygen, nitrogen, and water, either singly or in combination, with a silver catalyst at temperatures up to 973 K resulted in the appearance of a reconstructed silver surface. More spectacular was the effect an oxygen/methanol mixture had on the silver morphology. At a temperature of ca. 713 K pinholes were created in the vicinity of defects as a consequence of subsurface explosions. These holes gradually increased in size and large platelet features were created. Elevation of the catalyst temperature to 843 K facilitated the wholescale oxygen induced restructuring of the entire silver surface. Methanol reacted with subsurface oxygen to produce subsurface hydroxyl species which ultimately formed water in the subsurface layers of silver. The resultant hydrostatic pressure forced the silver surface to adopt a "hill and valley" conformation in order to minimise the surface free energy. Upon approaching typical industrial operating conditions widespread explosions occurred on the catalyst and it was also apparent that the silver surface was extremely mobile under the applied conditions. The interaction of methanol alone with silver resulted in the initial formation of pinholes primarily in the vicinity of defects, due to reaction with oxygen species incorporated in the catalyst during electrochemical synthesis. However, dramatic reduction in the hole concentration with time occurred as all the available oxygen became consumed. A remarkable correlation between formaldehyde production and hole concentration was found.

Characterisation of combustion related acoustic emission sources for diesel engine condition monitoring

This thesis represents a major step forward in understanding the link between the development of combustion related faults in diesel engines and the generation of acoustic emissions. The findings presented throughout the thesis provide a foundation so that future diesel engine monitoring systems are able to more effectively detect and monitor developing faults. In undertaking this research knowledge concerning engine function and relevant failure mechanisms was combined with different modelling methods to generate a framework that was used to effectively identify fault related activity within acoustic emissions recorded from different engines.

Confirmation of the assignment of vibrations of goethite: An ATR and IES study of goethite structure

Transmission electron microscopy (TEM), field emission scanning electron microscopy (FE-SEM/EDS) and X-ray diffraction (XRD) were used to characterize the morphology of synthetic goethite. The behavior of the hydroxyl/water molecular units of goethite and its thermally treated products were characterized using Fourier transform-infrared emission spectroscopy (FT-IES) and attenuated total reflectance–Fourier transform infrared (ATR–FTIR) spectroscopy. The results showed that all the expected vibrational bands between 4000 and 650 cm−1 including the resolved bands (3800–2200 cm−1) were confirmed. A band attributed to a new type of hydroxyl unit was found at 3708 cm−1 and assigned to the FeO–H stretching vibration without hydrogen bonding. This hydroxyl unit was retained up to the thermal treatment temperature of 500 °C. On the whole, seven kinds of hydroxyl units, involving three surface hydroxyls, a bulk hydroxyl, a FeO–H without hydrogen bonding, a nonstoichiometric hydroxyl and a reversed hydroxyl were observed, and three kinds of adsorbed water were found in/on goethite.

Removal of phosphorus using NZVI derived from reducing natural goethite

Nano zero valent iron (NZVI) prepared by reducing natural goethite in hydrogen at 550 °C was used to remove phosphate. The effect of particle size, reaction time, NZVI dose, pH, initial phosphorus concentration, and oxygen amount in reaction system on phosphorus removal was investigated. The characterization of X-ray fluorescence (XRF), X-ray diffraction (XRD), N2 adsorption and desorption (BET analysis), transmission electron microscope (TEM), field emission scanning electron microscope with a energy dispersive X-ray detector (FESEM/EDS) and X-ray photoelectron spectroscopy (XPS) indicated that nanoscale of iron (around 80–150 nm length and 5–30 nm width) was prepared successfully with high dispersion and relative large surface area around 22 m2/g. The results of batch experiments and XPS analysis suggested that this kind of NZVI had a good performance on removal of phosphate (over 99%) despite in slightly acidic media as the initial concentration of P was 5 mg/L. The reason was ascribed to the effective corrosion of this NZVI under the function of proton and dissolved oxygen in spite of the existence of thin passive films.

Excursion guidance criteria to guide control of peak emission and exposure to airborne engineered particles

The overall aim of our research was to characterize airborne particles from selected nanotechnology processes and to utilize the data to develop and test quantitative particle concentration-based criteria that can be used to trigger an assessment of particle emission controls. We investigated particle number concentration (PNC), particle mass (PM) concentration, count median diameter (CMD), alveolar deposited surface area, elemental composition, and morphology from sampling of aerosols arising from six nanotechnology processes. These included fibrous and non-fibrous particles, including carbon nanotubes (CNTs). We adopted standard occupational hygiene principles in relation to controlling peak emission and exposures, as outlined by both Safe Work Australia, (1) and the American Conference of Governmental Industrial Hygienists (ACGIH®). (2) The results from the study were used to analyses peak and 30-minute averaged particle number and mass concentration values measured during the operation of the nanotechnology processes. Analysis of peak (highest value recorded) and 30-minute averaged particle number and mass concentration values revealed: Peak PNC20–1000 nm emitted from the nanotechnology processes were up to three orders of magnitude greater than the local background particle concentration (LBPC). Peak PNC300–3000 nm was up to an order of magnitude greater, and PM2.5 concentrations up to four orders of magnitude greater. For three of these nanotechnology processes, the 30-minute average particle number and mass concentrations were also significantly different from the LBPC (p-value < 0.001). We propose emission or exposure controls may need to be implemented or modified, or further assessment of the controls be undertaken, if concentrations exceed three times the LBPC, which is also used as the local particle reference value, for more than a total of 30 minutes during a workday, and/or if a single short-term measurement exceeds five times the local particle reference value. The use of these quantitative criteria, which we are terming the universal excursion guidance criteria, will account for the typical variation in LBPC and inaccuracy of instruments, while precautionary enough to highlight peaks in particle concentration likely to be associated with particle emission from the nanotechnology process. Recommendations on when to utilize local excursion guidance criteria are also provided.

Validating SMPS-measured size distribution of double-mode spherical Silica nanoparticles by Transmission Electron Microscope (TEM)

A nanoparticles size is one of their key physical characteristics that can affect their fate in a human’s respiratory tract (in case of inhalation) and also in the environment. Hence, measuring the size distribution of nanoparticles is absolutely essential and contributes greatly to their characterization. For years, Scanning Mobility Particle Sizers (SMPS), which rely on measuring the electrical mobility diameter of particles, have been used as one of the most reliable real-time instruments for the size distribution measurement of nanoparticles. Despite its benefits, this instrument has some drawbacks, including equivalency problems for non-spherical particles (i.e. assuming a non-spherical particle is equal to a spherical particle of diameter d due to the same electrical mobility), as well as limitations in terms of its use in workplaces, because of its large size and the complexity of its operation...

Remote measurement of diesel locomotive emission factors and particle size distributions

A technique for analysing exhaust emission plumes from unmodified locomotives under real world conditions is described and applied to the task of characterizing plumes from railway trains servicing an Australian shipping port. The method utilizes the simultaneous measurement, downwind of the railway line, of the following pollutants; particle number, PM2.5 mass fraction, SO2, NOx and CO2, with the last of these being used as an indicator of fuel combustion. Emission factors are then derived, in terms of number of particles and mass of pollutant emitted per unit mass of fuel consumed. Particle number size distributions are also presented. The practical advantages of the method are discussed including the capacity to routinely collect emission factor data for passing trains and to thereby build up a comprehensive real world database for a wide range of pollutants. Samples from 56 train movements were collected, analyzed and presented. The quantitative results for emission factors are: EF(N)=(1.7±1)×1016 kg-1, EF(PM2.5)= (1.1±0.5) g·kg-1, EF(NOx)= (28±14) g·kg-1, and EF(SO2 )= (1.4±0.4) g·kg-1. The findings are compared with comparable previously published work. Statistically significant (p<α, α=0.05) correlations within the group of locomotives sampled were found between the emission factors for particle number and both SO2 and NOx.