Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Scattering from phase-separated vesicles. I. An analytical form factor for multiple static domains

This is the first in a series of studies considering elastic scattering from laterally heterogeneous lipid vesicles containing multiple domains. Unique among biophysical tools, small-angle neutron scattering can in principle give detailed information about the size, shape and spatial arrangement of domains. A general theory for scattering from laterally heterogeneous vesicles is presented, and the analytical form factor for static domains with arbitrary spatial configuration is derived, including a simplification for uniformly sized round domains. The validity of the model, including series truncation effects, is assessed by comparison with simulated data obtained from a Monte Carlo method. Several aspects of the analytical solution for scattering intensity are discussed in the context of small-angle neutron scattering data, including the effect of varying domain size and number, as well as solvent contrast. Finally, the analysis indicates that effects of domain formation are most pronounced when the vesicle's average scattering length density matches that of the surrounding solvent.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Model-free reconstruction of magnetic correlations in frustrated magnets

Frustrated magnetic systems exhibit extraordinary physical properties, but quantification of their magnetic correlations poses a serious challenge to experiment and theory. Current insight into frustrated magnetic correlations relies on modelling techniques such as reverse Monte-Carlo methods, which require knowledge about the exact ordered atomic structure. Here, we present a method for direct reconstruction of magnetic correlations in frustrated magnets by three-dimensional difference pair distribution function analysis of neutron total scattering data. The methodology is applied to the disordered frustrated magnet bixbyite, (Mn 1– x Fe x ) 2 O 3 , which reveals nearest-neighbor antiferromagnetic correlations for the metal sites up to a range of approximately 15 Å. Pertinently, this technique allows for magnetic correlations to be determined directly from the experimental data without any assumption about the atomic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Progress Towards Prognostic Health Management of Passive Components in Advanced Reactors – Model Selection and Evaluation

This paper presents recent progress towards developing a prognostic health management framework for passive components of advanced reactors (AR). The focus of this paper is on lifecycle prognostics for passive components using a Bayesian prognostic algorithm that provides a natural framework for incorporating different sources of variability, and uncertainties inherent in the operations of AR. High-temperature creep damage, a prototypic failure mechanism in AR materials, is used as the context for this research. A Bayesian model selection approach is implemented to select the appropriate creep degradation model at any given time, using relevant sensor measurements reflecting the material degradation state. The model selection approach, based on reversible jump Markov chain Monte Carlo methods, is integrated with Bayesian particle filter-based prognostic framework. The proposed approach is evaluated using strain measurements obtained from accelerated creep testing of stainless steel specimens. Results indicate feasibility of the proposed approach in accurately identifying the creep degradation stage from the available measurements at a given time. Effect of uncertainties in material degradation model and measurement noise on the performance of the prognostic algorithm is also investigated.

advanced reactor prognostics↗
DOE OSTI2019

A Multi-Model Adaptive Kalman Filtering Approach to Power System Dynamic State Estimation

Accurate information about dynamic states (such as rotor angle and speed of a synchronous machine) is important for monitoring and controlling power system rotor-angle stability. In this paper, a multi-model adaptive Kalman filtering (MMAKF) approach is proposed to accurately and robustly estimate power system dynamic states. This approach consists of three major steps: (i) multiple Kalman filtering approaches, i.e., the extended Kalman filter (EKF), unscented Kalman filter (UKF), ensemble Kalman filter (EnKF), and cubature Kalman filter (CKF), are run concurrently in parallel to estimate the dynamic states of a synchronous generator using phasor measurement unit data; (ii) probability indexes, which quantify the likelihood of each estimation model, are determined at each time step using hypothesis testing based on the measurement innovation; (iii) the a posteriori estimate of states is obtained using the best-fix approach. The two-area four-machine system is used to evaluate the effectiveness of the proposed MMAKF approach. It is shown through the Monte-Carlo method that the estimation accuracy and robustness of the proposed approach is better than those from any individual filtering algorithm.

multi-model adaptive Kalman Filtering, Power Syste↗
DOE OSTI2017

Characterizing Quality Factor of Niobium Resonators Using a Markov Chain Monte Carlo Approach

