Macro Monte Carlo for dose calculation of proton beams

Although the Monte Carlo (MC) method allows accurate dose calculation for proton radiotherapy, its usage is limited due to long computing time. In order to gain efficiency, a new macro MC (MMC) technique for proton dose calculations has been developed. The basic principle of the MMC transport is a local to global MC approach. The local simulations using GEANT4 consist of mono-energetic proton pencil beams impinging perpendicularly on slabs of different thicknesses and different materials (water, air, lung, adipose, muscle, spongiosa, cortical bone). During the local simulation multiple scattering, ionization as well as elastic and inelastic interactions have been taken into account and the physical characteristics such as lateral displacement, direction distributions and energy loss have been scored for primary and secondary particles. The scored data from appropriate slabs is then used for the stepwise transport of the protons in the MMC simulation while calculating the energy loss along the path between entrance and exit position. Additionally, based on local simulations the radiation transport of neutrons and the generated ions are included into the MMC simulations for the dose calculations. In order to validate the MMC transport, calculated dose distributions using the MMC transport and GEANT4 have been compared for different mono-energetic proton pencil beams impinging on different phantoms including homogeneous and inhomogeneous situations as well as on a patient CT scan. The agreement of calculated integral depth dose curves is better than 1% or 1 mm for all pencil beams and phantoms considered. For the dose profiles the agreement is within 1% or 1 mm in all phantoms for all energies and depths. The comparison of the dose distribution calculated using either GEANT4 or MMC in the patient also shows an agreement of within 1% or 1 mm. The efficiency of MMC is up to 200 times higher than for GEANT4. The very good level of agreement in the dose comparisons demonstrate that the newly developed MMC transport results in very accurate and efficient dose calculations for proton beams.

X-ray phase contrast imaging of metallized and non-metallized biological samples

Imaging of biological samples has been performed with a variety of techniques for example electromagnetic waves, electrons, neutrons, ultrasound and X-rays. Also conventional X-ray imaging represents the basis of medical diagnostic imaging, it remains of limited use in this application because it is based solely on the differential absorption of X-rays by tissues. Coherent and bright photon beams, such as those produced by third-generation synchrotron X-ray sources, provide further information on subtle X-ray phase changes at matter interfaces. This complements conventional X-ray absorption by edge enhancement phenomena. Thus, phase contrast imaging has the potential to improve the detection of structures on images by detecting those structures that are invisible with X-ray absorption imaging. Images of a weakly absorbing nylon fibre were recorded in in-line holography geometry using a high resolution low-noise CCD camera at the ESRF in Grenoble. The method was also applied to improve image contrast for images of biological tissues. This paper presents phase contrast microradiographs of vascular tree casts and images of a housefly. These reveal very fine structures, that remain invisible with conventional absorption contrast only.

Coupled W–Os–Pt isotope systematics in IVB iron meteorites: In situ neutron dosimetry for W isotope chronology

Tungsten isotope compositions of magmatic iron meteorites yield ages of differentiation that are within ±2 Ma of the formation of CAIs, with the exception of IVB irons that plot to systematically less radiogenic compositions yielding erroneously old ages. Secondary neutron capture due to galactic cosmic ray (GCR) irradiation is known to lower the ε182W of iron meteorites, adequate correction of which requires a measure of neutron dosage which has not been available, thus far. The W, Os and Pt isotope systematics of 12 of the 13 known IVB iron meteorites were determined by MC-ICP-MS (W, Os, Pt) and TIMS (Os). On the same dissolutions that yield precise ε182W, stable Os and Pt isotopes were determined as in situ neutron dosimeters for empirical correction of the ubiquitous cosmic-ray induced burn-out of 182W in iron meteorites. The W isotope data reveal a main cluster with ε182W of ∼−3.6, but a much larger range than observed in previous studies including irons (Weaver Mountains and Warburton Range) that show essentially no cosmogenic effect on their ε182W. The IVB data exhibits resolvable negative anomalies in ε189Os (−0.6ε) and complementary ε190Os anomalies (+0.4ε) in Tlacotepec due to neutron capture on 189Os which has approximately the same neutron capture cross section as 182W, and captures neutrons to produce 190Os. The least irradiated IVB iron, Warburton Range, has ε189Os and ε190Os identical to terrestrial values. Similarly, Pt isotopes, which are presented as ε192Pt, ε194Pt and ε196Pt range from +4.4ε to +53ε, +1.54ε to −0.32ε and +0.73ε to −0.20ε, respectively, also identify Tlacotepec and Dumont as the most GCR-damaged samples. In W–Os and W–Pt isotope space, the correlated isotope data back-project toward a 0-epsilon value of ε192Pt, ε189Os and ε190Os from which a pre-GCR irradiation ε182W of −3.42±0.09 (2σ) is derived. This pre-GCR irradiation ε182W is within uncertainty of the currently accepted CAI initial ε182W. The Pt and Os isotope correlations in the IVB irons are in good agreement with a nuclear model for spherical irons undergoing GCR spallation, although this model over-predicts the change of ε182W by ∼2×, indicating a need for better W neutron capture cross section determinations. A nucleosynthetic effect in ε184W in these irons of −0.14±0.08 is confirmed, consistent with the presence of Mo and Ru isotope anomalies in IVB irons. The lack of a non-GCR Os isotope anomaly in these irons requires more complex explanations for the production of W, Ru and Mo anomalies than nebular heterogeneity in the distribution of s-process to r-process nuclides.

