Generalized eMC implementation for Monte Carlo dose calculation of electron beams from different machine types

The electron Monte Carlo (eMC) dose calculation algorithm available in the Eclipse treatment planning system (Varian Medical Systems) is based on the macro MC method and uses a beam model applicable to Varian linear accelerators. This leads to limitations in accuracy if eMC is applied to non-Varian machines. In this work eMC is generalized to also allow accurate dose calculations for electron beams from Elekta and Siemens accelerators. First, changes made in the previous study to use eMC for low electron beam energies of Varian accelerators are applied. Then, a generalized beam model is developed using a main electron source and a main photon source representing electrons and photons from the scattering foil, respectively, an edge source of electrons, a transmission source of photons and a line source of electrons and photons representing the particles from the scrapers or inserts and head scatter radiation. Regarding the macro MC dose calculation algorithm, the transport code of the secondary particles is improved. The macro MC dose calculations are validated with corresponding dose calculations using EGSnrc in homogeneous and inhomogeneous phantoms. The validation of the generalized eMC is carried out by comparing calculated and measured dose distributions in water for Varian, Elekta and Siemens machines for a variety of beam energies, applicator sizes and SSDs. The comparisons are performed in units of cGy per MU. Overall, a general agreement between calculated and measured dose distributions for all machine types and all combinations of parameters investigated is found to be within 2% or 2 mm. The results of the dose comparisons suggest that the generalized eMC is now suitable to calculate dose distributions for Varian, Elekta and Siemens linear accelerators with sufficient accuracy in the range of the investigated combinations of beam energies, applicator sizes and SSDs.

A critical evaluation of secondary cancer risk models applied to Monte Carlo dose distributions of 2-dimensional, 3-dimensional conformal and hybrid intensity-modulated radiation therapy for breast cancer

The comparison of radiotherapy techniques regarding secondary cancer risk has yielded contradictory results possibly stemming from the many different approaches used to estimate risk. The purpose of this study was to make a comprehensive evaluation of different available risk models applied to detailed whole-body dose distributions computed by Monte Carlo for various breast radiotherapy techniques including conventional open tangents, 3D conformal wedged tangents and hybrid intensity modulated radiation therapy (IMRT). First, organ-specific linear risk models developed by the International Commission on Radiological Protection (ICRP) and the Biological Effects of Ionizing Radiation (BEIR) VII committee were applied to mean doses for remote organs only and all solid organs. Then, different general non-linear risk models were applied to the whole body dose distribution. Finally, organ-specific non-linear risk models for the lung and breast were used to assess the secondary cancer risk for these two specific organs. A total of 32 different calculated absolute risks resulted in a broad range of values (between 0.1% and 48.5%) underlying the large uncertainties in absolute risk calculation. The ratio of risk between two techniques has often been proposed as a more robust assessment of risk than the absolute risk. We found that the ratio of risk between two techniques could also vary substantially considering the different approaches to risk estimation. Sometimes the ratio of risk between two techniques would range between values smaller and larger than one, which then translates into inconsistent results on the potential higher risk of one technique compared to another. We found however that the hybrid IMRT technique resulted in a systematic reduction of risk compared to the other techniques investigated even though the magnitude of this reduction varied substantially with the different approaches investigated. Based on the epidemiological data available, a reasonable approach to risk estimation would be to use organ-specific non-linear risk models applied to the dose distributions of organs within or near the treatment fields (lungs and contralateral breast in the case of breast radiotherapy) as the majority of radiation-induced secondary cancers are found in the beam-bordering regions.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

Saddlepoint approximations to sensitivities of tail probabilities of random sums and comparisons with Monte Carlo estimators

This article proposes computing sensitivities of upper tail probabilities of random sums by the saddlepoint approximation. The considered sensitivity is the derivative of the upper tail probability with respect to the parameter of the summation index distribution. Random sums with Poisson or Geometric distributed summation indices and Gamma or Weibull distributed summands are considered. The score method with importance sampling is considered as an alternative approximation. Numerical studies show that the saddlepoint approximation and the method of score with importance sampling are very accurate. But the saddlepoint approximation is substantially faster than the score method with importance sampling. Thus, the suggested saddlepoint approximation can be conveniently used in various scientific problems.

Vegetation and fire history of coastal north-eastern Sardinia (Italy) under changing Holocene climates and land use

Little is known about the vegetation and fire history of Sardinia, and especially the long-term history of the thermo-Mediterranean belt that encompasses its entire coastal lowlands. A new sedimentary record from a coastal lake based on pollen, spores, macrofossils and microscopic charcoal analysis is used to reconstruct the vegetation and fire history in north-eastern Sardinia. During the mid-Holocene (c. 8,100–5,300 cal bp), the vegetation around Stagno di Sa Curcurica was characterised by dense Erica scoparia and E. arborea stands, which were favoured by high fire activity. Fire incidence declined and evergreen broadleaved forests of Quercus ilex expanded at the beginning of the late Holocene. We relate the observed vegetation and fire dynamics to climatic change, specifically moister and cooler summers and drier and milder winters after 5,300 cal bp. Agricultural activities occurred since the Neolithic and intensified after c. 7,000 cal bp. Around 2,750 cal bp, a further decline of fire incidence and Erica communities occurred, while Quercus ilex expanded and open-land communities became more abundant. This vegetation shift coincided with the historically documented beginning of Phoenician period, which was followed by Punic and Roman civilizations in Sardinia. The vegetational change at around 2,750 cal bp was possibly advantaged by a further shift to moister and cooler summers and drier and milder winters. Triggers for climate changes at 5,300 and 2,750 cal bp may have been gradual, orbitally-induced changes in summer and winter insolation, as well as centennial-scale atmospheric reorganizations. Open evergreen broadleaved forests persisted until the twentieth century, when they were partly substituted by widespread artificial pine plantations. Our results imply that highly flammable Erica vegetation, as reconstructed for the mid-Holocene, could re-emerge as a dominant vegetation type due to increasing drought and fire, as anticipated under global change conditions.

