Comments on the primary and scatter dose-spread kernels used for convolution methods

The basic primary and scatter dose-spread kernels used for convolution methods are usually produced by Monte Carlo simulations with the interaction point forced to the center of a large water phantom. However, it is still not clear whether such Monte Carlo based kernels allow accurate dose calculations with a wide range of field sizes and depths, especially in thorax phantoms. Using the differential primary and scatter concept, this paper proposes another type of basic kernel, with which perfectly accurate primary and scatter absorbed dose calculations can be performed under conditions that the beam is parallel, the incident beam intensity is uniform within and zero outside the field, and the primary beam attenuation coefficient along raylines is not a function of depth and off-axis distance.     

Monte Carlo analysis of the electron thermalization process in the afterglow of a microsecond dc pulsed glow discharge

A Monte Carlo model is utilized for studying the behavior of electrons in the afterglow of an analytical microsecond dc pulsed glow discharge. This model uses several quantities as input data, such as electric field and potential, ion flux at the cathode, the fast argon ion and atom impact ionization rates, slow electron density, the electrical characterization of the pulse (voltage and current profiles) and temperature profile. These quantities were obtained by earlier Monte Carlo — fluid calculations for a pulsed discharge. Our goal is to study the behavior of the so-called Monte Carlo electrons (i.e., those electrons created at the cathode or by ionization collisions in the plasma which are followed by using the Monte Carlo model) from their origin to the moment when they are absorbed at the cell walls or when they have lost their energy by collisions (being transferred to the group of slow electrons) in the afterglow of the pulsed discharge. The thermalization of the electrons is a phenomenon where the electron-electron Coulomb collisions acquire a special importance. Indeed, in the afterglow the cross sections of the other electron reactions taken into account in the model are very low, because of the very low electron energy. We study the electron energy distributions at several times during and after the pulse and at several positions in the plasma cell, focusing on the thermalization and on the behavior of the electrons in the afterglow. Also, the time evolution of the rates of the various collision processes, the average electron energy, the densities of Monte Carlo and slow electrons and the ionization degree are investigated.     

Optimization of electron beam crosslinking of wire and cable insulation

The computer simulations based on Monte Carlo (MC) method and the ModeCEB software were carried out in connection with electron beam (EB) radiation set-up for crosslinking of electric wire and cable insulation. The theoretical predictions for absorbed dose distribution in irradiated electric insulation induced by scanned EB were compared to the experimental results of irradiation that was carried out in the experimental set-up based on ILU 6 electron accelerator with electron energy 0.5–2.0 MeV.The computer simulation of the dose distributions in two-sided irradiation system by a scanned electron beam in multilayer circular objects was performed for various process parameters, namely electric wire and cable geometry (thickness of insulation layers and copper wire diameter), type of polymer insulation, electron energy, energy spread and geometry of electron beam, electric wire and cable layout in irradiation zone. The geometry of electron beam distribution in the irradiation zone was measured using CTA and PVC foil dosimeters for available electron energy range. The temperature rise of the irradiated electric wire and irradiation homogeneity were evaluated for different experimental conditions to optimize technological process parameters. The results of computer simulation are consistent with the experimental data of dose distribution evaluated by gel-fraction measurements. Such conformity indicates that ModeCEB computer simulation is reliable and sufficient for optimization absorbed dose distribution in the multi-layer circular objects irradiated with scanned electron beams.     

Effects of physics change in Monte Carlo code on electron pencil beam dose distributions

Pencil beam algorithms used in computerized electron beam dose planning are usually described using the small angle multiple scattering theory. Alternatively, the pencil beams can be generated by Monte Carlo simulation of electron transport. In a previous work, the 4th version of the Electron Gamma Shower (EGS) Monte Carlo code was used to obtain dose distributions from monoenergetic electron pencil beam, with incident energy between 1 MeV and 50 MeV, interacting at the surface of a large cylindrical homogeneous water phantom. In 2000, a new version of this Monte Carlo code has been made available by the National Research Council of Canada (NRC), which includes various improvements in its electron-transport algorithms. In the present work, we were interested to see if the new physics in this version produces pencil beam dose distributions very different from those calculated with oldest one. The purpose of this study is to quantify as well as to understand these differences. We have compared a series of pencil beam dose distributions scored in cylindrical geometry, for electron energies between 1 MeV and 50 MeV calculated with two versions of the Electron Gamma Shower Monte Carlo Code. Data calculated and compared include isodose distributions, radial dose distributions and fractions of energy deposition. Our results for radial dose distributions show agreement within 10% between doses calculated by the two codes for voxels closer to the pencil beam central axis, while the differences are up to 30% for longer distances. For fractions of energy deposition, the results of the EGS4 are in good agreement (within 2%) with those calculated by EGSnrc at shallow depths for all energies, whereas a slightly worse agreement (15%) is observed at deeper distances. These differences may be mainly attributed to the different multiple scattering for electron transport adopted in these two codes and the inclusion of spin effect, which produces an increase of the effective range of electrons.     

