A peripheral dose calculation model for estimating second cancer risk after breast radiotherapy

AbstractBreast cancer is one of the most common cancers affecting one in eight women during their lives. Survival rates have increased steadily thanks to early diagnosis with mammography screening and more efficient treatment strategies. Post-operative radiation therapy is a standard of care in the management of breast cancer and has been shown to reduce efficiently both local recurrence rate and breast cancer mortality. Radiation therapy is however associated with some late effects for long-term survivors. Radiation-induced secondary cancer is a relatively rare but severe late effect of radiation therapy. Currently, radiotherapy plans are essentially optimized to maximize tumor control and minimize late deterministic effects (tissue reactions) that are mainly associated with high doses (» 1 Gy). With improved cure rates and new radiation therapy technologies, it is also important to evaluate and minimize secondary cancer risks for different treatment techniques. This is a particularly challenging task due to the large uncertainties in the dose-response relationship.In contrast with late deterministic effects, secondary cancers may be associated with much lower doses and therefore out-of-field doses (also called peripheral doses) that are typically inferior to 1 Gy need to be determined accurately. Out-of-field doses result from patient scatter and head scatter from the treatment unit. These doses are particularly challenging to compute and we characterized it by Monte Carlo (MC) calculation. A detailed MC model of the Siemens Primus linear accelerator has been thoroughly validated with measurements. We investigated the accuracy of such a model for retrospective dosimetry in epidemiological studies on secondary cancers. Considering that patients in such large studies could be treated on a variety of machines, we assessed the uncertainty in reconstructed peripheral dose due to the variability of peripheral dose among various linac geometries. For large open fields (> 10x10 cm2), the uncertainty would be less than 50%, but for small fields and wedged fields the uncertainty in reconstructed dose could rise up to a factor of 10. It was concluded that such a model could be used for conventional treatments using large open fields only.The MC model of the Siemens Primus linac was then used to compare out-of-field doses for different treatment techniques in a female whole-body CT-based phantom. Current techniques such as conformai wedged-based radiotherapy and hybrid IMRT were investigated and compared to older two-dimensional radiotherapy techniques. MC doses were also compared to those of a commercial Treatment Planning System (TPS). While the TPS is routinely used to determine the dose to the contralateral breast and the ipsilateral lung which are mostly out of the treatment fields, we have shown that these doses may be highly inaccurate depending on the treatment technique investigated. MC shows that hybrid IMRT is dosimetrically similar to three-dimensional wedge-based radiotherapy within the field, but offers substantially reduced doses to out-of-field healthy organs.Finally, many different approaches to risk estimations extracted from the literature were applied to the calculated MC dose distribution. Absolute risks varied substantially as did the ratio of risk between two treatment techniques, reflecting the large uncertainties involved with current risk models. Despite all these uncertainties, the hybrid IMRT investigated resulted in systematically lower cancer risks than any of the other treatment techniques. More epidemiological studies with accurate dosimetry are required in the future to construct robust risk models. In the meantime, any treatment strategy that reduces out-of-field doses to healthy organs should be investigated. Electron radiotherapy might offer interesting possibilities with this regard.RésuméLe cancer du sein affecte une femme sur huit au cours de sa vie. Grâce au dépistage précoce et à des thérapies de plus en plus efficaces, le taux de guérison a augmenté au cours du temps. La radiothérapie postopératoire joue un rôle important dans le traitement du cancer du sein en réduisant le taux de récidive et la mortalité. Malheureusement, la radiothérapie peut aussi induire des toxicités tardives chez les patients guéris. En particulier, les cancers secondaires radio-induits sont une complication rare mais sévère de la radiothérapie. En routine clinique, les plans de radiothérapie sont essentiellement optimisées pour un contrôle local le plus élevé possible tout en minimisant les réactions tissulaires tardives qui sont essentiellement associées avec des hautes doses (» 1 Gy). Toutefois, avec l'introduction de différentes nouvelles techniques et avec l'augmentation des taux de survie, il devient impératif d'évaluer et de minimiser les risques de cancer secondaire pour différentes techniques de traitement. Une telle évaluation du risque est une tâche ardue étant donné les nombreuses incertitudes liées à la relation dose-risque.Contrairement aux effets tissulaires, les cancers secondaires peuvent aussi être induits par des basses doses dans des organes qui se trouvent hors des champs d'irradiation. Ces organes reçoivent des doses périphériques typiquement inférieures à 1 Gy qui résultent du diffusé du patient et du diffusé de l'accélérateur. Ces doses sont difficiles à calculer précisément, mais les algorithmes Monte Carlo (MC) permettent de les estimer avec une bonne précision. Un modèle MC détaillé de l'accélérateur Primus de Siemens a été élaboré et validé avec des mesures. La précision de ce modèle a également été déterminée pour la reconstruction de dose en épidémiologie. Si on considère que les patients inclus dans de larges cohortes sont traités sur une variété de machines, l'incertitude dans la reconstruction de dose périphérique a été étudiée en fonction de la variabilité de la dose périphérique pour différents types d'accélérateurs. Pour de grands champs (> 10x10 cm ), l'incertitude est inférieure à 50%, mais pour de petits champs et des champs filtrés, l'incertitude de la dose peut monter jusqu'à un facteur 10. En conclusion, un tel modèle ne peut être utilisé que pour les traitements conventionnels utilisant des grands champs.Le modèle MC de l'accélérateur Primus a été utilisé ensuite pour déterminer la dose périphérique pour différentes techniques dans un fantôme corps entier basé sur des coupes CT d'une patiente. Les techniques actuelles utilisant des champs filtrés ou encore l'IMRT hybride ont été étudiées et comparées par rapport aux techniques plus anciennes. Les doses calculées par MC ont été comparées à celles obtenues d'un logiciel de planification commercial (TPS). Alors que le TPS est utilisé en routine pour déterminer la dose au sein contralatéral et au poumon ipsilatéral qui sont principalement hors des faisceaux, nous avons montré que ces doses peuvent être plus ou moins précises selon la technTque étudiée. Les calculs MC montrent que la technique IMRT est dosimétriquement équivalente à celle basée sur des champs filtrés à l'intérieur des champs de traitement, mais offre une réduction importante de la dose aux organes périphériques.Finalement différents modèles de risque ont été étudiés sur la base des distributions de dose calculées par MC. Les risques absolus et le rapport des risques entre deux techniques de traitement varient grandement, ce qui reflète les grandes incertitudes liées aux différents modèles de risque. Malgré ces incertitudes, on a pu montrer que la technique IMRT offrait une réduction du risque systématique par rapport aux autres techniques. En attendant des données épidémiologiques supplémentaires sur la relation dose-risque, toute technique offrant une réduction des doses périphériques aux organes sains mérite d'être étudiée. La radiothérapie avec des électrons offre à ce titre des possibilités intéressantes.

