Monte Carlo calculation of Varian 2300C/D Linac photon beam characteristics: a comparison between MCNP4C, GEANT3 and measurements

Different codes are used for Monte Carlo (MC) calculations in radiation therapy. In this research, MCNP4C and GEANT3 codes have been compared in calculations of dosimetric characteristics of Varian Clinac 2300C/D. The parameters of influence in the differences seen in dosimetric features were discussed. This study emphasizes that both MCNP4C and GEANT3 MC can be used in radiation therapy computations and their differences in photon spectra calculations have a negligible effect on percentage depth dose computations in radiation therapy.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

Computer algebra for x-ray spectral reconstruction between 6 and 25 MV

A previously presented algorithm for the reconstruction of bremsstrahlung spectra from transmission data has been implemented into MATHEMATICA. Spectra vectorial algebra has been used to solve the matrix system A * F = T. The new implementation has been tested by reconstructing photon spectra from transmission data acquired in narrow beam conditions, for nominal energies of 6, 15, and 25 MV. The results were in excellent agreement with the original calculations. Our implementation has the advantage to be based on a well-tested mathematical kernel. Furthermore it offers a comfortable user interface.

Vibronic Spectra of Jet-Cooled 2-Aminopurine center dot H2O Clusters Studied by UV Resonant Two-Photon Ionization Spectroscopy and Quantum Chemical Calculations

For understanding the major- and minor-groove hydration patterns of DNAs and RNAs, it is important to understand the local solvation of individual nucleobases at the molecular level. We have investigated the 2-aminopurine center dot H2O. monohydrate by two-color resonant two-photon ionization and UV/UV hole-burning spectroscopies, which reveal two isomers, denoted A and B. The electronic spectral shift delta nu of the S-1 <- S-0 transition relative to bare 9H-2-aminopurine (9H-2AP) is small for isomer A (-70 cm(-1)), while that of isomer B is much larger (delta nu = 889 cm(-1)). B3LYP geometry optimizations with the TZVP basis set predict four cluster isomers, of which three are doubly H-bonded, with H2O acting as an acceptor to a N-H or -NH2 group and as a donor to either of the pyrimidine N sites. The "sugar-edge" isomer A is calculated to be the most stable form with binding energy D-e = 56.4 kJ/mol. Isomers B and C are H-bonded between the -NH2 group and pyrimidine moieties and are 2.5 and 6.9 kJ/mol less stable, respectively. Time-dependent (TD) B3LYP/TZVP calculations predict the adiabatic energies of the lowest (1)pi pi* states of A and B in excellent agreement with the observed 0(0)(0) bands; also, the relative intensities of the A and B origin bands agree well with the calculated S-0 state relative energies. This allows unequivocal identification of the isomers. The R2PI spectra of 9H-2AP and of isomer A exhibit intense low-frequency out-of-plane overtone and combination bands, which is interpreted as a coupling of the optically excited (1)pi pi* state to the lower-lying (1)n pi* dark state. In contrast, these overtone and combination bands are much weaker for isomer B, implying that the (1)pi pi* state of B is planar and decoupled from the (1)n pi* state. These observations agree with the calculations, which predict the (1)n pi* above the (1)pi pi* state for isomer B but below the (1)pi pi* for both 9H-2AP and isomer A.

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.

Watson–Crick and Sugar-Edge Base Pairing of Cytosine in the Gas Phase: UV and Infrared Spectra of Cytosine·2-Pyridone

While keto-amino cytosine is the dominant species in aqueous solution, spectroscopic studies in molecular beams and in noble gas matrices show that other cytosine tautomers prevail in apolar environments. Each of these offers two or three H-bonding sites (Watson–Crick, wobble, sugar-edge). The mass- and isomer-specific S1 ← S0 vibronic spectra of cytosine·2-pyridone (Cyt·2PY) and 1-methylcytosine·2PY are measured using UV laser resonant two-photon ionization (R2PI), UV/UV depletion, and IR depletion spectroscopy. The UV spectra of the Watson–Crick and sugar-edge isomers of Cyt·2PY are separated using UV/UV spectral hole-burning. Five different isomers of Cyt·2PY are observed in a supersonic beam. We show that the Watson–Crick and sugar-edge dimers of keto-amino cytosine with 2PY are the most abundant in the beam, although keto-amino-cytosine is only the third most abundant tautomer in the gas phase. We identify the different isomers by combining three different diagnostic tools: (1) methylation of the cytosine N1–H group prevents formation of both the sugar-edge and wobble isomers and gives the Watson–Crick isomer exclusively. (2) The calculated ground state binding and dissociation energies, relative gas-phase abundances, excitation and the ionization energies are in agreement with the assignment of the dominant Cyt·2PY isomers to the Watson–Crick and sugar-edge complexes of keto-amino cytosine. (3) The comparison of calculated ground state vibrational frequencies to the experimental IR spectra in the carbonyl stretch and NH/OH/CH stretch ranges strengthen this identification.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.