Monte Carlo simulations of dynamic radiotherapy treatments

The effects of tumour motion during radiation therapy delivery have been widely investigated. Motion effects have become increasingly important with the introduction of dynamic radiotherapy delivery modalities such as enhanced dynamic wedges (EDWs) and intensity modulated radiation therapy (IMRT) where a dynamically collimated radiation beam is delivered to the moving target, resulting in dose blurring and interplay effects which are a consequence of the combined tumor and beam motion. Prior to this work, reported studies on the EDW based interplay effects have been restricted to the use of experimental methods for assessing single-field non-fractionated treatments. In this work, the interplay effects have been investigated for EDW treatments. Single and multiple field treatments have been studied using experimental and Monte Carlo (MC) methods. Initially this work experimentally studies interplay effects for single-field non-fractionated EDW treatments, using radiation dosimetry systems placed on a sinusoidaly moving platform. A number of wedge angles (60º, 45º and 15º), field sizes (20 × 20, 10 × 10 and 5 × 5 cm2), amplitudes (10-40 mm in step of 10 mm) and periods (2 s, 3 s, 4.5 s and 6 s) of tumor motion are analysed (using gamma analysis) for parallel and perpendicular motions (where the tumor and jaw motions are either parallel or perpendicular to each other). For parallel motion it was found that both the amplitude and period of tumor motion affect the interplay, this becomes more prominent where the collimator tumor speeds become identical. For perpendicular motion the amplitude of tumor motion is the dominant factor where as varying the period of tumor motion has no observable effect on the dose distribution. The wedge angle results suggest that the use of a large wedge angle generates greater dose variation for both parallel and perpendicular motions. The use of small field size with a large tumor motion results in the loss of wedged dose distribution for both parallel and perpendicular motion. From these single field measurements a motion amplitude and period have been identified which show the poorest agreement between the target motion and dynamic delivery and these are used as the „worst case motion parameters.. The experimental work is then extended to multiple-field fractionated treatments. Here a number of pre-existing, multiple–field, wedged lung plans are delivered to the radiation dosimetry systems, employing the worst case motion parameters. Moreover a four field EDW lung plan (using a 4D CT data set) is delivered to the IMRT quality control phantom with dummy tumor insert over four fractions using the worst case parameters i.e. 40 mm amplitude and 6 s period values. The analysis of the film doses using gamma analysis at 3%-3mm indicate the non averaging of the interplay effects for this particular study with a gamma pass rate of 49%. To enable Monte Carlo modelling of the problem, the DYNJAWS component module (CM) of the BEAMnrc user code is validated and automated. DYNJAWS has been recently introduced to model the dynamic wedges. DYNJAWS is therefore commissioned for 6 MV and 10 MV photon energies. It is shown that this CM can accurately model the EDWs for a number of wedge angles and field sizes. The dynamic and step and shoot modes of the CM are compared for their accuracy in modelling the EDW. It is shown that dynamic mode is more accurate. An automation of the DYNJAWS specific input file has been carried out. This file specifies the probability of selection of a subfield and the respective jaw coordinates. This automation simplifies the generation of the BEAMnrc input files for DYNJAWS. The DYNJAWS commissioned model is then used to study multiple field EDW treatments using MC methods. The 4D CT data of an IMRT phantom with the dummy tumor is used to produce a set of Monte Carlo simulation phantoms, onto which the delivery of single field and multiple field EDW treatments is simulated. A number of static and motion multiple field EDW plans have been simulated. The comparison of dose volume histograms (DVHs) and gamma volume histograms (GVHs) for four field EDW treatments (where the collimator and patient motion is in the same direction) using small (15º) and large wedge angles (60º) indicates a greater mismatch between the static and motion cases for the large wedge angle. Finally, to use gel dosimetry as a validation tool, a new technique called the „zero-scan method. is developed for reading the gel dosimeters with x-ray computed tomography (CT). It has been shown that multiple scans of a gel dosimeter (in this case 360 scans) can be used to reconstruct a zero scan image. This zero scan image has a similar precision to an image obtained by averaging the CT images, without the additional dose delivered by the CT scans. In this investigation the interplay effects have been studied for single and multiple field fractionated EDW treatments using experimental and Monte Carlo methods. For using the Monte Carlo methods the DYNJAWS component module of the BEAMnrc code has been validated and automated and further used to study the interplay for multiple field EDW treatments. Zero-scan method, a new gel dosimetry readout technique has been developed for reading the gel images using x-ray CT without losing the precision and accuracy.

