Monte Carlo dose calculation improvements for low energy electron beams using eMC

The electron Monte Carlo (eMC) dose calculation algorithm in Eclipse (Varian Medical Systems) is based on the macro MC method and is able to predict dose distributions for high energy electron beams with high accuracy. However, there are limitations for low energy electron beams. This work aims to improve the accuracy of the dose calculation using eMC for 4 and 6 MeV electron beams of Varian linear accelerators. Improvements implemented into the eMC include (1) improved determination of the initial electron energy spectrum by increased resolution of mono-energetic depth dose curves used during beam configuration; (2) inclusion of all the scrapers of the applicator in the beam model; (3) reduction of the maximum size of the sphere to be selected within the macro MC transport when the energy of the incident electron is below certain thresholds. The impact of these changes in eMC is investigated by comparing calculated dose distributions for 4 and 6 MeV electron beams at source to surface distance (SSD) of 100 and 110 cm with applicators ranging from 6 x 6 to 25 x 25 cm(2) of a Varian Clinac 2300C/D with the corresponding measurements. Dose differences between calculated and measured absolute depth dose curves are reduced from 6% to less than 1.5% for both energies and all applicators considered at SSD of 100 cm. Using the original eMC implementation, absolute dose profiles at depths of 1 cm, d(max) and R50 in water lead to dose differences of up to 8% for applicators larger than 15 x 15 cm(2) at SSD 100 cm. Those differences are now reduced to less than 2% for all dose profiles investigated when the improved version of eMC is used. At SSD of 110 cm the dose difference for the original eMC version is even more pronounced and can be larger than 10%. Those differences are reduced to within 2% or 2 mm with the improved version of eMC. In this work several enhancements were made in the eMC algorithm leading to significant improvements in the accuracy of the dose calculation for 4 and 6 MeV electron beams of Varian linear accelerators.

Low Energy Neutral Atoms from the Heliosheath

In the heliosheath beyond the termination shock, low energy (<0.5 keV) neutral atoms are created by charge exchange with interstellar neutrals. Detecting these neutrals from Earth's orbit is difficult because their flux is reduced substantially by ionization losses as they propagate from about 100 to 1 AU and because there are a variety of other signals and backgrounds that compete with this weak signal. Observations from IBEX-Lo and -Hi from two opposing vantage points in Earth's orbit established a lower energy limit of about 0.1 keV on measurements of energetic neutral atoms (ENAs) from the heliosphere and the form of the energy spectrum from about 0.1 to 6 keV in two directions in the sky. Below 0.1 keV, the detailed ENA spectrum is not known, and IBEX provides only upper limits on the fluxes. However, using some assumptions and taking constraints on the spectrum into account, we find indications that the spectrum turns over at an energy between 0.1 and 0.2 keV.

Measurement of the muon neutrino inclusive charged-current cross section in the energy range of 1–3 GeV with the T2K INGRID detector

We report a measurement of the νµ-nucleus inclusive charged current cross section (=σ cc) on ironusing data from exposed to the J-PARC neutrino beam. The detector consists of 14 modules in total, which are spread over a range of off-axis angles from 0◦ to 1.1◦. The variation in the neutrino energy spectrum as a function of the off-axis angle, combined with event topology information, is used to calculate this cross section as a function of neutrino energy. The cross section is measured to be σcc(1.1 GeV) = 1.10±0.15 (10^−38cm^2/nucleon), σcc(2.0 GeV) = 2.07±0.27 (10^−38cm^2/nucleon), and σcc(3.3 GeV) = 2.29 ± 0.45 (10^−38cm^2/nucleon), at energies of 1.1, 2.0, and 3.3 GeV, respectively. These results are consistent with the cross section calculated by the neutrino interaction generators currently used by T2K. More importantly, the method described here opens up a new way to determine the energy dependence of neutrino-nucleus cross sections.

Photon-beam subsource sensitivity to the initial electron-beam parameters

