PubMed Health⌕ Search

Biomedical subjects

Jan Seuntjens

Publications and source records attributed to Jan Seuntjens.

16 recordsLinked to original sources
PubMed2006

A direct voxel tracking method for four-dimensional Monte Carlo dose calculations in deforming anatomy.

In this work we present a method of calculating dose in deforming anatomy where the position and shape of each dose voxel is tracked as the anatomy changes. The EGSnrc/DOSXYZnrc Monte Carlo code was modified to calculate dose in voxels that are deformed according to deformation vectors obtained from a nonlinear image registration algorithm. The defDOSXYZ code was validated by consistency checks and by comparing calculations against DOSXYZnrc calculations. Calculations in deforming phantoms were compared with a dose remapping method employing trilinear interpolation. Dose calculations with the deforming voxels agree with DOSXYZnrc calculations within 1%. In simple deforming rectangular phantoms the trilinear dose remapping method was found to underestimate the dose by up to 29% for a 1.0 cm voxel size within the field, with larger discrepancies in regions of steep dose gradients. The agreement between the two calculation methods improved with smaller voxel size and deformation magnitude. A comparison of dose remapping from Inhale to Exhale in an anatomical breathing phantom demonstrated that dose deformations are underestimated by up to 16% in the penumbra and 8% near the surface with trilinear interpolation.

Algorithms↗
PubMed2005

Dosimetric and microdosimetric study of contrast-enhanced radiotherapy with kilovolt x-rays.

Kilovolt x-rays are clearly suboptimal compared to MV photon beams for radiotherapy of deep-seated tumours because of the increased attenuation in tissue, causing a rapid dose fall-off. This picture could change drastically when tumours can be labelled with contrast medium, containing high atomic number elements. This causes a significant dose enhancement to the tumour by exploiting the high cross sections for the photo-electric effect for kV x-rays. In this work, we have investigated the dosimetric and microdosimetric characteristics of kV contrast-enhanced radiation therapy (CERT) for different photon energies, contrast-medium concentrations and types (I and Gd). Two idealized patient treatment plans (head and lung) for irradiation with CT-arc beams were calculated. It is shown that the dose enhancement in tumours can be highly significant (up to about sixfold for realistic 80-120 kVp x-ray spectra and an iodine concentration of 50 mg ml-1) but that dose homogeneity in the tumour depends on photon energy, contrast-medium concentration and type, and irradiation scheme. An attempt to optimize the irradiation scheme is discussed. The microdosimetric study of the dose mean lineal energy shows that radiation quality changes in the contrast-medium-labelled region compared to homogeneous tissue are fairly small and limited to 10%. It is concluded that kV-CERT is a promising radiotherapy technique, provided contrast medium can be delivered reliably to tumours.

Body Burden↗
PubMed2005

Consistency test of the electron transport algorithm in the GEANT4 Monte Carlo code.

In this work, the condensed history algorithm in GEANT4 (version 4.6.2.p01) is examined. We performed simulations of an ionization chamber composed of water for 1.25 MeV incident photon beams under Fano conditions, and evaluated the consistency of the cavity response for several combinations of electron transport parameters. GEANT4 permits electrons to reach geometric boundaries in large steps, and underestimates lateral displacement near interfaces. Step size artefacts due to distortions in electron fluence and angular distributions reduce the cavity dose by up to 39%. Accurate cavity response can be achieved using severe user-imposed step size restrictions. We suggest that improvements in the electron transport algorithm in GEANT4 should address the handling of boundary crossing.

Algorithms↗
PubMed2005

Experimental determination of electron source parameters for accurate Monte Carlo calculation of large field electron therapy.

