A method of determining the bremsstrahlung flux-weighted average photonuclear cross section

Based on a developed analytical model, a method is proposed for measuring the photonuclear cross section averaged over bremsstrahlung flux without application of additional target-monitor of photon flux. The method involves the use of a thin isotopic target, that completely overlaps the photon beam (a photonuclear converter), as well as an algorithm for processing the data on the yield of a reaction under study in such a target. The novel technique was validated on the reactions 100Mo(γ,n)99Mo and 58Ni(γ,n)57Ni in the range of photon end-point energy of 40.7-93.9 MeV. The photon flux-weighted average cross sections of the reactions measured experimentally are in good agreement with Monte Carlo simulations and TALYS predictions on their excitation functions.

Extending deterministic transport capabilities for very-high and ultra-high energy electron beams

Focused Very-High Energy Electron (VHEE, 50-300 MeV) and Ultra-High Energy Electron (UHEE, > 300 MeV) beams can accurately target both large and deeply seated human tumors with high sparing properties, while avoiding the spatial requirements and cost of proton and heavy ion facilities. Advanced testing phases are underway at the CLEAR facilities at CERN (Switzerland), NLCTA at Stanford (USA), and SPARC at INFN (Italy), aiming to accelerate the transition to clinical application. Currently, Monte Carlo (MC) transport is the sole paradigm supporting preclinical trials and imminent clinical deployment. In this paper, we propose an alternative: the first extension of the nuclear-reactor deterministic chain NJOY-DRAGON for VHEE and UHEE applications. We have extended the Boltzmann-Fokker-Planck (BFP) multigroup formalism and validated it using standard radio-oncology benchmarks, complex assemblies with a wide range of atomic numbers, and comprehensive irradiation of the entire periodic table. We report that [Formula: see text] of water voxels exhibit a BFP-MC deviation below [Formula: see text] for electron energies under [Formula: see text]. Additionally, we demonstrate that at least [Formula: see text] of voxels of bone, lung, adipose tissue, muscle, soft tissue, tumor, steel, and aluminum meet the same criterion between [Formula: see text] and [Formula: see text]. For water, the thorax, and the breast intra-operative benchmark, typical average BFP-MC deviations of [Formula: see text] and [Formula: see text] were observed at [Formula: see text] and [Formula: see text], respectively. By irradiating the entire periodic table, we observed similar performance between lithium ([Formula: see text]) and cerium ([Formula: see text]). Deficiencies observed between praseodymium ([Formula: see text]) and einsteinium ([Formula: see text]) have been reported, analyzed, and quantified, offering critical insights for the ongoing development of the Evaluated Nuclear Data File mode in NJOY.

Feasibility of a multigroup Boltzmann-Fokker-Planck solution for electron beam dose calculations

Legacy nuclear-reactor Boltzmann solvers start clinical deployment as an alternative to Monte Carlo (MC) codes and Fermi-Eyges semiemprical models in radiation oncology treatment planning. Today's certified clinical solvers are limited to photon beams. In this paper, ELECTR, a state-of-the-art multigroup electron cross sections generation module in NJOY is presented and validated against Lockwood's calorimetric measurements, EGS-nrc and GEANT-4 for 1-20 MeV unidirectional electron beams. The nuclear-reactor DRAGON-5 solver is upgraded to access the library and solve the Boltzmann-Fokker-Planck (BFP) equation. A variety of heterogeneous radiotherapy and radiosurgery phantom configurations were used for validation purpose. Case studies include a thorax benchmark, that of a typical breast Intra-Operative Radiotherapy and a high-heterogeneity patient-like benchmark. For all beams, [Formula: see text] of the water voxels satisfied the American Association of Physicists in Medicine accuracy criterion for a BFP-MC dose error below [Formula: see text]. At least, [Formula: see text] of adipose, muscle, bone, lung, tumor and breast voxels satisfied the [Formula: see text] criterion. The average BFP-MC relative error was about [Formula: see text] for all voxels, beams and materials combined. By irradiating homogeneous slabs from [Formula: see text] (hydrogen) to [Formula: see text] (einsteinium), we reported performance and defects of the CEPXS mode [US. Sandia National Lab., SAND-89-1685] in ELECTR for the entire periodic table. For all Lockwood's benchmarks, NJOY-DRAGON dose predictions are within the experimental data precision for [Formula: see text] of voxels.