Monte Carlo Model of Stricture Formation in Photodynamic Therapy of Normal Pig Esophagus

Photodynamic therapy (PDT) is FDA-approved for use in patients with Barrett's esophagus using porfimer sodium (2 mg per kg) and a recommended light dose of 130 J cm⁻¹ for high grade dysplasia. Despite uniform drug and light doses, the clinical outcome of PDT is variable. A significant number of PDT cases result in esophageal strictures, a side effect related to excessive energy absorption. The purpose of this project was to model esophageal stricture formation with a Monte Carlo simulation. An original multilayer Monte Carlo computer simulation was developed for esophageal PDT. Optical absorption and scattering coefficients were derived for mucosal and muscle layers of normal porcine esophagus. Porfimer sodium was added to each layer by increasing the absorption coefficient by the appropriate amount. A threshold-absorbed light dose was assumed to be required for stricture formation and ablation. The simulation predicted irreversible damage to the mucosa with a 160 J cm⁻¹ light dose and damage to the muscle layer with an additional 160 J cm⁻¹ light dose for a tissue porfimer sodium content of 3.5 mg kg⁻¹. The simulation accurately modeled photodynamic stricture formation in normal pig in vivo esophageal tissue. This preliminary work suggests that the absorbed light threshold for stricture formation may be between 2 and 4 J per gram of tissue.     

Slot scintillation detector modeling for digital radiography

A scanning slot X-ray digital detector operating at relatively high photon energies, mainly 70 kVp for dental radiography, has been designed and modeled. In this energy range, further improvements are required in order to attain optimal image quality for the lowest achievable mean absorbed dose to the patient. This paper concerns the investigation of the physics of imaging using a structured screen composed of a scintillating fiber optic (LKH-6/CCD system) and its performance by Monte Carlo (MC) simulation. The A_s factor, X-ray cross-talk, quantum noise effect on resultant values of Detective Quantum Efficiency (DQE) and Modulation Transfer Function (MTF) are discussed. The results of this paper should serve as a useful guide in the engineering design of future structured screen/CCD based systems.     

Super-Monte Carlo: A photon/electron dose calculation algorithm for radiotherapy

Super-Monte Carlo (SMC) is a method of dose calculation for radiotherapy which combines both analytical calculations and Monte Carlo electron transport. Analytical calculations are used where possible, such as the determination of photon interaction density, to decrease computation time. A Monte Carlo method is used for the electron transport in order to obtain high accuracy of results. To further speed computation, Monte Carlo is used once only, to form an electron track kernel (etk). The etk is a dataset containing the lengths and energy deposition of each step of a number of electron tracks. The etk is transported from each incident particle interaction site, from which the dose is calculated. Dose distributions calculated in heterogeneous media show SMC results similar to those of Monte Carlo. For the same statistical uncertainty, SMC takes an order of magnitude less computation time than a full Monte Carlo simulation. SMC has only been implemented for photons and electrons, however the same basic method could be used for the transport of other particles. Current development includes the optimisation of the etks and the code in order to decrease computation time, and also the inclusion of SMC onto a clinical planning system.     

Advantages of mesh tallying in MCNPX for 3D dose calculations in radiotherapy

The energy deposition mesh tally option of MCNPX Monte Carlo code is very useful for 3-Dimentional (3D) dose calculations. In this study, the 3D dose calculation was done for CT-based Monte Carlo treatment planning in which the energy deposition mesh tally were superimposed on merged voxel model. The results were compared with those of obtained from the common energy deposition (*F8) tally method for all cells of non-merged voxel model. The results of these two tallies and their respective computational times are compared, and the advantages of the proposed method are discussed. For this purpose, a graphical user interface (GUI) application was developed for reading CT slice data of patient, creating voxelized model of patient, optionally merging adjacent cells with the same material to reduce the total number of cells, reading beam configuration from commercial treatment planning system transferred in DICOM-RT format, and showing the isodose distribution on the CT images. To compare the results of Monte Carlo calculated and TiGRT planning system (LinaTech LLC, USA), treatment head of the Siemens ONCOR Impression accelerator was also simulated and the phase-space data on the scoring plane just above the Y-jaws was created and used. The results for a real prostate intensity-modulated radiation therapy (IMRT) plan showed that the proposed method was fivefold faster while the precision was almost the same.     

