Monte Carlo simulation for scanning technique with scattering foil free electron beam: A proof of concept study.
Ontology highlight
ABSTRACT:
SUBMITTER: Sung W
PROVIDER: S-EPMC5426680 | biostudies-literature | 2017
REPOSITORIES: biostudies-literature
Publications
Monte Carlo simulation for scanning technique with scattering foil free electron beam: A proof of concept study.
Sung Wonmo W Park Jong In JI Kim Jung-In JI Carlson Joel J Ye Sung-Joon SJ Park Jong Min JM
PloS one 20170511 5
This study investigated the potential of a newly proposed scattering foil free (SFF) electron beam scanning technique for the treatment of skin cancer on the irregular patient surfaces using Monte Carlo (MC) simulation. After benchmarking of the MC simulations, we removed the scattering foil to generate SFF electron beams. Cylindrical and spherical phantoms with 1 cm boluses were generated and the target volume was defined from the surface to 5 mm depth. The SFF scanning technique with 6 MeV ele ...[more]
Similar Datasets
Project description:The novel coronavirus pneumonia triggered by COVID-19 is now raging the whole world. As a rapid and reliable killing COVID-19 method in industry, electron beam irradiation can interact with virus molecules and destroy their activity. With the unexpected appearance and quickly spreading of the virus, it is urgently necessary to figure out the mechanism of electron beam irradiation on COVID-19. In this study, we establish a virus structure and molecule model based on the detected gene sequence of Wuhan patient, and calculate irradiated electron interaction with virus atoms via a Monte Carlo simulation that track each elastic and inelastic collision of all electrons. The characteristics of irradiation damage on COVID-19, atoms' ionizations and electron energy losses are calculated and analyzed with regions. We simulate the different situations of incident electron energy for evaluating the influence of incident energy on virus damage. It is found that under the major protecting of an envelope protein layer, the inner RNA suffers the minimal damage. The damage for a ∼100-nm-diameter virus molecule is not always enhanced by irradiation energy monotonicity, for COVID-19, the irradiation electron energy of the strongest energy loss damage is 2 keV.
| S-EPMC7351116 | biostudies-literature
Project description:Recent technical and strategical developments have increased the survival chances for avalanche victims. Still hundreds of people, primarily recreationists, get caught and buried by snow avalanches every year. About 100 die each year in the European Alps-and many more worldwide. Refining concepts for avalanche rescue means to optimize the procedures such that the survival chances are maximized in order to save the greatest possible number of lives. Avalanche rescue includes several parameters related to terrain, natural hazards, the people affected by the event, the rescuers, and the applied search and rescue equipment. The numerous parameters and their complex interaction make it unrealistic for a rescuer to take, in the urgency of the situation, the best possible decisions without clearly structured, easily applicable decision support systems. In order to analyse which measures lead to the best possible survival outcome in the complex environment of an avalanche accident, we present a numerical approach, namely a Monte Carlo simulation. We demonstrate the application of Monte Carlo simulations for two typical, yet tricky questions in avalanche rescue: (1) calculating how deep one should probe in the first passage of a probe line depending on search area, and (2) determining for how long resuscitation should be performed on a specific patient while others are still buried. In both cases, we demonstrate that optimized strategies can be calculated with the Monte Carlo method, provided that the necessary input data are available. Our Monte Carlo simulations also suggest that with a strict focus on the "greatest good for the greatest number", today's rescue strategies can be further optimized in the best interest of patients involved in an avalanche accident.
| S-EPMC5414947 | biostudies-literature
Project description:Conventional mammography can suffer from poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter slot scan imaging is an imaging technique which provides additional information and is compatible with conventional mammography. A Monte Carlo simulation of coherent scatter slot scan imaging was performed to assess its performance and provide system optimization. Coherent scatter could be exploited using a system similar to conventional slot scan mammography system with antiscatter grids tilted at the characteristic angle of cancerous tissues. System optimization was performed across several parameters, including source voltage, tilt angle, grid distances, grid ratio, and shielding geometry. The simulated carcinomas were detectable for tumors as small as 5 mm in diameter, so coherent scatter analysis using a wide-slot setup could be promising as an enhancement for screening mammography. Employing coherent scatter information simultaneously with conventional mammography could yield a conventional high spatial resolution image with additional coherent scatter information.
| S-EPMC4999523 | biostudies-other
Project description:Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4-167.1 cc) and motion amplitude (2.9-30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (?es), time required for magnet settling (?ss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior-inferior motion amplitude alone. Larger spot sizes (? ~ 9-16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (? ~ 2-4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes respectively. The interplay effect is highly patient specific, depending on the motion amplitude, tumor location and the delivery parameters. Large degradations of the dose distribution in a single fraction were observed, but improved significantly using conventional fractionation.
| S-EPMC3752993 | biostudies-literature
Project description:Background:Proton pencil beam (PB) dose calculation algorithms have limited accuracy within heterogeneous tissues of lung cancer patients, which may be addressed by modern commercial Monte Carlo (MC) algorithms. We investigated clinical pencil beam scanning (PBS) dose differences between PB and MC-based treatment planning for lung cancer patients. Methods:With IRB approval, a comparative dosimetric analysis between RayStation MC and PB dose engines was performed on ten patient plans. PBS gantry plans were generated using single-field optimization technique to maintain target coverage under range and setup uncertainties. Dose differences between PB-optimized (PBopt), MC-recalculated (MCrecalc), and MC-optimized (MCopt) plans were recorded for the following region-of-interest metrics: clinical target volume (CTV) V95, CTV homogeneity index (HI), total lung V20, total lung VRX (relative lung volume receiving prescribed dose or higher), and global maximum dose. The impact of PB-based and MC-based planning on robustness to systematic perturbation of range (±3% density) and setup (±3 mm isotropic) was assessed. Pairwise differences in dose parameters were evaluated through non-parametric Friedman and Wilcoxon sign-rank testing. Results:In this ten-patient sample, CTV V95 decreased significantly from 99-100% for PBopt to 77-94% for MCrecalc and recovered to 99-100% for MCopt (P<10-5). The median CTV HI (D95/D5) decreased from 0.98 for PBopt to 0.91 for MCrecalc and increased to 0.95 for MCopt (P<10-3). CTV D95 robustness to range and setup errors improved under MCopt (?D95 =-1%) compared to MCrecalc (?D95 =-6%, P=0.006). No changes in lung dosimetry were observed for large volumes receiving low to intermediate doses (e.g., V20), while differences between PB-based and MC-based planning were noted for small volumes receiving high doses (e.g., VRX). Global maximum patient dose increased from 106% for PBopt to 109% for MCrecalc and 112% for MCopt (P<10-3). Conclusions:MC dosimetry revealed a reduction in target dose coverage under PB-based planning that was regained under MC-based planning along with improved plan robustness. MC-based optimization and dose calculation should be integrated into clinical planning workflows of lung cancer patients receiving actively scanned proton therapy.
| S-EPMC5960654 | biostudies-literature
Project description:The electron Monte Carlo (eMC) dose calculation algorithm of the Eclipse treatment planning system is based heavily upon Monte Carlo simulation of the linac head and modeling of the linac beam characteristics with minimal measurement of beam data. Commissioning of the eMC algorithm on multiple identical linacs provided a unique opportunity to systematically evaluate the algorithm with actual measurements of clinically relevant beam and dose parameters. In this study, measured and eMC calculated dose distributions were compared both along and perpendicular to electron beam direction for electron energy/applicator/depth combination using measurement data from four Varian 21EX CLINAC linear accelerator (Varian Medical System, Palo Alto, CA). Cutout factors for sizes down to 3 x 3 cm were also compared. Comparisons between the measurement and the eMC calculated values show that the R90, R80, R50, and R10 values mostly agree within 3 mm. Measure and Calculated bremsstrahlung dose Dx correlates well statistically although eMC calculated Dx values are consistently smaller than the measured, with maximum discrepancy of 1% for the 20 MeV electron beams. Surface dose agrees mostly within 2%. Field width and penumbra agree mostly within 3mm. Calculation grid size is found to have a significant effect on the dose calculation. A grid size of 5 mm can produce erroneous dose distributions. Using a grid size of 2.5 mm and a 3% accuracy specified for the eMC to stop calculation iteration, the absolute output agrees with measurements within 3% for field sizes of 5 x 5 cm or larger. For cutout of 3 x 3 cm, however, the output disagreement can reach 8%. Our result indicate that eMC algorithm in Eclipse provides acceptable agreement with measurement data for most clinical situations. Calculation grid size of 2.5 mm or smaller is recommended.
| S-EPMC5722292 | biostudies-literature
Project description:Traditional multiple-group confirmatory factor analysis (multiple-group CFA) is usually criticized for having too restrictive model assumption, namely the scalar measurement invariance. The new multiple-group analysis methodology, alignment, has become an effective alternative. The alignment evaluates measurement invariance and more importantly, permits factor mean comparisons without requiring scalar invariance which is usually required in traditional multiple-group CFA. Some simulation studies and empirical studies have investigated the applicability of alignment under different conditions, but some areas remain unexplored. Based on the simulation studies of Asparouhov and Muthén and of Flake and McCoach, this current simulation study is broken into two sections. The first study investigates the minimal group sizes required for alignment in three-factor models. The second study compares the performance of multiple-group CFA, multiple-group exploratory structural equation model (multiple-group ESEM), and alignment by including different proportions and magnitudes of cross-loadings in the models. Study 1 shows that when the model has no noninvariant parameters, the alignment requires relatively lower group sizes. Explicitly, the minimal group size required for alignment was 250 when the amount of groups was three, the minimal group size was 150 when the amount of groups was nine, and 200 when the amount of groups was 15. When there are noninvariant parameters in the model and the amount of groups is low, a group size of 350 is a safe rule of thumb. When there are noninvariant parameters in the model and the amount of groups is high, a group size of 250 is required for trustworthy results. The magnitude of noninvariance and the noninvariance rate do not affect the minimal group size required for alignment. Study 2 shows that multiple-group CFA provides accurate factor mean estimates when each factor had 20% factor loading (1 factor loading) with small-sized cross-loading. Multiple-group ESEM provides accurate factor mean estimates when the magnitude of cross-loading is small or when each factor had 20% factor loading (1 factor loading) with medium-sized cross-loading. Alignment provides accurate factor mean estimates when there are only small-sized cross-loadings in the model. The parameter estimates, coverage rates and ratios of average standard error to standard deviation for each methodology are not influenced by the amount of groups. Recommendations are concluded for using multiple-group CFA, multiple-group ESEM, traditional alignment and aligned ESEM (AESEM) based on the results. Multiple-group CFA is more suitable for use when scalar invariance is established. Multiple-group ESEM works best when there are small-sized or only a few medium-sized cross-loadings in the model. Traditional alignment allows for small-sized cross-loadings and a few noninvariant parameters in the model. AESEM integrates the advantages of alignment and ESEM, can provide accurate estimates when noninvariant parameters and cross-loadings both exist in the model. Compared to multiple-group CFA, multiple-group ESEM, the alignment methodology performs well in more situations.
| S-EPMC9263979 | biostudies-literature
Project description:A computational method is developed to carry out explicit solvent simulations of complex molecular systems under conditions of constant pH. In constant-pH simulations, preidentified ionizable sites are allowed to spontaneously protonate and deprotonate as a function of time in response to the environment and the imposed pH. The method, based on a hybrid scheme originally proposed by H. A. Stern (J. Chem. Phys. 2007, 126, 164112), consists of carrying out short nonequilibrium molecular dynamics (neMD) switching trajectories to generate physically plausible configurations with changed protonation states that are subsequently accepted or rejected according to a Metropolis Monte Carlo (MC) criterion. To ensure microscopic detailed balance arising from such nonequilibrium switches, the atomic momenta are altered according to the symmetric two-ends momentum reversal prescription. To achieve higher efficiency, the original neMD-MC scheme is separated into two steps, reducing the need for generating a large number of unproductive and costly nonequilibrium trajectories. In the first step, the protonation state of a site is randomly attributed via a Metropolis MC process on the basis of an intrinsic pKa; an attempted nonequilibrium switch is generated only if this change in protonation state is accepted. This hybrid two-step inherent pKa neMD-MC simulation method is tested with single amino acids in solution (Asp, Glu, and His) and then applied to turkey ovomucoid third domain and hen egg-white lysozyme. Because of the simple linear increase in the computational cost relative to the number of titratable sites, the present method is naturally able to treat extremely large systems.
| S-EPMC4535364 | biostudies-literature
Project description:BACKGROUND: In the biological sciences the TCID50 (median tissue culture infective dose) assay is often used to determine the strength of a virus. METHODS: When the so-called Spearman-Kaerber calculation is used, the ratio between the pfu (the number of plaque forming units, the effective number of virus particles) and the TCID50, theoretically approaches a simple function of Eulers constant. Further, the standard deviation of the logarithm of the TCID50 approaches a simple function of the dilution factor and the number of wells used for determining the ratios in the assay. However, these theoretical calculations assume that the dilutions of the assay are independent, and in practice this is not completely correct. The assay was simulated using Monte Carlo techniques. RESULTS: Our simulation studies show that the theoretical results actually hold true for practical implementations of the assay. Furthermore, the simulation studies show that the distribution of the (the log of) TCID50, although discrete in nature, has a close relationship to the normal distribution. CONCLUSION: The pfu is proportional to the TCID50 titre with a factor of about 0.56 when using the Spearman-Kaerber calculation method. The normal distribution can be used for statistical inferences and ANOVA on the (the log of) TCID50 values is meaningful with group sizes of 5 and above.
| S-EPMC3331834 | biostudies-other
Project description:Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10(7) - 10(9) detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS.The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels.On the order of only 10-100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10(7) reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10(6) CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU.The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.
| S-EPMC3978362 | biostudies-other