Calculation of dose kernels for radiotherapy applications using the GEANT 4 Monte Carlo toolkit

Dose kernels may be used to calculate dose distributions in radiotherapy (as described by Ahnesjo et al., 1999). Their calculation requires use of Monte Carlo methods, usually by forcing interactions to occur at a point. The Geant4 Monte Carlo toolkit provides a capability to force interactions to occur in a particular volume. We have modified this capability and created a Geant4 application to calculate dose kernels in cartesian, cylindrical, and spherical scoring systems. The simulation considers monoenergetic photons incident at the origin of a 3 m x 3 x 9 3 m water volume. Photons interact via compton, photo-electric, pair production, and rayleigh scattering. By default, Geant4 models photon interactions by sampling a physical interaction length (PIL) for each process. The process returning the smallest PIL is then considered to occur. In order to force the interaction to occur within a given length, L_FIL, we scale each PIL according to the formula: PIL_forced = L_FIL 9 (1 - exp(-PIL/PILo)) where PILo is a constant. This ensures that the process occurs within L_FIL, whilst correctly modelling the relative probability of each process. Dose kernels were produced for an incident photon energy of 0.1, 1.0, and 10.0 MeV. In order to benchmark the code, dose kernels were also calculated using the EGSnrc Edknrc user code. Identical scoring systems were used; namely, the collapsed cone approach of the Edknrc code. Relative dose difference images were then produced. Preliminary results demonstrate the ability of the Geant4 application to reproduce the shape of the dose kernels; median relative dose differences of 12.6, 5.75, and 12.6 % were found for an incident photon energy of 0.1, 1.0, and 10.0 MeV respectively.

The influence of Monte Carlo source parameters on detector design and dose perturbation in small field dosimetry

To obtain accurate Monte Carlo simulations of small radiation fields, it is important model the initial source parameters (electron energy and spot size) accurately. However recent studies have shown that small field dosimetry correction factors are insensitive to these parameters. The aim of this work is to extend this concept to test if these parameters affect dose perturbations in general, which is important for detector design and calculating perturbation correction factors. The EGSnrc C++ user code cavity was used for all simulations. Varying amounts of air between 0 and 2 mm were deliberately introduced upstream to a diode and the dose perturbation caused by the air was quantified. These simulations were then repeated using a range of initial electron energies (5.5 to 7.0 MeV) and electron spot sizes (0.7 to 2.2 FWHM). The resultant dose perturbations were large. For example 2 mm of air caused a dose reduction of up to 31% when simulated with a 6 mm field size. However these values did not vary by more than 2 % when simulated across the full range of source parameters tested. If a detector is modified by the introduction of air, one can be confident that the response of the detector will be the same across all similar linear accelerators and the Monte Carlo modelling of each machine is not required.

The lowest singlet and triplet states of the oxyallyl diradical

Small S-T splitting : The photoelectron spectrum of the oxyallyl radical anion (see picture) reveals that the electronic ground state of oxyallyl is singlet, and the lowest triplet state is separated from the singlet state by only (55 ± 2) meV in adiabatic energy.

Surface dosimetry for very small X-ray beams