The next generation of radiation detectors in high precision Cosmology, Astronomy, and particle-astrophysics experiments will rely heavily on superconducting microwave resonators and kinetic inductance devices. Understanding the physics of energy loss in these devices, in particular at low temperatures and powers, is vital. We present a comprehensive analysis framework, using Markov Chain Monte Carlo methods, to characterize loss due to two-level system in concert with quasi-particle dynamics in thin-film Nb resonators in the GHz range.

MKID↗
DOE OSTI2011

Particle-Filter-Based Multisensor Fusion For Solving Low-Frequency Electromagnetic NDE Inverse Problems

Flaw profile characterization from NDE measurements is a typical inverse problem. A novel transformation of this inverse problem into a tracking problem, and subsequent application of a sequential Monte Carlo method called particle filtering, has been proposed by the authors in an earlier publication [1]. In this study, the problem of flaw characterization from multi-sensor data is considered. The NDE inverse problem is posed as a statistical inverse problem and particle filtering is modified to handle data from multiple measurement modes. The measurement modes are assumed to be independent of each other with principal component analysis (PCA) used to legitimize the assumption of independence. The proposed particle filter based data fusion algorithm is applied to experimental NDE data to investigate its feasibility.

Data Fusion↗
DOE OSTI2018

A Multi-Step Adaptive Interpolation Approach to Mitigating the Impact of Nonlinearity on Dynamic State Estimation

Accurate estimation of dynamic states (such as synchronous machine rotor angle and speed) is critical for monitoring and controlling transient stability. Extended Kalman filter (EKF)--based approaches have been developed for estimating dynamic states. In order to improve the EKF's performance in estimating dynamic states of a synchronous machine, this paper proposes a multi-step adaptive interpolation (MSAI) approach to achieve balance between estimation accuracy and computational time. This approach consists of three major steps. First, two indexes are calculated to quantify the nonlinearity of the state transition function and measurement function, respectively. Second, based on the nonlinearity indexes, the interpolation factor is determined using a finite state machine model. And finally, to mitigate the negative impact of nonlinearity on the estimation accuracy, pseudo-measurements are added between consecutive measurements through linear interpolation. A simple example is used to validate the proposed nonlinear indexes. The two-area four-machine system and 16-machine 68-bus system are used to evaluate the effectiveness of the proposed MSAI approach. It is shown through the Monte-Carlo method that the estimation accuracy can be improved through interpolation. In addition, a good trade-off between estimation accuracy and computational time can be achieved effectively through the proposed MSAI approach.

