Cosmogenic production rates and recoil loss effects in micrometeorites and interplanetary dust particles

We present a purely physical model to determine cosmogenic production rates for noble gases and radionuclides in micrometeorites (MMs) and interplanetary dust particles (IDPs) by solar cosmic-rays (SCR) and galactic cosmic-rays (GCR) fully considering recoil loss effects. Our model is based on various nuclear model codes to calculate recoil cross sections, recoil ranges, and finally the percentages of the cosmogenic nuclides that are lost as a function of grain size, chemical composition of the grain, and the spectral distribution of the projectiles. The main advantage of our new model compared with earlier approaches is that we consider the entire SCR particle spectrum up to 240 MeV and not only single energy points. Recoil losses for GCR-produced nuclides are assumed to be equal to recoil losses for SCR-produced nuclides. Combining the model predictions with Poynting-Robertson orbital lifetimes, we calculate cosmic-ray exposure ages for recently studied MMs, cosmic spherules, and IDPs. The ages for MMs and the cosmic-spherule are in the range <2.2–233 Ma, which corresponds, according to the Poynting-Robertson drag, to orbital distances in the range 4.0–34 AU. For two IDPs, we determine exposure ages of longer than 900 Ma, which corresponds to orbital distances larger than 150 AU. The orbital distance in the range 4–6 AU for one MM and the cosmic spherule indicate an origin either in the asteroid belt or release from comets coming either from the Kuiper Belt or the Oort Cloud. Three of the studied MMs have orbital distances in the range 23–34 AU, clearly indicating a cometary origin, either from short-period comets from the Kuiper Belt or from the Oort Cloud. The two IDPs have orbital distances of more than 150 AU, indicating an origin from Oort Cloud comets.

Calibration of cosmogenic noble gas production in ordinary chondrites based on 36 Cl- 36 Ar ages. Part 1: Refined produced rates for cosmogenic 21 Ne and 38 Ar

We measured the concentrations and isotopic compositions of He, Ne, and Ar in bulk samples and metal separates of 14 ordinary chondrite falls with long exposure ages and high metamorphic grades. In addition, we measured concentrations of the cosmogenic radionuclides 10Be, 26Al, and 36Cl in metal separates and in the nonmagnetic fractions of the selected meteorites. Using cosmogenic 36Cl and 36Ar measured in the metal separates, we determined 36Cl-36Ar cosmic-ray exposure (CRE) ages, which are shielding-independent and therefore particularly reliable. Using the cosmogenic noble gases and radionuclides, we are able to decipher the CRE history for the studied objects. Based on the correlation 3He/21Ne versus 22Ne/21Ne, we demonstrate that, among the meteorites studied, only one suffered significant diffusive losses (about 35%). The data confirm that the linear correlation 3He/21Ne versus 22Ne/21Ne breaks down at high shielding. Using 36Cl-36Ar exposure ages and measured noble gas concentrations, we determine 21Ne and 38Ar production rates as a function of 22Ne/21Ne. The new data agree with recent model calculations for the relationship between 21Ne and 38Ar production rates and the 22Ne/21Ne ratio, which does not always provide unique shielding information. Based on the model calculations, we determine a new correlation line for 21Ne and 38Ar production rates as a function of the shielding indicator 22Ne/21Ne for H, L, and LL chondrites with preatmospheric radii less than about 65 cm. We also calculated the 10Be/21Ne and 26Al/21Ne production rate ratios for the investigated samples, which show good agreement with recent model calculations.

Search for long-lived stopped R-hadrons decaying out-of-time with pp collisions using the ATLAS detector

An updated search is performed for gluino, top squark, or bottom squark R-hadrons that have come to rest within the ATLAS calorimeter, and decay at some later time to hadronic jets and a neutralino, using 5.0 and 22.9 fb(-1) of pp collisions at 7 and 8 TeV, respectively. Candidate decay events are triggered in selected empty bunch crossings of the LHC in order to remove pp collision backgrounds. Selections based on jet shape and muon system activity are applied to discriminate signal events from cosmic ray and beam-halo muon backgrounds. In the absence of an excess of events, improved limits are set on gluino, stop, and sbottom masses for different decays, lifetimes, and neutralino masses. With a neutralino of mass 100 GeV, the analysis excludes gluinos with mass below 832 GeV (with an expected lower limit of 731 GeV), for a gluino lifetime between 10 mu s and 1000 s in the generic R-hadron model with equal branching ratios for decays to q (q) over bar(chi) over tilde (0) and g (chi) over tilde (0). Under the same assumptions for the neutralino mass and squark lifetime, top squarks and bottom squarks in the Regge R-hadron model are excluded with masses below 379 and 344 GeV, respectively.

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.

