On-the-fly memory compression for multibody algorithms.

Memory and bandwidth demands challenge developers of particle-based codes that have to scale on new architectures, as the growth of concurrency outperforms improvements in memory access facilities, as the memory per core tends to stagnate, and as communication networks cannot increase bandwidth arbitrary. We propose to analyse each particle of such a code to find out whether a hierarchical data representation storing data with reduced precision caps the memory demands without exceeding given error bounds. For admissible candidates, we perform this compression and thus reduce the pressure on the memory subsystem, lower the total memory footprint and reduce the data to be exchanged via MPI. Notably, our analysis and transformation changes the data compression dynamically, i.e. the choice of data format follows the solution characteristics, and it does not require us to alter the core simulation code

Simulating ice core 10Be on the glacial–interglacial timescale

10Be ice core measurements are an important tool for paleoclimate research, e.g., allowing for the reconstruction of past solar activity or changes in the geomagnetic dipole field. However, especially on multi-millennial timescales, the share of production and climate-induced variations of respective 10Be ice core records is still up for debate. Here we present the first quantitative climatological model of the 10Be ice concentration up to the glacial–interglacial timescale. The model approach is composed of (i) a coarse resolution global atmospheric transport model and (ii) a local 10Be air–firn transfer model. Extensive global-scale observational data of short-lived radionuclides as well as new polar 10Be snow-pit measurements are used for model calibration and validation. Being specifically configured for 10Be in polar ice, this tool thus allows for a straightforward investigation of production- and non-production-related modulation of this nuclide. We find that the polar 10Be ice concentration does not immediately record the globally mixed cosmogenic production signal. Using geomagnetic modulation and revised Greenland snow accumulation rate changes as model input, we simulate the observed Greenland Summit (GRIP and GISP2) 10Be ice core records over the last 75 kyr (on the GICC05modelext timescale). We show that our basic model is capable of reproducing the largest portion of the observed 10Be changes. However, model–measurement differences exhibit multi-millennial trends (differences up to 87% in case of normalized to the Holocene records) which call for closer investigation. Focusing on the (12–37) b2k (before the year AD 2000) period, mean model–measurement differences of 30% cannot be attributed to production changes. However, unconsidered climate-induced changes could likely explain the model–measurement mismatch. In fact, the 10Be ice concentration is very sensitive to snow accumulation changes. Here the reconstructed Greenland Summit (GRIP) snow accumulation rate record would require revision of +28% to solely account for the (12–37) b2k model–measurement differences

Settling the half-life of ⁶⁰Fe: fundamental for a versatile astrophysical chronometer

In order to resolve a recent discrepancy in the half-life of ⁶⁰Fe, we performed an independent measurement with a new method that determines the ⁶⁰Fe content of a material relative to Fe55 (t1/2=2.744yr) with accelerator mass spectrometry. Our result of (2.50±0.12)×10⁶yr clearly favors the recently reported value (2.62±0.04)×10⁶yr, and rules out the older result of (1.49±0.27)×10⁶yr. The present weighted mean half-life value of (2.60±0.05)×10⁶yr substantially improves the reliability as an important chronometer for astrophysical applications in the million-year time range. This includes its use as a sensitive probe for studying recent chemical evolution of our Galaxy, the formation of the early Solar System, nucleosynthesis processes in massive stars, and as an indicator of a recent nearby supernova.Part of this work was funded by the Austrian Science Fund (FWF) Projects No. AP20434 and AI00428 (FWF and CoDustMas, Eurogenesis via ESF)

The Search for Supernova-produced Radionuclides in Terrestrial Deep-sea Archives

