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.

An investigation into factors affecting electron density calibration for a megavoltage cone-beam CT system

There is a growing interest in the use of megavoltage cone-beam computed tomography (MV CBCT) data for radiotherapy treatment planning. To calculate accurate dose distributions, knowledge of the electron density (ED) of the tissues being irradiated is required. In the case of MV CBCT, it is necessary to determine a calibration-relating CT number to ED, utilizing the photon beam produced for MV CBCT. A number of different parameters can affect this calibration. This study was undertaken on the Siemens MV CBCT system, MVision, to evaluate the effect of the following parameters on the reconstructed CT pixel value to ED calibration: the number of monitor units (MUs) used (5, 8, 15 and 60 MUs), the image reconstruction filter (head and neck, and pelvis), reconstruction matrix size (256 by 256 and 512 by 512), and the addition of extra solid water surrounding the ED phantom. A Gammex electron density CT phantom containing EDs from 0.292 to 1.707 was imaged under each of these conditions. The linear relationship between MV CBCT pixel value and ED was demonstrated for all MU settings and over the range of EDs. Changes in MU number did not dramatically alter the MV CBCT ED calibration. The use of different reconstruction filters was found to affect the MV CBCT ED calibration, as was the addition of solid water surrounding the phantom. Dose distributions from treatment plans calculated with simulated image data from a 15 MU head and neck reconstruction filter MV CBCT image and a MV CBCT ED calibration curve from the image data parameters and a 15 MU pelvis reconstruction filter showed small and clinically insignificant differences. Thus, the use of a single MV CBCT ED calibration curve is unlikely to result in any clinical differences. However, to ensure minimal uncertainties in dose reporting, MV CBCT ED calibration measurements could be carried out using parameter-specific calibration measurements.

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.

Distribution of ion-current density on substrate between carbon nanotubes grown from low-temperature plasma

Using Monte Carlo simulation technique, we have calculated the distribution of ion current extracted from low-temperature plasmas and deposited onto the substrate covered with a nanotube array. We have shown that a free-standing carbon nanotube is enclosed in a circular bead of the ion current, whereas in square and hexagonal nanotube patterns, the ion current is mainly concentrated along the lines connecting the nearest nanotubes. In a very dense array (with the distance between nanotubes/nanotube-height ratio less than 0.05), the ions do not penetrate to the substrate surface and deposit on side surfaces of the nanotubes.

Nanopowder management and control of plasma parameters in electronegative SiH4 plasmas

The nanopowder management and control of plasma parameters in electronegative SiH4 plasmas were discussed. The spatial profiles of electron and positive/negative ion number densities, electron temperature and charge of the fine particles were obtained. It was found that management of powder charge distribution is also possible through control of the external parameters.

Behavior of the electron temperature in nonuniform complex plasmas

The response of complex ionized gas systems to the presence of nonuniform distribution of charged grains is investigated using a kinetic model. Contrary to an existing view that the electron temperature inevitably increases in the grain-occupied region because of enhanced ionization to compensate for the electrons lost to the grains, it is shown that this happens only when the ionizing electric field increases in the electron depleted region. The results for two typical plasma systems suggest that when the ionizing electric field depends on the spatially averaged electron density, the electron temperature in the grain containing region can actually decrease.

Dispersion properties of surface waves in a two-layer plasma structure bounded by a metal

The effect of the nonuniformity of the electron density on the dispersion properties of surface waves propagating in a direction transverse to an external magnetic field is studied for the model of a two-layer plasma structure bounded by a metal. It is shown that the spectra of the waves can be effectively controlled by varying the degree of nonuniformity of the density and the dimensions of the layers.

The dependence of computed tomography number to relative electron density conversion on phantom geometry and its impact on planned dose

A computed tomography number to relative electron density (CT-RED) calibration is performed when commissioning a radiotherapy CT scanner by imaging a calibration phantom with inserts of specified RED and recording the CT number displayed. In this work, CT-RED calibrations were generated using several commercially available phantoms to observe the effect of phantom geometry on conversion to electron density and, ultimately, the dose calculation in a treatment planning system. Using an anthropomorphic phantom as a gold standard, the CT number of a material was found to depend strongly on the amount and type of scattering material surrounding the volume of interest, with the largest variation observed for the highest density material tested, cortical bone. Cortical bone gave a maximum CT number difference of 1,110 when a cylindrical insert of diameter 28 mm scanned free in air was compared to that in the form of a 30 × 30 cm2 slab. The effect of using each CT-RED calibration on planned dose to a patient was quantified using a commercially available treatment planning system. When all calibrations were compared to the anthropomorphic calibration, the largest percentage dose difference was 4.2 % which occurred when the CT-RED calibration curve was acquired with heterogeneity inserts removed from the phantom and scanned free in air. The maximum dose difference observed between two dedicated CT-RED phantoms was ±2.1 %. A phantom that is to be used for CT-RED calibrations must have sufficient water equivalent scattering material surrounding the heterogeneous objects that are to be used for calibration.

