Tsi-Chian Chao

Vip-man: An Image-based Whole-body Adult Male Model Constructed From Color Photographs Of The Visible Human Project For Multi-particle Monte Carlo Calculations

Human anatomical models have been indispensable to radiation protection dosimetry using Monte Carlo calculations. Existing MIRD-based mathematical models are easy to compute and standardize, but they are simplified and crude compared to human anatomy. This article describes the development of an image-based whole-body model, called VIP-Man, using transversal color photographic images obtained from the National Library of Medicine’s Visible Human Project for Monte Carlo organ dose calculations involving photons, electron, neutrons, and protons. As the first of a series of papers on dose calculations based on VIP-Man, this article provides detailed information about how to construct an image-based model, as well as how to adopt it into well-tested Monte Carlo codes, EGS4, MCNP4B, and MCNPX.

Specific absorbed fractions from the image-based VIP-Man body model and EGS4-VLSI Monte Carlo code : internal electron emitters

VIP-Man is a whole-body anatomical model newly developed at Rensselaer from the high-resolution colour images of the National Library of Medicines Visible Human Project. This paper summarizes the use of VIP-Man and the Monte Carlo method to calculate specific absorbed fractions from internal electron emitters. A specially designed EGS4 user code, named EGS4-VLSI, was developed to use the extremely large number of image data contained in the VIP-Man. Monoenergetic and isotropic electron emitters with energies from 100 keV to 4 MeV are considered to be uniformly distributed in 26 organs. This paper presents, for the first time, results of internal electron exposures based on a realistic whole-body tomographic model. Because VIP-Man has many organs and tissues that were previously not well defined (or not available) in other models, the efforts at Rensselaer and elsewhere bring an unprecedented opportunity to significantly improve the internal dosimetry.

Fluence-to-dose conversion coefficients from monoenergetic neutrons below 20 MeV based on the VIP-Man anatomical model

A new set of fluence-to-absorbed dose and fluence-to-effective dose conversion coefficients have been calculated for neutrons below 20 MeV using a whole-body anatomical model, VIP-Man, developed from the high-resolution transverse colour photographic images of the National Library of Medicines Visible Human Project. Organ dose calculations were performed using the Monte Carlo code MCNP for 20 monoenergetic neutron beams between 1 x 10(-9) MeV and 20 MeV under six different irradiation geometries: anterior-posterior, posterior-anterior, right lateral, left lateral, rotational and isotropic. The absorbed dose for 24 major organs and effective dose results based on the realistic VIP-Man are presented and compared with those based on the simplified MIRD-based phantoms reported in the literature. Effective doses from VIP-Man are not significantly different from earlier results for neutrons in the energy range studied. There are, however, remarkable deviations in organ doses due to the anatomical differences between the image-based and the earlier mathematical models.

Fluence-to-dose conversion coefficients based on the VIP-Man anatomical model and MCNPX code for monoenergetic neutrons above 20 MeV.

A new set of fluence-to-absorbed dose and fluence-to-effective dose conversion coefficients has been calculated for high-energy neutrons using a whole-body anatomical model, VIP-Man, developed from the high-resolution transversal color photographic images of the National Library of Medicines Visible Human Project. Organ dose calculations were performed using the Monte Carlo code MCNPX for 20 monoenergetic neutron beams between 20 MeV and 10,000 MeV under 6 different irradiation geometries: anterior-posterior, posterior-anterior, left lateral, right lateral, isotropic, and rotational. For neutron Monte Carlo calculations, results based on an image-based whole-body model were not available in the literature. The absorbed dose results for 24 major organs of VIP-Man are presented in the form of tables and selected figures that compare with those based on simplified mathematical phantoms reported in the literature. VIP-Man yields up to 40% larger values of effective dose and many organ doses, thus suggesting that the results reported in the past may not be conservative.

Organ dose conversion coefficients for 0.1-10 MeV electrons calculated for the VIP-man tomographic model

A whole-body tomographic model, called VIP-Man, was recently developed at Rensselaer Polytechnic Institute from the high-resolution color photographic images of the National Library of Medicines Visible Human Project. An EGS4-based Monte Carlo user code, named EGS4-VLSI, was developed to efficiently transport electrons using the large image data set for VIP-Man. VIP-Man has been used to calculate doses for neutrons and photons. This paper presents a new set of fluence-to-absorbed-dose conversion coefficients for monoenergetic electron beams between 100 keV and 10 MeV for VIP-Man. Irradiation conditions include anterior-posterior, posterior-anterior, right lateral, left lateral, rotational, and isotropic source geometries. Comparisons between organ doses from VIP-Man, which is taller and heavier than the Reference Man, and existing data from mathematical models show significant discrepancies. It appears that even slight differences between body models can cause dramatic dosimetric deviations for low penetrating electron irradiation. This suggests that a single standard body model may poorly represent a large population and may not be acceptable for electron dosimetry.

Use of the VIP-Man model to calculate energy imparted and effective dose for x-ray examinations.

