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.

Electrical property and characterization of nano-SnO2/wollastonite composite materials

Nano-tin oxide was deposited on the surface of wollastonite using the mixed solution including stannic chloride pentahydrate precursor and wollastonite by a hydrolysis precipitation process. The antistatic properties of the wollastonite materials under different calcined conditions and composite materials (nano-SnO2/wollastonite, SW) were measured by rubber sheeter and four-point probe (FPP) sheet resistance measurement. Effects of hydrolysis temperature and time, calcination temperature and time, pH value and nano-SnO2 coating amount on the resistivity of SW powders were studied, and the optimum experimental conditions were obtained. The microstructure and surface properties of wollastonite, precipitate and SW were characterized by transmission electron microscopy (TEM), scanning electron microscopy (SEM), energy-dispersive X-ray spectrometry (EDS), specific surface area analyzer (BET), thermogravimetry (TG), X-ray photoelectron spectroscopy (XPS), X-ray diffraction (XRD), and Fourier translation infrared spectroscopy (FTIR) respectively. The results showed that the nano-SnO2/wollastonite composite materials under optimum preparation conditions showed better antistatic properties, the resistivity of which was reduced from 1.068 × 104 Ω cm to 2.533 × 103 Ω cm. From TG and XRD analysis, the possible mechanism for coating of SnO2 nanoparticles on the surface of wollastonite was proposed. The infrared spectrum indicated that there were a large number of the hydroxyl groups on the surface of wollastonite. This is beneficial to the heterogeneous nucleation reaction. Through morphology, EDS and XPS analysis, the surface of wollastonite fiber was coated with a layer of 10–15 nm thickness of tin oxide grains the distribution of which was uniform.

Fabrication and characterisation of graphene oxide-epoxy nanocomposite

Adequate amount of graphene oxide (GO) was firstly prepared by oxidation of graphite and GO/epoxy nanocomposites were subsequently prepared by typical solution mixing technique. X-ray diffraction (XRD) pattern, X-ray photoelectron (XPS), Raman and Fourier transform infrared (FTIR) spectroscopy indicated the successful preparation of GO. Scanning electron microscopy (SEM) and Transmission electron microscopy (TEM) images of the graphite oxide showed that they consist of a large amount of graphene oxide platelets with a curled morphology containing of a thin wrinkled sheet like structure. AFM image of the exfoliated GO signified that the average thickness of GO sheets is ~1.0 nm which is very similar to GO monolayer. Mechanical properties of as prepared GO/epoxy nanocomposites were investigated. Significant improvements in both Young’s modulus and tensile strength were observed for the nanocomposites at very low level of GO loading. The Young’s modulus of the nanocomposites containing 0.5 wt% GO was 1.72 GPa, which was 35 % higher than that of the pure epoxy resin (1.28 GPa). The effective reinforcement of the GO based epoxy nanocomposites can be attributed to the good dispersion and the strong interfacial interactions between the GO sheets and the epoxy resin matrices.

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'.

Characterisation of SiO2-supported nickel catalysts for carbon dioxide reforming of methane

The CO2-methane reformation reaction over Ni/SiO2 catalysts has been extensively studied using a range of temperature-programmed techniques and characterisation of the catalysts by thermogravimetry (TG), X-ray diffraction (XRD) and electron microscopy (TEM). The results indicate a strong correlation between the microstructure of the catalyst and its performance. The role of both CO2 and CH4 in the reaction has been investigated and the role of methyl radicals in the reaction mechanism highlighted. A reaction mechanism involving dissociatively adsorbed CO2 and methyl radicals has been proposed.

Characterization of precursors to methanol synthesis catalysts Cu/ZnO system

The composition of a series of hydroxycarbonate precursors to copper/zinc oxide methanol synthesis catalysts prepared under conditions reported as optimum for catalytic activity has been studied. Techniques employed included thermogravimetry (TG), temperature-programmed decomposition (TPD), X-ray diffraction (XRD), high-resolution transmission electron microscopy (HRTEM), and Raman and FTIR spectroscopies. Evidence was obtained for various structural phases including hydrozincite, copper hydrozincite, aurichalcite, zincian malachite and malachite (the concentrations of which depended upon the exact Cu/Zn ratio used). Significantly, previously reported phases such as gerhardite and rosasite were not identified when catalysts were synthesized at optimum solution pH and temperature values, and after appropriate aging periods. Calcination of the hydroxycarbonate precursors resulted in the formation of catalysts containing an intimate mixture of copper and zinc oxides. Temperature-programmed reduction (TPR) revealed that a number of discrete copper oxide species were present in the catalyst, the precise concentrations of which were determined to be related to the structure of the catalyst precursor. Copper hydrozincite decomposed to give zinc oxide particles decorated by highly dispersed, small copper oxide species. Aurichalcite appeared to result ultimately in the most intimately mixed catalyst structure whereas zincian malachite decomposed to produce larger copper oxide and zinc oxide grains. The reason for the stabilization of small copper oxide and zinc oxide clusters by aurichalcite was investigated by using carefully selected calcination temperatures. It was concluded that the unique formation of an 'anion-modified' oxide resulting from the initial decomposition stage of aurichalcite was responsible for the 'binding' of copper species to zinc moieties.

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.

Fabrication and mechanical and thermal behaviour of graphene oxide/epoxy nanocomposites

Bulk amount of graphite oxide was prepared by oxidation of graphite using the modified Hummers method and its ultrasonication in organic solvents yielded graphene oxide (GO). X-ray diffraction (XRD) pattern, X-ray photoelectron (XPS), Raman and Fourier transform infrared (FTIR) spectroscopy indicated the successful preparation of GO. XPS survey spectrum of GO revealed the presence of 66.6 at% C and 30.4 at% O. Scanning electron microscopy (SEM) and Transmission electron microscopy (TEM) images of the graphene oxide showed that they consist of a large amount of graphene oxide platelets with a curled morphology containing of a thin wrinkled sheet like structure. AFM image of the exfoliated GO signified that the average thickness of GO sheets is ~1.0 nm which is very similar to GO monolayer. GO/epoxy nanocomposites were prepared by typical solution mixing technique and influence of GO on mechanical and thermal properties of nanocomposites were investigated. As for the mechanical behaviour of GO/epoxy nanocomposites, 0.5 wt% GO in the nanocomposite achieved the maximum increase in the elastic modulus (~35%) and tensile strength (~7%). The TEM analysis provided clear image of microstructure with homogeneous dispersion of GO in the polymer matrix. The improved strength properties of GO/epoxy nanocomposites can be attributed to inherent strength of GO, the good dispersion and the strong interfacial interactions between the GO sheets and the polymer matrix. However, incorporation of GO showed significant negative effect on composite glass transition temperature (Tg). This may arise due to the interference of GO on curing reaction of epoxy.

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.