Fabrication and characterization of plasma polymer thin films from monoterpene alcohols for applications in organic electronics and biotechnology

After more than twenty years of basic and applied research, the use of nanotechnology in the design and manufacture of nanoscale materials is rapidly increasing, particularly in commercial applications that span from electronics across renewable energy areas, and biomedical devices. Novel polymers are attracting significant attention for they promise to provide a low−cost high−performance alternative to existing materials. Furthermore, these polymers have the potential to overcome limitations imposed by currently available materials thus enabling the development of new technologies and applications that are currently beyond our reach. This work focuses on the development of a range of new low−cost environmentally−friendly polymer materials for applications in areas of organic (flexible) electronics, optics, and biomaterials. The choice of the monomer reflects the environmentally−conscious focus of this project. Terpinen−4−ol is a major constituent of Australian grown Melaleuca alternifolia (tea tree) oil, attributed with the oil's antimicrobial and anti−inflammatory properties. Plasma polymerisation was chosen as a deposition technique for it requires minimal use of harmful chemicals and produces no hazardous by−products. Polymer thin films were fabricated under varied process conditions to attain materials with distinct physico−chemical, optoelectrical, biological and degradation characteristics. The resultant materials, named polyterpenol, were extensively characterised using a number of well−accepted and novel techniques, and their fundamental properties were defined. Polyterpenol films were demonstrated to be hydrocarbon rich, with variable content of oxygen moieties, primarily in the form of hydroxyl and carboxyl functionalities. The level of preservation of original monomer functionality was shown to be strongly dependent on the deposition energy, with higher applied power increasing the molecular fragmentation and substrate temperature. Polyterpenol water contact angle contact angle increased from 62.7° for the 10 W samples to 76.3° for the films deposited at 100 W. Polymers were determined to resist solubilisation by water, due to the extensive intermolecular and intramolecular hydrogen bonds present, and other solvents commonly employed in electronics and biomedical processing. Independent of deposition power, the surface topography of the polymers was shown to be smooth (Rq <0.5 nm), uniform and defect free. Hardness of polyterpenol coatings increased from 0.33 GPa for 10 W to 0.51 GPa for 100 W (at 500 μN load). Coatings deposited at higher input RF powers showed less mechanical deformation during nanoscratch testing, with no considerable damage, cracking or delamination observed. Independent of the substrate, the quality of film adhesion improved with RF power, suggesting these coatings are likely to be more stable and less susceptible to wear. Independent of fabrication conditions, polyterpenol thin films were optically transparent, with refractive index approximating that of glass. Refractive index increased slightly with deposition power, from 1.54 (10 W) to 1.56 (100 W) at 500 nm. The optical band gap values declined with increasing power, from 2.95 eV to 2.64 eV, placing the material within the range for semiconductors. Introduction of iodine impurity reduced the band gap of polyterpenol, from 2.8 eV to 1.64 eV, by extending the density of states more into the visible region of the electromagnetic spectrum. Doping decreased the transparency and increased the refractive index from 1.54 to 1.70 (at 500 nm). At optical frequencies, the real part of permittivity (k) was determined to be between 2.34 and 2.65, indicating a potential low-k material. These permittivity values were confirmed at microwave frequencies, where permittivity increased with input RF energy – from 2.32 to 2.53 (at 10 GHz ) and from 2.65 to 2.83 (at 20 GHz). At low frequencies, the dielectric constant was determined from current−voltage characteristics of Al−polyterpenol−Al devices. At frequencies below 100 kHz, the dielectric constant varied with RF power, from 3.86 to 4.42 at 1 kHz. For all samples, the resistivity was in order of 10⁸−10⁹ _m (at 6 V), confirming the insulating nature of polyterpenol material. In situ iodine doping was demonstrated to increase the conductivity of polyterpenol, from 5.05 × 10⁻⁸ S/cm to 1.20 × 10⁻⁶ S/cm (at 20 V). Exposed to ambient conditions over extended period of time, polyterpenol thin films were demonstrated to be optically, physically and chemically stable. The bulk of ageing occurred within first 150 h after deposition and was attributed to oxidation and volumetric relaxation. Thermal ageing studies indicated thermal stability increased for the films manufactured at higher RF powers, with degradation onset temperature associated with weight loss shifting from 150 ºC to 205 ºC for 10 W and 100 W polyterpenol, respectively. Annealing the films to 405 °C resulted in full dissociation of the polymer, with minimal residue. Given the outcomes of the fundamental characterisation, a number of potential applications for polyterpenol have been identified. Flexibility, tunable permittivity and loss tangent properties of polyterpenol suggest the material can be used as an insulating layer in plastic electronics. Implementation of polyterpenol as a surface modification of the gate insulator in pentacene-based Field Effect Transistor resulted in significant improvements, shifting the threshold voltage from + 20 V to –3 V, enhancing the effective mobility from 0.012 to 0.021 cm²/Vs, and improving the switching property of the device from 10⁷ to 10⁴. Polyterpenol was demonstrated to have a hole transport electron blocking property, with potential applications in many organic devices, such as organic light emitting diodes. Encapsulation of biomedical devices is also proposed, given that under favourable conditions, the original chemical and biological functionality of terpinen−4−ol molecule can be preserved. Films deposited at low RF power were shown to successfully prevent adhesion and retention of several important human pathogens, including P. aeruginosa, S. aureus, and S. epidermidis, whereas films deposited at higher RF power promoted bacterial cell adhesion and biofilm formation. Preliminary investigations into in vitro biocompatibility of polyterpenol demonstrated the coating to be non−toxic for several types of eukaryotic cells, including Balb/c mice macrophage and human monocyte type (HTP−1 non-adherent) cells. Applied to magnesium substrates, polyterpenol encapsulating layer significantly slowed down in vitro biodegradation of the metal, thus increasing the viability and growth of HTP−1 cells. Recently, applied to varied nanostructured titanium surfaces, polyterpenol thin films successfully reduced attachment, growth, and viability of P. aeruginosa and S. aureus.