Abstract— A male human tomographic model was used to calculate values of energy imparted (&egr;) and effective dose (E) for monoenergetic photons (30–150 keV) in radiographic examinations. Energy deposition in the organs and tissues of the human phantom were obtained using Monte Carlo simulations. Values of E/&egr; were obtained for three common projections [anterior-posterior (AP), posterior-anterior (PA), and lateral (LAT)] of the head, cervical spine, chest, and abdomen, respectively. For head radiographs, all three projections yielded similar E/&egr; values. At 30 keV, the value of E/&egr; was ∼1.6 mSv J−1, which is increased to ∼7 mSv J−1 for 150 keV photons. The AP cervical spine was the only projection investigated where the value of E/&egr; decreased with increasing photon energy. Above 70 keV, cervical spine E/&egr; values showed little energy dependence and ranged between ∼8.5 mSv J−1 for PA projections and ∼17 mSv J−1 for AP projections. The values of E/&egr; for AP chest examinations showed very little variation with photon energy, and had values of ∼23 mSv J−1. Values of E/&egr; for PA and LAT chest projections were substantially lower than the AP projections and increased with increasing photon energy. For abdominal radiographs, differences between the PA and LAT projections were very small. All abdomen projections showed an increase in the E/&egr; ratio with increasing photon energy, and reached a maximum value of ∼13.5 mSv J−1 for AP projections, and ∼9.5 mSv J−1 for PA/lateral projections. These monoenergetic E/&egr; values can generate values of E/&egr; for any x-ray spectrum, and can be used to convert values of energy imparted into effective dose for patients undergoing common head and body radiological examinations.

S-values calculated from a tomographic head/brain model for brain imaging

A tomographic head/brain model was developed from the Visible Human images and used to calculate S-values for brain imaging procedures. This model contains 15 segmented sub-regions including caudate nucleus, cerebellum, cerebral cortex, cerebral white matter, corpus callosum, eyes, lateral ventricles, lenses, lentiform nucleus, optic chiasma, optic nerve, pons and middle cerebellar peduncle, skull CSF, thalamus and thyroid. S-values for C-11, O-15, F-18, Tc-99m and I-123 have been calculated using this model and a Monte Carlo code, EGS4. Comparison of the calculated S-values with those calculated from the MIRD (1999) stylized head/brain model shows significant differences. In many cases, the stylized head/brain model resulted in smaller S-values (as much as 88%), suggesting that the doses to a specific patient similar to the Visible Man could have been underestimated using the existing clinical dosimetry.

Measurement‐based Monte Carlo dose calculation system for IMRT pretreatment and on‐line transit dose verifications

The aim of this study was to develop a dose simulation system based on portal dosimetry measurements and the BEAM Monte Carlo code for intensity-modulated (IM) radiotherapy dose verification. This measurement-based Monte Carlo (MBMC) system can perform, within one systematic calculation, both pretreatment and on-line transit dose verifications. BEAMnrc and DOSXYZnrc 2006 were used to simulate radiation transport from the treatment head, through the patient, to the plane of the aS500 electronic portal imaging device (EPID). In order to represent the nonuniform fluence distribution of an IM field within the MBMC simulation, an EPID-measured efficiency map was used to redistribute particle weightings of the simulated phase space distribution of an open field at a plane above a patient/phantom. This efficiency map was obtained by dividing the measured energy fluence distribution of an IM field to that of an open field at the EPID plane. The simulated dose distribution at the midplane of a homogeneous polystyrene phantom was compared to the corresponding distribution obtained from the Eclipse treatment planning system (TPS) for pretreatment verification. It also generated a simulated transit dose distribution to serve as the on-line verification reference for comparison to that measured by the EPID. Two head-and-neck (NPC1 and NPC2) and one prostate cancer fields were tested in this study. To validate the accuracy of the MBMC system, film dosimetry was performed and served as the dosimetry reference. Excellent agreement between the film dosimetry and the MBMC simulation was obtained for pretreatment verification. For all three cases tested, gamma evaluation with 3%/3 mm criteria showed a high pass percentage (> 99.7%) within the area in which the dose was greater than 30% of the maximum dose. In contrast to the TPS, the MBMC system was able to preserve multileaf collimator delivery effects such as the tongue-and-groove effect and interleaf leakage. In the NPC1 field, the TPS showed 16.5% overdose due to the tongue-and-groove effect and 14.6% overdose due to improper leaf stepping. Similarly, in the NPC2 field, the TPS showed 14.1% overdose due to the tongue-and-groove effect and 8.9% overdose due to improper leaf stepping. In the prostate cancer field, the TPS showed 6.8% overdose due to improper leaf stepping. No tongue-and-groove effect was observed for this field. For transit dose verification, agreements among the EPID measurement, the film dosimetry, and the MBMC system were also excellent with a minimum gamma pass percentage of 99.6%.

Novel diffusion anisotropy indices: An evaluation

To systematically evaluate diffusion anisotropy (DA) using newly defined indices based on the diffusion deviation and mean diffusivity approach.

Monte Carlo simulations of the relative biological effectiveness for DNA double strand breaks from 300 MeV u(-1) carbon-ion beams.

Monte Carlo simulations are used to calculate the relative biological effectiveness (RBE) of 300 MeV u(-1) carbon-ion beams at different depths in a cylindrical water phantom of 10 cm radius and 30 cm long. RBE values for the induction of DNA double strand breaks (DSB), a biological endpoint closely related to cell inactivation, are estimated for monoenergetic and energy-modulated carbon ion beams. Individual contributions to the RBE from primary ions and secondary nuclear fragments are simulated separately. These simulations are based on a multi-scale modelling approach by first applying the FLUKA (version 2011.2.17) transport code to estimate the absorbed doses and fluence energy spectra, then using the MCDS (version 3.10A) damage code for DSB yields. The approach is efficient since it separates the non-stochastic dosimetry problem from the stochastic DNA damage problem. The MCDS code predicts the major trends of the DSB yields from detailed track structure simulations. It is found that, as depth is increasing, RBE values increase slowly from the entrance depth to the plateau region and change substantially in the Bragg peak region. RBE values reach their maxima at the distal edge of the Bragg peak. Beyond this edge, contributions to RBE are entirely from nuclear fragments. Maximum RBE values at the distal edges of the Bragg peak and the spread-out Bragg peak are, respectively, 3.0 and 2.8. The present approach has the flexibility to weight RBE contributions from different DSB classes, i.e. DSB0, DSB+ and DSB++.