An implementation of radiative transfer in the cosmological simulation code gadget

We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code gadget . It is based on a fast, robust and photon-conserving integration scheme where the radiation transport problem is approximated in terms of moments of the transfer equation and by using a variable Eddington tensor as a closure relation, following the Optically Thin Variable Eddington Tensor suggestion of Gnedin & Abel. We derive a suitable anisotropic diffusion operator for use in the SPH discretization of the local photon transport, and we combine this with an implicit solver that guarantees robustness and photon conservation. This entails a matrix inversion problem of a huge, sparsely populated matrix that is distributed in memory in our parallel code. We solve this task iteratively with a conjugate gradient scheme. Finally, to model photon sink processes we consider ionization and recombination processes of hydrogen, which is represented with a chemical network that is evolved with an implicit time integration scheme. We present several tests of our implementation, including single and multiple sources in static uniform density fields with and without temperature evolution, shadowing by a dense clump and multiple sources in a static cosmological density field. All tests agree quite well with analytical computations or with predictions from other radiative transfer codes, except for shadowing. However, unlike most other radiative transfer codes presently in use for studying re-ionization, our new method can be used on-the-fly during dynamical cosmological simulation, allowing simultaneous treatments of galaxy formation and the re-ionization process of the Universe.     

INTEGRAL/XMM views on the MeV source GRO J1411-64

The COMPTEL unidentified source GRO J 1411-64 was observed by INTEGRAL and XMM-Newton in 2005. The Circinus Galaxy is the only source detected within the 4σ location error of GRO J1411-64, but in here excluded as the possible counterpart. At soft X-rays, 22 reliable and statistically significant sources (likelihood >10) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6-635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. At the best location of the source, detections at hard X-rays show only upper limits, which, together with MeV results obtained by COMPTEL suggest the existence of a peak in power output located somewhere between 300–700 keV for the so-called low state. Such a spectrum resembles those in blazars or microquasars, and might suggest at work by the models accordingly. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons, shows that it is hard to comply with all observational constrains. This fact and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, what is discussed here.     

The fine-grained phase-space structure of cold dark matter haloes

We present a new and completely general technique for calculating the fine-grained phase-space structure of dark matter (DM) throughout the Galactic halo. Our goal is to understand this structure on the scales relevant for direct and indirect detection experiments. Our method is based on evaluating the geodesic deviation equation along the trajectories of individual DM particles. It requires no assumptions about the symmetry or stationarity of the halo formation process. In this paper we study general static potentials which exhibit more complex behaviour than the separable potentials studied previously. For ellipsoidal logarithmic potentials with a core, phase mixing is sensitive to the resonance structure, as indicated by the number of independent orbital frequencies. Regions of chaotic mixing can be identified by the very rapid decrease in the real-space density of the associated DM streams. We also study the evolution of stream-density in ellipsoidal NFW haloes with radially varying isopotential shape, showing that if such a model is applied to the Galactic halo, at least 10⁵ streams are expected near the Sun. The most novel aspect of our approach is that general non-static systems can be studied through implementation in a cosmological N -body code. Such an implementation allows a robust and accurate evaluation of the enhancements in annihilation radiation due to fine-scale structure such as caustics. We embed the scheme in the current state-of-the-art code gadget -3 and present tests which demonstrate that N -body discreteness effects can be kept under control in realistic configurations.     

Spatial variability of phytoplankton, nutrients and new production estimates in the waters around Svalbard