Extensive work has been performed to validate Monte Carlo models for both photon and electron beams under standard conditions. However, for large field electron beam therapy, Monte Carlo simulations have not been able to provide good agreement when compared to the measured dose distributions. Since the accuracy of the calculation relies heavily on the geometry parameters of the linear accelerator and the characteristics of the incident electron beam, it is crucial to have a complete comprehension of these independent factors. In this work, the electron focal spot size for a CL21EX linac with various energies (6, 9, 12 and 16 MeV) was measured with a slit camera composed of alternating lead and paper sheets. For all the energies investigated, the electron focal spot is found to be elliptical and has a full width at half maximum (FWHM) ranging from 1.69 mm to 2.24 mm. A shift with respect to the crosshair was associated with each measured focal spot. In addition, we present an improved result for the large field in-air profile by utilizing a proposed divergent beam model in conjunction with the experimental focal spot dimension. This model can potentially provide a solution to the Monte Carlo validation of large field electron beam therapy.

Electrons↗
PubMed2005

Monte Carlo based modulated electron beam treatment planning using a few-leaf electron collimator--feasibility study.

Energy modulated electron beam therapy with conventional clinical accelerators has lagged behind photon IMRT despite its potential to achieve highly conformal dose distributions in superficial targets. One of the reasons for this is the absence of an automated collimating device that allows for the flexible delivery of a series of variable field openings. Electron-specific multileaf collimators attached to the bottom of the applicator require the use of a large number of motors and suffer from being relatively bulky and impractical for head and neck sites. In this work, we investigate the treatment planning aspects of a proposed 'few-leaf' electron collimator (FLEC) that consists of four motor-driven trimmer bars at the end of the applicator. The device is designed to serve as an accessory to standard equipment and allows for the shaping of any irregular field by combination of rectangular fieldlets. Using a Monte Carlo model of the FLEC, dose distributions are optimized using a simulated annealing (SA) inverse planning algorithm based on a limited number of Monte Carlo pre-generated, realistic phantom-specific dose kernels and user-specified dose-volume constraints. Using a phantom setup with an artificial target enclosed by organs at risk (OAR) as well as using a realistic patient case, we demonstrate that highly conformal distributions can be generated. Estimates of delivery times are made and show that a full treatment fraction can be kept to 15 min or less.

Algorithms↗
PubMed2005

Precise radiochromic film dosimetry using a flat-bed document scanner.

In this study, a measurement protocol is presented that improves the precision of dose measurements using a flat-bed document scanner in conjunction with two new GafChromic film models, HS and Prototype A EBT exposed to 6 MV photon beams. We established two sources of uncertainties in dose measurements, governed by measurement and calibration curve fit parameters contributions. We have quantitatively assessed the influence of different steps in the protocol on the overall dose measurement uncertainty. Applying the protocol described in this paper on the Agfa Arcus II flat-bed document scanner, the overall one-sigma dose measurement uncertainty for an uniform field amounts to 2% or less for doses above around 0.4 Gy in the case of the EBT (Prototype A), and for doses above 5 Gy in the case of the HS model GafChromic film using a region of interest 2 X 2 mm2 in size.

Algorithms↗
PubMed2005

Energy modulated electron therapy using a few leaf electron collimator in combination with IMRT and 3D-CRT: Monte Carlo-based planning and dosimetric evaluation.

Energy modulated electron therapy (EMET) based on Monte Carlo dose calculation is a promising technique that enhances the treatment planning and delivery of superficially located tumors. This study investigated the application of EMET using a novel few-leaf electron collimator (FLEC) in head and neck and breast sites in comparison with three-dimensional conventional radiation therapy (3D-CRT) and intensity modulated radiation therapy (IMRT) techniques. Treatment planning was performed for two parotid cases and one breast case. Four plans were compared for each case: 3D-CRT, IMRT, 3D-CRT in conjunction with EMET (EMET-CRT), and IMRT in conjunction with EMET (EMET-IMRT), all of which were performed and calculated with Monte Carlo techniques. For all patients, dose volume histograms (DVHs) were obtained for all organs of interest and the DVHs were used as a means of comparing the plans. Homogeneity and conformity of dose distributions were calculated, as well as a sparing index that compares the effect of the low isodose lines. In addition, the whole-body dose equivalent (WBDE) was estimated for each plan. Adding EMET delivered with the FLEC to 3D-CRT improves sparing of normal tissues. For the two head and neck cases, the mean dose to the contralateral parotid and brain stem was reduced relative to IMRT by 43% and 84%, and by 57% and 71%, respectively. Improved normal tissue sparing was quantified as an increase in sparing index of 47% and 30% for the head and neck and the breast cases, respectively. Adding EMET to either 3D-CRT or IMRT results in preservation of target conformity and dose homogeneity. When adding EMET to the treatment plan, the WBDE was reduced by between 6% and 19% for 3D-CRT and by between 21% and 33% for IMRT, while WBDE for EMET-CRT was reduced by up to 72% when compared with IMRT. FLEC offers a practical means of delivering modulated electron therapy. Although adding EMET delivered using the FLEC results in perturbation of target conformity when compared to IMRT, it significantly improves normal tissue sparing while offering enhanced target conformity to the 3D-CRT planning. The addition of EMET systematically leads to a reduction in WBDE especially when compared with IMRT.