Dynamic State Estimation, extended Kalman filter (↗
DOE OSTI2021

An Iterative Response-Surface-Based Approach for Chance-Constrained AC Optimal Power Flow Considering Dependent Uncertainty

A modern power system is characterized by a stochastic variation of the loads and an increasing penetration of renewable energy generation, which results in large uncertainties in its states. These uncertainties bring formidable challenges to the power system planning and operation process. To address these challenges, we propose a cost-effective, iterative response-surface-based approach for the chance-constrained AC optimal power-flow problem that aims to ensure the secure operation of the power systems considering dependent uncertainties. Starting from a stochastic-sampling-based framework, we first utilize the copula theory to simulate the dependence among multivariate uncertain inputs. Then, to reduce the prohibitive computational time required in the traditional Monte-Carlo method, we propose, instead of using the original complicated power-system model, to rely on a polynomial-chaos-based response surface. This response surface allows us to efficiently evaluate the time-consuming power-system model at arbitrary distributed sampled values with a negligible computational cost. This further enables us to efficiently conduct an online stochastic testing for the system states that not only screens out the statistical active constraints, but also assists in a better design of the tightened bounds without using any Gaussian or symmetric assumption. Finally, an iterative procedure is executed to fine-tune the optimal solution that better satisfies a predefined probability. The simulations conducted in multiple test systems demonstrate the excellent performance of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A Bayesian Approach for Estimating Uncertainty in Stochastic Economic Dispatch considering Wind Power Penetration

The increasing penetration of renewable energy resources in power systems, represented as random processes, converts the traditional deterministic economic dispatch problem into a stochastic one. To estimate the uncertainty in this stochastic economic dispatch problem for forecasting purposes, the conventional Monte-Carlo method is prohibitively time-consuming for practical applications. To overcome this problem, here we propose a novel Gaussian-process-emulator-based approach to quantify the uncertainty in the stochastic economic dispatch considering wind power penetration. Facing high-dimensional real-world data representing the correlated uncertainties from wind generation, a manifold-learning-based Isomap algorithm is proposed to efficiently represent the low-dimensional hidden probabilistic structure of the data. In this low-dimensional latent space, with Latin hypercube sampling as the computer experimental design, a Gaussian-process emulator is used, for the first time, to serve as a nonparametric, surrogate model for the original complicated stochastic economic dispatch model. This reduced-order representative allows us to evaluate the economic dispatch solver at sampled values with a negligible computational cost while maintaining a desirable accuracy. Simulation results conducted on the IEEE 118-bus test system reveal the impressive performance of the proposed method.

17 WIND ENERGY↗
DOE OSTI2026

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian↗
DOE OSTI2024

Development of Modified Perturbation Solutions to the One-Phase Stefan Problems With a Convective Boundary

The classical Stefan problem is used to track the moving solid-liquid interface during the freezing process. Perturbation theory has often been applied to find an approximate analytical solution due to the nonlinearity of the moving interface. However, the Stefan number (i.e., the sensible over latent heat) must be small and usually less than 0.01 to assume the perturbation expansion, which in turn limits the thermal engineering applications. In this study, a modified perturbation solution is developed by adding a correction term after the leading-order solution to be valid for a much wider range of Stefan numbers (i.e., 0.01 less than or equal to Ste less than or equal to 1). Specifically, a one-phase Stefan problem is first formulated subjected to a convective boundary in the Cartesian, cylindrical, and spherical coordinate systems for inward solidification. The leading-order solution is calculated based on the regular perturbation theory, while the correction term is obtained using the Monte-Carlo method and a multi-variant regression. Results show that the correction term has a linear relationship with the Stefan number and is not significantly influenced by the Biot number. The proposed modified perturbation solution can accurately and rapidly predict the nonlinear moving interface motion for the freezing process.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Electron contamination modeling and skin dose in 6 MV longitudinal field MRIgRT: Impact of the MRI and MRI fringe field

Purpose: In recent times, longitudinal field MRI-linac systems have been proposed for 6 MV MRI-guided radiotherapy (MRIgRT). The magnetic field is parallel with the beam axis and so will alter the transport properties of any electron contamination particles. The purpose of this work is to provide a first investigation into the potential effects of the MR and fringe magnetic fields on the electron contamination as it is transported toward a phantom, in turn, providing an estimate of the expected patient skin dose changes in such a modality. Methods: Geant4 Monte Carlo simulations of a water phantom exposed to a 6 MV x-ray beam were performed. Longitudinal magnetic fields of strengths between 0 and 3 T were applied to a 30 x 30 x 20 cm{sup 3} phantom. Surrounding the phantom there is a region where the magnetic field is at full MRI strength, consistent with clinical MRI systems. Beyond this the fringe magnetic field entering the collimation system is also modeled. The MRI-coil thickness, fringe field properties, and isocentric distance are varied and investigated. Beam field sizes of 5 x 5, 10 x 10, 15 x 15 and 20 x 20 cm{sup 2} were simulated. Central axis dose, 2D virtual entry skin dose films, and 70 {mu}m skin depth doses were calculated using high resolution scoring voxels. Results: In the presence of a longitudinal magnetic field, electron contamination from the linear accelerator is encouraged to travel almost directly toward the patient surface with minimal lateral spread. This results in a concentration of electron contamination within the x-ray beam outline. This concentration is particularly encouraged if the fringe field encompasses the collimation system. Skin dose increases of up to 1000% were observed for certain configurations and increases above Dmax were common. In nonmagnetically shielded cases, electron contamination generated from the jaw faces and air column is trapped and propagated almost directly to the phantom entry region, giving rise to intense dose hot spots inside the x-ray treatment field. These range up to 1000% or more of Dmax at the CAX, depending on field size, isocenter, and coil thickness. In the case of a fully magnetically shielded collimation system and the lowest MRI field of 0.25 T, the entry skin dose is expected to increase to at least 40%, 50%, 65%, and 80% of Dmax for 5 x 5, 10 x 10, 15 x 15, and 20 x 20 cm{sup 2}, respectively. Conclusions: Electron contamination from the linac head and air column may cause considerable skin dose increases or hot spots at the beam central axis on the entry side of a phantom or patient in longitudinal field 6 MV MRIgRT. This depends heavily on the properties of the magnetic fringe field entering the linac beam collimation system. The skin dose increase is also related to the MRI-coil thickness, the fringe field, and the isocenter distance of the linac. The results of this work indicate that the properties of the MRI fringe field, electron contamination production, and transport must be considered carefully during the design stage of a longitudinal MRI-linac system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Dedicated breast CT: radiation dose for circle-plus-line trajectory

Purpose: Dedicated breast CT prototypes used in clinical investigations utilize single circular source trajectory and cone-beam geometry with flat-panel detectors that do not satisfy data-sufficiency conditions and could lead to cone beam artifacts. Hence, this work investigated the glandular dose characteristics of a circle-plus-line trajectory that fulfills data-sufficiency conditions for image reconstruction in dedicated breast CT. Methods: Monte Carlo-based computer simulations were performed using the GEANT4 toolkit and was validated with previously reported normalized glandular dose coefficients for one prototype breast CT system. Upon validation, Monte Carlo simulations were performed to determine the normalized glandular dose coefficients as a function of x-ray source position along the line scan. The source-to-axis of rotation distance and the source-to-detector distance were maintained constant at 65 and 100 cm, respectively, in all simulations. The ratio of the normalized glandular dose coefficient at each source position along the line scan to that for the circular scan, defined as relative normalized glandular dose coefficient (RD{sub g}N), was studied by varying the diameter of the breast at the chest wall, chest-wall to nipple distance, skin thickness, x-ray beam energy, and glandular fraction of the breast. Results: The RD{sub g}N metric when stated as a function of source position along the line scan, relative to the maximum length of line scan needed for data sufficiency, was found to be minimally dependent on breast diameter, chest-wall to nipple distance, skin thickness, glandular fraction, and x-ray photon energy. This observation facilitates easy estimation of the average glandular dose of the line scan. Polynomial fit equations for computing the RD{sub g}N and hence the average glandular dose are provided. Conclusions: For a breast CT system that acquires 300-500 projections over 2{pi} for the circular scan, the addition of a line trajectory with equal source spacing and constant x-ray beam quality (kVp and HVL) and mAs matched to the circular scan, will result in less than 0.18% increase in average glandular dose to the breast per projection along the line scan.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Biological effect of dose distortion by fiducial markers in spot-scanning proton therapy with a limited number of fields: A simulation study

Purpose: In accurate proton spot-scanning therapy, continuous target tracking by fluoroscopic x ray during irradiation is beneficial not only for respiratory moving tumors of lung and liver but also for relatively stationary tumors of prostate. Implanted gold markers have been used with great effect for positioning the target volume by a fluoroscopy, especially for the cases of liver and prostate with the targets surrounded by water-equivalent tissues. However, recent studies have revealed that gold markers can cause a significant underdose in proton therapy. This paper focuses on prostate cancer and explores the possibility that multiple-field irradiation improves the underdose effect by markers on tumor-control probability (TCP). Methods: A Monte Carlo simulation was performed to evaluate the dose distortion effect. A spherical gold marker was placed at several characteristic points in a water phantom. The markers were with two different diameters of 2 and 1.5 mm, both visible on fluoroscopy. Three beam arrangements of single-field uniform dose (SFUD) were examined: one lateral field, two opposite lateral fields, and three fields (two opposite lateral fields + anterior field). The relative biological effectiveness (RBE) was set to 1.1 and a dose of 74 Gy (RBE) was delivered to the target of a typical prostate size in 37 fractions. The ratios of TCP to that without the marker (TCP{sub r}) were compared with the parameters of the marker sizes, number of fields, and marker positions. To take into account the dependence of biological parameters in TCP model, {alpha}/{beta} values of 1.5, 3, and 10 Gy (RBE) were considered. Results: It was found that the marker of 1.5 mm diameter does not affect the TCPs with all {alpha}/{beta} values when two or more fields are used. On the other hand, if the marker diameter is 2 mm, more than two irradiation fields are required to suppress the decrease in TCP from TCP{sub r} by less than 3%. This is especially true when multiple (two or three) markers are used for alignment of a patient. Conclusions: It is recommended that 1.5-mm markers be used to avoid the reduction of TCP as well as to spare the surrounding critical organs, as long as the markers are visible on x-ray fluoroscopy. When 2-mm markers are implanted, more than two fields should be used and the markers should not be placed close to the distal edge of any of the beams.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Dose response of selected ion chambers in applied homogeneous transverse and longitudinal magnetic fields

Purpose: The magnetic fields of an integrated MR-Linac system will alter the paths of electrons that produce ions in the ionization chambers. The dose response of selected ion chambers is evaluated in the presence of varying transverse and longitudinal magnetic fields. The investigation is useful in calibration of therapeutic x-ray beams associated with MR-Linac systems. Methods: The Monte Carlo code PENELOPE was used to model the irradiation of NE2571, and PR06C ionization chambers in the presence of a transverse and longitudinal (with respect to the photon beam) magnetic fields of varying magnitude. The long axis of each chamber was simulated both parallel and perpendicular to the incident photon beam for each magnetic field case. The dose deposited in each chamber for each case was compared to the case with zero magnetic field by means of a ratio. The PR06C chamber's response was measured in the presence of a transverse magnetic field with field strengths ranging from 0.0 to 0.2 T to compare to simulated results. Results: The simulations and measured data show that in the presence of a transverse magnetic field there is a considerable dose response (maximum of 11% near 1.0 T in the ion chambers investigated, which depends on the magnitude of magnetic field, and relative orientation of the magnetic field, radiation beam, and ion chamber. Measurements made with the PR06C chamber verify these results in the region of measurement. In contrast, a longitudinal magnetic field produces only a slight increase in dose response (2% at 1.5 T) that rises slowly with increasing magnetic field and is seemingly independent of chamber orientation. Response trends were similar for the two ion chambers and relative orientations considered, but slight variations are present from chamber to chamber. Conclusions: Care must be taken when making ion chamber measurements in a transverse magnetic field. Ion chamber responses vary not only with transverse field strength, but with chamber orientation and type, and can be considerable. Longitudinal magnetic fields influence ion chamber responses relatively little (2% at 1.5 T), and only at field strengths in excess of 1.0 T.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Out-of-field organ doses and associated radiogenic risks from para-aortic radiotherapy for testicular seminoma

Purpose: The aims of this study were to (a) calculate the radiation dose to out-of-field organs from radiotherapy for stage I testicular seminoma and (b) estimate the associated radiogenic risks. Methods: Monte Carlo methodology was employed to model radiation therapy with typical anteroposterior and posteroanterior para-aortic fields on an anthropomorphic phantom simulating an average adult. The radiation dose received by all main and remaining organs that defined by the ICRP publication 103 and excluded from the treatment volume was calculated. The effect of field dimensions on each organ dose was determined. Additional therapy simulations were generated by introducing shielding blocks to protect the kidneys from primary radiation. The gonadal dose was employed to assess the risk of heritable effects for irradiated male patients of reproductive potential. The lifetime attributable risks (LAR) of radiotherapy-induced cancer were estimated using gender- and organ-specific risk coefficients for patient ages of 20, 30, 40, and 50 years old. The risk values were compared with the respective nominal risks. Results: Para-aortic irradiation to 20 Gy resulted in out-of-field organ doses of 5.0–538.6 mGy. Blocked field treatment led to a dose change up to 28%. The mean organ dose variation by increasing or decreasing the applied field dimensions was 18.7% ± 3.9% and 20.8% ± 4.5%, respectively. The out-of-field photon doses increased the lifetime intrinsic risk of developing thyroid, lung, bladder, prostate, and esophageal cancer by (0.1–1.4)%, (0.4–1.1)%, (2.5–5.4)%, (0.2–0.4)%, and (6.4–9.2)%, respectively, depending upon the patient age at exposure and the field size employed. A low risk for heritable effects of less than 0.029% was found compared with the natural incidence of these defects. Conclusions: Testicular cancer survivors are subjected to an increased risk for the induction of bladder and esophageal cancer following para-aortic radiotherapy. The probability for the appearance of any other malignant disease to out-of-field organs was slightly elevated in respect to the nominal cancer incidence rates.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

MO-A-BRD-07: Feasibility of X-Ray Acoustic Computed Tomography as a Tool for Calibration and In Vivo Dosimetry of Radiotherapy Electron and Photon Beams

Purpose: This work simulates radiation-induced acoustic waves to assess the feasibility of x-ray acoustic computed tomography (XACT) as a dosimeter. XACT exploits the phenomenon that acoustic waves with amplitude proportional to the dose deposited are induced following a radiation pulse. After detecting these acoustic waves with an ultrasound transducer, an image of the dose distribution can be reconstructed in realtime. Methods: Monte Carlo was used to simulate the dose distribution for monoenergetic 6 MeV photon and 9 MeV electron beams incident on a water tank. The dose distribution for a prostate patient planned with a photon 4-field box technique was calculated using clinical treatment planning software. All three dose distributions were converted into initial pressure distributions, and transportation of the induced acoustic waves was simulated using an open-source toolkit. Ideal transducers were placed around the circumference of the target to detect the acoustic waves, and a time reversal reconstruction algorithm was used to obtain an XACT image of the dose for each radiation pulse. Results: For the photon water tank relative dosimetry case, it was found that the normalized acoustic signal amplitude agreed with the normalized dose at depths from 0 cm to 10 cm, with an average percent difference of 0.5%. For the reconstructed in-plane dose distribution of an electron water tank irradiation, all pixels passed a 3%–3 mm 2D gamma test. The reconstructed prostate dose distribution closely resembled the plan, with 89% of pixels passing a 3%–3 mm 2D gamma test. For all situations, the amplitude of the induced acoustic waves ranged from 0.01 Pa to 1 Pa. Conclusion: Based on the amplitude of the radiation-induced acoustic waves and accuracy of the reconstructed dose distributions, XACT is a feasible technique for dosimetry in both calibration and in vivo environments for photon and electron beams and merits further investigation. Funding from NSERC, CIHR and McGill University. SH acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-G-BRE-02: Biological Modeling of Gold Nanoparticle Radiosensitization for Proton Therapy

Purpose: The aim of this work is to investigate the radiosensitization effect of gold nanoparticles (GNP) in a proton beam. A computational model was built using the Local Effect Model (LEM) to predict the biological outcome of gold nanoparticle (GNP) sensitization. We present the results using a clinical proton beam, 6MV photon beam and two kilovoltage photon beams. Methods: First, Monte Carlo simulations were carried out using TOPAS (TOol for PArticle Simulation) to obtain the spatial dose distribution in the vicinity of GNPs. The dose distribution was then used as an input for LEM, which predicts dose-response curves for high linear energy transfer radiation using the track structure. The cell survival curves were evaluated for three particle sources (proton beam, MV photon beam and kV photon beam), various treatment depths for each particle source, various GNP uptakes and two different GNP sizes. Results: For proton therapy, the GNP sensitization effect is highly dependent on the treatment depth due to the energy-dependent interaction probability. We predict that if GNPs can be taken up by the cell nucleus, proton therapy can be significantly enhanced. If GNPs are only internalized into the cytoplasm, proton therapy can still be enhanced by GNPs and if GNPs are not internalized into cells, there will be no direct damage to the nucleus. For the same GNP uptake and concentration, the cell survival at 2Gy is reduced by 80% using kilovoltage photons, 50% using protons and only 2% using clinical MV photons. Finally, for the same weight of GNPs taken up by the cells, 10 nm GNPs causes 3 times more damage than 50 nm GNPs. Conclusion: We showed that GNPs have potential to be used to enhance radiation therapy for clinical proton beams.

07 ISOTOPE AND RADIATION SOURCES↗