Intensity-modulated brachytherapy techniques, compensator-based intensity modulated brachytherapy (CBT) and interstitial rotating shield brachytherapy (I-RSBT), are two novel conceptual radiation therapies for treating cervical and prostate cancer, respectively. Compared to conventional brachytherapy techniques for treating cervical cancer, CBT can potentially improve the dose conformity to the high-risk clinical target volume (CTV) of the cervix in a less invasive approach. I-RSBT can reduce the dose delivered to the prostate organ at risks (OARs) with the same radiation dose delivered to the prostate CTV. In this work, concepts and prototypes for CBT and I-RSBT were introduced and developed. Preliminary dosimetric measurements were performed for CBT and I-RSBT, respectively. A CBT prototype system was constructed and experimentally validated. A prototype cylindrical compensator with eight octants, each with different thicknesses, was designed. Direct metal laser sintering (DMLS) was used to construct CoCr and Ti compensator prototypes, and a 4-D milling technique was used to construct a Ti compensator prototype. Gafchromic EBT2 films, held by an acrylic quality assurance (QA) phantom, were irradiated to approximately 125 cGy with an electronic brachytherapy (eBT) source for both shielded and unshielded cases. The dose at each point on the films were calculated using a TG-43 calculation model that was modified to account for the presence of a compensator prototype by ray-tracing. With I-RSBT, a multi-pass dose delivery mechanism with prototypes was developed. Dosimetric measurements for a Gd-153 radioisotope was performed to demonstrate that using multiple partially shielded Gd-153 sources for I-RSBT is feasible. A treatment planning model was developed for applying I-RSBT clinically. A custom-built, stainless steel encapsulated 150 mCi Gd-153 capsule with an outer length of 12.8 mm, outer diameter of 2.10 mm, active length of 9.98 mm, and active diameter of 1.53 mm was used. A partially shielded catheter was constructed with a 500 micron platinum shield and a 500 micron aluminum emission window, both with 180° azimuthal coverage. An acrylic phantom was constructed to measure the dose distributions from the shielded catheter in the transverse plane using Gafchromic EBT3 films. Film calibration curves were generated from 50, 70, and 100 kVp x-ray beams with NIST-traceable air kerma values to account for energy variation. In conclusion, CBT, which is a non-invasive alternative to supplementary interstitial brachytherapy, is expected to improve dose conformity to bulky cervical tumors relative to conventional intracavitary brachytherapy. However, at the current stage, it would be time-consuming to construct a patient-specific compensator using DMLS, and the quality assurance of the compensator would be difficult. I-RSBT is a promising approach to reducing radiation dose delivered to prostate OARs. The next step in making Gd-153 based I-RSBT feasible in clinic is developing a Gd-153 source that is small enough such that the source, shield, and catheter all fit within a 16 guage needle, which has a 1.65 mm diameter.

Ce travail de thèse a pour but d'améliorer la précision de mesure de la dosimétrie basée sur le transducteur en Nitrure de Gallium (GaN) et de développer son application en radiothérapie. L'étude comprend des phases de modélisation, de simulation et de caractérisation de cette réponse pour la radiothérapie externe et la curiethérapie. En modélisation, nous avons proposé deux approches pour modéliser la réponse du transducteur GaN en radiothérapie externe. Dans la première approche, un modèle a été construit à partir de données expérimentales et en séparant les composantes primaires et diffusées du faisceau. Pour la deuxième approche, nous avons adapté un modèle initialement proposé pour les diodes silicium pour l'adapter au transducteur radioluminescent GaN. Nous avons également proposé un concept original de dosimétrie bi-média qui permet à partir des réponses mesurées des deux média de remonter à la dose dans les tissus, sans connaissance à priori des conditions d'irradiation. Ce concept a été démontré par des simulations Monte Carlo. Par ailleurs GaN pour la curiethérapie à Haut Débit de Dose, la réponse du transducteur GaN sous irradiation des sources d'iridium 192 et de cobalt 60 a été évaluée par simulation Monte Carlo et confirmée par des mesures. Des études de caractérisation des propriétés du transducteur radioluminescent GaN ont été réalisées avec ces sources. Un prototype de fantôme instrumenté avec des sondes GaN a été développé pour le contrôle qualité en curiethérapie HDR. Il permet de vérifier en temps réel les paramètres physiques du traitement (position de la source, le temps d'exposition, activité de la source) / The work in this thesis has the objective to increase the measurement precision of the dosimetry based on the Gallium Nitride (GaN) transducer and develop its applications on radiotherapy. The study includes the aspects of modeling, simulation and characterization of this response in external radiotherapy and brachytherapy. In modeling, we have proposed two approaches to model the GaN transducer’s response in external radiotherapy. For the first approach, a model has been built based on experimental data, while separating the primary and scattering component of the beam. For the second approach, we have adopted a response model initially developed for the silicon diodes for the GaN radioluminescent transducer. We have also proposed an original concept of bi-media dosimetry which evaluates the dose in tissue according to different responses from two media without prior information on the conditions of irradiation. This concept has been shown by Monte Carlo simulation. Moreover, for High Dose Rate brachytherapy, the response of GaN transducer irradiated by iridium 192 and cobalt 60 sources has been evaluated by Monte Carlo simulation and confirmed by the measurements. Studies on the property characterization of GaN radioluminescent transducer has been carried out with these sources as well. An instrumented phantom prototype with GaN probe has been developed for the HDR brachytherapy quality control. It allows a real-time verification of the physics parameters of a treatment (source dwell position, source dwell time, source activity)