Breast Neoplasms↗
PubMed2005

Absorbed dose to water reference dosimetry using solid phantoms in the context of absorbed-dose protocols.

For reasons of phantom material reproducibility, the absorbed dose protocols of the American Association of Physicists in Medicine (AAPM) (TG-51) and the International Atomic Energy Agency (IAEA) (TRS-398) have made the use of liquid water as a phantom material for reference dosimetry mandatory. In this work we provide a formal framework for the measurement of absorbed dose to water using ionization chambers calibrated in terms of absorbed dose to water but irradiated in solid phantoms. Such a framework is useful when there is a desire to put dose measurements using solid phantoms on an absolute basis. Putting solid phantom measurements on an absolute basis has distinct advantages in verification measurements and quality assurance. We introduce a phantom dose conversion factor that converts a measurement made in a solid phantom and analyzed using an absorbed dose calibration protocol into absorbed dose to water under reference conditions. We provide techniques to measure and calculate the dose transfer from solid phantom to water. For an Exradin A12 ionization chamber, we measured and calculated the phantom dose conversion factor for six Solid Water phantoms and for a single Lucite phantom for photon energies between 60Co and 18 MV photons. For Solid Water of certified grade, the difference between measured and calculated factors varied between 0.0% and 0.7% with the average dose conversion factor being low by 0.4% compared with the calculation whereas for Lucite, the agreement was within 0.2% for the one phantom examined. The composition of commercial plastic phantoms and their homogeneity may not always be reproducible and consistent with assumed composition. By comparing measured and calculated phantom conversion factors, our work provides methods to verify the consistency of a given plastic for the purpose of clinical reference dosimetry.

Algorithms↗
PubMed2004

Accurate determination of dose-point-kernel functions close to the origin using Monte Carlo simulations.

Dose-point-kernel (DPK) functions are used extensively for the dosimetry of gamma and beta emitters in many physical problems. These functions are usually obtained from Monte Carlo simulations where the energy deposited in concentric spherical shells around a point source is tallied. The energy scored in a spherical shell divided by the shell mass is taken as the dose at some effective radius R(eff) of the shell. The effective radius R(eff), defined as the distance of a hypothetical zero-thickness scoring region from the source, can be evaluated in different ways for a finite thickness scoring region. For a shell thickness that is very small compared to the distance from the origin, this exact evaluation method becomes unimportant and the arithmetic mean is usually an accurate estimator for R(eff). However, accurately determining R(eff) can be problematic for the innermost regions when the radial dose function D(r) varies considerably over the finite spherical shell thickness. In this work, a new method for determining R(eff) is introduced which yields consistent results for any shell thickness, thus improving on previous Monte Carlo calculations for DPKs at or near the origin. Dimensionless DPK functions for monoenergetic electrons were reevaluated using EGSnrc with an emphasis on accuracy and consistency near the origin using our new method for determining R(eff). These improved functions were implemented in a software code to calculate the DPKs for an exhaustive list of 546 beta emitters, thus extending the compilation from previous works.

Computer Simulation↗
PubMed2004

Dosimetric properties of improved GafChromic films for seven different digitizers.