Dark Matter Detection with Bound Nuclear Targets: The Poisson Phonon Tail

Dark matter (DM) scattering with nuclei in solid-state systems may produce elastic nuclear recoil at high energies and single-phonon excitation at low energies. When the DM momentum is comparable to the momentum spread of nuclei bound in a lattice, q_{0}=sqrt[2m_{N}ω_{0}] where m_{N} is the mass of the nucleus and ω_{0} is the optical phonon energy, an intermediate scattering regime characterized by multiphonon excitations emerges. We study a greatly simplified model of a single nucleus in a harmonic potential and show that, while the mean energy deposited for a given momentum transfer q is equal to the elastic value q^{2}/(2m_{N}), the phonon occupation number follows a Poisson distribution and thus the energy spread is ΔE=qsqrt[ω_{0}/(2m_{N})]. This observation suggests that low-threshold calorimetric detectors may have significantly increased sensitivity to sub-GeV DM compared to the expectation from elastic scattering, even when the energy threshold is above the single-phonon energy, by exploiting the tail of the Poisson distribution for phonons above the elastic energy. We use a simple model of electronic excitations to argue that this multiphonon signal will also accompany ionization signals induced from DM-electron scattering or the Migdal effect. In well-motivated models where DM couples to a heavy, kinetically mixed dark photon, we show that these signals can probe experimental milestones for cosmological DM production via thermal freeze-out, including the thermal target for Majorana fermion DM.

Intensities of the strongest γ-ray transitions originating from the 112gIn decay determined via photoactivation yield measurements

Analyzing the γ-ray spectra of the decay of the ^112m,gIn isomer pair produced by the (γ, n)-reaction the intensities of strongest 606.8, 617.5, 851.2, 1253.5, and 1468.8 keV radiation transitions following the ground state decay were correctly determined. They amounted to 0.87(9), 3.9(4), 0.12(2), 0.17(3), and 0.076(12), respectively, and differ from those given in the current databases.

Role of Chiral Two-Body Currents in ^{6}Li Magnetic Properties in Light of a New Precision Measurement with the Relative Self-Absorption Technique

A direct measurement of the decay width of the excited 0_{1}^{+} state of ^{6}Li using the relative self-absorption technique is reported. Our value of Γ_{γ,0_{1}^{+}→1_{1}^{+}}=8.17(14)_{stat.}(11)_{syst.}  eV provides sufficiently low experimental uncertainties to test modern theories of nuclear forces. The corresponding transition rate is compared to the results of ab initio calculations based on chiral effective field theory that take into account contributions to the magnetic dipole operator beyond leading order. This enables a precision test of the impact of two-body currents that enter at next-to-leading order.

Isomer ratios for products of photonuclear reactions on Rh

Over the past several years various preequilibrium models for nuclear reaction mechanisms description were developed. Diversified detailed experimental data in the medium excitation energy region for nuclei are needed for reasonable selection among these theoretical models. Lack of experimental data in this energy region does essentially limit the possibilities for analysis and comparison of different preequilibrium theoretical models. For photonuclear reactions this energy range covers 30-100 MeV. Experimental measurements and estimations of isomer ratios for products of photonuclear reactions with multiple particle escape on antimony were performed using bremsstrahlung spectrum as projectile with end-point energies 74,9 and 85,7 MeV. Method of the induced activity measurement was applied. For acquisition of gamma spectra we used HPGe spectrometer with 20% relative efficiency. Linear accelerator of electrons LU-40 was a source of bremsstrahlung. Energy resolution of electron beam was about 1% and a mean electron current varied within (3.8 - 5.3)μA.

Simulation of a medical linear accelerator for teaching purposes

Simulation software for medical linear accelerators that can be used in a teaching environment was developed. The components of linear accelerators were modeled to first order accuracy using analytical expressions taken from the literature. The expressions used constants that were empirically set such that realistic response could be expected. These expressions were programmed in a MATLAB environment with a graphical user interface in order to produce an environment similar to that of linear accelerator service mode. The program was evaluated in a systematic fashion, where parameters affecting the clinical properties of medical linear accelerator beams were adjusted independently, and the effects on beam energy and dose rate recorded. These results confirmed that beam tuning adjustments could be simulated in a simple environment. Further, adjustment of service parameters over a large range was possible, and this allows the demonstration of linear accelerator physics in an environment accessible to both medical physicists and linear accelerator service engineers. In conclusion, a software tool, named SIMAC, was developed to improve the teaching of linear accelerator physics in a simulated environment. SIMAC performed in a similar manner to medical linear accelerators. The authors hope that this tool will be valuable as a teaching tool for medical physicists and linear accelerator service engineers.