Neutron capture on Pt isotopes in iron meteorites and the Hf–W chronology of core formation in planetesimals

The short-lived 182Hf–182W isotope system can provide powerful constraints on the timescales of planetary core formation, but its application to iron meteorites is hampered by neutron capture reactions on W isotopes resulting from exposure to galactic cosmic rays. Here we show that Pt isotopes in magmatic iron meteorites are also affected by capture of (epi)thermal neutrons and that the Pt isotope variations are correlated with variations in 182W/184W. This makes Pt isotopes a sensitive neutron dosimeter for correcting cosmic ray-induced W isotope shifts. The pre-exposure 182W/184W derived from the Pt–W isotope correlations of the IID, IVA and IVB iron meteorites are higher than most previous estimates and are more radiogenic than the initial 182W/184W of Ca–Al-rich inclusions (CAI). The Hf–W model ages for core formation range from +1.6±1.0 million years (Ma; for the IVA irons) to +2.7±1.3 Ma after CAI formation (for the IID irons), indicating that there was a time gap of at least ∼1 Ma between CAI formation and metal segregation in the parent bodies of some iron meteorites. From the Hf–W ages a time limit of <1.5–2 Ma after CAI formation can be inferred for the accretion of the IID, IVA and IVB iron meteorite parent bodies, consistent with earlier conclusions that the accretion of differentiated planetesimals predated that of most chondrite parent bodies.

Thermal neutron capture effects in radioactive and stable nuclide systems

Neutron capture effects in meteorites and lunar surface samples have been successfully used in the past to study exposure histories and shielding conditions. In recent years, however, it turned out that neutron capture effects produce a nuisance for some of the short-lived radionuclide systems. The most prominent example is the 182Hf-182W system in iron meteorites, for which neutron capture effects lower the 182W/184W ratio, thereby producing too old apparent ages. Here, we present a thorough study of neutron capture effects in iron meteorites, ordinary chondrites, and carbonaceous chondrites, whereas the focus is on iron meteorites. We study in detail the effects responsible for neutron production, neutron transport, and neutron slowing down and find that neutron capture in all studied meteorite types is not, as usually expected, exclusively via thermal neutrons. In contrast, most of the neutron capture in iron meteorites is in the epithermal energy range and there is a significant contribution from epithermal neutron capture even in stony meteorites. Using sophisticated particle spectra and evaluated cross section data files for neutron capture reactions we calculate the neutron capture effects for Sm, Gd, Cd, Pd, Pt, and Os isotopes, which all can serve as neutron-dose proxies, either in stony or in iron meteorites. In addition, we model neutron capture effects in W and Ag isotopes. For W isotopes, the GCR-induced shifts perfectly correlate with Os and Pt isotope shifts, which therefore can be used as neutron-dose proxies and permit a reliable correction. We also found that GCR-induced effects for the 107Pd-107Ag system can be significant and need to be corrected, a result that is in contrast to earlier studies.

The abundance and isotopic composition of Cd in iron meteorites