Decreased local control following radiation therapy alone in early-stage glottic carcinoma with anterior commissure extension.

PURPOSE: To assess the patterns of failure in the treatment of early-stage squamous cell carcinoma of the glottic larynx. PATIENTS AND METHODS: Between 1983-2000, 122 consecutive patients treated for early laryngeal cancer (UICC T1N0 and T2N0) by radical radiation therapy (RT) were retrospectively studied. Male-to-female ratio was 106 : 16, and median age 62 years (35-92 years). There were 68 patients with T1a, 18 with T1b, and 36 with T2 tumors. Diagnosis was made by biopsy in 104 patients, and by laser vaporization or stripping in 18. Treatment planning consisted of three-dimensional (3-D) conformal RT in 49 (40%) patients including nine patients irradiated using arytenoid protection. A median dose of 70 Gy (60-74 Gy) was given (2 Gy/fraction) over a median period of 46 days (21-79 days). Median follow-up period was 85 months. RESULTS: The 5-year overall, cancer-specific, and disease-free survival amounted to 80%, 94%, and 70%, respectively. 5-year local control was 83%. Median time to local recurrence in 19 patients was 13 months (5-58 months). Salvage treatment consisted of surgery in 17 patients (one patient refused salvage and one was inoperable; total laryngectomy in eleven, and partial laryngectomy or cordectomy in six patients). Six patients died because of laryngeal cancer. Univariate analyses revealed that prognostic factors negatively influencing local control were anterior commissure extension, arytenoid protection, and total RT dose < 66 Gy. Among the factors analyzed, multivariate analysis (Cox model) demonstrated that anterior commissure extension, arytenoid protection, and male gender were the worst independent prognostic factors in terms of local control. CONCLUSION: For early-stage laryngeal cancer, outcome after RT is excellent. In case of anterior commissure extension, surgery or higher RT doses are warranted. Because of a high relapse risk, arytenoid protection should not be attempted.