One limitation to the widespread implementation of Monte Carlo (MC) patient dose-calculation algorithms for radiotherapy is the lack of a general and accurate source model of the accelerator radiation source. Our aim in this work is to investigate the sensitivity of the photon-beam subsource distributions in a MC source model (with target, primary collimator, and flattening filter photon subsources and an electron subsource) for 6- and 18-MV photon beams when the energy and radial distributions of initial electrons striking a linac target change. For this purpose, phase-space data (PSD) was calculated for various mean electron energies striking the target, various normally distributed electron energy spread, and various normally distributed electron radial intensity distributions. All PSD was analyzed in terms of energy, fluence, and energy fluence distributions, which were compared between the different parameter sets. The energy spread was found to have a negligible influence on the subsource distributions. The mean energy and radial intensity significantly changed the target subsource distribution shapes and intensities. For the primary collimator and flattening filter subsources, the distribution shapes of the fluence and energy fluence changed little for different mean electron energies striking the target, however, their relative intensity compared with the target subsource change, which can be accounted for by a scaling factor. This study indicates that adjustments to MC source models can likely be limited to adjusting the target subsource in conjunction with scaling the relative intensity and energy spectrum of the primary collimator, flattening filter, and electron subsources when the energy and radial distributions of the initial electron-beam change.

Monte Carlo simulation of a dynamic MLC based on a multiple source model

Detailed knowledge of the characteristics of the radiation field shaped by a multileaf collimator (MLC) is essential in intensity modulated radiotherapy (IMRT). A previously developed multiple source model (MSM) for a 6 MV beam was extended to a 15 MV beam and supplemented with an accurate model of an 80-leaf dynamic MLC. Using the supplemented MSM and the MC code GEANT, lateral dose distributions were calculated in a water phantom and a portal water phantom. A field which is normally used for the validation of the step and shoot technique and a field from a realistic IMRT treatment plan delivered with dynamic MLC are investigated. To assess possible spectral changes caused by the modulation of beam intensity by an MLC, the energy spectra in five portal planes were calculated for moving slits of different widths. The extension of the MSM to 15 MV was validated by analysing energy fluences, depth doses and dose profiles. In addition, the MC-calculated primary energy spectrum was verified with an energy spectrum which was reconstructed from transmission measurements. MC-calculated dose profiles using the MSM for the step and shoot case and for the dynamic MLC case are in very good agreement with the measured data from film dosimetry. The investigation of a 13 cm wide field shows an increase in mean photon energy of up to 16% for the 0.25 cm slit compared to the open beam for 6 MV and of up to 6% for 15 MV, respectively. In conclusion, the MSM supplemented with the dynamic MLC has proven to be a powerful tool for investigational and benchmarking purposes or even for dose calculations in IMRT.

A multiple source model for 6 MV photon beam dose calculations using Monte Carlo

A multiple source model (MSM) for the 6 MV beam of a Varian Clinac 2300 C/D was developed by simulating radiation transport through the accelerator head for a set of square fields using the GEANT Monte Carlo (MC) code. The corresponding phase space (PS) data enabled the characterization of 12 sources representing the main components of the beam defining system. By parametrizing the source characteristics and by evaluating the dependence of the parameters on field size, it was possible to extend the validity of the model to arbitrary rectangular fields which include the central 3 x 3 cm2 field without additional precalculated PS data. Finally, a sampling procedure was developed in order to reproduce the PS data. To validate the MSM, the fluence, energy fluence and mean energy distributions determined from the original and the reproduced PS data were compared and showed very good agreement. In addition, the MC calculated primary energy spectrum was verified by an energy spectrum derived from transmission measurements. Comparisons of MC calculated depth dose curves and profiles, using original and PS data reproduced by the MSM, agree within 1% and 1 mm. Deviations from measured dose distributions are within 1.5% and 1 mm. However, the real beam leads to some larger deviations outside the geometrical beam area for large fields. Calculated output factors in 10 cm water depth agree within 1.5% with experimentally determined data. In conclusion, the MSM produces accurate PS data for MC photon dose calculations for the rectangular fields specified.

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.

Calibration of a portal imaging device for high-precision dosimetry: a Monte Carlo study