Characteristics of epoxy resin/SiO2 nanocomposite insulation : effects of plasma surface treatment on the nanoparticles

The present study compares the effects of two different material processing techniques on modifying hydrophilic SiO2 nanoparticles. In one method, the nanoparticles undergo plasma treatment by using a custom-developed atmospheric-pressure non-equilibrium plasma reactor. With the other method, they undergo chemical treatment which grafts silane groups onto their surface and turns them into hydrophobic. The treated nanoparticles are then used to synthesize epoxy resin-based nanocomposites for electrical insulation applications. Their characteristics are investigated and compared with the pure epoxy resin and nanocomposite fabricated with unmodified nanofillers counterparts. The dispersion features of the nanoparticles in the epoxy resin matrix are examined through scanning electron microscopy (SEM) images. All samples show evidence that the agglomerations are smaller than 30 nm in their diameters. This indicates good dispersion uniformity. The Weibull plot of breakdown strength and the recorded partial discharge (PD) events of the epoxy resin/plasma-treated hydrophilic SiO2 nanocomposite (ER/PTI) suggest that the plasma-treated specimen yields higher breakdown strength and lower PD magnitude as compared to the untreated ones. In contrast, surprisingly, lower breakdown strength is found for the nanocomposite made by the chemically treated hydrophobic particles, whereas the PD magnitude and PD numbers remain at a similar level as the plasma-treated ones.

Deterministic control of structural and optical properties of plasma-grown vertical graphene nanosheet networks via nitrogen gas variation

The effect of nitrogen on the growth of vertically oriented graphene nanosheets on catalyst-free silicon and glass substrates in a plasma-assisted process is studied. Different concentrations of nitrogen were found to act as versatile control knobs that could be used to tailor the length, number density and structural properties of the nanosheets. Nanosheets with different structural characteristics exhibit markedly different optical properties. The nanosheet samples were treated with a bovine serum albumin protein solution to investigate the effects of this variation on the optical properties for biosensing through confocal micro-Raman spectroscopy and UV-Vis spectrophotometry. © 2012 Optical Society of America.

Generation of uniform plasmas by crossed internal oscillating current sheets : key concepts and experimental verification

The results of comprehensive experimental studies of the operation, stability, and plasma parameters of the low-frequency (0.46 MHz) inductively coupled plasmas sustained by the internal oscillating rf current are reported. The rf plasma is generated by using a custom-designed configuration of the internal rf coil that comprises two perpendicular sets of eight currents in each direction. Various diagnostic tools, such as magnetic probes, optical emission spectroscopy, and an rf-compensated Langmuir probe were used to investigate the electromagnetic, optical, and global properties of the argon plasma in wide ranges of the applied rf power and gas feedstock pressure. It is found that the uniformity of the electromagnetic field inside the plasma reactor is improved as compared to the conventional sources of inductively coupled plasmas with the external flat coil configuration. A reasonable agreement between the experimental data and computed electromagnetic field topography inside the chamber is reported. The Langmuir probe measurements reveal that the spatial profiles of the electron density, the effective electron temperature, plasma potential, and electron energy distribution/probability functions feature a high degree of the radial and axial uniformity and a weak azimuthal dependence, which is consistent with the earlier theoretical predictions. As the input rf power increases, the azimuthal dependence of the global plasma parameters vanishes. The obtained results demonstrate that by introducing the internal oscillated rf currents one can noticeably improve the uniformity of electromagnetic field topography, rf power deposition, and the plasma density in the reactor.

