Monte Carlo modeling and simulations of the High Definition (HD120) micro MLC and validation against measurements for a 6 MV beam

Purpose: The most recent Varian® micro multileaf collimator(MLC), the High Definition (HD120) MLC, was modeled using the BEAMNRCMonte Carlo code. This model was incorporated into a Varian medical linear accelerator, for a 6 MV beam, in static and dynamic mode. The model was validated by comparing simulated profiles with measurements. Methods: The Varian® Trilogy® (2300C/D) accelerator model was accurately implemented using the state-of-the-art Monte Carlo simulation program BEAMNRC and validated against off-axis and depth dose profiles measured using ionization chambers, by adjusting the energy and the full width at half maximum (FWHM) of the initial electron beam. The HD120 MLC was modeled by developing a new BEAMNRC component module (CM), designated HDMLC, adapting the available DYNVMLC CM and incorporating the specific characteristics of this new micro MLC. The leaf dimensions were provided by the manufacturer. The geometry was visualized by tracing particles through the CM and recording their position when a leaf boundary is crossed. The leaf material density and abutting air gap between leaves were adjusted in order to obtain a good agreement between the simulated leakage profiles and EBT2 film measurements performed in a solid water phantom. To validate the HDMLC implementation, additional MLC static patterns were also simulated and compared to additional measurements. Furthermore, the ability to simulate dynamic MLC fields was implemented in the HDMLC CM. The simulation results of these fields were compared with EBT2 film measurements performed in a solid water phantom. Results: Overall, the discrepancies, with and without MLC, between the opened field simulations and the measurements using ionization chambers in a water phantom, for the off-axis profiles are below 2% and in depth-dose profiles are below 2% after the maximum dose depth and below 4% in the build-up region. On the conditions of these simulations, this tungsten-based MLC has a density of 18.7 g cm− 3 and an overall leakage of about 1.1 ± 0.03%. The discrepancies between the film measured and simulated closed and blocked fields are below 2% and 8%, respectively. Other measurements were performed for alternated leaf patterns and the agreement is satisfactory (to within 4%). The dynamic mode for this MLC was implemented and the discrepancies between film measurements and simulations are within 4%. Conclusions: The Varian® Trilogy® (2300 C/D) linear accelerator including the HD120 MLC was successfully modeled and simulated using the Monte CarloBEAMNRC code by developing an independent CM, the HDMLC CM, either in static and dynamic modes.

Dosimetric effect by shallow air cavities in high energy electron beams

This study evaluates the dosimetric impact caused by an air cavity located at 2 mm depth from the top surface in a PMMA phantom irradiated by electron beams produced by a Siemens Primus linear accelerator. A systematic evaluation of the effect related to the cavity area and thickness as well as to the electron beam energy was performed by using Monte Carlo simulations (EGSnrc code), Pencil Beam algorithm and Gafchromic EBT2 films. A home-PMMA phantom with the same geometry as the simulated one was specifically constructed for the measurements. Our results indicate that the presence of the cavity causes an increase (up to 70%) of the dose maximum value as well as a shift forward of the position of the depthedose curve, compared to the homogeneous one. Pronounced dose discontinuities in the regions close to the lateral cavity edges are observed. The shape and magnitude of these discontinuities change with the dimension of the cavity. It is also found that the cavity effect is more pronounced (6%) for the 12 MeV electron beam and the presence of cavities with large thickness and small area introduces more significant variations (up to 70%) on the depthedose curves. Overall, the Gafchromic EBT2 film measurements were found in agreement within 3% with Monte Carlo calculations and predict well the fine details of the dosimetric change near the cavity interface. The Pencil Beam calculations underestimate the dose up to 40% compared to Monte Carlo simulations; in particular for the largest cavity thickness (2.8 cm).

Validação do uso de cone-beam computed tomography de kilovoltagem de um LINAC para a realização de estudos dosimétricos

Mestrado em Radioterapia

Radioterapia adaptativa: utilização do cone-beam computed tomography para o cálculo de dose em radioterapia

Mestrado em Radiações Aplicadas às Tecnologias da Saúde - Ramo de especialização: Terapia com Radiações

Femtosecond laser ablation of dentin