No presente trabalho, nós desenvolvemos um modelo de mecânica estatística para a dinâmica do cotransporte. Em particular, foram analisados o simporte de lactose e H+ catalisado pela LacY, e o antiporte de triose-P e Pi catalisado pelo TPT. Em ambos os casos, nosso ponto de partida foi determinar a curva de energia livre do ciclo de transporte com base na estrutura bioquímica de cada estado. Então, usando o algoritmo de Metropolis em um modelo do caminho aleatório não-homogêneo, nós determinamos a dinâmica do ciclo expressa na curva de energia livre obtida anteriormente. Tanto no caso da LacY quanto no caso do TPT, a concordância entre as predições do nosso modelo e os dados experimentais sugere que a nossa curva de energia livre é apropriada para descrever os processos de transporte através da membrana plasmática. / In the present work, we developed a statistical mechanical model for the dynamics of cotransport. In particular, the symport of lactose and H+ and the antiport of triose-P and Pi were analyzed. In both cases, our starting point was to determine the free-energy curve of the transport cycle based on the biochemical structure of each state. Then, using the Metropolis algorithm in a nonhomogeneous random walk model, we determined the dynamics of the cycle expressed in the free-energy curve obtained previously. Both in the case of LacY and TPT, the agreement between our model predictions and the experimental data suggests that our free-energy curve is appropriate for describing the transport process across the plasma membrane.

<p>The aim of this work is to characterize the beam of the 60Co therapy unit “Siemens Gammatron 1”, used at the Swedish Radiation Protection Authority (SSI) to calibrate therapy level ionization chambers. SSI wants to know the spectra in the laboratory’s reference points and a verified, virtual model of the 60Co unit to be able to compare current and future experiments to Monte Carlo simulations.</p><p>EGSnrc is a code for performing Monte Carlo simulations. By using BEAMnrc, which is an additional package that simplifies the building process of a geometry in the EGS-code, the whole Gammatron at SSI was defined virtually. In this work virtual models for two experimental setups were built: the Gammatron irradiating in air to simulate the air-kerma calibration geometry and the Gammatron irradiating a water phantom similar to that used for the absorbed dose to water calibrations.</p><p>The simulations are divided into two different substeps: one for the fixed part of the Gammatron and one for the variable part to be able to study different entities and to shorten simulation times.</p><p>The virtual geometries are verified by comparing Monte Carlo results with measurements. When it was verified that the virtual geometries were to be trusted, they were used to generate the Gammatron photon spectra in air and water with different field sizes and at different depths. The contributions to the photon spectra from different regions in the Gammatron were also collected. This is something that is easy to achieve with Monte Carlo calculations, but difficult to obtain with ordinary detectors in real life measurements.</p><p>The results from this work give SSI knowledge of the photon spectra in their reference points for calibrations in air and in water phantom. The first step of the virtual model (fixed part of Gammatron) can be used for future experimental setups at SSI.</p>