Traitement non chirurgical des rectites postactiniques ou radiques chroniques [Non-surgical treatment of acute and late radiation proctitis]

Pelvic external radiotherapy with or without brachytherapy plays an important role in the management of pelvic cancers. Despite recent technical innovations including conformal three-dimensional (3D) external beam radiotherapy and more recently intensity modulated radiotherapy (IMRT), local side effects can occur secondary to normal tissue damage caused by ionising radiation. Morbidity depends on the anatomic position of the rectum within the pelvis and the fast turnover rate of the mucosa, as well as the characteristics of the radiation treatment and patient co-morbidities. Medical management is sometimes complex and merits herein a short review.

Dummy run and conformity indices in the ongoing EORTC low-grade glioma trial 22033-26033: First evaluation of quality of radiotherapy planning.

PURPOSE: Early assessment of radiotherapy (RT) quality in the ongoing EORTC trial comparing primary temozolomide versus RT in low-grade gliomas. MATERIALS AND METHODS: RT plans provided for dummy cases were evaluated and compared against expert plans. We analysed: (1) tumour and organs-at-risk delineation, (2) geometric and dosimetric characteristics, (3) planning parameters, compliance with dose prescription and Dmax for OAR (4) indices: RTOG conformity index (CI), coverage factor (CF), tissue protection factor (PF); conformity number (CN = PF x CF); dose homogeneity in PTV (U). RESULTS: Forty-one RT plans were evaluated. Only two (5%) centres were requested to repeat CTV-PTV delineations. Three (7%) plans had a significant under-dosage and dose homogeneity in one deviated > 10%. Dose distribution was good with mean values of 1.5, 1, 0.68, and 0.68 (ideal values = 1) for CI, CF, PF, and CN, respectively. CI and CN strongly correlated with PF and they correlated with PTV. Planning with more beams seems to increase PTV(Dmin), improving CF. U correlated with PTV(Dmax). CONCLUSION: Preliminary results of the dummy run procedure indicate that most centres conformed to protocol requirements. To quantify plan quality we recommend systematic calculation of U and either CI or CN, both of which measure the amount of irradiated normal brain tissue.

Determination of Very Slow Diffusion of Paramagnetic Ions by NMR Spectroscopy

