Monte Carlo-based diode design for correction-less small field dosimetry

Due to their small collecting volume diodes are commonly used in small field dosimetry. However the relative sensitivity of a diode increases with decreasing small field size. Conversely, small air gaps have been shown to cause a significant decrease in the sensitivity of a detector as the field size is decreased. Therefore this study uses Monte Carlo simulations to look at introducing air upstream to diodes such that they measure with a constant sensitivity across all field sizes in small field dosimetry. Varying thicknesses of air were introduced onto the upstream end of two commercial diodes (PTW 60016 photon diode and PTW 60017 electron diode), as well as a theoretical unenclosed silicon chip using field sizes as small as 5 mm × 5 mm . The metric D_(w,Q)/D_(Det,Q) used in this study represents the ratio of the dose to a point of water to the dose to the diode active volume, for a particular field size and location. The optimal thickness of air required to provide a constant sensitivity across all small field sizes was found by plotting D_(w,Q)/D_(Det,Q) as a function of introduced air gap size for various field sizes, and finding the intersection point of these plots. That is, the point at which D_(w,Q)/D_(Det,Q) was constant for all field sizes was found. The optimal thickness of air was calculated to be 3.3 mm, 1.15 mm and 0.10 mm for the photon diode, electron diode and unenclosed silicon chip respectively. The variation in these results was due to the different design of each detector. When calculated with the new diode design incorporating the upstream air gap, k_(Q_clin 〖,Q〗_msr)^(f_clin 〖,f〗_msr ) was equal to unity to within statistical uncertainty (0.5 %) for all three diodes. Cross-axis profile measurements were also improved with the new detector design. The upstream air gap could be implanted on the commercial diodes via a cap consisting of the air cavity surrounded by water equivalent material. The results for the unclosed silicon chip show that an ideal small field dosimetry diode could be created by using a silicon chip with a small amount of air above it.

Monte Carlo dosimetry for stereotactic radiosurgery

Stereotactic radiosurgery (SRS) treatments for brain cancers require small and precisely shaped photon beams. These beams can be generated by fitting a linear accelerator with a micro-multileaf collimator (mMLC) such as the BrainLAB m3, which offers greater flexibility for field shaping than standard SRS cone collimators

Bayesian model comparison in cosmology with Population Monte Carlo

We use Bayesian model selection techniques to test extensions of the standard flat LambdaCDM paradigm. Dark-energy and curvature scenarios, and primordial perturbation models are considered. To that end, we calculate the Bayesian evidence in favour of each model using Population Monte Carlo (PMC), a new adaptive sampling technique which was recently applied in a cosmological context. The Bayesian evidence is immediately available from the PMC sample used for parameter estimation without further computational effort, and it comes with an associated error evaluation. Besides, it provides an unbiased estimator of the evidence after any fixed number of iterations and it is naturally parallelizable, in contrast with MCMC and nested sampling methods. By comparison with analytical predictions for simulated data, we show that our results obtained with PMC are reliable and robust. The variability in the evidence evaluation and the stability for various cases are estimated both from simulations and from data. For the cases we consider, the log-evidence is calculated with a precision of better than 0.08. Using a combined set of recent CMB, SNIa and BAO data, we find inconclusive evidence between flat LambdaCDM and simple dark-energy models. A curved Universe is moderately to strongly disfavoured with respect to a flat cosmology. Using physically well-motivated priors within the slow-roll approximation of inflation, we find a weak preference for a running spectral index. A Harrison-Zel'dovich spectrum is weakly disfavoured. With the current data, tensor modes are not detected; the large prior volume on the tensor-to-scalar ratio r results in moderate evidence in favour of r=0.

Monte Carlo filtering and smoothing with application to time-varying spectral estimation

We develop methods for performing filtering and smoothing in non-linear non-Gaussian dynamical models. The methods rely on a particle cloud representation of the filtering distribution which evolves through time using importance sampling and resampling ideas. In particular, novel techniques are presented for generation of random realisations from the joint smoothing distribution and for MAP estimation of the state sequence. Realisations of the smoothing distribution are generated in a forward-backward procedure, while the MAP estimation procedure can be performed in a single forward pass of the Viterbi algorithm applied to a discretised version of the state space. An application to spectral estimation for time-varying autoregressions is described.

An overview of sequential Monte Carlo methods for parameter estimation in general state-space models