PACS number: 87.55Gh, 87.56bd

Accuracy evaluation of a Compton X-ray spectrometer with bremsstrahlung X-rays generated by a 6 MeV electron bunch

A Compton-scattering-based X-ray spectrometer is developed to obtain the energy distribution of fast electrons produced by intense laser and matter interactions. Bremsstrahlung X-rays generated by fast electrons in a material are used to measure fast electrons' energy distribution in matter. In the Compton X-ray spectrometer, X-rays are converted into recoil electrons by Compton scattering in a converter made from fused silica glass, and a magnet-based electron energy analyzer is used to measure the energy distribution of the electrons that recoil in the direction of the incident X-rays. The spectrum of the incident X-rays is reconstructed from the energy distribution of the recoil electrons. The accuracy of this spectrometer is evaluated using a quasi-monoenergetic 6 MeV electron bunch that emanates from a linear accelerator. An electron bunch is injected into a 1.5 mm thick tungsten plate to produce bremsstrahlung X-rays. The spectrum of these bremsstrahlung X-rays is obtained in the range from 1 to 9 MeV. The energy of the electrons in the bunch is estimated using a Monte Carlo simulation of particle-matter interactions. The result shows that the spectrometer's energy accuracy is ±0.5 MeV for 6.0 MeV electrons.

[Presumption of the energy-spectrum of high-energy electron beam based on the beta-distribution model]

The energy spectra of high-energy electron beams used in radiotherapy are the most important data for evaluating absorbed doses and/or dose distributions in the body of a patient. However, it is impossible to measure the actual spectra of a high-energy electron beam. In this study, we suggest a method to presume the spectra of high-energy electron beams by use of the beta distribution model. The procedure of this method is as follows: (1) the spectrum of the high-energy electron beam was assumed to have a maximum energy Emax, and α, β parameters of the beta probability density function. (2) The percentage depth dose (PDD) based on the assumed spectrum was calculated by a Monte Carlo simulation. (3) The best matching energy spectrum was searched in comparison with the experimental PDD curves. Finally, the optimal energy spectrum of the electron beam was estimated after reiterating the process from (1) to (3). With our method, the measured PDD curves were optimally simulated following the experimental data. It appeared that the assumed spectra approximated well to the actual spectra. However, the error between the assumed and experimental data was observed in the region under the incident surface. We believe this was due to the influence of low-energy electrons scattered at installed collimators, etc. In order to simulate PDDs in this region accurately, a further correction process is required for a spectrum based on the beta distribution model.

Dose profile measurement using an imaging plate: Evaluation of filters using Monte Carlo simulation of 4 MV x-rays

Computed radiography (CR) is gradually replacing film. The application of CR for two-dimensional profiles and off-axis ratio (OAR) measurement using an imaging plate (IP) in a CR system is currently under discussion. However, a well known problem for IPs in dosimetry is that they use high atomic number (Z) materials, such as Ba, which have an energy dependency in a photon interaction. Although there are some reports that it is possible to compensate for the energy dependency with metal filters, the appropriate thicknesses of these filters and where they should be located have not been investigated. The purpose of this study is to find the most suitable filter for use with an IP as a dosimetric tool. Monte Carlo simulation (Geant4 8.1) was used to determine the filter to minimize the measurement error in OAR measurements of 4 MV x-rays. In this simulation, the material and thickness of the filter and distance between the IP and the filter were varied to determine most suitable filter conditions that gave the best fit to the MC calculated OAR in water. With regard to changing the filter material, we found that using higher Z and higher density material increased the effectiveness of the filter. Also, increasing the distance between the filter and the IP reduced the effectiveness, whereas increasing the thickness of the filter increased the effectiveness. The result of this study showed that the most appropriate filter conditions consistent with the calculated OAR in water were the ones with the IP sandwiched between two 2 mm thick lead filters at a distance of 5 mm from the IP or the IP sandwiched directly between two 1 mm lead filters. Using these filters, we measured the OAR at 10 cm depth with 100 cm source-to-surface distance and surface 10x10 cm(2) field size. The results of this measurement represented that it is possible to achieve measurements with less than within 2.0% and 2.0% in the field and with less than 1.1% and 0.6% out of the field by using 2 and 1 mm lead filters, respectively.