Introduction The dose to skin surface is an important factor for many radiotherapy treatment techniques. It is known that TPS predicted surface doses can be significantly different from actual ICRP skin doses as defined at 70 lm. A number of methods have been implemented for the accurate determination of surface dose including use of specific dosimeters such as TLDs and radiochromic film as well as Monte Carlo calculations. Stereotactic radiosurgery involves delivering very high doses per treatment fraction using small X-ray fields. To date, there has been limited data on surface doses for these very small field sizes. The purpose of this work is to evaluate surface doses by both measurements and Monte Carlo calculations for very small field sizes. Methods All measurements were performed on a Novalis Tx linear accelerator which has a 6 MV SRS X-ray beam mode which uses a specially thin flattening filter. Beam collimation was achieved by circular cones with apertures that gave field sizes ranging from 4 to 30 mm at the isocentre. The relative surface doses were measured using Gafchromic EBT3 film which has the active layer at a depth similar to the ICRP skin dose depth. Monte Carlo calculations were performed using the BEAMnrc/EGSnrc Monte Carlo codes (V4 r225). The specifications of the linear accelerator, including the collimator, were provided by the manufacturer. Optimisation of the incident X-ray beam was achieved by an iterative adjustment of the energy, spatial distribution and radial spread of the incident electron beam striking the target. The energy cutoff parameters were PCUT = 0.01 MeV and ECUT = 0.700 - MeV. Directional bremsstrahlung splitting was switched on for all BEAMnrc calculations. Relative surface doses were determined in a layer defined in a water phantom of the same thickness and depth as compared to the active later in the film. Results Measured surface doses using the EBT3 film varied between 13 and 16 % for the different cones with an uncertainty of 3 %. Monte Carlo calculated surface doses were in agreement to better than 2 % to the measured doses for all the treatment cones. Discussion and conclusions This work has shown the consistency of surface dose measurements using EBT3 film with Monte Carlo predicted values within the uncertainty of the measurements. As such, EBT3 film is recommended for in vivo surface dose measurements.

Photoelectron spectroscopic study of the oxyallyl diradical

