Biomedical subjects

E Fourkal

Publications and source records attributed to E Fourkal.

9 recordsLinked to original sources

PubMed2005

Coulomb explosion effect and the maximum energy of protons accelerated by high-power lasers.

The acceleration of light ions (protons) through the interaction of a high-power laser pulse with a double-layer target is theoretically studied by means of two-dimensional particle-in-cell simulations and a one-dimensional analytical model. It is shown that the maximum energy acquired by the accelerated light ions (protons) depends on the physical characteristics of a heavy-ion layer (electron-ion mass ratio and effective charge state of the ions). In our theoretical model, the hydrodynamic equations for both electron and heavy-ion species are solved and the test-particle approximation for the light ions (protons) is applied. The heavy-ion motion is found to modify the longitudinal electric field distribution, thus changing the acceleration conditions for the protons.

Biophysics↗

PubMed2005

A particle track-repeating algorithm for proton beam dose calculation.

A particle track-repeating algorithm has been developed for proton beam dose calculation for radiotherapy. Monoenergetic protons with 250 MeV kinetic energy were simulated in an infinite water phantom using the GEANT3 Monte Carlo code. The changes in location, angle and energy for every transport step and the energy deposition along the track were recorded for the primary protons and all secondary particles. When calculating dose for a patient with a realistic proton beam, the pre-generated particle tracks were repeated in the patient geometry consisting of air, soft tissue and bone. The medium and density for each dose scoring voxel in the patient geometry were derived from patient CT data. The starting point, at which a proton track was repeated, was determined according to the incident proton energy. Thus, any protons with kinetic energy less than 250 MeV can be simulated. Based on the direction of the incident proton, the tracks were first rotated and for the subsequent steps, the scattering angles were simply repeated for air and soft tissue but adjusted properly based on the scattering power for bone. The particle step lengths were adjusted based on the density for air and soft tissue and also on the stopping powers for bone while keeping the energy deposition unchanged in each step. The difference in nuclear interactions and secondary particle generation between water and these materials was ignored. The algorithm has been validated by comparing the dose distributions in uniform water and layered heterogeneous phantoms with those calculated using the GEANT3 code for 120, 150, 180 and 250 MeV proton beams. The differences between them were within 2%. The new algorithm was about 13 times faster than the GEANT3 Monte Carlo code for a uniform phantom geometry and over 700 times faster for a heterogeneous phantom geometry.

Algorithms↗

PubMed2004

Monitor unit calculation for Monte Carlo treatment planning.

In this work, we investigate a formalism for monitor unit (MU) calculation in Monte Carlo based treatment planning. By relating MU to dose measured under reference calibration conditions (central axis, depth of dose maximum in water, 10 cm x 10 cm field defined at 100 cm source-to-surface distance) our formalism determines the MU required for a treatment plan based on the prescription dose and Monte Carlo calculated dose distribution. Detailed descriptions and formulae are given for various clinical situations including conventional treatments and advanced techniques such as intensity-modulated radiotherapy (IMRT) and modulated electron radiotherapy (MERT). Analysis is made of the effects of source modelling, beam modifier simulation and patient dose calculation accuracy, all of which are important factors for absolute dose calculations using Monte Carlo simulations. We have tested the formalism through phantom measurements and the predicted MU values were consistent with measured values to within 2%. The formalism has been used for MU calculation and plan comparison for advanced treatment techniques such as MERT, extracranial stereotactic IMRT, MRI-based treatment planning and intensity-modulated laser-proton therapy studies. It is also used for absolute dose calculations using Monte Carlo simulations for treatment verification, which has become part of our comprehensive IMRT quality assurance programme.

Dose-Response Relationship, Radiation↗

PubMed2004

Stereotactic IMRT for prostate cancer: dosimetric impact of multileaf collimator leaf width in the treatment of prostate cancer with IMRT.