E and H regimes of plasma enhanced chemical vapor deposition of diamond-like carbon film in low frequency inductively coupled plasma reactor

This paper reports on the efficient deposition of hydrogenated diamond-like carbon (DLC) film in a plasma reactor that features both the capacitively and inductively coupled operation regimes. The hydrogenated DLC films have been prepared on silicon wafers using a low-frequency (500 kHz) inductively coupled plasma (ICP) chemical vapor deposition (CVD) system. At low RF powers, the system operates as an asymmetric capacitively coupled plasma source, and the film deposition process is undertaken in the electrostatic (E) discharge regime. Above the mode transition threshold, the high-density inductively coupled plasma is produced in the electromagnetic (H) discharge regime. It has been shown that the deposition rate and hardness of the DLC film are much higher in the H-mode deposition regime. For a 2.66-Pa H-mode CH4 + Ar gas mixture discharge, the deposited DLC film exhibits a mechanical hardness of 18 GPa, Young's modulus of 170 GPa, and compressive stress of 1.3 GPa.

Effect of ion current density on the properties of vacuum arc-deposited TiN coatings

The influence of ion current density on the thickness of coatings deposited in a vacuum arc setup has been investigated to optimize the coating porosity. A planar probe was used to measure the ion current density distribution across plasma flux. A current density from 20 to 50 A/m2 was obtained, depending on the probe position relative to the substrate center. TiN coatings were deposited onto the cutting inserts placed at different locations on the substrate, and SEM was used to characterize the surfaces of the coatings. It was found that lowdensity coatings were formed at the decreased ion current density. A quantitative dependence of the coating thickness on the ion current density in the range of 20-50 A/m2 were obtained for the films deposited at substrate bias of 200 V and nitrogen pressure 0.1 Pa, and the coating porosity was calculated. The coated cutting inserts were tested by lathe machining of the martensitic stainless steel AISI 431. The results may be useful for controlling ion flux distribution over large industrial-scale substrates.

Infrared spectroscopic investigation of the effects of titania photocatalyst on the degradation of linear low density polyethylene film for commercial applications

There is a need in industry for a commodity polyethylene film with controllable degradation properties that will degrade in an environmentally neutral way, for applications such as shopping bags and packaging film. Additives such as starch have been shown to accelerate the degradation of plastic films, however control of degradation is required so that the film will retain its mechanical properties during storage and use, and then degrade when no longer required. By the addition of a photocatalyst it is hoped that polymer film will breakdown with exposure to sunlight. Furthermore, it is desired that the polymer film will degrade in the dark, after a short initial exposure to sunlight. Research has been undertaken into the photo- and thermo-oxidative degradation processes of 25 ìm thick LLDPE (linear low density polyethylene) film containing titania from different manufacturers. Films were aged in a suntest or in an oven at 50 °C, and the oxidation product formation was followed using IR spectroscopy. Degussa P25, Kronos 1002, and various organic-modified and doped titanias of the types Satchleben Hombitan and Hunstsman Tioxide incorporated into LLDPE films were assessed for photoactivity. Degussa P25 was found to be the most photoactive with UVA and UVC exposure. Surface modification of titania was found to reduce photoactivity. Crystal phase is thought to be among the most important factors when assessing the photoactivity of titania as a photocatalyst for degradation. Pre-irradiation with UVA or UVC for 24 hours of the film containing 3% Degussa P25 titania prior to aging in an oven resulted in embrittlement in ca. 200 days. The multivariate data analysis technique PCA (principal component analysis) was used as an exploratory tool to investigate the IR spectral data. Oxidation products formed in similar relative concentrations across all samples, confirming that titania was catalysing the oxidation of the LLDPE film without changing the oxidation pathway. PCA was also employed to compare rates of degradation in different films. PCA enabled the discovery of water vapour trapped inside cavities formed by oxidation by titania particles. Imaging ATR/FTIR spectroscopy with high lateral resolution was used in a novel experiment to examine the heterogeneous nature of oxidation of a model polymer compound caused by the presence of titania particles. A model polymer containing Degussa P25 titania was solvent cast onto the internal reflection element of the imaging ATR/FTIR and the oxidation under UVC was examined over time. Sensitisation of 5 ìm domains by titania resulted in areas of relatively high oxidation product concentration. The suitability of transmission IR with a synchrotron light source to the study of polymer film oxidation was assessed as the Australian Synchrotron in Melbourne, Australia. Challenges such as interference fringes and poor signal-to-noise ratio need to be addressed before this can become a routine technique.