Cadmium is a highly volatile element and its abundance in meteorites may help better understand volatility-controlled processes in the solar nebula and on meteorite parent bodies. The large thermal neutron capture cross section of 113Cd suggests that Cd isotopes might be well suited to quantify neutron fluences in extraterrestrial materials. The aims of this study were (1) to evaluate the range and magnitude of Cd concentrations in magmatic iron meteorites, and (2) to assess the potential of Cd isotopes as a neutron dosimeter for iron meteorites. Our new Cd concentration data determined by isotope dilution demonstrate that Cd concentrations in iron meteorites are significantly lower than in some previous studies. In contrast to large systematic variations in the concentration of moderately volatile elements like Ga and Ge, there is neither systematic variation in Cd concentration amongst troilites, nor amongst metal phases of different iron meteorite groups. Instead, Cd is strongly depleted in all iron meteorite groups, implying that the parent bodies accreted well above the condensation temperature of Cd (i.e., ≈650 K) and thus incorporated only minimal amounts of highly volatile elements. No Cd isotope anomalies were found, whereas Pt and W isotope anomalies for the same iron meteorite samples indicate a significant fluence of epithermal and higher energetic neutrons. This observation demonstrates that owing to the high Fe concentrations in iron meteorites, neutron capture mainly occurs at epithermal and higher energies. The combined Cd-Pt-W isotope results from this study thus demonstrate that the relative magnitude of neutron capture-induced isotope anomalies is strongly affected by the chemical composition of the irradiated material. The resulting low fluence of thermal neutrons in iron meteorites and their very low Cd concentrations make Cd isotopes unsuitable as a neutron dosimeter for iron meteorites.

The neutron background of the XENON100 dark matter experiment

TheXENON100 experiment, installed underground at the LaboratoriNazionali del Gran Sasso, aims to directly detect dark matter in the form of weakly interacting massive particles (WIMPs) via their elastic scattering off xenon nuclei. This paper presents a study on the nuclear recoil background of the experiment, taking into account neutron backgrounds from (alpha, n) reactions and spontaneous fission due to natural radioactivity in the detector and shield materials, as well as muon-induced neutrons. Based on MonteCarlo simulations and using measured radioactive contaminations of all detector components, we predict the nuclear recoil backgrounds for the WIMP search results published by theXENON100 experiment in 2011 and 2012, 0.11(-0.04)(+0.08) events and 0.17(-0.07)(+0.12) events, respectively, and conclude that they do not limit the sensitivity of the experiment.

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

Conceptual design and simulation of a water Cherenkov muon veto for the XENON1T experiment

XENON is a dark matter direct detection project, consisting of a time projection chamber (TPC) filled with liquid xenon as detection medium. The construction of the next generation detector, XENON1T, is presently taking place at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy. It aims at a sensitivity to spin-independent cross sections of 2 10-47 c 2 for WIMP masses around 50 GeV2, which requires a background reduction by two orders of magnitude compared to XENON100, the current generation detector. An active system that is able to tag muons and muon-induced backgrounds is critical for this goal. A water Cherenkov detector of ~ 10 m height and diameter has been therefore developed, equipped with 8 inch photomultipliers and cladded by a reflective foil. We present the design and optimization study for this detector, which has been carried out with a series of Monte Carlo simulations. The muon veto will reach very high detection efficiencies for muons (>99.5%) and showers of secondary particles from muon interactions in the rock (>70%). Similar efficiencies will be obtained for XENONnT, the upgrade of XENON1T, which will later improve the WIMP sensitivity by another order of magnitude. With the Cherenkov water shield studied here, the background from muon-induced neutrons in XENON1T is negligible.

Liulin silicon semiconductor spectrometers as cosmic ray monitors at the high mountain observatories Jungfraujoch and Lomnický stit

Currently, most cosmic ray data are obtained by detectors on satellites, aircraft, high-altitude balloons and ground (neutron monitors). In our work, we examined whether Liulin semiconductor spectrometers (simple silicon planar diode detectors with spectrometric properties) located at high mountain observatories could contribute new information to the monitoring of cosmic rays by analyzing data from selected solar events between 2005 and 2013. The decision thresholds and detection limits of these detectors placed at Jungfraujoch (Switzerland; 3475 m a.s.l.; vertical cut-off rigidity 4.5 GV) and Lomnicky stıt (Slovakia; 2633 m a.s.l.; vertical cut-off rigidity 3.84 GV) highmountain observatories were determined. The data showed that only the strongest variations of the cosmic ray flux in this period were detectable. The main limitation in the performance of these detectors is their small sensitive volume and low sensitivity of the PIN photodiode to neutrons.