Nonlinear non-Gaussian state-space models arise in numerous applications in control and signal processing. Sequential Monte Carlo (SMC) methods, also known as Particle Filters, are numerical techniques based on Importance Sampling for solving the optimal state estimation problem. The task of calibrating the state-space model is an important problem frequently faced by practitioners and the observed data may be used to estimate the parameters of the model. The aim of this paper is to present a comprehensive overview of SMC methods that have been proposed for this task accompanied with a discussion of their advantages and limitations.

The application of a Trous wave filtering and Monte Carlo analysis on secis 2001 solar eclipse observations

Eight thousand images of the solar corona were captured during the June 2001 total solar eclipse. New software for the alignment of the images and an automated technique for detecting intensity oscillations using multi-scale wavelet analysis were developed. Large areas of the images covered by the Moon and the upper corona were scanned for oscillations and the statistical properties of the atmospheric effects were determined. The a Trous wavelet transform was used for noise reduction and Monte Carlo analysis as a significance test of the detections. The effectiveness of those techniques is discussed in detail.

Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Body shadowing mitigation using differentiated LOS/NLOS channel models for RSSI-based Monte Carlo personnel localization

Research into localization has produced a wealth of algorithms and techniques to estimate the location of wireless network nodes, however the majority of these schemes do not explicitly account for non-line of sight conditions. Disregarding this common situation reduces their accuracy and their potential for exploitation in real world applications. This is a particular problem for personnel tracking where the user's body itself will inherently cause time-varying blocking according to their movements. Using empirical data, this paper demonstrates that, by accounting for non-line of sight conditions and using received signal strength based Monte Carlo localization, meter scale accuracy can be achieved for a wrist-worn personnel tracking tag in a 120 m indoor office environment. © 2012 IEEE.

Root Cause Analysis by a Combined Sparse Classification and Monte Carlo Approach

Classification methods with embedded feature selection capability are very appealing for the analysis of complex processes since they allow the analysis of root causes even when the number of input variables is high. In this work, we investigate the performance of three techniques for classification within a Monte Carlo strategy with the aim of root cause analysis. We consider the naive bayes classifier and the logistic regression model with two different implementations for controlling model complexity, namely, a LASSO-like implementation with a L1 norm regularization and a fully Bayesian implementation of the logistic model, the so called relevance vector machine. Several challenges can arise when estimating such models mainly linked to the characteristics of the data: a large number of input variables, high correlation among subsets of variables, the situation where the number of variables is higher than the number of available data points and the case of unbalanced datasets. Using an ecological and a semiconductor manufacturing dataset, we show advantages and drawbacks of each method, highlighting the superior performance in term of classification accuracy for the relevance vector machine with respect to the other classifiers. Moreover, we show how the combination of the proposed techniques and the Monte Carlo approach can be used to get more robust insights into the problem under analysis when faced with challenging modelling conditions.

Otimização de esquemas terapêuticos do carcinoma da próstata envolvendo braquiterapia de baixa taxa de dose e radioterapia externa : abordagens física e radiobiológica