Measurement of the TeV atmospheric muon charge ratio with the complete OPERA data set

The OPERA detector, designed to search for νμ → ντ oscillations in the CNGS beam, is located in the underground Gran Sasso laboratory, a privileged location to study TeV-scale cosmic rays. For the analysis here presented, the detector was used to measure the atmospheric muon charge ratio in the TeV region. OPERA collected chargeseparated cosmic ray data between 2008 and 2012. More than 3 million atmospheric muon events were detected and reconstructed, among which about 110000 multiple muon bundles. The charge ratio Rμ ≡ Nμ+/Nμ− was measured separately for single and for multiple muon events. The analysis exploited the inversion of the magnet polarity which was performed on purpose during the 2012 Run. The combination of the two data sets with opposite magnet polarities allowedminimizing systematic uncertainties and reaching an accurate determination of the muon charge ratio. Data were fitted to obtain relevant parameters on the composition of primary cosmic rays and the associated kaon production in the forward fragmentation region. In the surface energy range 1–20 TeV investigated by OPERA, Rμ is well described by a parametric model including only pion and kaon contributions to themuon flux, showing no significant contribution of the prompt component. The energy independence supports the validity of Feynman scaling in the fragmentation region up to 200 TeV/nucleon primary energy.

A steerable UV laser system for the calibration of liquid argon time projection chambers

A number of liquid argon time projection chambers (LAr TPCs) are being built or are proposed for neutrino experiments on long- and short baseline beams. For these detectors, a distortion in the drift field due to geometrical or physics reasons can affect the reconstruction of the events. Depending on the TPC geometry and electric drift field intensity, this distortion could be of the same magnitude as the drift field itself. Recently, we presented a method to calibrate the drift field and correct for these possible distortions. While straight cosmic ray muon tracks could be used for calibration, multiple coulomb scattering and momentum uncertainties allow only a limited resolution. A UV laser instead can create straight ionization tracks in liquid argon, and allows one to map the drift field along different paths in the TPC inner volume. Here we present a UV laser feed-through design with a steerable UV mirror immersed in liquid argon that can point the laser beam at many locations through the TPC. The straight ionization paths are sensitive to drift field distortions, a fit of these distortion to the linear optical path allows to extract the drift field, by using these laser tracks along the whole TPC volume one can obtain a 3D drift field map. The UV laser feed-through assembly is a prototype of the system that will be used for the MicroBooNE experiment at the Fermi National Accelerator Laboratory (FNAL).

Cosmogenic 180W variations in meteorites and re-assessment of a possible 184Os–180W decay system

We measured tungsten (W) isotopes in 23 iron meteorites and the metal phase of the CB chondrite Gujba in order to ascertain if there is evidence for a large-scale nucleosynthetic heterogeneity in the p-process isotope 180W in the solar nebula as recently suggested by Schulz et al. (2013). We observed large excesses in 180W (up to ≈ 6 ε) in some irons. However, significant within-group variations in magmatic IIAB and IVB irons are not consistent with a nucleosynthetic origin, and the collateral effects on 180W from an s-deficit in IVB irons cannot explain the total variation. We present a new model for the combined effects of spallation and neutron capture reactions on 180W in iron meteorites and show that at least some of the observed within-group variability is explained by cosmic ray effects. Neutron capture causes burnout of 180W, whereas spallation reactions lead to positive shifts in 180W. These effects depend on the target composition and cosmic-ray exposure duration; spallation effects increase with Re/W and Os/W ratios in the target and with exposure age. The correlation of 180W/184W with Os/W ratios in iron meteorites results in part from spallogenic production of 180W rather than from 184Os decay, contrary to a recent study by Peters et al. (2014). Residual ε180W excesses after correction for an s-deficit and for cosmic ray effects may be due to ingrowth of 180W from 184Os decay, but the magnitude of this ingrowth is at least a factor of ≈2 smaller than previously suggested. These much smaller effects strongly limit the applicability of the putative 184Os-180W system to investigate geological problems.