O—酰基—α—酮肟光分解反应的Monte Carlo处理

本文首次对光化学反应体系用Monte Carlo方法进行模拟处理。通过5个O-酰基-α-酮肟光分解反应的Monte Carlo模拟,可避免解析解中由于对吸收光强须采用一级近似求解动力学微分方程组,而造成拟合反应在后期产生与实验结果的偏差。     

The operator dose assessment of landmine detection systems using the neutron backscattering method

Personal safety issue is one of the major limitations in landmine detection by nuclear techniques. In this paper, the dose rate received by the operator of different hand-held landmine detection systems using the neutron backscattering method was investigated. Monte Carlo method was used to simulate a computational model of the body, instruments with different shielding configurations and soil with various moisture contents. The effective dose received by the different organs of the operator has been estimated by using two different approximation methods including average equivalent dose and dose equivalent assessment. The results obtained by these two methods were compared. The common results showed that the operator dose rate depends on the facility shielding, soil moisture level, and source-to-operator distance. Also, although the absorbed dose received by most organs generally decreases as a function of source to operator distance, for some organs such as kidneys and lungs it firstly increases when source-to-operator distance increases up to 0.7–1.2 m and after that the absorbed dose decreases. Furthermore, the results showed that the effective dose received by the operator has its maximum value when the source to operator distance is 0.60 m.     

Monte Carlo treatment planning for photon and electron beams

During the last few decades, accuracy in photon and electron radiotherapy has increased substantially. This is partly due to enhanced linear accelerator technology, providing more flexibility in field definition (e.g. the usage of computer-controlled dynamic multileaf collimators), which led to intensity modulated radiotherapy (IMRT). Important improvements have also been made in the treatment planning process, more specifically in the dose calculations. Originally, dose calculations relied heavily on analytic, semi-analytic and empirical algorithms. The more accurate convolution/superposition codes use pre-calculated Monte Carlo dose “kernels” partly accounting for tissue density heterogeneities. It is generally recognized that the Monte Carlo method is able to increase accuracy even further. Since the second half of the 1990s, several Monte Carlo dose engines for radiotherapy treatment planning have been introduced. To enable the use of a Monte Carlo treatment planning (MCTP) dose engine in clinical circumstances, approximations have been introduced to limit the calculation time. In this paper, the literature on MCTP is reviewed, focussing on patient modeling, approximations in linear accelerator modeling and variance reduction techniques. An overview of published comparisons between MC dose engines and conventional dose calculations is provided for phantom studies and clinical examples, evaluating the added value of MCTP in the clinic. An overview of existing Monte Carlo dose engines and commercial MCTP systems is presented and some specific issues concerning the commissioning of a MCTP system are discussed.     

An increase of utilization efficiency of X-ray beam

The features of absorbed dose field formation in objects irradiated with scanned X-ray beams at double—and four-sided irradiation were investigated both analytically and by Monte Carlo methods. An analytical approach uses an angular/spectrum X-ray characteristics calculated with PENELOPE, JEANT 4 and ModeXR codes. It was shown that the special angular orientation of electron beam incidence on the X-ray converter leads to X-ray dose smoothing on the surface of the irradiated object. At the same time, a double-sided irradiation can provide high X-ray beam efficiency at dose uniformity ration (DUR) <1.5 for sizeable object thickness. At four-sided irradiation, the angular orientation of electron beam incidence on the X-ray converter should be changed so as to focus the electrons to the center of the converter. At this mode X-ray beam efficiency is more than 60%.     

Radial dose profiles for pencil electron beams

The relationship between the Boltzmann and Fermi-Eyges-Yang equations governing electron transport is examined. Radial dose profiles for a pencil beam obtained by numerical solution of the Boltzmann equation in the small angle approximation are compared with both the Gaussian approximation and with Monte Carlo simulations for a carbon medium. For energies ranging from 5 to 20 MeV and penetration depths up to 75% of the range the numerical results are within 10% of the Monte Carlo results for the radial distance encompassing 63% of the energy deposition.     

Monte Carlo sampling algorithm for searching a scale-transformed energy space of polypeptides