Giant electromagnetic vortex and MeV monoenergetic electrons generated by short laser pulses in underdense plasma near quarter critical density region

Very efficient generation of monoenergetic, about 1MeV , electrons from underdense plasma with its electron density close to the critical, when irradiated by an intense femtosecond laser pulse, is found via two dimensional particle-in-cell simulation. The stimulated Raman scattering of a laser pulse with frequency omega< or =2omega(pl max) gives rise to a giant electromagnetic vortex. In contrast to electron acceleration by the well-known laser pulse wake, injected plasma electrons are accelerated up to vortex ponderomotive potential forming a quite monoenergetic distribution. A relatively high charge of such an electron source makes very efficient generation of soft gamma rays with homega>300 keV .

Examination of beam profile measurement using an imaging plate by the light exposure fading method for quality assurance of external radiation therapy

High-resolution film dosimetry has been used for several decades to check and to measure two-dimensional dose distributions. However, in recent years, the automatic processor has been replaced by the spread of computed radiography, or has been little used hospitals. In this study, we measured the off-center ratio (OCR) of the open field, after an irradiating radiation beam was delivered to the imaging plate (IP) under conditions in which the IP was exposed to a fixed amount of light with fading, and compared these data with the OCR measured by an ionization-chamber dosimeter, which is the standard method used for measuring radiation dose. Profile measurement using IP could be achieved by performing light fading, even at a range of more than 100 MU. Further, by using a metallic filter, we succeeded in demonstrating that the profile measurement of IP in an open irradiation field could approximate the values of those obtained by an ionization chamber dosimeter. This method can serve as a simple, easy-to-use method for evaluating the QA of dose distribution in radiation therapy.

[Differential dose albedo for high-energy X-rays on concrete slab]

We computed the differential dose albedo (alpha(D)) for high-energy X-rays on a concrete slab when the incident angle, reflection angle, and azimuth angle were changed, by means of Monte Carlo simulation. We found that alpha(D) changed with incident, reflection, and azimuth angles to the concrete slab. On the whole, the larger the incident angle, the larger alpha(D) tended to become. If the incident angle and reflection angle were the same, the larger the azimuth angle, the smaller alpha(D) tended to become. When the incident, reflection, and azimuth angles were the same, the smaller the X-ray energy was, the larger alpha(D) became, in the order of 10 MV, 6 MV, and 4 MV X-rays.

[Air kerma transmission factors of Scattered X-rays in the maze of a Linac room for lead shield]

Spectra of scattered X-rays in the maze of a Linac (X-ray energies of 4, 6, and 10 MV) room were estimated by means of the Monte Carlo simulation, and air kerma transmission factors of the X-rays scattered through a lead shield were evaluated based on those spectra. Spectra of scattered X-rays showed a maximum in the energy area below 200 keV. The higher the accelerated electron energy, also, the smaller the scattering angle that tended to spread to the higher energy area of the distribution of spectra. The air kerma transmission factor of 120 degrees scattered X-rays of 4 MV X-rays obtained in this study was larger than the transmission factors of 124 degrees scattered photons of (60)Co gamma rays through a lead shield given in ICRP. The air kerma transmission factors of 120 degrees scattered X-rays of 6 MV X-rays were smaller than the transmission factors of 90 degrees scattered photons of (60)Co gamma rays. The air kerma transmission factors of 120 degrees scattered X-rays of 10 MV X-rays was slightly larger than transmission factors of 90 degrees scattered photons of (60)Co gamma rays. Therefore, in the case of a 4 MV X-ray Linac room, the calculation method given in the "Manual of Practical Shield Calculation of Radiation Facilities (2000)" causes underestimation of leakage doses.

An MCNP-based model of a medical linear accelerator x-ray photon beam

