Simulations of multiphase turbulence in jet cocoons

The interaction of optically emitting clouds with warm X-ray gas and hot, tenuous radio plasma in radio jet cocoons is modelled by 2D compressible hydrodynamic simulations. The initial setup is the Kelvin–Helmholtz instability at a contact surface of density contrast 10⁴. The denser medium contains clouds of higher density. Optically thin radiation is realized via a cooling source term. The cool phase effectively extracts energy from the other gas which is both, radiated away and used for acceleration of the cold phase. This increases the system's cooling rate substantially and leads to a massively amplified cold mass dropout. We show that it is feasible, given small seed clouds of the order of  100 M_⊙  , that all of the optically emitting gas in a radio jet cocoon may be produced by this mechanism on the propagation time-scale of the jet. The mass is generally distributed as   T ^−1/2  with temperature, with a prominent peak at 14 000 K. This peak is likely to be related to the counteracting effects of shock heating and a strong rise in the cooling function. The volume filling factor of cold gas in this peak is of the order of  10⁻⁵–10⁻³  and generally increases during the simulation time.
The simulations tend towards an isotropic scale-free Kolmogorov-type energy spectrum over the simulation time-scale. We find the same Mach-number density relation as Kritsuk & Norman and show that this relation may explain the velocity widths of emission lines associated with high-redshift radio galaxies, if the environmental temperature is lower, or the jet-ambient density ratio is less extreme than in their low-redshift counterparts.     

26Al‐26Mg systematics of Ca‐Al‐rich inclusions,amoeboid olivine aggregates,and chondrules from the ungrouped carbonaceous chondrite Acfer 094

Abstract— We report in situ magnesium isotope measurements of 7 porphyritic magnesium‐rich (type I) chondrules, 1 aluminum‐rich chondrule, and 16 refractory inclusions (14 Ca‐Al‐rich inclusions [CAIs] and 2 amoeboid olivine aggregates [AOAs]) from the ungrouped carbonaceous chondrite Acfer 094 using a Cameca IMS 6f ion microprobe. Both AOAs and 9 CAIs show radiogenic ²⁶Mg excesses corresponding to initial ²⁶Al/²⁷Al ratios between ～5 × 10^?5 ～7 × 10^?5 suggesting that formation of the Acfer 094 CAIs may have lasted for ～300,000 years. Four CAIs show no evidence for radiogenic ²⁶Mg; three of these inclusions (a corundum‐rich, a grossite‐rich, and a pyroxene‐hibonite spherule CAI) are very refractory objects and show deficits in ²⁶Mg, suggesting that they probably never contained ²⁶Al. The fourth object without evidence for radiogenic ²⁶Mg is an anorthite‐rich, igneous (type C) CAI that could have experienced late‐stage melting that reset its Al‐Mg systematics. Significant excesses in ²⁶Mg were observed in two chondrules. The inferred ²⁶Al/²⁷Al ratios in these two chondrules are (10.3 ± 7.4) × 10^?6 (6.0 ± 3.8) × 10^?6 (errors are 2σ), suggesting formation 1.6^+1.2_‐0.6 and 2.2^+0.4_‐0.3 Myr after CAIs with the canonical ²⁶Al/²⁷Al ratio of 5 × 10^?5. These age differences are consistent with the inferred age differences between CAIs and chondrules in primitive ordinary (LL3.0–LL3.1) and carbonaceous (CO3.0) chondrites.     

Evolution of a Keplerian disk of colliding and fragmenting particles: a kinetic model with application to the Edgeworth-Kuiper belt

We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune.     

Remelting of refractory inclusions in the chondrule‐forming regions: Evidence from chondrule‐bearing type C calcium‐aluminum‐rich inclusions from Allende