Today electronic portal imaging devices (EPID's) are used primarily to verify patient positioning. They have, however, also the potential as 2D-dosimeters and could be used as such for transit dosimetry or dose reconstruction. It has been proven that such devices, especially liquid filled ionization chambers, have a stable dose response relationship which can be described in terms of the physical properties of the EPID and the pulsed linac radiation. For absolute dosimetry however, an accurate method of calibration to an absolute dose is needed. In this work, we concentrate on calibration against dose in a homogeneous water phantom. Using a Monte Carlo model of the detector we calculated dose spread kernels in units of absolute dose per incident energy fluence and compared them to calculated dose spread kernels in water at different depths. The energy of the incident pencil beams varied between 0.5 and 18 MeV. At the depth of dose maximum in water for a 6 MV beam (1.5 cm) and for a 18 MV beam (3.0 cm) we observed large absolute differences between water and detector dose above an incident energy of 4 MeV but only small relative differences in the most frequent energy range of the beam energy spectra. It is shown that for a 6 MV beam the absolute reference dose measured at 1.5 cm water depth differs from the absolute detector dose by 3.8%. At depth 1.2 cm in water, however, the relative dose differences are almost constant between 2 and 6 MeV. The effects of changes in the energy spectrum of the beam on the dose responses in water and in the detector are also investigated. We show that differences larger than 2% can occur for different beam qualities of the incident photon beam behind water slabs of different thicknesses. It is therefore concluded that for high-precision dosimetry such effects have to be taken into account. Nevertheless, the precise information about the dose response of the detector provided in this Monte Carlo study forms the basis of extracting directly the basic radiometric quantities photon fluence and photon energy fluence from the detector's signal using a deconvolution algorithm. The results are therefore promising for future application in absolute transit dosimetry and absolute dose reconstruction.

Tree-level contribution to B[over ¯]→X_{d}γ using fragmentation functions

We evaluate the most important tree-level contributions connected with the b→uu ¯ ¯ dγ transition to the inclusive radiative decay B ¯ ¯ ¯ →X d γ using fragmentation functions. In this framework the singularities arising from collinear photon emission from the light quarks (u , u ¯ ¯ , and d ) can be absorbed into the (bare) quark-to-photon fragmentation function. We use as input the fragmentation function extracted by the ALEPH group from the two-jet cross section measured at the large electron positron (LEP) collider, where one of the jets is required to contain a photon. To get the quark-to-photon fragmentation function at the fragmentation scale μ F ∼m b , we use the evolution equation, which we solve numerically. We then calculate the (integrated) photon energy spectrum for b→uu ¯ ¯ dγ related to the operators P u 1,2 . For comparison, we also give the corresponding results when using nonzero (constituent) masses for the light quarks.

Harmonic Oscillator in a 1D or 2D Cavity with General Perfectly Reflecting Walls

We investigate the simple harmonic oscillator in a 1-d box, and the 2-d isotropic harmonic oscillator problem in a circular cavity with perfectly reflecting boundary conditions. The energy spectrum has been calculated as a function of the self-adjoint extension parameter. For sufficiently negative values of the self-adjoint extension parameter, there are bound states localized at the wall of the box or the cavity that resonate with the standard bound states of the simple harmonic oscillator or the isotropic oscillator. A free particle in a circular cavity has been studied for the sake of comparison. This work represents an application of the recent generalization of the Heisenberg uncertainty relation related to the theory of self-adjoint extensions in a finite volume.

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-Axis Beam

The T2K collaboration reports a precision measurement of muon neutrino disappearance with an off-axis neutrino beam with a peak energy of 0.6 GeV. Near detector measurements are used to constrain the neutrino flux and cross section parameters. The Super-Kamiokande far detector, which is 295 km downstream of the neutrino production target, collected data corresponding to 3.01×1020 protons on target. In the absence of neutrino oscillations, 205±17 (syst.) events are expected to be detected and only 58 muon neutrino event candidates are observed. A fit to the neutrino rate and energy spectrum assuming three neutrino flavors, normal mass hierarchy and θ23≤π/4 yields a best-fit mixing angle sin2(2θ23)=1.000 and mass splitting |Δm232|=2.44×10−3 eV2/c4. If θ23≥π/4 is assumed, the best-fit mixing angle changes to sin2(2θ23)=0.999 and the mass splitting remains unchanged.

Shaping frequency entangled qudits

We demonstrate the creation, characterization, and manipulation of frequency-entangled qudits by shaping the energy spectrum of entangled photons. The generation of maximally entangled qudit states is verified up to dimension d=4 through tomographic quantum-state reconstruction. Subsequently, we measure Bell parameters for qubits and qutrits as a function of their degree of entanglement. In agreement with theoretical predictions, we observe that for qutrits the Bell parameter is less sensitive to a varying degree of entanglement than for qubits. For frequency-entangled photons, the dimensionality of a qudit is ultimately limited by the bandwidth of the pump laser and can be on the order of a few millions.