The major components in the x-ray photon beam path of the treatment head of the VARIAN Clinac 2300 EX medical linear accelerator were modeled and simulated using the Monte Carlo N-Particle radiation transport computer code (MCNP). Simulated components include x-ray target, primary conical collimator, x-ray beam flattening filter and secondary collimators. X-ray photon energy spectra and angular distributions were calculated using the model. The x-ray beam emerging from the secondary collimators were scored by considering the total x-ray spectra from the target as the source of x-rays at the target position. The depth dose distribution and dose profiles at different depths and field sizes have been calculated at a nominal operating potential of 6 MV and found to be within acceptable limits. It is concluded that accurate specification of the component dimensions, composition and nominal accelerating potential gives a good assessment of the x-ray energy spectra.

A simple method for bremsstrahlung spectra reconstruction from transmission measurements

A new method for evaluation of bremsstrahlung spectra from transmission measurements has been developed. In this method some very well known facts relating to thick target bremsstrahlung spectra are a priori included in the calculation procedure. Some characteristics of the method are preliminarily illustrated on a 6 MV therapy linear accelerator.

The energy-dependent electron loss model for pencil beam dose kernels

The 'monoenergetic' electron loss model was derived in a previous work to account for pathlength straggling in the Fermi-Eyges pencil beam problem. In this paper, we extend this model to account for energy-loss straggling and secondary knock-on electron transport in order to adequately predict a depth dose curve. To model energy-loss straggling, we use a weighted superposition of a discrete number of monoenergetic pencil beams with different initial energies where electrons travel along the depth-energy characteristics in the continuous slowing down approximation (CSDA). The energy straggling spectrum at depth determines the weighting assigned to each monoenergetic pencil beam. Supplemented by a simple transport model for the secondary knock-on electrons, the 'energy-dependent' electron loss model predicts both lateral and depth dose distributions from the electron pencil beams in good agreement with Monte Carlo calculations and measurements. The calculation of dose distribution from a pencil beam takes 0.2 s on a Pentium III 500 MHz computer. Being computationally fast, the 'energy-dependent' electron loss model can be used for the calculation of 3D energy deposition kernels in dose optimization schemes without using precalculated or measured data.

Reconstruction of megavoltage photon spectra from electronic portal imager derived transmission measurements

Three variants of the Schiff equation are investigated to model the spectra produced by megavoltage linear accelerators. These models are tested against well-validated Monte Carlo (MC) generated spectra on the central axis of large-area fields, and show excellent agreement. Numerical reconstructions of 6 and 10 MV spectra using the same models are then presented, using experimental attenuation data derived from an electronic portal imager. The process of deriving spectra from experimental attenuation data is shown to be inherently badly constrained mathematically, with the derived spectrum being highly sensitive to noise in the source data, and non-unique. By placing a priori constraints on the Schiff model from both physical knowledge of the construction of the accelerator and MC data, physically useful results are gained and presented for both the energy dependence and off-axis behaviour of photon spectra.

An MCNP-based model of a linear accelerator x-ray beam

The Monte Carlo N-Particle radiation transport computer code (MCNP) has been employed on a personal computer to develop a simple model simulating the major components within the beam path of a linear accelerator radiation head, namely the electron target, primary conical collimator, beam flattening filter, wedge filter and the secondary collimators. The model was initially used to calculate the energy spectra and angular distributions of the x-ray beam for the Philips SL 75/5 linear accelerator, in a plane immediately beneath the flattening filter. These data were subsequently used as a 'source' of x-rays at the target position, to assess the emergent beam from the secondary collimators. The depth dose distributions and dose profiles at constant depth for various field sizes have been calculated for a nominal operating potential of 4 MV and found to be within acceptable limits. It is concluded that the technique may be used to calculate the energy spectra of any linear accelerator upon specification of the component dimensions, materials and nominal accelerating potential. It is anticipated that this work will serve as the basis of a quality control tool for linear accelerators and treatment planning systems.

115In as a probe for the characterization of therapy bremsstrahlung beams and the detection of photoneutrons

Variation of the photoactivation rate across radiation fields of three different bremsstrahlung beams of two medical accelerators has been measured, making use of the photonuclear reactions in natural indium probes: 115In(y,y')115mIn and 115In(y,n)114mIn. The third nuclear reaction, 115In(n,y)116mIn, was used to detect the presence of neutrons in the photon beam and to estimate the spatial distribution of thermal and fast neutrons in the patient plane as a function of collimator opening.