Abstract— We describe the mineralogy, petrology, oxygen, and magnesium isotope compositions of three coarse‐grained, igneous, anorthite‐rich (type C) Ca‐Al‐rich inclusions (CAIs) (ABC, TS26, and 93) that are associated with ferromagnesian chondrule‐like silicate materials from the CV carbonaceous chondrite Allende. The CAIs consist of lath‐shaped anorthite (An₉₉), Cr‐bearing Al‐Ti‐diopside (Al and Ti contents are highly variable), spinel, and highly åkermanitic and Na‐rich melilite (Åk_63–74, 0.4–0.6 wt% Na₂O). TS26 and 93 lack Wark‐Lovering rim layers; ABC is a CAI fragment missing the outermost part. The peripheral portions of TS26 and ABC are enriched in SiO₂ and depleted in TiO₂ and Al₂O₃ compared to their cores and contain relict ferromagnesian chondrule fragments composed of forsteritic olivine (Fa_6–8) and low‐Ca pyroxene/pigeonite (Fs₁Wo_1–9). The relict grains are corroded by Al‐Ti‐diopside of the host CAIs and surrounded by haloes of augite (Fs_0.5Wo_30–42). The outer portion of CAI 93 enriched in spinel is overgrown by coarse‐grained pigeonite (Fs_0.5–2Wo_5–17), augite (Fs_0.5Wo_38–42), and anorthitic plagioclase (An₈₄). Relict olivine and low‐Ca pyroxene/pigeonite in ABC and TS26, and the pigeonite‐augite rim around 93 are ¹⁶O‐poor (Δ17O ～ ?1‰ to ?8‰). Spinel and Al‐Ti‐diopside in cores of CAIs ABC, TS26, and 93 are ¹⁶O‐enriched (Δ17O down to ?20‰), whereas Al‐Ti‐diopside in the outer zones, as well as melilite and anorthite, are ¹⁶O‐depleted to various degrees (Δ17O = ?11‰ to 2‰). In contrast to typical Allende CAIs that have the canonical initial ²⁶Al/²⁷Al ratio of ～5 × 10^?5 ABC, 93, and TS26 are ²⁶Al‐poor with (²⁶Al/²⁷Al)₀ ratios of (4.7 ± 1.4) × 10^?6 (1.5 ± 1.8) × 10^?6 <1.2 × 10^?6 respectively. We conclude that ABC, TS26, and 93 experienced remelting with addition of ferromagnesian chondrule silicates and incomplete oxygen isotopic exchange in an ¹⁶O‐poor gaseous reservoir, probably in the chondrule‐forming region. This melting episode could have reset the ²⁶Al‐²⁶Mg systematics of the host CAIs, suggesting it occurred ～2 Myr after formation of most CAIs. These observations and the common presence of relict CAIs inside chondrules suggest that CAIs predated formation of chondrules.     

Grenvillian and Caledonian evolution of eastern Svalbard – a tale of two orogenies

Svalbard is located in the north-west corner of the Barents Sea shelf and the Eurasian Plate, in a key area for interpreting Caledonian and older orogens in the Arctic region. Recent U–Pb dating in the Nordaustlandet Terrane of eastern Svalbard shows this terrane to consist of a Grenville-age basement, overlain by Neoproterozoic to early Palaeozoic platformal sediments, and intruded by Caledonian anatectic granites. Deformation, metamorphism and crustal anatectic magmatism occurred both during the Grenvillian (960–940 Ma) and Caledonian (450–410 Ma) orogenies. This evolution shows great similarities with that of eastern Greenland. In the classical model, eastern Svalbard is placed outboard of central east Greenland in pre-Caledonian time. Alternatively, it may have been located north-east of Greenland and transferred west and rotated anticlockwise during Caledonian continent–continent collision. In the Neoproterozoic, easternmost Svalbard may have been part of a wider area of Grenville-age crust, now fragmented and dispersed around the Arctic.     

Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

Preliminary studies concerning subsurface probes for the exploration of icy planetary bodies

It is well-known that the permanent terrestrial ice sheets (glaciers and polar caps) contain a lot of information about the recent geological history and in particular about climatic changes. Extrapolating this fact to other ice sheets in the solar system (e.g. the Mars polar regions, the icy moons of the outer planets, etc.), we may expect a similar wealth of information. To obtain this information it is possible to drill holes or melt the ice by a heated probe, which in this way is able to penetrate the surface and investigate the deeper layers in situ. In the latter case the driving agent is the heating power and the weight of the probe. In this paper we consider the application of such “melting probes” for exploring the structure of ice sheets in extraterrestrial environments. We describe several laboratory experiments with simple melting probes performed under cryo-vacuum conditions and compare the results with tests in a terrestrial environment. The experiments revealed that under space conditions the downward motion of a heated probe in an ice sheet is characterized by intermittent periods of sublimation and melting of the surrounding ice, sometimes interrupted by periods where a part of the probe's outer surface is frozen to the surrounding ice. This leads to a temporary blocking of the probe's downward motion. A similar situation can occur when the trailing tether is frozen in behind the probe. During the periods of ice sublimation the penetration process is significantly more power consuming, due to the large difference between the latent heat of sublimation and the latent heat of melting for water ice.