In this paper, we propose a new method of measuring the very slow paramagnetic ion diffusion coefficient using a commercial high-resolution spectrometer. If there are distinct paramagnetic ions influencing the hydrogen nuclear magnetic relaxation time differently, their diffusion coefficients can be measured separately. A cylindrical phantom filled with Fricke xylenol gel solution and irradiated with gamma rays was used to validate the method. The Fricke xylenol gel solution was prepared with 270 Bloom porcine gelatin, the phantom was irradiated with gamma rays originated from a (60)Co source and a high-resolution 200 MHz nuclear magnetic resonance (NMR) spectrometer was used to obtain the phantom (1)H profile in the presence of a linear magnetic field gradient. By observing the temporal evolution of the phantom NMR profile, an apparent ferric ion diffusion coefficient of 0.50 mu m(2)/ms due to ferric ions diffusion was obtained. In any medical process where the ionizing radiation is used, the dose planning and the dose delivery are the key elements for the patient safety and success of treatment. These points become even more important in modern conformal radio therapy techniques, such as stereotactic radiosurgery, where the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Several methods have been proposed to obtain the three-dimensional (3-D) dose distribution. Recently, we proposed an alternative method for the 3-D radiation dose mapping, where the ionizing radiation modifies the local relative concentration of Fe(2+)/Fe(3+) in a phantom containing Fricke gel and this variation is associated to the MR image intensity. The smearing of the intensity gradient is proportional to the diffusion coefficient of the Fe(3+) and Fe(2+) in the phantom. There are several methods for measurement of the ionic diffusion using NMR, however, they are applicable when the diffusion is not very slow.

Um estudo algorítmico para otimização do plano de tratamento da radioterapia conformal

This work performs an algorithmic study of optimization of a conformal radiotherapy plan treatment. Initially we show: an overview about cancer, radiotherapy and the physics of interaction of ionizing radiation with matery. A proposal for optimization of a plan of treatment in radiotherapy is developed in a systematic way. We show the paradigm of multicriteria problem, the concept of Pareto optimum and Pareto dominance. A generic optimization model for radioterapic treatment is proposed. We construct the input of the model, estimate the dose given by the radiation using the dose matrix, and show the objective function for the model. The complexity of optimization models in radiotherapy treatment is typically NP which justifyis the use of heuristic methods. We propose three distinct methods: MOGA, MOSA e MOTS. The project of these three metaheuristic procedures is shown. For each procedures follows: a brief motivation, the algorithm itself and the method for tuning its parameters. The three method are applied to a concrete case and we confront their performances. Finally it is analyzed for each method: the quality of the Pareto sets, some solutions and the respective Pareto curves

Programação linear no planejamento do tratamento de câncer por radiocirurgia

Medical Physics has been reaching an important role among several lines in Science, providing means for the improvement of several theories and procedures. Currently, its main application is related with the use of ionizing radiations, specially, in treatment procedures such as Radiotherapy. Radiosurgery is a Radiotherapy technique which consists in administering a single tumoricidal dose of radiation exclusively to the tumorous lesion. It becomes then an interesting alternative to surgical treatment, mainly in cerebral metastases, which are the most frequent cerebral tumors in the central nervous system. The radio neurosurgical team works out a planning for the Radiosurgery treatment, aiming for obtaining an appropriate ideal treatment for each case. For the working out of this treatment planning, Computed Tomography images of the region to be treated are obtained, digitalized and later, fused with nuclear magnetic resonance images. Through these images, critical structures, organs at risk and lesions are localized. After this, calculations are made to determine three-dimensional positions of isocenters, isodose curves, prescribed dose, collimators sizes, position, numbers and respective weight of isocentric conformal fields, and others. The treatment planning is commonly based in desired levels of dose for specific types of tumors and organs at risk concerning the irradiated region. Theses levels of dose are chosen in a way that a high probability of cure may be achieved and meanwhile, that the probability of complications, in whichever organ at risk, may be minimal. Thus, many researches have been carried out, showing that mathematical techniques may help to obtain an optimal planning for the treatment of cerebral metastases. Among the methods of optimization in the study...(Complete abstract click electronic access below)

Tomotherapy dose distribution verification using MAGIC-f polymer gel dosimetry