Analysis of Cone-Beam CT using prior information

Treatment plans for conformal radiotherapy are based on an initial CT scan. The aim is to deliver the prescribed dose to the tumour, while minimising exposure to nearby organs. Recent advances make it possible to also obtain a Cone-Beam CT (CBCT) scan, once the patient has been positioned for treatment. A statistical model will be developed to compare these CBCT scans with the initial CT scan. Changes in the size, shape and position of the tumour and organs will be detected and quantified. Some progress has already been made in segmentation of prostate CBCT scans [1],[2],[3]. However, none of the existing approaches have taken full advantage of the prior information that is available. The planning CT scan is expertly annotated with contours of the tumour and nearby sensitive objects. This data is specific to the individual patient and can be viewed as a snapshot of spatial information at a point in time. There is an abundance of studies in the radiotherapy literature that describe the amount of variation in the relevant organs between treatments. The findings from these studies can form a basis for estimating the degree of uncertainty. All of this information can be incorporated as an informative prior into a Bayesian statistical model. This model will be developed using scans of CT phantoms, which are objects with known geometry. Thus, the accuracy of the model can be evaluated objectively. This will also enable comparison between alternative models.

Structure and deformation correlation of closed-cell aluminium foam subject to uniaxial compression

We report the results of an experimental and numerical study conducted on a closed-cell aluminium foam that was subjected to uniaxial compression with lateral constraint. X-ray computed tomography was utilized to gain access into the three-dimensional (3-D) structure of the foam and some aspects of the deformation mechanisms. A series of advanced 3-D image analyses are conducted on the 3-D images aimed at characterizing the strain localization regions. We identify the morphological/geometrical features that are responsible for the collapse of the cells and the strain localization. A novel mathematical approach based on a Minkowski tensor analysis along with the mean intercept length technique were utilized to search for signatures of anisotropy across the foam sample and its evolution as a function of loading. Our results show that regions with higher degrees of anisotropy in the undeformed foam have a tendency to initiate the onset of cell collapse. Furthermore, we show that strain hardening occurs predominantly in regions with large cells and high anisotropy. We combine the finite element method with the tomographic images to simulate the mechanical response of the foam. We predict further deformation in regions where the foam is already deformed. Crown Copyright (C) 2012 Published by Elsevier Ltd. on behalf of Acta Materialia Inc. All rights reserved.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

A Plasmonic Gold Nanostar Theranostic Probe for In Vivo Tumor Imaging and Photothermal Therapy.

Nanomedicine has attracted increasing attention in recent years, because it offers great promise to provide personalized diagnostics and therapy with improved treatment efficacy and specificity. In this study, we developed a gold nanostar (GNS) probe for multi-modality theranostics including surface-enhanced Raman scattering (SERS) detection, x-ray computed tomography (CT), two-photon luminescence (TPL) imaging, and photothermal therapy (PTT). We performed radiolabeling, as well as CT and optical imaging, to investigate the GNS probe's biodistribution and intratumoral uptake at both macroscopic and microscopic scales. We also characterized the performance of the GNS nanoprobe for in vitro photothermal heating and in vivo photothermal ablation of primary sarcomas in mice. The results showed that 30-nm GNS have higher tumor uptake, as well as deeper penetration into tumor interstitial space compared to 60-nm GNS. In addition, we found that a higher injection dose of GNS can increase the percentage of tumor uptake. We also demonstrated the GNS probe's superior photothermal conversion efficiency with a highly concentrated heating effect due to a tip-enhanced plasmonic effect. In vivo photothermal therapy with a near-infrared (NIR) laser under the maximum permissible exposure (MPE) led to ablation of aggressive tumors containing GNS, but had no effect in the absence of GNS. This multifunctional GNS probe has the potential to be used for in vivo biosensing, preoperative CT imaging, intraoperative detection with optical methods (SERS and TPL), as well as image-guided photothermal therapy.