Two recently introduced GafChromic film models, HS and XR-T, have been developed as more sensitive and uniform alternatives to GafChromic MD-55-2 film. The HS model has been specifically designed for measurement of absorbed dose in high-energy photon beams (above 1 MeV), while the XR-T model has been introduced for dose measurements of low energy (0.1 MeV) photons. The goal of this study is to compare the sensitometric curves and estimated dosimetric uncertainties associated with seven different GafChromic film dosimetry systems for the two new film models. The densitometers tested are: LKB Pharmacia UltroScan XL, Molecular Dynamics Personal Densitometer, Nuclear Associates Radiochromic Densitometer Model 37-443, Photoelectron Corporation CMR-604, Laser Pro 16, Vidar VXR-16, and AGFA Arcus II document scanner. Pieces of film were exposed to different doses in a dose range from 0.5 to 50 Gy using 6 MV photon beam. Functional forms for dose vs net optical density have been determined for each of the GafChromic film-dosimetry systems used in this comparison. Two sources of uncertainties in dose measurements, governed by the experimental measurement and calibration curve fit procedure, have been compared for the densitometers used. Among the densitometers tested, it is found that for the HS film type the uncertainty caused by the experimental measurement varies from 1% to 3% while the calibration fit uncertainty ranges from 2% to 4% for doses above 5 Gy. Corresponding uncertainties for XR-T film model are somewhat higher and range from 1% to 5% for experimental and from 2% to 7% for the fit uncertainty estimates. Notwithstanding the significant variations in sensitivity, the studied densitometers exhibit very similar precision for GafChromic film based dose measurements above 5 Gy.

Equipment Failure Analysis↗
PubMed2004

Ionization chamber-based reference dosimetry of intensity modulated radiation beams.

The present paper addresses reference dose measurements using thimble ionization chambers for quality assurance in IMRT fields. In these radiation fields, detector fluence perturbation effects invalidate the application of open-field dosimetry protocol data for the derivation of absorbed dose to water from ionization chamber measurements. We define a correction factor C(Q)IMRT to correct the absorbed dose to water calibration coefficient N(D, w)Q for fluence perturbation effects in individual segments of an IMRT delivery and developed a calculation method to evaluate the factor. The method consists of precalculating, using accurate Monte Carlo techniques, ionization chamber, type-dependent cavity air dose, and in-phantom dose to water at the reference point for zero-width pencil beams as a function of position of the pencil beams impinging on the phantom surface. These precalculated kernels are convolved with the IMRT fluence distribution to arrive at the dose-to-water-dose-to-cavity air ratio [D(a)w (IMRT)] for IMRT fields and with a 10x10 cm2 open-field fluence to arrive at the same ratio D(a)w (Q) for the 10x10 cm2 reference field. The correction factor C(Q)IMRT is then calculated as the ratio of D(a)w (IMRT) and D(a)w (Q). The calculation method was experimentally validated and the magnitude of chamber correction factors in reference dose measurements in single static and dynamic IMRT fields was studied. The results show that, for thimble-type ionization chambers the correction factor in a single, realistic dynamic IMRT field can be of the order of 10% or more. We therefore propose that for accurate reference dosimetry of complete n-beam IMRT deliveries, ionization chamber fluence perturbation correction factors must explicitly be taken into account.

Algorithms↗
PubMed2004

Dosimetric evaluation of the clinical implementation of the first commercial IMRT Monte Carlo treatment planning system at 6 MV.

In this work we dosimetrically evaluated the clinical implementation of a commercial Monte Carlo treatment planning software (PEREGRINE, North American Scientific, Cranberry Township, PA) intended for quality assurance (QA) of intensity modulated radiation therapy treatment plans. Dose profiles calculated in homogeneous and heterogeneous phantoms using this system were compared to both measurements and simulations using the EGSnrc Monte Carlo code for the 6 MV beam of a Varian CL21EX linear accelerator. For simple jaw-defined fields, calculations agree within 2% of the dose at d(max) with measurements in homogeneous phantoms with the exception of the buildup region where the calculations overestimate the dose by up to 8%. In heterogeneous lung and bone phantoms the agreement is within 3%, on average, up to 5% for a 1 x 1 cm2 field. We tested two consecutive implementations of the MLC model. After matching the calculated and measured MLC leakage, simulations of static and dynamic MLC-defined fields using the most recent MLC model agreed to within 2% with measurements.