<p>Due to its close tissue-equivalence, high radiation sensitivity, dose and dose-rate linearity, diamond is a very promising detector for radiation therapy applications. The present thesis focuses on the development of a chemical vapour deposited (CVD) diamond detector with special attention on the arrangement of the electrodes and encapsulation having minimal influence on the measured signal. Several prototype detectors were designed by using CVD-diamond substrates with attached silver electrodes.</p><p>Interface effects in the electrode-diamond-electrode structure are investigated using the Monte Carlo (MC) code PENELOPE. The studies cover a wide range of electrode and diamond thicknesses, electrode materials and photon beam energies. An appreciable enhancement of the absorbed dose to diamond was found for high-Z electrodes. The influence of the electrodes diminishes with decreasing atomic number difference and layer thickness, so that from this point of view thin graphite electrodes would be ideal. The effect of encapsulation, cable and electrical connections on the detector response is also addressed employing MC techniques. For Co-60, 6 and 18 MV photon beam qualities it is shown that the prototypes exhibit energy and directional dependence of about 3% and 2%, respectively. By modifying the geometry and using graphite electrodes the dependencies are reduced to 1%.</p><p>Although experimental studies disclose some limitations of the prototypes (high leakage current, priming effect and slow signal stabilisation), diamonds of higher quality, suitable for dosimetry, can be produced with better-controlled CVD process. With good crystals and a well-designed encapsulation, the CVD-diamond detector could become competitive for routine dosimetry. It is then important for correct dose determination to use a collision stopping power for diamond incorporating proper mean excitation energy and density-effect corrections. A new mean excitation energy of 88 eV has been calculated.</p>

<p>The purpose of this thesis is to provide a detailed description of recent methods and scientific basis for characterizing earthquake sources within a certain region with distinct tectonic environments. The focus will be on those characteristics that are most significant to the ground-shaking hazard and on how we can incorporate our current knowledge into hazard analyses for engineering design purposes. I treat two particular geographical areas where I think current hazard analysis methods are in need of significant improvement, and suggest some approaches that have proven to be effective in past applications elsewhere. A combined hazard procedure is used to estimate seismicity in <i>northern Central America</i>, where there appear to be four tectonic environments for modeling the seismogenic sources and in <i>Iran</i>, where the large earthquakes usually occur on known faults. A preferred seismic hazard model for northern Central America and the western Caribbean plate based on earthquake catalogs, geodetic measurements, and geological information is presented. I used the widely practiced method of relating seismicity data to geological data to assess the various seismic hazard parameters and test parameter sensitivities. </p><p>The sensitivity and overall uncertainty in peak ground acceleration (PGA) estimates are calculated for northwestern Iran by using a <i>specific randomized blocks design</i>. A Monte Carlo approach is utilized to evaluate the ground motion hazard and its uncertainties in northern Central America. A set of new seismic hazard maps, exhibiting probabilistic values of peak ground acceleration (PGA) with 50%, 10%, and 5% probabilities of exceedance (PE) in 50 years, is presented for the area of relevance. <i>Disaggregation of seismic hazard</i> is carried out for cities of San Salvador and Guatemala by using a spatial distribution of epicenters around these sites to select design ground motion for seismic risk decisions. </p><p>In conclusion, consideration of the effect of parameters such as seismic moment, fault rupture, rupture directivity and stress drop are strongly recommended in estimating the near field ground motions. The rupture process of the 2002 Changureh earthquake (<i>M</i><i>w</i> = 6.5), Iran, was analyzed by using the<i> empirical Green’s function (EGF) method</i>. This method simulates strong ground motions for future large earthquakes at particular sites where no empirical data are available.</p>

Computation methods based on the Wiener chaos expansion have been developed to study the behaviors of the aeroelastic system with randomparameters. It is proven that the discrete wavelet transformation is one ofthe most accurate and efficient numerical schemes for this uncertainty quantizationproblem. In this thesis, we propose the stochastic collocation methods(SCM), whichis a type of sampling method combining the strength of the MonteCarlo simulation and the stochastic Galerkin method. The convergence with respect to the number of the nodal points is investigated, and simulation results to aeroelastic models in the presence of uncertainty in the system parameter and due to the initial condition are reported. It is demonstrated that the accuracy of the SCM is comparable to those achieved by using the wavelet chaos expansion. However, the SCM is more straightforward, efficient and easy to implement. / Applied Mathematics