The surface morphology, structure and composition of human dentin treated with a femtosecond infrared laser (pulse duration 500 fs, wavelength 1030 nm, fluences ranging from 1 to 3 J cm(-2)) was studied by scanning electron microscopy, x-ray diffraction, x-ray photoelectron spectroscopy and Fourier transform infrared spectroscopy. The average dentin ablation threshold under these conditions was 0.6 +/- 0.2 J cm(-2) and the ablation rate achieved in the range 1 to 2 mu m/pulse for an average fluence of 3 J cm(-2). The ablation surfaces present an irregular and rugged appearance, with no significant traces of melting, deformation, cracking or carbonization. The smear layer was entirely removed by the laser treatment. For fluences only slightly higher than the ablation threshold the morphology of the laser-treated surfaces was very similar to the dentin fracture surfaces and the dentinal tubules remained open. For higher fluences, the surface was more porous and the dentin structure was partially concealed by ablation debris and a few resolidified droplets. Independently on the laser processing parameters and laser processing method used no sub-superficial cracking was observed. The dentin constitution and chemical composition was not significantly modified by the laser treatment in the processing parameter range used. In particular, the organic matter is not preferentially removed from the surface and no traces of high temperature phosphates, such as the beta-tricalcium phosphate, were observed. The achieved results are compatible with an electrostatic ablation mechanism. In conclusion, the high beam quality and short pulse duration of the ultrafast laser used should allow the accurate preparation of cavities, with negligible damage of the underlying material.

Structural and Microanalytical Studies of CrO2 Thin Films on c-Sapphire by High Resolution Electron Microscopy Methods

Chromium dioxide (CrO2) has been extensively used in the magnetic recording industry. However, it is its ferromagnetic half-metallic nature that has more recently attracted much attention, primarily for the development of spintronic devices. CrO2 is the only stoichiometric binary oxide theoretically predicted to be fully spin polarized at the Fermi level. It presents a Curie temperature of ∼ 396 K, i.e. well above room temperature, and a magnetic moment of 2 mB per formula unit. However an antiferromagnetic native insulating layer of Cr2O3 is always present on the CrO2 surface which enhances the CrO2 magnetoresistance and might be used as a barrier in magnetic tunnel junctions.

Desvio de posicionamento em radioterapia para patologias de cabeça e pescoço e próstata: revisão de literatura

Introdução – Numa era em que os tratamentos de Radioterapia Externa (RTE) exigem cada vez mais precisão, a utilização de imagem médica permitirá medir, quantificar e avaliar o impacto do erro provocado pela execução do tratamento ou pelos movimentos dos órgãos. Objetivo – Analisar os dados existentes na literatura acerca de desvios de posicionamento (DP) em patologias de cabeça e pescoço (CP) e próstata, medidos com Cone Beam Computed Tomography (CBCT) ou Electronic Portal Image Device (EPID). Metodologia – Para esta revisão da literatura foram pesquisados artigos recorrendo às bases de dados MEDLINE/PubMed e b-on. Foram incluídos artigos que reportassem DP em patologias CP e próstata medidos através de CBCT e EPID. Seguidamente foram aplicados critérios de validação, que permitiram a seleção dos estudos. Resultados – Após a análise de 35 artigos foram incluídos 13 estudos e validados 9 estudos. Para tumores CP, a média (μ) dos DP encontra-se entre 0,0 e 1,2mm, com um desvio padrão (σ) máximo de 1,3mm. Para patologias de próstata observa-se μDP compreendido entre 0,0 e 7,1mm, com σ máximo de 7,5mm. Discussão/Conclusão – Os DP em patologias CP são atribuídos, maioritariamente, aos efeitos secundários da RTE, como mucosite e dor, que afetam a deglutição e conduzem ao emagrecimento, contribuindo para a instabilidade da posição do doente durante o tratamento, aumentando as incertezas de posicionamento. Os movimentos da próstata devem-se principalmente às variações de preenchimento vesical, retal e gás intestinal. O desconhecimento dos DP afeta negativamente a precisão da RTE. É importante detetá-los e quantificá-los para calcular margens adequadas e a magnitude dos erros, aumentando a precisão da administração de RTE, incluindo o aumento da segurança do doente. - ABSTRACT - Background and Purpose – In an era where precision is an increasing necessity in external radiotherapy (RT), modern medical imaging techniques provide means for measuring, quantifying and evaluating the impact of treatment execution and movement error. The aim of this paper is to review the current literature on the quantification of setup deviations (SD) in patients with head and neck (H&N) or prostate tumors, using Cone Beam Computed Tomography (CBCT) or Electronic Portal Image Device (EPID). Methods – According to the study protocol, MEDLINE/PubMed and b-on databases were searched for trials, which were analyzed using selection criteria based on the quality of the articles. Results – After assessment of 35 papers, 13 studies were included in this analysis and nine were authenticated (6 for prostate and 3 for H&N tumors). The SD in the treatment of H&N cancer patients is in the interval of 0.1 to 1.2mm, whereas in prostate cancer this interval is 0.0 to 7.1mm. Discussion – The reproducibility of patient positioning is the biggest barrier for higher precision in RT, which is affected by geometrical uncertainty, positioning errors and inter or intra-fraction organ movement. There are random and systematic errors associated to patient positioning, introduced since the treatment planning phase or through physiological organ movement. Conclusion – The H&N SD are mostly assigned to the Radiotherapy adverse effects, like mucositis and pain, which affect swallowing and decrease secretions, contributing for the instability of patient positioning during RT treatment and increasing positioning uncertainties. Prostate motion is mainly related to the variation in bladder and rectal filling. Ignoring SD affects negatively the accuracy of RT. Therefore, detection and quantification of SD is crucial in order to calculate appropriate margins, the magnitude of error and to improve accuracy in RTE and patient safety.