Differences in published characteristics of GLE60 and their consequences on computed radiation dose rates along selected flight paths

The radiation dose rates at flight altitudes can increase by orders of magnitude for a short time during energetic solar cosmic ray events, so called ground level enhancements (GLEs). Especially at high latitudes and flight altitudes, solar energetic particles superposed on galactic cosmic rays may cause radiation that exceeds the maximum allowed dosage limit for the general public. Therefore the determination of the radiation dose rate during GLEs should be as reliable as possible. Radiation dose rates along flight paths are typically determined by computer models that are based on cosmic ray flux and anisotropy parameters derived from neutron monitor and/or satellite measurements. The characteristics of the GLE on 15 April 2001 (GLE60) were determined and published by various authors. In this work we compare these results and investigate the consequences on the computed radiation dose rates along selected flight paths. In addition, we compare the computed radiation dose rates with measurements that were made during GLE60 on board two transatlantic flights.

Cosmogenic activation of xenon and copper

Rare event search experiments using liquid xenon as target and detection medium require ultra-low background levels to fully exploit their physics potential. Cosmic ray induced activation of the detector components and, even more importantly, of the xenon itself during production, transportation and storage at the Earth's surface, might result in the production of radioactive isotopes with long half-lives, with a possible impact on the expected background. We present the first dedicated study on the cosmogenic activation of xenon after 345 days of exposure to cosmic rays at the Jungfraujoch research station at 3470m above sea level, complemented by a study of copper which has been activated simultaneously. We have directly observed the production of 7Be, 101Rh, 125Sb, 126I and 127Xe in xenon, out of which only 125Sb could potentially lead to background for a multi-ton scale dark matter search. The production rates for five out of eight studied radioactive isotopes in copper are in agreement with the only existing dedicated activation measurement, while we observe lower rates for the remaining ones. The specific saturation activities for both samples are also compared to predictions obtained with commonly used software packages, where we observe some underpredictions, especially for xenon activation.

Correlated cosmogenic W and Os isotopic variations in Carbo and implications for Hf-W chronology

An obstacle for establishing the chronology of iron meteorite formation using 182Hf-182W systematics (t1/2 = 8.9 Myr) is to find proper neutron fluence monitors to correct for cosmic ray modification of W isotopic composition. Recent studies showed that siderophile elements such as Pt and Os could serve such a purpose. To test and calibrate these neutron dosimeters, the isotopic compositions of W and Os were measured in a slab of the IID iron meteorite Carbo. This slab has a well-characterized noble gas depth profile reflecting different degrees of shielding to cosmic rays. The results show that W and Os isotopic ratios correlate with distance from the pre-atmospheric center. Negative correlations, barely resolved within error, were found between epsilo190Os-epsilo189Os and epsilo186Os-epsilo189Os with slopes of -0.64 ± 0.45 and -1.8(+1.9/-2.1), respectively. These Os isotope correlations broadly agree with model predictions for capture of secondary neutrons produced by cosmic ray irradiation and results reported previously for other groups of iron meteorites. Correlations were also found between epsilo182W-epsilo189Os (slope = 1.02 ± 0.37) and epsilo182W-epsilo190Os (slope = -1.38 ± 0.58). Intercepts of these two correlations yield pre-exposure epsilo182W values of -3.32 ± 0.51 and -3.62 ± 0.23, respectively (weighted average epsilo182W = -3.57 ± 0.21). This value relies on a large extrapolation leading to a large uncertainty but gives a metal-silicate segregation age of -0.5 ± 2.4 Myr after formation of the solar system. Combining the iron meteorite measurements with simulations of cosmogenic effects in iron meteorites, equations are presented to calculate and correct for cosmogenic effects on 182W using Os isotopes.

Radiocarbon age determination of 16 sites