RESUMO: Este trabalho teve como objetivo a determinação de esquemas de tratamento alternativos para o carcinoma da próstata com radioterapia externa (EBRT) e braquiterapia de baixa taxa de dose (LDRBT) com implantes permanentes de Iodo-125, biologicamente equivalentes aos convencionalmente usados na prática clínica, com recurso a modelos teóricos e a métodos de Monte Carlo (MC). Os conceitos de dose biológica efetiva (BED) e de dose uniforme equivalente (EUD) foram utilizados, com o modelo linear-quadrático (LQ), para a determinação de regimes de tratamento equivalentes. Numa primeira abordagem, utilizou-se a BED para determinar: 1) esquemas hipofracionados de EBRT mantendo as complicações retais tardias de regimes convencionais com doses totais de 75,6 Gy, 77,4 Gy, 79,2 Gy e 81,0 Gy; e 2) a relação entre as doses totais de EBRT e LDRBT de modo a manter a BED do regime convencional de 45 Gy de EBRT e 110 Gy de LDRBT. Numa segunda abordagem, recorreu-se ao código de MC MCNPX para a simulação de distribuições de dose de EBRT e LDRBT em dois fantomas de voxel segmentados a partir das imagens de tomografia computorizada de pacientes com carcinoma da próstata. Os resultados das simulações de EBRT e LDRBT foram somados e determinada uma EUD total de forma a obterem-se: 1) esquemas equivalentes ao tratamento convencional de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT; e 2) esquemas equivalentes a EUD na próstata de 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy e 110 Gy. Em todos os resultados nota-se um ganho terapêutico teórico na utilização de esquemas hipofracionados de EBRT. Para uma BED no reto equivalente ao esquema convencional, tem-se um aumento de 2% na BED da próstata com menos 5 frações. Este incremento dá-se de forma cada vez mais visível à medida que se reduz o número de frações, sendo da ordem dos 10-11% com menos 20 frações e dos 35-45% com menos 40 frações. Considerando os resultados das simulações de EBRT, obteve-se uma EUD média de 107 Gy para a próstata e de 42 Gy para o reto, com o esquema convencional de 110 Gy de LDRBT, seguidos de 25 frações de 1,8 Gy de EBRT. Em termos de probabilidade de controlo tumoral (igual EUD), é equivalente a este tratamento a administração de EBRT em 66 frações de 1,8 Gy, 56 de 2 Gy, 40 de 2,5 Gy, 31 de 3 Gy, 20 de 4 Gy ou 13 de 5 Gy. Relativamente à administração de 66 frações de 1,8 Gy, a EUD generalizada no reto reduz em 6% com o recurso a frações de 2,5 Gy e em 10% com frações de 4 Gy. Determinou-se uma BED total de 162 Gy para a administração de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT. Variando-se a dose total de LDRBT (TDLDRBT) em função da dose total de EBRT (TDEBRT), de modo a garantir uma BED de 162 Gy, obteve-se a seguinte relação:.......... Os resultados das simulações mostram que a EUD no reto diminui com o aumento da dose total de LDRBT para dose por fração de EBRT (dEBRT) inferiores a 2, Gy e aumenta para dEBRT a partir dos 3 Gy. Para quantidades de TDLDRBT mais baixas (<50 Gy), o reto beneficia de frações maiores de EBRT. À medida que se aumenta a TDLDRBT, a EUD generalizada no reto torna-se menos dependente da dEBRT. Este trabalho mostra que é possível a utilização de diferentes regimes de tratamento para o carcinoma da próstata com radioterapia que possibilitem um ganho terapêutico, quer seja administrando uma maior dose biológica com efeitos tardios constantes, quer mantendo a dose no tumor e diminuindo a toxicidade retal. A utilização com precaução de esquemas hipofracionados de EBRT, para além do benefício terapêutico, pode trazer vantagens ao nível da conveniência para o paciente e economia de custos. Os resultados das simulações deste estudo e conversão para doses de efeito biológico para o tratamento do carcinoma da próstata apresentam linhas de orientação teórica de interesse para novos ensaios clínicos. --------------------------------------------------ABSTRACT: The purpose of this work was to determine alternative radiotherapy regimens for the treatment of prostate cancer using external beam radiotherapy (EBRT) and low dose-rate brachytherapy (LDRBT) with Iodine-125 permanent implants which are biologically equivalent to conventional clinical treatments, by the use of theoretical models and Monte Carlo techniques. The concepts of biological effective dose (BED) and equivalent uniform dose (EUD), together with the linear-quadratic model (LQ), were used for determining equivalent treatment regimens. In a first approach, the BED concept was used to determine: 1) hypofractionated schemes of EBRT maintaining late rectal complications as with the conventional regimens with total doses of 75.6 Gy, 77.4 Gy, 79.2 Gy and 81.0 Gy; and 2) the relationship between total doses of EBRT and LDRBT in order to keep the BED of the conventional treatment of 45 Gy of EBRT and 110 Gy of LDRBT. In a second approach, the MC code MCNPX was used for simulating dose distributions of EBRT and LDRBT in two voxel phantoms segmented from the computed tomography of patients with prostate cancer. The results of the simulations of EBRT and LDRBT were added up and given an overall EUD in order to obtain: 1) equivalent to conventional treatment regimens of 25 fraction of 1.8 Gy of EBRT in combination with 110Gy of LDRBT; and 2) equivalent schemes of EUD of 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy, and 110Gy to the prostate. In all the results it is noted a therapeutic gain using hypofractionated EBRT schemes. For a rectal BED equivalent to the conventional regimen, an increment of 2% in the prostate BED was achieved with less 5 fractions. This increase is visibly higher as the number of fractions decrease, amounting 10-11% with less 20 fractions and 35-45% with less 20 fractions. Considering the results of the EBRT simulations an average EUD of 107 Gy was achieved for the prostate and of 42 Gy for the rectum with the conventional scheme of 110 Gy of LDRBT followed by 25 fractions of 1.8 Gy of EBRT. In terms of tumor control probability (same EUD) it is equivalent to this treatment, for example, delivering the EBRT in 66 fractions of 1.8 Gy, 56 fractions of 2 Gy, 40 fractions of 2.5 Gy, 31 fractions of 3 Gy, 20 fractions of 4 Gy or 13 fractions of 5 Gy. Regarding the use of 66 fractions of 1.8 Gy, the rectum EUD is reduced to 6% with 2.5 Gy per fraction and to 10% with 4 Gy. A total BED of 162 Gy was achieved for the delivery of 25 fractions of 1.8 Gy of EBRT in combination with 110 Gy of LDRBT. By varying the total dose of LDRBT (TDLDRBT) with the total dose of EBRT (TDEBRT) so as to ensure a BED of 162 Gy, the following relationship was obtained: ....... The simulation results show that the rectum EUD decreases with the increase of the TDLDRBT, for EBRT dose per fracion (dEBRT) less than 2.5 Gy and increases for dEBRT above 3 Gy. For lower amounts of TDLDRBT (< 50Gy), the rectum benefits of larger EBRT fractions. As the TDLDRBT increases, the rectum gEUD becomes less dependent on the dEBRT. The use of different regimens which enable a therapeutic gain, whether deivering a higher dose with the same late biological effects or maintaining the dose to the tumor and reducing rectal toxicity is possible. The use with precaution of hypofractionated regimens, in addition to the therapeutic benefit, can bring advantages in terms of convenience for the patient and cost savings. The simulation results of this study together with the biological dose conversion for the treatment of prostate cancer serve as guidelines of interest for new clinical trials.