A Monte Carlo sampling algorithm for searching a scale-transformed conformational energy space of polypeptides is presented. This algorithm is based on the assumption that energy barriers can be overcome by a uniform sampling of the logarithmically transformed energy space. This algorithm is tested with Met-enkephalin. For comparison, the entropy sampling Monte Carlo (ESMC) simulation is performed. First, the global minimum is easily found by the optimization of a scale-transformed energy space. With a new Monte Carlo sampling, energy barriers of 3000 kcal/mol are frequently overcome, and low-energy conformations are sampled more efficiently than with ESMC simulations. Several thermodynamic quantities are calculated with good accuracy.     

X-ray treatment at 5 MeV and above

There is agreement in the scientific community that X-ray treatment of food at 7.5 MeV can be safe. Possible process improvements for treating at higher than 5 MeV X-rays have been re-visited. Monte Carlo methods have been applied to simulate the X-ray conversion process and to calculate dose distributions in homogeneous phantoms. Experimental data obtained using X-rays produced with a Rhodotron TT200 at 5 and 10 MeV verifies a representative set of data which is calculated with the presented method.

With this qualified Monte Carlo tool, calculations at 7.5 MeV incident electron energy were performed. The analysis gives special attention to higher photon yield, improved product penetration, as well as surface and edge effects.     

Photon-electron energy deposition in CANDU reactor channels: Simulation and modelling

Simulation of photon-electron transport in a CANDU reactor fuel channel using the Monte Carlo method for calculating the energy deposition in the coolant is studied. The geometry of the CANDU fuel channel is very complex so methods that make such simulations more practicable, without adversely affecting the results, are introduced. In this regard, the use of simplifying assumptions and simplified geometrical models on the performance of two different Monte Carlo codes has been compared. An ETRAN-based code (SANDYL), and the code EGS4 produced comparable results, although the former performs faster in accounting for low energy electrons. A simplified computational model is also introduced. This model is based on decoupling photon-electron transport simulations by the use of electron-energy-transfer functions. The results obtained using the model are successfully validated using the EGS4 and SANDYL codes. A significant computational speedup (about a factor of seven compared to Monte Carlo simulations) is achieved with this model.     

Monte Carlo simulation of electrodeposition of copper: a multistep free energy calculation

Electrodeposition of copper (Cu) involves length scales of a micrometer or even less. Several theoretical techniques such as continuum Monte Carlo, kinetic Monte Carlo (KMC), and molecular dynamics have been used for simulating this problem. However the multiphenomena characteristics of the problem pose a challenge for an efficient simulation algorithm. Traditional KMC methods are slow, especially when modeling surface diffusion with large number of particles and frequent particle jumps. Parameter estimation involving thousands of KMC runs is very time-consuming. Thus a less time-consuming and novel multistep continuum Monte Carlo simulation is carried out to evaluate the step wise free energy change in the process of electrochemical copper deposition. The procedure involves separate Monte Carlo codes employing different random number criterion (using hydrated radii, bare radii, hydration number of the species, redox potentials, etc.) to obtain the number of species (CuCl(2) or CuSO(4) or Cu as the case may be) and in turn the free energy. The effect of concentration of electrolyte, influence of electric field and presence of chloride ions on the free energy change for the processes is studied. The rate determining step for the process of electrodeposition of copper from CuCl(2) and CuSO(4) is also determined.     

Free energies and phase equilibria of chain molecules

In this Article, a review is presented of recent developments in Monte Carlo simulations of chain molecules. The Rosenbluth chain insertion technique is used to calculate the free energy of the chain molecules. Furthermore, this insertion method is used to generate biased Monte Carlo moves. It is shown that this bias can be removed by adjusting the acceptance rules such that configurations are generated with their correct Boltzmann weight. This configurational-bias Monte Carlo method can be combined with the Gibbs-ensemble technique which results in an efficient method to simulate phase equilibria of chain molecules.     

Quantum Monte Carlo simulations of selected ammonia clusters (n = 2-5): isotope effects on the ground state of typical hydrogen bonded systems

Variational Monte Carlo, diffusion Monte Carlo, and stereographic projection path integral simulations are performed on eight selected species from the (NH(3))(n), (ND(3))(n), (NH(2)D)(n), and (NH(3))(n-1)(ND(3)) clusters. Each monomer is treated as a rigid body with the rotation spaces mapped by the stereographic projection coordinates. We compare the energy obtained from path integral simulations at several low temperatures with those obtained by diffusion Monte Carlo, for two dimers, and we find that at 4 K, the fully deuterated dimer energy is in excellent agreement with the ground state energy of the same. The ground state wavefunction for the (NH(3))(2-5) clusters is predominantly localized in the global minimum of the potential energy. In all simulations of mixed isotopic substitutions, we find that the heavier isotope is almost exclusively the participant in the hydrogen bond.