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.

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.

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.

Maximum-entropy calculation of the electron density at 4 A resolution of Pf1 filamentous bacteriophage

A 4 A electron-density map of Pf1 filamentous bacterial virus has been calculated from x-ray fiber diffraction data by using the maximum-entropy method. This method produces a map that is free of features due to noise in the data and enables incomplete isomorphous-derivative phase information to be supplemented by information about the nature of the solution. The map shows gently curved (banana-shaped) rods of density about 70 A long, oriented roughly parallel to the virion axis but slewing by about 1/6th turn while running from a radius of 28 A to one of 13 A. Within these rods, there is a helical periodicity with a pitch of 5 to 6 A. We interpret these rods to be the helical subunits of the virion. The position of strongly diffracted intensity on the x-ray fiber pattern shows that the basic helix of the virion is right handed and that neighboring nearly parallel protein helices cross one another in an unusual negative sense.

Low-frequency variability and the galactic distribution of electron-density fluctuation

A sample of 96 compact flat-spectrum extragalactic sources, spread evenly over all galactic latitudes, has been studied at 327 MHz for variability over a time interval of about 15 yr. The variability shows a dependence on galactic latitude being less both at low and high latitudes and peaking around absolute value of b approximately 15-degrees. The latitude dependence is surprisingly similar in both the galactic centre and anticentre directions. Assuming various single and multi-component distributions for the ionized, irregular interstellar plasma, we have tried to generate the observed dependence using a semi-qualitative treatment of refractive interstellar scintillations. We find that it is difficult to fit our data with any single or double component cylindrical distribution. Our data suggests that the observed variability could be influenced by the spiral structure of our Galaxy.

Topological electron density analysis and electrostatic properties of aspirin: an experimental and theoretical study

The topological and the electrostatic properties of the aspirin drug molecule were determined from high-resolution X-ray diffraction data at 90 K, and the corresponding results are compared with the theoretical calculations. The electron density at the bond critical point of all chemical bonds induding the intermolecular interactions of aspirin has been quantitatively described using Bader's quantum theory of ``Atoms in Molecules''. The electrostatic potential of the molecule emphasizes the preferable binding sites of the drug and the interaction features of the molecule, which are crucial for drug-receptor recognition. The topological analysis of hydrogen bonds reveals the strength of intermolecular interactions.

Effect of analyte concentration on the laser-induced plasma temperature and electron density in liquid matrix

In the present work, we report spectroscopic studies of laser-induced plasmas produced by focusing the second harmonic (532nm) of a Nd:YAG laser onto the laminar flow of a liquid containing chromium. The plasma temperature is determined from the coupled Saha-Boltzmann plot and the electron density is evaluated from the Stark broadening of an ionic line of chromium Cr(II)] at 267.7nm. Our results reveal a decrease in plasma temperature with an increase in Cr concentration up to a certain concentration level; after that, it becomes approximately constant, while the electron density increases with an increase in analyte (Cr) concentration in liquid matrix.

Topological analysis of electron density and the electrostatic properties of isoniazid: an experimental and theoretical study

Isoniazid (isonicotinohydrazide) is an important first-line antitubercular drug that targets the InhA enzyme which synthesizes the critical component of the mycobacterial cell wall. An experimental charge-density analysis of isoniazid has been performed to understand its structural and electronic properties in the solid state. A high-resolution single-crystal X-ray intensity data has been collected at 90 K. An aspherical multipole refinement was carried out to explore the topological and electrostatic properties of the isoniazid molecule. The experimental results were compared with the theoretical charge-density calculations performed using CRYSTAL09 with the B3LYP/6-31G** method. A topological analysis of the electron density reveals that the Laplacian of electron density of the N-N bond is significantly less negative, which indicates that the charges at the b.c.p. (bond-critical point) of the bond are least accumulated, and so the bond is considered to be weak. As expected, a strong negative electrostatic potential region is present in the vicinity of the O1, N1 and N3 atoms, which are the reactive locations of the molecule. The C-H center dot center dot center dot N, C-H center dot center dot center dot O and N-H center dot center dot center dot N types of intermolecular hydrogen-bonding interactions stabilize the crystal structure. The topological analysis of the electron density on hydrogen bonding shows the strength of intermolecular interactions.

Measurement of Electron Density Profile Behind Strong Shock Waves with a Langmuir Probe

At the shock velocity range of 7~9km/s, the variations of electron density behind strong normal shock waves are measured in a low-density shock tube by using the Langmuir electrostatic probe technique. The electron temperature, calculated based on Park’s three-temperature model, is used in interpreting the probe current data. The peak electron densities determined in the present experiment are shown to be in a good agreement with those predicted by Lin’s calculation. The experimentally obtained ratios of the characteristic ionization distance to the mean free path of freestream ahead of the shock wave are found to be in a good agreement with the existing experiments and Park’s calculation.

Comparison of the electron density measurements using Thomson scattering and emission spectroscopy for laser induced breakdown in one atmosphere of helium

Thomson scattering from laser-induced plasma in atmospheric helium was used to obtain temporally and spatially resolved electron temperature and density profiles. Electron density measurements at 5 s after breakdown are compared with those derived from the separation of the allowed and forbidden components of the 447.1 nm He I line. Plasma is created using 9 ns, 140 mJ pulses from Nd:YAG laser at 1064 nm. Electron densities of ~5 × 10 cm are in good agreement with Thomson scattering measurements, benchmarking this emission line as a useful diagnostic for high density plasmas. © 2011 American Institute of Physics.

Temporal and Spatial Behavior of Electron Density and Temperature in a Laser-Produced Plasma from YBa2Cu3O7

Spectroscopic studies of laser -induced plasma from a high-temperature superconducting material, viz., YBa2Cu3O7 (YBCO), have been carried out. Electron temperature and electron density measurements were made from spectral data. The Stark broad ening of emission lines was used to determine the electron density, and the ratio of line in tensities was exploited for the determination of electron temperature. An initial electron temperature of 2.35 eV and electron density of 2.5 3 1017 cm2 3 were observed. The dependence on electron temperature and density on different experimental parameters such as distance from the target, delay time after the in itiation of the plasm a, and laser irradiance is also discussed in detail. Index Headings: Laser -plasma spectroscopy; Plasma diagnostics; Emission spectroscop y; YBa2Cu3O7.

Time evolution of the electron density and temperature in laser-produced plasmas from YBa2Cu3O7

Laser radiation at 1.06 µm from a pulsed Nd:YAG laser was focused onto a multielement YBa2Cu3O7 target in vacuum and the plasma thus generated was studied using time-resolved spectroscopic techniques. Line broadening of the Ba I emission line at 553.5 nm was monitored as a function of time elapsed after the incidence of a laser pulse on the target. Measured line profiles of barium species were used to infer the electron density and temperature, and the time evolution of these important plasma parameters has been worked out.

Electron density and temperature measurements in a laser produced carbon plasma

Plasma generated by fundamental radiation from a Nd:YAG laser focused onto a graphite target is studied spectroscopically. Measured line profiles of several ionic species were used to infer electron temperature and density at several sections located in front of the target surface. Line intensities of successive ionization states of carbon were used for electron temperature calculations. Stark broadened profiles of singly ionized species have been utilized for electron density measurements. Electron density as well as electron temperature were studied as functions of laser irradiance and time elapsed after the incidence of laser pulse. The validity of the assumption of local thermodynamic equilibrium is discussed in light of the results obtained.