Regulation of leukocyte recruitment by the long pentraxin PTX3

Pentraxins are a superfamily of conserved proteins involved in the acute-phase response and innate immunity. Pentraxin 3 (PTX3), a prototypical member of the long pentraxin subfamily, is a key component of the humoral arm of innate immunity that is essential for resistance to certain pathogens. A regulatory role for pentraxins in inflammation has long been recognized, but the underlying mechanisms remain unclear. Here we report that PTX3 bound P-selectin and attenuated neutrophil recruitment at sites of inflammation. PTX3 released from activated leukocytes functioned locally to dampen neutrophil recruitment and regulate inflammation. Antibodies have glycosylation-dependent regulatory effect on inflammation. Therefore, PTX3, which is an essential component of humoral innate immunity, and immunoglobulins share functional outputs, including complement activation, opsonization and, as shown here, glycosylation-dependent regulation of inflammation.

Contamination of Alpine snow and ice at Colle Gnifetti, Swiss/Italian Alps, from nuclear weapons tests

Plutonium is present in the environment as a consequence of atmospheric nuclear tests, nuclear weapons production and industrial releases over the past 50 years. To study temporal trends, a high resolution Pu record was obtained by analyzing 52 discrete samples of an alpine firn/ice core from Colle Gnifetti (Monte Rosa, 4450 m a.s.l.), dating from 1945 to 1990. The 239Pu signal was recorded directly, without decontamination or preconcentration steps, using an Inductively Coupled Plasma - Sector Field Mass Spectrometer (ICP-SFMS) equipped with an high efficiency sample introduction system, thus requiring much less sample preparation than previously reported methods. The 239Pu profile reflects the three main periods of atmospheric nuclear weapons testing: the earliest peak lasted from 1954/55 to 1958 and was caused by the first testing period reaching a maximum in 1958. Despite a temporary halt of testing in 1959/60, the Pu concentration decreased only by half with respect to the 1958 peak due to long atmospheric residence times. In 1961/62 Pu concentrations rapidly increased reaching a maximum in 1963, which was about 40% more intense than the 1958 peak. After the signing of the "Limited Test Ban Treaty" between USA and USSR in 1964, Pu deposition decreased very sharply reaching a minimum in 1967. The third period (1967-1975) is characterized by irregular Pu concentrations with smaller peaks (about 20-30% of the 1964 peak) which might be related to the deposition of Saharan dust contaminated by the French nuclear tests of the 1960s. The data presented are in very good agreement with Pu profiles previously obtained from the Col du Dome ice core (by multi-collector ICP-MS) and Belukha ice core (by Accelerator Mass Spectrometry, AMS). Although a semi-quantitative method was employed here, the results are quantitatively comparable to previously published results.

Measuring the gravitational free-fall of antihydrogen

Antihydrogen holds the promise to test, for the first time, the universality of freefall with a system composed entirely of antiparticles. The AEgIS experiment at CERN’s antiproton decelerator aims to measure the gravitational interaction between matter and antimatter by measuring the deflection of a beam of antihydrogen in the Earths gravitational field (g). The principle of the experiment is as follows: cold antihydrogen atoms are synthesized in a Penning-Malberg trap and are Stark accelerated towards a moir´e deflectometer, the classical counterpart of an atom interferometer, and annihilate on a position sensitive detector. Crucial to the success of the experiment is the spatial precision of the position sensitive detector.We propose a novel free-fall detector based on a hybrid of two technologies: emulsion detectors, which have an intrinsic spatial resolution of 50 nm but no temporal information, and a silicon strip / scintillating fiber tracker to provide timing and positional information. In 2012 we tested emulsion films in vacuum with antiprotons from CERN’s antiproton decelerator. The annihilation vertices could be observed directly on the emulsion surface using the microscope facility available at the University of Bern. The annihilation vertices were successfully reconstructed with a resolution of 1–2 μmon the impact parameter. If such a precision can be realized in the final detector, Monte Carlo simulations suggest of order 500 antihydrogen annihilations will be sufficient to determine gwith a 1 % accuracy. This paper presents current research towards the development of this technology for use in the AEgIS apparatus and prospects for the realization of the final detector.

Pharao Haremhab und die nördliche Levante. Bemerkungen zu einem Skarabäenabdruck aus Tall Mishrife/Qatna

This paper deals with a clay sealing which seems to feature one scarab impression with the personal name of Pharaoh Haremhab, last king of the 18th Dynasty (ca. 1319/1305?–1292 BCE). The clay sealing was found in a pit within the central area of the Royal Palace of Qaṭna. The impression thus could shed new light on the chronology of the palace, as well as Haremhab´s political involvement in the northern Levant.