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.

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.

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.

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.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

Personalized X-Ray Reconstruction of the Proximal Femur via Intensity-Based Non-rigid 2D-3D Registration

This paper presents a new approach for reconstructing a patient-specific shape model and internal relative intensity distribution of the proximal femur from a limited number (e.g., 2) of calibrated C-arm images or X-ray radiographs. Our approach uses independent shape and appearance models that are learned from a set of training data to encode the a priori information about the proximal femur. An intensity-based non-rigid 2D-3D registration algorithm is then proposed to deformably fit the learned models to the input images. The fitting is conducted iteratively by minimizing the dissimilarity between the input images and the associated digitally reconstructed radiographs of the learned models together with regularization terms encoding the strain energy of the forward deformation and the smoothness of the inverse deformation. Comprehensive experiments conducted on images of cadaveric femurs and on clinical datasets demonstrate the efficacy of the present approach.

Hardware accelerated ray cast of volume data and volume gradient for an optimized splines-based multi-resolution 2D-3D registration

This paper describes a method for DRR generation as well as for volume gradients projection using hardware accelerated 2D texture mapping and accumulation buffering and demonstrates its application in 2D-3D registration of X-ray fluoroscopy to CT images. The robustness of the present registration scheme are guaranteed by taking advantage of a coarse-to-fine processing of the volume/image pyramids based on cubic B-splines. A human cadaveric spine specimen together with its ground truth was used to compare the present scheme with a purely software-based scheme in three aspects: accuracy, speed, and capture ranges. Our experiments revealed an equivalent accuracy and capture ranges but with much shorter registration time with the present scheme. More specifically, the results showed 0.8 mm average target registration error, 55 second average execution time per registration, and 10 mm and 10° capture ranges for the present scheme when tested on a 3.0 GHz Pentium 4 computer.

In situ dehydration behavior of zeolite-like cavansite: A single-crystal X-ray study

To track dehydration behavior of cavansite, Ca(VO)(Si4O10)·4H2O space group Pnma, a = 9.6329(2), b = 13.6606(2), c = 9.7949(2) Å, V = 1288.92(4) Å3 single-crystal X-ray diffraction data on a crystal from Wagholi quarry, Poona district (India) were collected up to 400 °C in steps of 25 °C up to 250 °C and in steps of 50 °C between 250 and 400 °C. The structure of cavansite is characterized by layers of silicate tetrahedra connected by V4+O5 square pyramids. This way a porous framework structure is formed with Ca and H2O as extraframework occupants. At room temperature, the hydrogen bond system was analyzed. Ca is eightfold coordinated by four bonds to O of the framework structure and four bonds to H2O molecules. H2O linked to Ca is hydrogen bonded to the framework and also to adjacent H2O molecules. The dehydration in cavansite proceeds in four steps.At 75 °C, H2O at O9 was completely expelled leading to 3 H2O pfu with only minor impact on framework distortion and contraction V = 1282.73(3) Å3. The Ca coordination declined from originally eightfold to sevenfold and H2O at O7 displayed positional disorder.At 175 °C, the split O7 sites approached the former O9 position. In addition, the sum of the three split positions O7, O7a, and O7b decreased to 50% occupancy yielding 2 H2O pfu accompanied by a strong decrease in volume V = 1206.89(8) Å3. The Ca coordination was further reduced from sevenfold to sixfold.At 350 °C, H2O at O8 was released leading to a formula with 1 H2O pfu causing additional structural contraction (V = 1156(11) Å3). At this temperature, Ca adopted fivefold coordination and O7 rearranged to disordered positions closer to the original O9 H2O site.At 400 °C, cavansite lost crystallinity but the VO2+ characteristic blue color was preserved. Stepwise removal of water is discussed on the basis of literature data reporting differential thermal analyses, differential thermo-gravimetry experiments and temperature dependent IR spectra in the range of OH stretching vibrations.