Spectrotemporal CT data acquisition and reconstruction at low dose.

PURPOSE: X-ray computed tomography (CT) is widely used, both clinically and preclinically, for fast, high-resolution anatomic imaging; however, compelling opportunities exist to expand its use in functional imaging applications. For instance, spectral information combined with nanoparticle contrast agents enables quantification of tissue perfusion levels, while temporal information details cardiac and respiratory dynamics. The authors propose and demonstrate a projection acquisition and reconstruction strategy for 5D CT (3D+dual energy+time) which recovers spectral and temporal information without substantially increasing radiation dose or sampling time relative to anatomic imaging protocols. METHODS: The authors approach the 5D reconstruction problem within the framework of low-rank and sparse matrix decomposition. Unlike previous work on rank-sparsity constrained CT reconstruction, the authors establish an explicit rank-sparse signal model to describe the spectral and temporal dimensions. The spectral dimension is represented as a well-sampled time and energy averaged image plus regularly undersampled principal components describing the spectral contrast. The temporal dimension is represented as the same time and energy averaged reconstruction plus contiguous, spatially sparse, and irregularly sampled temporal contrast images. Using a nonlinear, image domain filtration approach, the authors refer to as rank-sparse kernel regression, the authors transfer image structure from the well-sampled time and energy averaged reconstruction to the spectral and temporal contrast images. This regularization strategy strictly constrains the reconstruction problem while approximately separating the temporal and spectral dimensions. Separability results in a highly compressed representation for the 5D data in which projections are shared between the temporal and spectral reconstruction subproblems, enabling substantial undersampling. The authors solved the 5D reconstruction problem using the split Bregman method and GPU-based implementations of backprojection, reprojection, and kernel regression. Using a preclinical mouse model, the authors apply the proposed algorithm to study myocardial injury following radiation treatment of breast cancer. RESULTS: Quantitative 5D simulations are performed using the MOBY mouse phantom. Twenty data sets (ten cardiac phases, two energies) are reconstructed with 88 μm, isotropic voxels from 450 total projections acquired over a single 360° rotation. In vivo 5D myocardial injury data sets acquired in two mice injected with gold and iodine nanoparticles are also reconstructed with 20 data sets per mouse using the same acquisition parameters (dose: ∼60 mGy). For both the simulations and the in vivo data, the reconstruction quality is sufficient to perform material decomposition into gold and iodine maps to localize the extent of myocardial injury (gold accumulation) and to measure cardiac functional metrics (vascular iodine). Their 5D CT imaging protocol represents a 95% reduction in radiation dose per cardiac phase and energy and a 40-fold decrease in projection sampling time relative to their standard imaging protocol. CONCLUSIONS: Their 5D CT data acquisition and reconstruction protocol efficiently exploits the rank-sparse nature of spectral and temporal CT data to provide high-fidelity reconstruction results without increased radiation dose or sampling time.

Training improves interobserver reliability for the diagnosis of scaphoid fracture displacement

Background The diagnosis of displacement in scaphoid fractures is notorious for poor interobserver reliability.

Questions/purposes We tested whether training can improve interobserver reliability and sensitivity, specificity, and accuracy for the diagnosis of scaphoid fracture displacement on radiographs and CT scans.