Quantum Monte Carlo study of electrostatic polarizabilities of H and He atoms /

The infinitesimal differential quantum Monte Carlo (QMC) technique is used to estimate electrostatic polarizabilities of the H and He atoms up to the sixth order in the electric field perturbation. All 542 different QMC estimators of the nonzero atomic polarizabilities are derived and used in order to decrease the statistical error and to obtain the maximum efficiency of the simulations. We are confident that the estimates are "exact" (free of systematic error): the two atoms are nodeless systems, hence no fixed-node error is introduced. Furthermore, we develope and use techniques which eliminate systematic error inherent when extrapolating our results to zero time-step and large stack-size. The QMC results are consistent with published accurate values obtained using perturbation methods. The precision is found to be related to the number of perturbations, varying from 2 to 4 significant digits.

Spectrum simulation of rough and nanostructured targets from their 2D and 3D image by Monte Carlo methods

Corteo is a program that implements Monte Carlo (MC) method to simulate ion beam analysis (IBA) spectra of several techniques by following the ions trajectory until a sufficiently large fraction of them reach the detector to generate a spectrum. Hence, it fully accounts for effects such as multiple scattering (MS). Here, a version of Corteo is presented where the target can be a 2D or 3D image. This image can be derived from micrographs where the different compounds are identified, therefore bringing extra information into the solution of an IBA spectrum, and potentially significantly constraining the solution. The image intrinsically includes many details such as the actual surface or interfacial roughness, or actual nanostructures shape and distribution. This can for example lead to the unambiguous identification of structures stoichiometry in a layer, or at least to better constraints on their composition. Because MC computes in details the trajectory of the ions, it simulates accurately many of its aspects such as ions coming back into the target after leaving it (re-entry), as well as going through a variety of nanostructures shapes and orientations. We show how, for example, as the ions angle of incidence becomes shallower than the inclination distribution of a rough surface, this process tends to make the effective roughness smaller in a comparable 1D simulation (i.e. narrower thickness distribution in a comparable slab simulation). Also, in ordered nanostructures, target re-entry can lead to replications of a peak in a spectrum. In addition, bitmap description of the target can be used to simulate depth profiles such as those resulting from ion implantation, diffusion, and intermixing. Other improvements to Corteo include the possibility to interpolate the cross-section in angle-energy tables, and the generation of energy-depth maps.

Chronological Monte Carlo-based assessment of distribution system reliability

Regulatory authorities in many countries, in order to maintain an acceptable balance between appropriate customer service qualities and costs, are introducing a performance-based regulation. These regulations impose penalties, and in some cases rewards, which introduce a component of financial risk to an electric power utility due to the uncertainty associated with preserving a specific level of system reliability. In Brazil, for instance, one of the reliability indices receiving special attention by the utilities is the Maximum Continuous Interruption Duration per customer (MCID). This paper describes a chronological Monte Carlo simulation approach to evaluate probability distributions of reliability indices, including the MCID, and the corresponding penalties. In order to get the desired efficiency, modern computational techniques are used for modeling (UML -Unified Modeling Language) as well as for programming (Object- Oriented Programming). Case studies on a simple distribution network and on real Brazilian distribution systems are presented and discussed. © Copyright KTH 2006.