The focus of this work is the dosimetric impact of multileaf collimator (MLC) leaf width on the treatment of prostate cancer with intensity-modulated radiation therapy (IMRT). Ten patients with prostate cancer were planned for IMRT delivery using two different MLC leaf widths--4mm and 10mm--representing the Radionics micro-multileaf collimator (mMLC) and Siemens MLC, respectively. Treatment planning was performed on the XKnifeRT2 treatment-planning system (Radionics, Burlington, MA). All beams and optimization parameters were identical for the mMLC and MLC plans. All the plans were normalized to ensure that 95% of the planning target volume (PTV) received 100% of the prescribed dose. The differences in dose distribution between the two different plans were assessed by dose-volume histogram (DVH) analysis of the target and critical organs. We specifically compared the volume of rectum receiving 40 Gy (V40), 50 Gy (V50), 60 Gy (V60), the dose received by 17% and 35% of rectum (D17 and D35), and the maximum dose to 1 cm3 of the rectum for a prescription dose of 74 Gy. For the urinary bladder, the dose received by 25% of bladder (D25), V40, and the maximum dose to 1 cm3 of the organ were recorded. For PTV we compared the maximum dose to the "hottest" 1 cm3 (Dmax1 cm3) and the dose to 99% of the PTV (D99). The dose inhomogeneity in the target, defined as the ratio of the difference in Dmax1 cm3 and D99 to the prescribed dose, was also compared between the two plans. In all cases studied, significant reductions in the volume of rectum receiving doses less than 65 Gy were seen using the mMLC. The average decrease in the volume of the rectum receiving 40 Gy, 50 Gy, and 60 Gy using the mMLC plans was 40.2%, 33.4%, and 17.7%, respectively, with p < 0.0001 for V40 and V50 and p < 0.012 for V60. The mean dose reductions for D17 and D35 for the rectum using the mMLC were 20.4% (p < 0.0001) and 18.3% (p < 0.0002), respectively. There were consistent reductions in all dose indices studied for the bladder. The target dose inhomogeneity was improved in the mMLC plans by an average of 29%. In the high-dose range, there was no significant difference in the dose deposited in the "hottest" 1 cm3 of the rectum between the two plans for all cases (p > 0.78). In conclusion, the use of the mMLC for IMRT of the prostate resulted in significant improvement in the DVH parameters of the prostate and critical organs, which may improve the therapeutic ratio.

Equipment Design↗

PubMed2003

Intensity modulated radiation therapy using laser-accelerated protons: a Monte Carlo dosimetric study.

In this paper we present Monte Carlo studies of intensity modulated radiation therapy using laser-accelerated proton beams. Laser-accelerated protons coming out of a solid high-density target have broad energy and angular spectra leading to dose distributions that cannot be directly used for therapeutic applications. Through the introduction of a spectrometer-like particle selection system that delivers small pencil beams of protons with desired energy spectra it is feasible to use laser-accelerated protons for intensity modulated radiotherapy. The method presented in this paper is a three-dimensional modulation in which the proton energy spectrum and intensity of each individual beamlet are modulated to yield a homogeneous dose in both the longitudinal and lateral directions. As an evaluation of the efficacy of this method, it has been applied to two prostate cases using a variety of beam arrangements. We have performed a comparison study between intensity modulated photon plans and those for laser-accelerated protons. For identical beam arrangements and the same optimization parameters, proton plans exhibit superior coverage of the target and sparing of neighbouring critical structures. Dose-volume histogram analysis of the resulting dose distributions shows up to 50% reduction of dose to the critical structures. As the number of fields is decreased, the proton modality exhibits a better preservation of the optimization requirements on the target and critical structures. It is shown that for a two-beam arrangement (parallel-opposed) it is possible to achieve both superior target coverage with 5% dose inhomogeneity within the target and excellent sparing of surrounding tissue.

Algorithms↗

PubMed2003

Particle selection for laser-accelerated proton therapy feasibility study.