Methods Sixty-four orthopaedic surgeons rated a set of radiographs and CT scans of 10 displaced and 10 nondisplaced scaphoid fractures for the presence of displacement, using a web-based rating application. Before rating, observers were randomized to a training group (34 observers) and a nontraining group (30 observers). The training group received an online training module before the rating session, and the nontraining group did not. Interobserver reliability for training and nontraining was assessed by Siegel’s multirater kappa and the Z-test was used to test for significance.

Results There was a small, but significant difference in the interobserver reliability for displacement ratings in favor of the training group compared with the nontraining group. Ratings of radiographs and CT scans combined resulted in moderate agreement for both groups. The average sensitivity, specificity, and accuracy of diagnosing displacement of scaphoid fractures were, respectively, 83%, 85%, and 84% for the nontraining group and 87%, 86%, and 87% for the training group. Assuming a 5% prevalence of fracture displacement, the positive predictive value was 0.23 in the nontraining group and 0.25 in the training group. The negative predictive value was 0.99 in both groups.

Conclusions Our results suggest training can improve interobserver reliability and sensitivity, specificity and accuracy for the diagnosis of scaphoid fracture displacement, but the improvements are slight. These findings are encouraging for future research regarding interobserver variation and how to reduce it further.

Interobserver reliability of radial head fracture classification : two-dimensional compared with three-dimensional CT

Background: The Broberg and Morrey modification of the Mason classification of radial head fractures has substantial interobserver variation. This study used a large web-based collaborative of experienced orthopaedic surgeons to test the hypothesis that three-dimensional reconstructions of computed tomography (CT) scans improve the interobserver reliability of the classification of radial head fractures according to the Broberg and Morrey modification of the Mason classification.

Methods: Eighty-five orthopaedic surgeons evaluated twelve radial head fractures. They were randomly assigned to review either radiographs and two-dimensional CT scans or radiographs and three-dimensional CT images to determine the fracture classification, fracture characteristics, and treatment recommendations. The kappa multirater measure (κ) was calculated to estimate agreement between observers.

Results: Three-dimensional CT had moderate agreement and two-dimensional CT had fair agreement among observers for the Broberg and Morrey modification of the Mason classification, a difference that was significant. Observers assessed seven fracture characteristics, including fracture line, comminution, articular surface involvement, articular step or gap of ≥2 mm, central impaction, recognition of more than three fracture fragments, and fracture fragments too small to repair. There was a significant difference in kappa values between three-dimensional CT and two-dimensional CT for fracture fragments too small to repair, recognition of three fracture fragments, and central impaction. The difference between the other four fracture characteristics was not significant. Among treatment recommendations, there was fair agreement for both three-dimensional CT and two-dimensional CT.

Conclusions: Although three-dimensional CT led to some small but significant decreases in interobserver variation, there is still considerable disagreement regarding classification and characterization of radial head fractures. Three-dimensional CT may be insufficient to optimize interobserver agreement.

Comportamento da sutura palatina mediana em crianças submetidas à expansão rápida da maxila: avaliação mediante imagem de tomografia computadorizada

OBJETIVO: o presente estudo, de caráter prospectivo, objetivou avaliar a imagem da sutura palatina mediana em crianças submetidas à expansão rápida da maxila por meio de tomografia computadorizada, após a fase de contenção. METODOLOGIA: a amostra constou de 17 crianças de ambos os gêneros, na faixa etária compreendida entre 5 anos e 2 meses e 10 anos e 5 meses no início do tratamento, provenientes da Clínica de Ortodontia Interceptiva da Sociedade de Promoção Social do Fissurado Lábio-Palatal da Universidade de São Paulo (PROFIS-USP), Bauru/SP. Tomografias computadorizadas foram implementadas para avaliar o comportamento da sutura palatina mediana em diferentes momentos do tratamento. RESULTADOS E CONCLUSÕES: constatou-se que, após um período médio de 8 a 9 meses de contenção com o aparelho expansor, a sutura palatina mediana mostrou-se completamente ossificada, desde a região da espinha nasal anterior até a espinha nasal posterior. Tal informação esclarece o comportamento da sutura palatina mediana frente à expansão rápida da maxila e reitera o senso comum quanto ao caráter biológico do procedimento.