Positive solutions of fourth order problems with clamped beam boundary conditions

n this paper we make an exhaustive study of the fourth order linear operator u((4)) + M u coupled with the clamped beam conditions u(0) = u(1) = u'(0) = u'(1) = 0. We obtain the exact values on the real parameter M for which this operator satisfies an anti-maximum principle. Such a property is equivalent to the fact that the related Green's function is nonnegative in [0, 1] x [0, 1]. When M < 0 we obtain the best estimate by means of the spectral theory and for M > 0 we attain the optimal value by studying the oscillation properties of the solutions of the homogeneous equation u((4)) + M u = 0. By using the method of lower and upper solutions we deduce the existence of solutions for nonlinear problems coupled with this boundary conditions. (C) 2011 Elsevier Ltd. All rights reserved.

Otimização da imagem em tomografia computorizada no exame de tórax

Mestrado em Radiações Aplicadas às Tecnologias da Saúde. Área de especialização: Imagem Digital.

Avaliação da distribuição de dose no CTV-PTV através de imagem guiada em doentes com carcinoma da próstata

Mestrado em Radioterapia.

IGRT com CBCT: o impacto na precisão em radioterapia

Mestrado em Radiações Aplicadas às Tecnologias da Saúde.

Comparação do cálculo de dose com os algoritmos Analytical Anisotropic Algorithm e Pencil Beam Convolution na patologia de pulmão

Mestrado em Radioterapia.

Modelling a Rotating Shaft as an Elastically Restrained Bernoulli-Euler Beam

Industrial rotating machines may be exposed to severe dynamic excitations due to resonant working regimes. Dealing with the bending vibration, problem of a machine rotor, the shaft - and attached discs - can be simply modelled using the Bernoulli-Euler beam theory, as a continuous beam subjected to a specific set of boundary conditions. In this study, the authors recall Rayleigh's method to propose an iterative strategy, which allows for the determination of natural frequencies and mode shapes of continuous beams taking into account the effect of attached concentrated masses and rotational inertias, including different stiffness coefficients at the right and the left end sides. The algorithm starts with the exact solutions from Bernoulli-Euler's beam theory, which are then updated through Rayleigh's quotient parameters. Several loading cases are examined in comparison with the experimental data and examples are presented to illustrate the validity of the model and the accuracy of the obtained values.

Comparação entre o pencil beam convolution algorithm e o analytical anisotropic algorithm em tumores de mama

Mestrado em Radioterapia

Monte Carlo simulation of the basic features of the GE Millennium MG single photon emission computed tomography gamma camera

Objective - To describe and validate the simulation of the basic features of GE Millennium MG gamma camera using the GATE Monte Carlo platform. Material and methods - Crystal size and thickness, parallel-hole collimation and a realistic energy acquisition window were simulated in the GATE platform. GATE results were compared to experimental data in the following imaging conditions: a point source of 99mTc at different positions during static imaging and tomographic acquisitions using two different energy windows. The accuracy between the events expected and detected by simulation was obtained with the Mann–Whitney–Wilcoxon test. Comparisons were made regarding the measurement of sensitivity and spatial resolution, static and tomographic. Simulated and experimental spatial resolutions for tomographic data were compared with the Kruskal–Wallis test to assess simulation accuracy for this parameter. Results - There was good agreement between simulated and experimental data. The number of decays expected when compared with the number of decays registered, showed small deviation (≤0.007%). The sensitivity comparisons between static acquisitions for different distances from source to collimator (1, 5, 10, 20, 30cm) with energy windows of 126–154 keV and 130–158 keV showed differences of 4.4%, 5.5%, 4.2%, 5.5%, 4.5% and 5.4%, 6.3%, 6.3%, 5.8%, 5.3%, respectively. For the tomographic acquisitions, the mean differences were 7.5% and 9.8% for the energy window 126–154 keV and 130–158 keV. Comparison of simulated and experimental spatial resolutions for tomographic data showed no statistically significant differences with 95% confidence interval. Conclusions - Adequate simulation of the system basic features using GATE Monte Carlo simulation platform was achieved and validated.