Simulation of the interaction of galactic cosmic‐ray protons with meteoroids: On the production of radionuclides in thick gabbro and iron targets irradiated isotropically with 1.6 GeV protons

Abstract— Thick spherical targets made of gabbro (R = 25 cm) and of steel (R = 10 cm) were irradiated isotropically with 1.6 GeV protons at the Saturne synchrotron at Laboratoire National Saturne (LNS)/CEN Saclay in order to simulate the interaction in space of galactic cosmic‐ray (GCR) protons with stony and iron meteoroids. Proton fluences of 1.32 × 10¹⁴ cm^?2 and 2.45 × 10¹⁴ cm^?2 were received by the gabbro and iron sphere, respectively, which corresponds to cosmic‐ray exposure ages of about 1.6 and 3.0 Ma. Both artificial meteoroids contained large numbers of high‐purity target foils of up to 28 elements at different depths. In these individual target foils, elementary production rates of radionuclides and rare gas isotopes were measured by x‐ and γ‐spectrometry, by low‐level counting, accelerator mass spectrometry (AMS), and by conventional rare gas mass spectrometry. Also samples of the gabbro itself were analyzed. Up to now, for each of the experiments, ～500 target‐product combinations were investigated of which the results for radionuclides are presented here. The experimental production rates show a wide range of depth profiles reflecting the differences between low‐, medium‐, and high‐energy products. The influence of the stony and iron matrices on the production of secondary particles and on particle transport, in general, and consequently on the production rates is clearly exhibited by the phenomenology of the production rates as well as by a detailed theoretical analysis. Theoretical production rates were calculated in an a priori way by folding depth‐dependent spectra of primary and secondary protons and secondary neutrons calculated by Monte Carlo techniques with the excitation functions of the underlying nuclear reactions. Discrepancies of up to a factor of 2 between the experimental and a priori calculated depth profiles are attributed to the poor quality of the mostly theoretical neutron excitation functions. Improved neutron excitation functions were obtained by least‐squares deconvolution techniques from experimental thick‐target production rates of up to five thick‐target experiments in which isotropic irradiations were performed. A posteriori calculations using the adjusted neutron cross sections describe the measured depth profiles of all these simulation experiments within 9%. The thus validated model calculations provide a basis for reliable physical model calculations of the production rates of cosmogenic nuclides in stony and iron meteorites as well as in lunar samples and terrestrial materials.     

Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging, and studied the effect of single- and multiple-sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionization structures including appropriate microphysics. At low densities, the grains do not play any significant rôle in the shock dynamics. At high densities, the ionization fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allow a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.     

Structure and kinematics of the interstellar medium around WR 134 and WR 135

Based on our H _α interferometry and 21-cm and CO observations, we analyze the structure and kinematics of the interstellar medium around the stars WR 134 and WR 135. We conclude that the HI bubble found here previously is associated with WR 135, not with WR 134. High-velocity motions of ionized gas that can be interpreted as expansion of the gas swept up by the stellar wind with a velocity up to 50–80 km s^?1 are observed around both stars. The line-of-sight velocity field of the ionized hydrogen in the Cygnus arm is shown to agree with the large-scale line-of-sight velocity distribution of the CO emission.     

Early evolution of newly born magnetars with a strong toroidal field

We present a state-of-the-art scenario for newly born magnetars as strong sources of gravitational waves (GWs) in the early days after formation. We address several aspects of the astrophysics of rapidly rotating, ultra-magnetized neutron stars (NSs), including early cooling before transition to superfluidity, the effects of the magnetic field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, oblique rotator and the strength of the electromagnetic torque on the newly born NS. We show that our scenario is consistent with recent studies of supernova remnant surrounding Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) in the Galaxy that constrains the electromagnetic energy input from the central NS to be  ≤ 10⁵¹  erg. We further show that if this condition is met, then the GW signal from such sources is potentially detectable with the forthcoming generation of GW detectors up to Virgo cluster distances where an event rate ∼1 yr⁻¹ can be estimated. Finally, we point out that the decay of an internal magnetic field in the 10¹⁶ G range couples strongly with the NS cooling at very early stages, thus significantly slowing down both processes: the field can remain this strong for at least 10³ yr, during which the core temperature stays higher than several times 10⁸ K.     

Trace constituents of Europa's atmosphere

Europa is bombarded by intense radiation that erodes the surface, launching molecules into a thin “atmosphere” representative of surface composition. In addition to atoms and molecules created in the mostly water ice surface such as H₂O, O₂, H₂, the atmosphere is known to have species representative of trace surface materials. These trace species are carried off with the 10-10⁴ H₂O molecules ejected by each energetic heavy ion, a process we have simulated using molecular dynamics. Using the results of those simulations, we found that a neutral mass spectrometer orbiting ∼100 km above the surface could detect species with surface concentrations above ∼0.03%. We have also modeled the atmospheric spatial structure of the volatile species CO₂ and SO₂ under a variety of assumptions. Detections of these species with moderate time and space resolution would allow us to constrain surface composition, chemistry and to study space weathering processes.     

The interior structure of Mercury and its core sulfur content

Depth-dependent interior structure models of Mercury are calculated for several plausible chemical compositions of the core and of the mantle. For those models, we compute the associated libration amplitude, obliquity, tidal deformation, and tidal changes in the external potential. In particular we study the relation between the interior structure parameters for five different mantle mineralogies and two different temperature profiles together with two extreme crust density values. We investigate the influence of the core light element concentration, temperature, and melting law on core state and inner and outer core size. We show that a sulfur concentration above 10 wt% is unlikely if the temperature at the core-mantle boundary is above 1850 K and the silicate shell at least 240 km thick. The interior models can only have an inner core if the sulfur weight fraction is below 5 wt% for core-mantle boundary temperature in the 1850-2200 K range. Within our modeling hypotheses, we show that with the expected precision on the moment of inertia the core size can be estimated to a precision of about 50 km and the core sulfur concentration with an error of about 2 wt%. This uncertainty can only be reduced when more information on the mantle mineralogy of Mercury becomes available. However, we show that the uncertainty on the core size estimation can be greatly reduced, to about 25 km, if tidal surface displacements and tidal variations in the external potential are considered.