Monte Carlo and analytical calculation of lateral deflection of proton beams in homogeneous targets

Proton radiation therapy is a precise form of radiation therapy, but the avoidance of damage to critical normal tissues and the prevention of geographical tumor misses require accurate knowledge of the dose delivered to the patient and the verification of his position demand a precise imaging technique. In proton therapy facilities, the X-ray Computed Tomography (xCT) is the preferred technique for the planning treatment of patients. This situation has been changing nowadays with the development of proton accelerators for health care and the increase in the number of treated patients. In fact, protons could be more efficient than xCT for this task. One essential difficulty in pCT image reconstruction systems came from the scattering of the protons inside the target due to the numerous small-angle deflections by nuclear Coulomb fields. The purpose of this study is the comparison of an analytical formulation for the determination of beam lateral deflection, based on Molière's theory and Rutherford scattering with Monte Carlo calculations by SRIM 2008 and MCNPX codes. © 2010 American Institute of Physics.

Análise tomográfica Cone Beam e histológica do preparo do canal radicular com três sistemas rotatórios

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Monitoramento do processo de densificação de compósitos termoestruturais de carbono reforçado com fibras de carbono por tomografia computadorizada e resistividade elétrica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Application of two segmentation protocols during the processing of virtual images in rapid prototyping: ex vivo study with human dry mandibles

The aim of this study was to evaluate the accuracy of virtual three-dimensional (3D) reconstructions of human dry mandibles, produced from two segmentation protocols (outline only and all-boundary lines).Twenty virtual three-dimensional (3D) images were built from computed tomography exam (CT) of 10 dry mandibles, in which linear measurements between anatomical landmarks were obtained and compared to an error probability of 5 %.The results showed no statistically significant difference among the dry mandibles and the virtual 3D reconstructions produced from segmentation protocols tested (p = 0,24).During the designing of a virtual 3D reconstruction, both outline only and all-boundary lines segmentation protocols can be used.Virtual processing of CT images is the most complex stage during the manufacture of the biomodel. Establishing a better protocol during this phase allows the construction of a biomodel with characteristics that are closer to the original anatomical structures. This is essential to ensure a correct preoperative planning and a suitable treatment.

Monitoring the process of densification of 3D and 4D (tridirectional and tetradirectional) carbon/carbon composites

Carbon fiber reinforced carbon composites can be made by iterative liquid impregnation or gas phase carbon deposition routes. In both cases, at the final processing stage the carbon fiber is embedded in carbon matrix which results in unique properties such as low density, high thermal conductivity and thermal shock resistance, low thermal expansion and high modulus, in relation to other refractory materials. In the present study assembled three-directional and four-directional preforms, having 50% volume of pores, were densified by iterative cycles of thermoset resin impregnation followed by pyrolysis under inert atmosphere, until appropriate densities were achieved. The thermoset resin is converted in a carbon matrix during pyrolysis. The iterative manufacturing process of the carbon fiber reinforced carbon composites is evaluated by means of nondestructive techniques based on X-ray computed tomography and electrical resistivity. X-ray computed tomography gives a general mapping view of the filling pores of the preforms which impacts results of the electrical resistivity. After six processing cycles and heat treatments up to 2000?, the final densities of the three-directional and four-directional carbon fiber reinforced carbon composites were 1.16g/cm(3) and an electrical resistivity of approximate to 0.07m. The configuration of preforms, three-directional or four-directional, did not alter the densification profile, in terms of increasing density and reducing porosity during the processing cycles.

Melatonin decreases breast cancer metastasis by modulating Rho-associated kinase protein-1 expression

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)