Net pay (NP) and net-to-gross ratio (NGR) are often crucial quantities to characterize a reservoir and assess the amount of hydrocarbons in place. Numerous methods in the industry have been developed to evaluate NP and NGR, depending on the intended purposes. These methods usually involve the use of cut-off values of one or more surrogate variables to discriminate non-reservoir from reservoir rocks. This study investigates statistical issues related to the selection of such cut-off values by considering the specific case of using porosity () as the surrogate. Four methods are applied to permeability-porosity datasets to estimate porosity cut-off values. All the methods assume that a permeability cut-off value has been previously determined and each method is based on minimizing the prediction error when particular assumptions are satisfied. The results show that delineating NP and evaluating NGR require different porosity cut-off values. In the case where porosity and the logarithm of permeability are joint normally distributed, NP delineation requires the use of the Y-on-X regression line to estimate the optimal porosity cut-off while the reduced major axis (RMA) line provides the optimal porosity cut-off value to evaluate NGR. Alternatives to RMA and regression lines are also investigated, such as discriminant analysis and a data-oriented method using a probabilistic analysis of the porosity-permeability crossplots. Joint normal datasets are generated to test the ability of the methods to predict accurately the optimal porosity cut-off value for sampled sub datasets. These different methods have been compared to one another on the basis of the bias, standard error and robustness of the estimates. A set of field data has been used from the Travis Peak formation to test the performance of the methods. The conclusions of the study have been confirmed when applied to field data: as long as the initial assumptions concerning the distribution of data are verified, it is recommended to use the Y-on-X regression line to delineate NP while either the RMA line or discriminant analysis should be used for evaluating NGR. In the case where the assumptions on data distribution are not verified, the quadrant method should be used.

This thesis addresses the challenges encountered when dealing with signal decomposition problems with an unknown number of components in a Bayesian framework. Particularly, we focus on the issue of summarizing the variable-dimensional posterior distributions that typically arise in such problems. Such posterior distributions are defined over union of subspaces of differing dimensionality, and can be sampled from using modern Monte Carlo techniques, for instance the increasingly popular Reversible-Jump MCMC (RJ-MCMC) sampler. No generic approach is available, however, to summarize the resulting variable-dimensional samples and extract from them component-specific parameters. One of the main challenges that needs to be addressed to this end is the label-switching issue, which is caused by the invariance of the posterior distribution to the permutation of the components. We propose a novel approach to this problem, which consists in approximating the complex posterior of interest by a "simple"--but still variable-dimensional parametric distribution. We develop stochastic EM-type algorithms, driven by the RJ-MCMC sampler, to estimate the parameters of the model through the minimization of a divergence measure between the two distributions. Two signal decomposition problems are considered, to show the capability of the proposed approach both for relabeling and for summarizing variable dimensional posterior distributions: the classical problem of detecting and estimating sinusoids in white Gaussian noise on the one hand, and a particle counting problem motivated by the Pierre Auger project in astrophysics on the other hand.

Due to its close tissue-equivalence, high radiation sensitivity, dose and dose-rate linearity, diamond is a very promising detector for radiation therapy applications. The present thesis focuses on the development of a chemical vapour deposited (CVD) diamond detector with special attention on the arrangement of the electrodes and encapsulation having minimal influence on the measured signal. Several prototype detectors were designed by using CVD-diamond substrates with attached silver electrodes. Interface effects in the electrode-diamond-electrode structure are investigated using the Monte Carlo (MC) code PENELOPE. The studies cover a wide range of electrode and diamond thicknesses, electrode materials and photon beam energies. An appreciable enhancement of the absorbed dose to diamond was found for high-Z electrodes. The influence of the electrodes diminishes with decreasing atomic number difference and layer thickness, so that from this point of view thin graphite electrodes would be ideal. The effect of encapsulation, cable and electrical connections on the detector response is also addressed employing MC techniques. For Co-60, 6 and 18 MV photon beam qualities it is shown that the prototypes exhibit energy and directional dependence of about 3% and 2%, respectively. By modifying the geometry and using graphite electrodes the dependencies are reduced to 1%. Although experimental studies disclose some limitations of the prototypes (high leakage current, priming effect and slow signal stabilisation), diamonds of higher quality, suitable for dosimetry, can be produced with better-controlled CVD process. With good crystals and a well-designed encapsulation, the CVD-diamond detector could become competitive for routine dosimetry. It is then important for correct dose determination to use a collision stopping power for diamond incorporating proper mean excitation energy and density-effect corrections. A new mean excitation energy of 88 eV has been calculated.