Purpose: This paper presents the application of MAGIC-f gel in a three-dimensional dose distribution measurement and its ability to accurately measure the dose distribution from a tomotherapy unit. Methods: A prostate intensity-modulated radiation therapy (IMRT) irradiation was simulated in the gel phantom and the treatment was delivered by a TomoTherapy equipment. Dose distribution was evaluated by the R2 distribution measured in magnetic resonance imaging. Results: A high similarity was found by overlapping of isodoses of the dose distribution measured with the gel and expected by the treatment planning system (TPS). Another analysis was done by comparing the relative absorbed dose profiles in the measured and in the expected dose distributions extracted along indicated lines of the volume and the results were also in agreement. The gamma index analysis was also applied to the data and a high pass rate was achieved (88.4% for analysis using 3%/3 mm and of 96.5% using 4%/4 mm). The real three-dimensional analysis compared the dose-volume histograms measured for the planning volumes and expected by the treatment planning, being the results also in good agreement by the overlapping of the curves. Conclusions: These results show that MAGIC-f gel is a promise for tridimensional dose distribution measurements. (C) 2012 American Association of Physicists in Medicine. [http://dx.doi.org/10.1118/1.4704496]

Monte Carlo source model for photon beam radiotherapy: photon source characteristics

A major barrier to widespread clinical implementation of Monte Carlo dose calculation is the difficulty in characterizing the radiation source within a generalized source model. This work aims to develop a generalized three-component source model (target, primary collimator, flattening filter) for 6- and 18-MV photon beams that match full phase-space data (PSD). Subsource by subsource comparison of dose distributions, using either source PSD or the source model as input, allows accurate source characterization and has the potential to ease the commissioning procedure, since it is possible to obtain information about which subsource needs to be tuned. This source model is unique in that, compared to previous source models, it retains additional correlations among PS variables, which improves accuracy at nonstandard source-to-surface distances (SSDs). In our study, three-dimensional (3D) dose calculations were performed for SSDs ranging from 50 to 200 cm and for field sizes from 1 x 1 to 30 x 30 cm2 as well as a 10 x 10 cm2 field 5 cm off axis in each direction. The 3D dose distributions, using either full PSD or the source model as input, were compared in terms of dose-difference and distance-to-agreement. With this model, over 99% of the voxels agreed within +/-1% or 1 mm for the target, within 2% or 2 mm for the primary collimator, and within +/-2.5% or 2 mm for the flattening filter in all cases studied. For the dose distributions, 99% of the dose voxels agreed within 1% or 1 mm when the combined source model-including a charged particle source and the full PSD as input-was used. The accurate and general characterization of each photon source and knowledge of the subsource dose distributions should facilitate source model commissioning procedures by allowing scaling the histogram distributions representing the subsources to be tuned.

The exchange of radiotherapy data as part of an electronic patient-referral system

PURPOSE: To describe the implementation and use of an electronic patient-referral system as an aid to the efficient referral of patients to a remote and specialized treatment center. METHODS AND MATERIALS: A system for the exchange of radiotherapy data between different commercial planning systems and a specially developed planning system for proton therapy has been developed through the use of the PAPYRUS diagnostic image standard as an intermediate format. To ensure the cooperation of the different TPS manufacturers, the number of data sets defined for transfer has been restricted to the three core data sets of CT, VOIs, and three-dimensional dose distributions. As a complement to the exchange of data, network-wide application-sharing (video-conferencing) technologies have been adopted to provide methods for the interactive discussion and assessment of treatments plans with one or more partner clinics. RESULTS: Through the use of evaluation plans based on the exchanged data, referring clinics can accurately assess the advantages offered by proton therapy on a patient-by-patient basis, while the practicality or otherwise of the proposed treatments can simultaneously be assessed by the proton therapy center. Such a system, along with the interactive capabilities provided by video-conferencing methods, has been found to be an efficient solution to the problem of patient assessment and selection at a specialized treatment center, and is a necessary first step toward the full electronic integration of such centers with their remotely situated referral centers.

Two dimensional analysis of inflatable structures by the positional FEM

This paper presents a, simple two dimensional frame formulation to deal with structures undergoing large motions due to dynamic actions including very thin inflatable structures, balloons. The proposed methodology is based on the minimum potential energy theorem written regarding nodal positions. Velocity, acceleration and strain are achieved directly from positions, not. displacements, characterizing the novelty of the proposed technique. A non-dimensional space is created and the deformation function (change of configuration) is written following two independent mappings from which the strain energy function is written. The classical New-mark equations are used to integrate time. Dumping and non-conservative forces are introduced into the mechanical system by a rheonomic energy function. The final formulation has the advantage of being simple and easy to teach, when compared to classical Counterparts. The behavior of a bench-mark problem (spin-up maneuver) is solved to prove the formulation regarding high circumferential speed applications. Other examples are dedicated to inflatable and very thin structures, in order to test the formulation for further analysis of three dimensional balloons.

Defect-induced effects on carrier migration through one-dimensional poly(para-phenylenevinylene) chains

Defects in one-dimensional (1D) systems can be intrinsically distinct from its three-dimensional counterparts, and polymer films are good candidates for showing both extremes that are difficult to individuate in the experimental data. We study theoretically the impact of simple hydrogen and oxygen defects on the electron transport properties of one-dimensional poly(para-phenylenevinylene) chains through a multiscale technique, starting from classical structural simulations for crystalline films to extensive ab initio calculations within density functional theory for the defects in single crystalline-constrained chains. The most disruptive effect on carrier transport comes from conjugation breaking imposed by the overcoordination of a carbon atom in the vinyl group independently from the chemical nature of the defect. The particular case of the [C=O] (keto-defect) shows in addition unexpected electron-hole separation, suggesting that the experimentally detected photoluminescence bleaching and photoconductivity enhancement could be due to exciton dissociation caused by the 1D characteristics of the defect.

Character of magnetic excitations in a quasi-one-dimensional antiferromagnet near the quantum critical points: Impact on magnetoacoustic properties

We report results of magnetoacoustic studies in the quantum spin-chain magnet NiCl(2)-4SC(NH(2))(2) (DTN) having a field-induced ordered antiferromagnetic (AF) phase. In the vicinity of the quantum critical points (QCPs) the acoustic c(33) mode manifests a pronounced softening accompanied by energy dissipation of the sound wave. The acoustic anomalies are traced up to T > T(N), where the thermodynamic properties are determined by fermionic magnetic excitations, the ""hallmark"" of one-dimensional (1D) spin chains. On the other hand, as established in earlier studies, the AF phase in DTN is governed by bosonic magnetic excitations. Our results suggest the presence of a crossover from a 1D fermionic to a three-dimensional bosonic character of the magnetic excitations in DTN in the vicinity of the QCPs.

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

Finite element and dimensional analysis algorithm for the prediction of mechanical properties of bulk materials and thin films

In this work, the applicability of a new algorithm for the estimation of mechanical properties from instrumented indentation data was studied for thin films. The applicability was analyzed with the aid of both three-dimensional finite element simulations and experimental indentation tests. The numerical approach allowed studying the effect of the substrate on the estimation of mechanical properties of the film, which was conducted based on the ratio h(max)/l between maximum indentation depth and film thickness. For the experimental analysis, indentation tests were conducted on AISI H13 tool steel specimens, plasma nitrated and coated with TiN thin films. Results have indicated that, for the conditions analyzed in this work, the elastic deformation of the substrate limited the extraction of mechanical properties of the film/substrate system. This limitation occurred even at low h(max)/l ratios and especially for the estimation of the values of yield strength and strain hardening exponent. At indentation depths lower than 4% of the film thickness, the proposed algorithm estimated the mechanical properties of the film with accuracy. Particularly for hardness, precise values were estimated at h(max)/l lower than 0.1, i.e. 10% of film thickness. (C) 2010 Published by Elsevier B.V.