Algorithms↗
PubMed2003

Development and validation of a BEAMnrc component module for accurate Monte Carlo modelling of the Varian dynamic Millennium multileaf collimator.

A new component module (CM), designated DYNVMLC, was developed to fully model the geometry of the Varian Millennium 120 leaf collimator using the BEAMnrc Monte Carlo code. The model includes details such as the leaf driving screw hole, support railing groove and leaf tips. Further modifications also allow sampling of leaf sequence files to simulate the movement of the multileaf collimator (MLC) leaves during an intensity modulated radiation therapy (IMRT) delivery. As an initial validation of the code, the individual leaf geometries were visualized by tracing particles through the component module and recording their position each time a leaf boundary was crossed. A model of the Varian CL21EX linear accelerator 6 MV photon beam incorporating the new CM was built with the BEAMnrc user code. The leaf material density and abutting leaf air gap were chosen to match simulated leaf leakage profiles with film measurements in a solid water phantom. Simulated depth dose and off-axis profiles for a variety of MLC defined static fields agreed to within 2% with ion chamber and diode measurements in a water phantom. Simulated dose distributions for IMRT intensity patterns delivered using both static and dynamic techniques were found to agree with film measurements to within 4%. A comparison of interleaf leakage profiles for the new CM and an equivalent leaf model using the existing VARMLC CM demonstrated that the simplified geometry of VARMLC is not able to accurately predict the details of the MLC leakage for the 120 leaf collimator.

Computer Simulation↗
PubMed2003

Monte Carlo modelling of external radiotherapy photon beams.

An essential requirement for successful radiation therapy is that the discrepancies between dose distributions calculated at the treatment planning stage and those delivered to the patient are minimized. An important component in the treatment planning process is the accurate calculation of dose distributions. The most accurate way to do this is by Monte Carlo calculation of particle transport, first in the geometry of the external or internal source followed by tracking the transport and energy deposition in the tissues of interest. Additionally, Monte Carlo simulations allow one to investigate the influence of source components on beams of a particular type and their contaminant particles. Since the mid 1990s, there has been an enormous increase in Monte Carlo studies dealing specifically with the subject of the present review, i.e., external photon beam Monte Carlo calculations, aided by the advent of new codes and fast computers. The foundations for this work were laid from the late 1970s until the early 1990s. In this paper we will review the progress made in this field over the last 25 years. The review will be focused mainly on Monte Carlo modelling of linear accelerator treatment heads but sections will also be devoted to kilovoltage x-ray units and 60Co teletherapy sources.

Animals↗
PubMed2003

Re-evaluation of the dose to the cyst wall in P-32 radiocolloid treatments of cystic brain tumors using the dose-point-kernel and Monte Carlo methods.

Intracavity instillation of beta-emitting colloid pharmaceuticals is a common technique used to treat cystic brain tumors. Most of the dosimetric calculations that have been reported in the literature for this problem are based on empirical formulas derived by Loevinger. Concentration of P-32 radiolabeled solution for the delivery of a prescribed dose (200 Gy to the cyst wall) has been published previously using this formalism in what we refer to as a standard nomogram. The calculations using the Loevinger formulas for calculating the P-32 activity necessary to achieve 200 Gy at the cyst wall is re-evaluated and compared to numerically computed results based on full Monte Carlo simulations (EGSnrc) and the dose-point-kernel (DPK) integration method. For cyst diameters greater than 1 cm, the new calculations agree well with previously published results (the standard nomogram) to within a few percents. However, for cyst diameters of less than 1 cm, it is shown that the standard nomogram results underestimate the therapeutic activity by a factor of approximately 3 for very small diameters (approximately 0.2 cm). New tables based on our calculations are presented and the sources of discrepancies are identified. It is concluded that the new set of data based on our calculations should replace the standard nomogram to administer accurately the target dose to the cyst wall for the smaller diameter cysts (< 1 cm).

Beta Particles↗