In this paper we present calculations for the design of a particle selection system for laser-accelerated proton therapy. Laser-accelerated protons coming from a thin high-density foil have broad energy and angular spectra leading to dose distributions that cannot be directly used for therapeutic applications. Our solution to this problem is a compact particle selection and collimation device that delivers small pencil beams of protons with desired energy spectra. We propose a spectrometer-like particle selection and beam modulation system in which the magnetic field will be used to spread the protons spatially according to their energies and emitting angles. Subsequently, an aperture will be used to select the protons within a therapeutic window of energy (energy modulation). It will be shown that for the effective proton spatial differentiation, the primary collimation device should be used, which will collimate protons to the desired angular distribution and limit the spatial mixing of different energy protons once they have traveled through the magnetic system. Due to the angular proton distribution, the spatial mixing of protons of different energies will always be present and it will result in a proton energy spread with the width depending on the energy. For 250 MeV protons, the width (from the maximum to the minimum energy) is found to be 50 MeV for the magnetic field configuration used in our calculations. As the proton energy decreases, its energy width decreases as well, and for 80 MeV protons it equals 9 MeV. The presence of the energy width in the proton energy distribution will modify the depth dose curves needed for the energy modulation calculation. The matching magnetic field setup will ensure the refocusing of the selected protons and the final beam will be collimated by the secondary collimator. The calculations presented in this article show that the dose rate that the selection system can yield is on the order of D=260 Gy/min for a field size of 1 x 1 cm2.

Equipment Design↗

PubMed2002

Particle in cell simulation of laser-accelerated proton beams for radiation therapy.

In this article we present the results of particle in cell (PIC) simulations of laser plasma interaction for proton acceleration for radiation therapy treatments. We show that under optimal interaction conditions protons can be accelerated up to relativistic energies of 300 MeV by a petawatt laser field. The proton acceleration is due to the dragging Coulomb force arising from charge separation induced by the ponderomotive pressure (light pressure) of high-intensity laser. The proton energy and phase space distribution functions obtained from the PIC simulations are used in the calculations of dose distributions using the GEANT Monte Carlo simulation code. Because of the broad energy and angular spectra of the protons, a compact particle selection and beam collimation system will be needed to generate small beams of polyenergetic protons for intensity modulated proton therapy.

Computer Simulation↗

PubMed2001

Detective quantum efficiency of a direct-detection active matrix flat panel imager at megavoltage energies.

The use of an amorphous selenium (a-Se) based direct-detection active matrix flat-panel imager (AMFPI) is studied for megavoltage imaging. The detector consists of a 1.2 mm copper front plate and 200 microm a-Se layer, and has a 85 microm pixel pitch. The Modulation Transfer Function (MTF), Noise Power Spectrum (NPS), and Detective Quantum Efficiency (DQE) are measured for 6 and 15 MV photon beams. A theoretical expression for the DQE is derived using a recently developed formalism for nonelementary cascade stages. A comparison of theory with experiment is good for the 6 and 15 MV beams. The model is used to explore the DQE for more typical pixel sizes. The results indicate that with proper modifications, such as a larger a-Se thickness, a direct flat-panel AMFPI is a very promising detector for megavoltage imaging.

Electrons↗

PubMed2001

Investigation into the physical characteristics of active matrix flat panel imagers for radiotherapy.

The effect of physical characteristics on active matrix flat panel imagers (AMFPIs) at megavoltage energies is studied. The detective quantum efficiency (DQE) of both direct and indirect AMFPIs is modeled using a modified cascade analysis combined with Monte Carlo simulations. It is found that for a given thickness of the sensitive layer less than about 1 mm, there should be no significant difference between the detection techniques, but for larger mass thicknesses (> or = 1 mm) there should be an advantage to using direct detection if such thick layers can be manufactured. The effect of the front plate on both direct and indirect techniques is also explored in terms of both the DQE and scatter rejection. It is found that for small sensitive layer thicknesses (< or =0.3 mm) a front plate thickness of about 1 mm Cu is optimal, whereas for larger mass thicknesses about 0.4 mm Cu should lead to better image quality.

Computer Simulation↗