Phytoplankton dynamics and carbon input into Arctic and sub-Arctic ecosystems were investigated around Svalbard, in summer 1991. Phytoplankton biomass, species composition and dissolved nutrient concentrations were analysed from water samples collected along seven transects. Phytoplankton biomass was low especially to the north (Chlorophyll-a mean 0.3 pg 1- '), where flagellates dominated the communities and only ice-diatoms were present. To the west, the phytoplankton composition was representative of a summer Atlantic community, in a post-bloom state. Zooplankton grazing, mainly by copepods, appeared to be the main control on biomass to the west and north of Svalbard.
In the Barents Sea (east of Svalbard), an ice edge bloom was observed (Chlorophyll-a max. 6.8 pgl-') and the ice edge receded at a rate of approximately 1 1 km day-'. The phytoplankton community was represented by marginal ice species, especially Phaeocystis poucherii and Chaeroceros socialis. South of the ice edge, Deep Chlorophyll Maxima (DCM) were observed, as surface waters became progressively nutrient-depleted. In these surface waters, the phytoplankton were predominantly auto- and heterotrophic flagellates.
Carbon production measurements revealed high net production (new and regenerated) to the north of the Barents Sea Polar Front (BSPF); it was especially high at the receding ice edge (reaching 1.44gC m-'day-'). To the south, a low level of production was maintained, mainly through regenerative processes.     

Stable carbon isotopic ratios of CH4-CO2-bearing fluid inclusions in fracture-fill mineralization from the Lower Saxony Basin (Germany) - A tool for tracing gas sources and maturity

The stable carbon isotopic ratios (δ¹³C) of methane (CH₄) and carbon dioxide (CO₂) of gas-rich fluid inclusions hosted in fracture-fill mineralization from the southern part of the Lower Saxony Basin, Germany have been measured online using a crushing device interfaced to an isotopic ratio mass spectrometer (IRMS). The data reveal that CH₄ trapped in inclusions seems to be derived from different source rocks with different organic matter types. The δ¹³C values of CH₄ in inclusions in quartz hosted by Carboniferous rocks range between −25 and −19‰, suggesting high-maturity coals as the source of methane. Methane in fluid inclusions in minerals hosted by Mesozoic strata has more negative carbon isotope ratios (−45 to −31‰) and appears to represent primary cracking products from type II kerogens, i.e., marine shales. There is a positive correlation between increasing homogenization temperatures of aqueous fluid inclusions and less negative δ¹³C(CH₄) values of in co-genetic gas inclusions probably indicating different mtaturity of the potential source rocks at the time the fluids were released. The CO₂ isotopic composition of CH₄-CO₂-bearing inclusions shows slight negative or even positive δ¹³C values indicating an inorganic source (e.g., water-rock interaction and dissolution of detrital, marine calcite) for CO₂ in inclusions. We conclude that the δ¹³C isotopic ratios of CH₄-CO₂-bearing fluid inclusions can be used to trace migration pathways, sources of gases, and alteration processes. Furthermore, the δ¹³C values of methane can be used to estimate the maturity of the rocks from which it was sourced. Results presented here are further supported by organic geochemical analysis of surface bitumens which coexist with the gas inclusion-rich fracture-fill mineralization and confirm the isotopic interpretations with respect to fluid source, type and maturity.     

Geoelektrische Untersuchungen zur Ermittlung der Blitzgefährdung

Zusammenfassung Aus den Untersuchungen, die vom Verf. in Oesterreich ausgeführt wurden, werden folgende Ergebnisse bewiesen: 1. Die elektrische Bodenleitfähigkeit beinflusst massgebend den erdnahen Teil der Blitzbahn; 2. Es gibt daher Zonen hoher Blitzgefährdung, die Blitznester genannt werden; 3. Die Scheitelstromstärke des Blitzes ist durch den spezifischen Bodenwiderstand mitbestimmt; 4. Es besteht ein Zusammenhang zwischen Blitzgefährdung und geoelektrischer Bodenbeschaffenheit.

---

Summary Experiments carried out by the Author in Austria show the following results: 1. The electrical conductivity of the grounds has a decisive influence on the course of the lightning: such zones are called «lightning nests». 3. The amplitude of the strength of lightning current is determined also by the specific resistivity of ground. 4. The agreement between the position of lightning nests and the geoelectric condition of ground is beyond dispute.

     

3D finite volume groundwater and heat transport modeling with non-orthogonal grids,using a coordinate transformation method

Many popular groundwater modeling codes are based on the finite differences or finite volume method for orthogonal grids. In cases of complex subsurface geometries this type of grid either leads to coarse geometric representations or to extremely fine meshes. We use a coordinate transformation method (CTM) to circumvent this shortcoming. In computational fluid dynamics (CFD), this method has been applied successfully to the general Navier–Stokes equation. The method is based on tensor analysis and performs a transformation of a curvilinear into a rectangular unit grid, on which a modified formulation of the differential equations is applied. Therefore, it is not necessary to reformulate the code in total. We applied the CTM to an existing three-dimensional code (SHEMAT), a simulator for heat conduction and advection in porous media. The finite volume discretization scheme for the non-orthogonal, structured, hexahedral grid leads to a 19-point stencil and a correspondingly banded system matrix. The implementation is straightforward and it is possible to use some existing routines without modification. The accuracy of the modified code is demonstrated for single phase flow on a two-dimensional analytical solution for flow and heat transport. Additionally, a simple case of potential flow is shown for a two-dimensional grid which is increasingly deformed. The result reveals that the corresponding error increases only slightly. Finally, a thermal free-convection benchmark is discussed. The result shows, that the solution obtained with the new code is in good agreement with the ones obtained by other codes.     

Direct dating of hydrothermal W mineralization: U-Pb age for hübnerite (MnWO4), Sweet Home Mine, Colorado

We have investigated the potential of hübnerite for U-Pb dating. Hübnerite forms typically at medium to low-temperatures in a wide range of pneumatolytic-hydrothermal mineral deposits, particularly porphyry molybdenum and Sn-specialized granites. Hübnerite from the Sweet Home Mine (Alma, Colorado) formed in a Pb-rich, U-poor environment, but still developed relatively radiogenic Pb isotopic compositions. The low Pb_common contents in hübnerite (0.075 to 0.155 ppm) demonstrate that Pb is efficiently excluded from the crystal lattice. In contrast, U may substitute for Mn. The U-Pb data of hübnerite scatter. Most of the scatter originates from samples with ²⁰⁶Pb/²⁰⁴Pb values below 50, where Pb_blank contributes up to 30% to Pb_total. Using the least radiogenic galena Pb, samples with ²⁰⁶Pb/²⁰⁴Pb values above 70 have overlapping ²⁰⁶Pb∗/²³⁸U and ²⁰⁷Pb∗/²³⁵U values and yield a ²⁰⁶Pb/²³⁸U age of 25.7 ± 0.3 Ma (2σ). Late stage apatite from the Sweet Home Mine yields a ²⁰⁶Pb/²⁰⁴Pb-²³⁸U/²⁰⁴Pb isochron corresponding to an age of 24.8 ± 0.5 Ma (2σ). A comparison of the U-Pb hübnerite ages with literature ⁴⁰Ar/³⁹Ar ages on earlier sericite and the U-Pb age on later apatite suggests that (i) hübnerite yields accurate U-Pb ages and (ii) the evolution of the Sweet Home mineralization from greisen-type mineralization to medium-temperature hydrothermal vein mineralization took place in a few hundred thousand years at most. Aqueous low-N₂-bearing and aqueous inclusions in the dated hübnerite have homogenization temperatures between 325 and 356 °C and moderate salinity (up to 6.7 wt% NaCl equiv.). Thus, hübnerite represents one of the rare examples of a mineral that can be dated accurately and carries petrological information.