An enhanced concentration of 60Fe was found in a deep ocean's crust in 2004 in a layer corresponding to an age of ~2 Myr. The confirmation of this signal in terrestrial archives as supernova-induced and detection of other supernova-produced radionuclides is of great interest. We have identified two suitable marine sediment cores from the South Australian Basin and estimated the intensity of a possible signal of the supernova-produced radionuclides 26Al, 53Mn, 60Fe and the pure r-process element 244Pu in these cores. A finding of these radionuclides in a sediment core might allow to improve the time resolution of the signal and thus to link the signal to a supernova event in the solar vicinity ~2 Myr ago. Furthermore, it gives an insight on nucleosynthesis scenarios in massive stars, the condensation into dust grains and transport mechanisms from the supernova shell into the solar system

Search for supernova-produced 60Fe in a marine sediment

An 60Fe peak in a deep-sea FeMn crust has been interpreted as due to the signature left by the ejecta of a supernova explosion close to the solar system 2.8 +/- 0.4 Myr ago [Knie et al., Phys. Rev. Lett. 93, 171103 (2004)]. To confirm this interpretation with better time resolution and obtain a more direct flux estimate, we measured 60Fe concentrations along a dated marine sediment. We find no 60Fe peak at the expected level from 1.7 to 3.2 Myr ago. However, applying the same chemistry used for the sediment, we confirm the 60Fe signal in the FeMn crust. The cause of the discrepancy is discussed.Comment: 15 pages, 5 figures, submitted to PR

From Gravitons to Gravity: Myths and Reality

There is a general belief, reinforced by statements in standard textbooks, that: (i) one can obtain the full non-linear Einstein's theory of gravity by coupling a massless, spin-2 field $h_{ab}$ self-consistently to the total energy momentum tensor, including its own; (ii) this procedure is unique and leads to Einstein-Hilbert action and (iii) it only uses standard concepts in Lorentz invariant field theory and does not involve any geometrical assumptions. After providing several reasons why such beliefs are suspect -- and critically re-examining several previous attempts -- we provide a detailed analysis aimed at clarifying the situation. First, we prove that it is \textit{impossible} to obtain the Einstein-Hilbert (EH) action, starting from the standard action for gravitons in linear theory and iterating repeatedly. Second, we use the Taylor series expansion of the action for Einstein's theory, to identify the tensor $\mathcal{S}^{ab}$, to which the graviton field $h_{ab}$ couples to the lowest order. We show that the second rank tensor $\mathcal{S}^{ab}$ is {\it not} the conventional energy momentum tensor $T^{ab}$ of the graviton and provide an explanation for this feature. Third, we construct the full nonlinear Einstein's theory with the source being spin-0 field, spin-1 field or relativistic particles by explicitly coupling the spin-2 field to this second rank tensor $\mathcal{S}^{ab}$ order by order and summing up the infinite series. Finally, we construct the theory obtained by self consistently coupling $h_{ab}$ to the conventional energy momentum tensor $T^{ab}$ order by order and show that this does {\it not} lead to Einstein's theory. (condensed).Comment: revtex; 19 pages; no figure

Sparling two-forms, the conformal factor and the gravitational energy density of the teleparallel equivalent of general relativity

It has been shown recently that within the framework of the teleparallel equivalent of general relativity (TEGR) it is possible to define the energy density of the gravitational field. The TEGR amounts to an alternative formulation of Einstein's general relativity, not to an alternative gravity theory. The localizability of the gravitational energy has been investigated in a number of space-times with distinct topologies, and the outcome of these analises agree with previously known results regarding the exact expression of the gravitational energy, and/or with the specific properties of the space-time manifold. In this article we establish a relationship between the expression for the gravitational energy density of the TEGR and the Sparling two-forms, which are known to be closely connected with the gravitational energy. We also show that our expression of energy yields the correct value of gravitational mass contained in the conformal factor of the metric field.Comment: 12 pages, Latex file, no figures, to be published in Gen. Rel. Gra

Hidden Consequence of Active Local Lorentz Invariance

In this paper we investigate a hidden consequence of the hypothesis that Lagrangians and field equations must be invariant under active local Lorentz transformations. We show that this hypothesis implies in an equivalence between spacetime structures with several curvature and torsion possibilities.Comment: Some misprints appearing in the published version have been correcte