The photoelectron spectrum of the oxyallyl (OXA) radical anion has been measured. The radical anion has been generated in the reaction of the atomic oxygen radical anion (O center dot-) with acetone. Three low-lying electronic states of OXA have been observed in the spectrum. Electronic structure calculations have been performed for the triplet states (B-3(2) and B-3(1)) of OXA and the ground doublet state ((2)A(2)) of the radical anion using density, functional theory (DFT). Spectral simulations have been carried out for the triplet statics based on the results of the DFT calculations. The simulation identifies a vibrational progression of the CCC bending mode of the B-3(2) state of OXA in the lower electron binding energy (eBE) portion of the spectrum. On top of the B-3(2) feature, however, the experimental spectrum exhibits additional photoelectron peaks whose angular distribution is distinct from that for the vibronic peaks of the B-3(2) state. Complete active space self-consistent field (CASSCF) method and second-order perturbation theory based on the CASSCF wave function (CASPT2) have been employed to study the lowest singlet state ((1)A(1)) of OXA. The simulation based on the results of these electronic structure calculations establishes that the overlapping peaks represent the vibrational ground level of the (1)A(1) state and its vibrational progression of the CO stretching mode. The A, state is the lowest electronic state of,OXA, and the electron affinity (EA) of OXA is 1.940 +/- 0.010 eV. The B-3(2) state is the first excited state with an electronic term energy of 55 +/- 2 meV. The widths of the vibronic peaks of the (X) over tilde (1)A(1) state are much broader than those of the (a) over tilde B-3(2) state, implying that the (1)A(1) state is indeed a transition state. The CASSCF and CASPT2 calculations suggest that the (1)A(1) state is at a potential maximum along the nuclear coordinate representing disrotatory motion of the two methylene groups, which leads to three-membered-ring formation, i.e., cydopropanone. The simulation of (b) over tilde B-3(1) OXA reproduces the higher eBE portion of the spectrum very well. The term energy of the B-3(1) state is 0.883 +/- 0.012 eV. Photoelectron spectroscopic measurements have also been conducted for the other ion products of the O center dot- reaction with acetone. The photoelectron imaging spectrum of the acetylcarbene (AC) radical anion exhibits a broad, structureless feature, which is assigned to the (X) over tilde (3)A '' state of AC. The ground ((2)A '') and first excited ((2)A') states of the 1-methylvinoxy (1-MVO) radical have been observed in the photoelectron spectrum of the 1-MVO ion, and their vibronic structure has been analyzed.

Gold nanoresistors with near-constant resistivity in the cryogenic-to-room temperature range

Using a multiple plasma deposition-annealing (MDA) technique, we have fabricated an Au nanoisland-based thin film nanoresistor with a very low temperature coefficient of electrical resistivity in a cryogenic-to-room temperature range of 10 to 300 K. The nanoislanded gold film was deposited on a SiO2/Si wafer (500 nm SiO2 thickness) between two 300 nm thick Au electrodes which were separated by 100 m. A sophisticated selection of the thickness of the nanoislanded gold film, the annealing temperature, as well as the number of deposition/annealing cycles resulted in the fabrication of a nanoresistor with a temperature coefficient of electrical resistivity of 2.1 × 10-3 K-1 and the resistivity deviation not exceeding 2% in a cryogenic-to-room temperature range. We have found that the constant resistivity regime of the nanoisland-based thin film nanoresistor corresponds to a minimized nanoisland activation energy (approximately 0.3 meV). This energy can be minimized by reducing the nearest neighbor distance and increasing the size of the Au nanoislands in the optimized nanoresistor structure. It is shown that the constant resistivity nanoresistor operates in the regime where the thermally activated electron tunneling is compensated by the negative temperature dependence of the metallic-type conductivity of nanoislands. Our results are relevant to the development of commercially viable methods of nanoresistor production for various nanoelectronics-based devices. The proposed MDA technique also provides the opportunity to fabricate large arrays of metallic nanoparticles with controllable size, shapes and inter-nanoparticle gaps.

Synthesis and study of conductivity behaviour of blended conducting polymer films irradiated with swift heavy ions of silicon

In recent times, blended polymers have shown a lot of promise in terms of easy processability in different shapes and forms. In the present work, polyaniline emeraldine base (PANi-EB) was doped with camphor sulfonic acid (CSA) and combined with the conducting polymer polyfluorene (PF) as well as the insulating polymer polyvinyl chloride (PVC) to synthesize CSA doped PANi-PF and PANi-PVC blended polymers. It is well known that PANi when doped with CSA becomes highly conducting. However, its poor mechanical properties, such as low tensile, compressive, and flexural strength render PANi a non-ideal material to be processed for its various practical applications, such as electromagnetic shielding, anti-corrosion shielding, photolithography and microelectronic devices etc. Thus the search for polymers which are easily processable and are capable of showing high conductivity still continues. PANi-PVC blend was prepared, which showed low conductivity which is limiting factor for certain applications. Therefore, another processable polymer PF was chosen as conducting matrix. Conducting PF can be easily processed into various shapes and forms. Therefore, a blend mixture was prepared by using PANi and PF through the use of CSA as a counter ion which forms a "bridge" between the two polymeric components of the inter-polymer complex. Two blended polymers have been synthesized and investigated for their conductivity behaviour. It was observed that the blended film of CSA doped PANi-PVC showed a room temperature electrical conductivity of 2.8 × 10-7 S/cm where as the blended film made by CSA doped PANi with conducting polymer PF showed a room temperature conductivity of 1.3 × 10-5 S/cm. Blended films were irradiated with 100 MeV silicon ions with a view to increase their conductivity with a fluence ranging from 1011 ions to 1013 per cm2 from 15 UD Pelletron accelerator at NSC, New Delhi.

Group 14 element-based non-centrosymmetric quantum spin hall insulators with large bulk gap

To date, a number of two-dimensional (2D) topological insulators (TIs) have been realized in Group 14 elemental honeycomb lattices, but all are inversionsymmetric. Here, based on first-principles calculations, we predict a new family of 2D inversion-asymmetric TIs with sizeable bulk gaps from 105 meV to 284 meV, in X2–GeSn (X = H, F, Cl, Br, I) monolayers, making them in principle suitable for room-temperature applications. The nontrivial topological characteristics of inverted band orders are identified in pristine X2–GeSn with X = (F, Cl, Br, I), whereas H2–GeSn undergoes a nontrivial band inversion at 8% lattice expansion. Topologically protected edge states are identified in X2–GeSn with X = (F, Cl, Br, I), as well as in strained H2–GeSn. More importantly, the edges of these systems, which exhibit single-Dirac-cone characteristics located exactly in the middle of their bulk band gaps, are ideal for dissipationless transport. Thus, Group 14 elemental honeycomb lattices provide a fascinating playground for the manipulation of quantum states.

On the interaction between radon progeny and particles generated by electronic and traditional cigarettes

During their entire lives, people are exposed to the pollutants present in indoor air. Recently, Electronic Nicotine Delivery Systems, mainly known as electronic cigarettes, have been widely commercialized: they deliver particles into the lungs of the users but a “second-hand smoke” has yet to be associated to this indoor source. On the other hand, the naturally-occurring radioactive gas, i.e. radon, represents a significant risk for lung cancer, and the cumulative action of these two agents could be worse than the agents separately would. In order to deepen the interaction between radon progeny and second-hand aerosol from different types of cigarettes, a designed experimental study was carried out by generating aerosol from e-cigarette vaping as well as from second-hand traditional smoke inside a walk-in radon chamber at the National Institute of Ionizing Radiation Metrology (INMRI) of Italy. In this chamber, the radon present in air comes naturally from the floor and ambient conditions are controlled. To characterize the sidestream smoke emitted by cigarettes, condensation particle counters and scanning mobility particle sizer were used. Radon concentration in the air was measured through an Alphaguard ionization chamber, whereas the measurement of radon decay product in the air was performed with the Tracelab BWLM Plus-2S Radon daughter Monitor. It was found an increase of the Potential Alpha-Energy Concentration (PAEC) due to the radon decay products attached to aerosol for higher particle number concentrations. This varied from 7.47 ± 0.34 MeV L−1 to 12.6 ± 0.26 MeV L−1 (69%) for the e-cigarette. In the case of traditional cigarette and at the same radon concentration, the increase was from 14.1 ± 0.43 MeV L−1 to 18.6 ± 0.19 MeV L−1 (31%). The equilibrium factor increases, varying from 23.4% ± 1.11% to 29.5% ± 0.26% and from 30.9% ± 1.0% to 38.1 ± 0.88 for the e-cigarette and traditional cigarette, respectively. These growths still continue for long time after the combustion, by increasing the exposure risk.

Use of electronic portal imaging devices for electron treatment verification

This study aims to help broaden the use of electronic portal imaging devices (EPIDs) for pre-treatment patient positioning verification, from photon-beam radiotherapy to photon- and electron-beam radiotherapy, by proposing and testing a method for acquiring clinicallyuseful EPID images of patient anatomy using electron beams, with a view to enabling and encouraging further research in this area. EPID images used in this study were acquired using all available beams from a linac configured to deliver electron beams with nominal energies of 6, 9, 12, 16 and 20 MeV, as well as photon beams with nominal energies of 6 and 10 MV. A widely-available heterogeneous, approximately-humanoid, thorax phantom was used, to provide an indication of the contrast and noise produced when imaging different types of tissue with comparatively realistic thicknesses. The acquired images were automatically calibrated, corrected for the effects of variations in the sensitivity of individual photodiodes, using a flood field image. For electron beam imaging, flood field EPID calibration images were acquired with and without the placement of blocks of water-equivalent plastic (with thicknesses approximately equal to the practical range of electrons in the plastic) placed upstream of the EPID, to filter out the primary electron beam, leaving only the bremsstrahlung photon signal. While the electron beam images acquired using a standard (unfiltered) flood field calibration were observed to be noisy and difficult to interpret, the electron beam images acquired using the filtered flood field calibration showed tissues and bony anatomy with levels of contrast and noise that were similar to the contrast and noise levels seen in the clinically acceptable photon beam EPID images. The best electron beam imaging results (highest contrast, signal-to-noise and contrast-to-noise ratios) were achieved when the images were acquired using the higher energy electron beams (16 and 20 MeV) when the EPID was calibrated using an intermediate (12 MeV) electron beam energy. These results demonstrate the feasibility of acquiring clinically-useful EPID images of patient anatomy using electron beams and suggest important avenues for future investigation, thus enabling and encouraging further research in this area. There is manifest potential for the EPID imaging method proposed in this work to lead to the clinical use of electron beam imaging for geometric verification of electron treatments in the future.