The thermal diffusion enrichment apparatus in use in Amsterdam before 1967, has been rebuilt in the Groningen Radiocarbon Dating Laboratory. It has been shown to operate reliably and reproducibly. A reasonable agreement exists between the theoretical calculations and the experimental results. The 14C enrichment of a CO sample is deduced from the simultaneous mass 30 enrichment, which is measured with a mass spectrometer. The relation between both enrichments follows from a series of calibration measurements. The over-all accuracy in the enrichment is a few percent, equivalent to a few hundred years in age. The main problem in dating very old samples is their possible contamination with recent carbon. Generally, careful sample selection and rigorous pretreatment reduce sample contamination to an acceptable value. Also, it has been established that laboratory contamination, due to a memory effect in the combustion system and to impurities in the oxygen and nitrogen gas used for combustion, can be eliminated. A detailed analysis shows that the counter background in our set-up is almost exclusively caused by cosmic ray muons. The measurement of 28 early glacial samples, mostly from North-west Europe, has yielded a consistent set of ages. These indicate the existence of three early glacial interstadials; using the Weichselian definitions: Amersfoort starting at 68 200 ± 1100, Brørup at 64 400 ± 800 and Odderade at 60 500 ± 600 years BP. This 14C chronology shows good agreement with the Camp Century chronology and the dated palaeo sea levels. The discrepancy in the age of the early part of the Last Glacial on the 14C time scale and on that adopted for the deep-sea d18 record, must probably be attributed to the use of a generalized d18 curve and a wrong interpretation of this curve in terms of three Barbados terraces.

A 18-years long (1993-2011) snow and meteorological dataset from a mid-altitude mountain site (Col de Porte, France, 1325 altitude)

A quality-controlled snow and meteorological dataset spanning the period 1 August 1993-31 July 2011 is presented, originating from the experimental station Col de Porte (1325 m altitude, Chartreuse range, France). Emphasis is placed on meteorological data relevant to the observation and modelling of the seasonal snowpack. In-situ driving data, at the hourly resolution, consist of measurements of air temperature, relative humidity, windspeed, incoming short-wave and long-wave radiation, precipitation rate partitioned between snow- and rainfall, with a focus on the snow-dominated season. Meteorological data for the three summer months (generally from 10 June to 20 September), when the continuity of the field record is not warranted, are taken from a local meteorological reanalysis (SAFRAN), in order to provide a continuous and consistent gap-free record. Data relevant to snowpack properties are provided at the daily (snow depth, snow water equivalent, runoff and albedo) and hourly (snow depth, albedo, runoff, surface temperature, soil temperature) time resolution. Internal snowpack information is provided from weekly manual snowpit observations (mostly consisting in penetration resistance, snow type, snow temperature and density profiles) and from a hourly record of temperature and height of vertically free ''settling'' disks. This dataset has been partially used in the past to assist in developing snowpack models and is presented here comprehensively for the purpose of multi-year model performance assessment. The data is placed on the PANGAEA repository (doi:10.1594/PANGAEA.774249) as well as on the public ftp server ftp://ftp-cnrm.meteo.fr/pub-cencdp/.

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, (sec theta)(max), sensitive to the mass composition of cosmic rays above 3 x 10(18) eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modeling that must be resolved before the mass composition can be inferred from (sec theta)(max).

Search for isotropic γ radiation in the cosmological window between 65 and 200 TeV

Electromagnetic energy injected into the universe above a few hundred TeV is expected to pile up as γ radiation in a relatively narrow energy interval below 100 TeV due to its interaction with the 2.7^°K background radiation. We present an upper limit (90% C.L.) on the ratio of primary γ to charged cosmic rays in the energy interval 65–160 TeV (80–200 TeV) of 10.3 • 10^−3 (7.8 • 10^−3). Data from the HEGRA cosmic-ray detector complex consisting of a wide angle Čerenkov array (AIROBICC) measuring the lateral distribution of air Čerenkov light and a scintillator array, were used with a novel method to discriminate γ-ray and hadron induced air showers. If the presently unmeasured universal far infrared background radiation is not too intense, the result rules out a topological-defect origin of ultrahigh energy cosmic rays for masses of the X particle released by the defects equal to or larger than about 10^16 GeV.