TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

The transport of cosmic rays in self-excited magnetic turbulence

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic-ray transport and acceleration in supernova remnants.

A cosmic ray current-driven instability in partially ionised media

Context. We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock.

Aims. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which temperature and ionisation fractions modify this picture.

Methods. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a streaming cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions.

Results. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic field, and determines the range of possible growth rates.

Conclusions. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth of the non-resonant mode. For typical supernova parameters in our Galaxy, thermal effects do not significantly alter the growth rates. For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes can be rapidly damped.

ENVIRONMENTAL LIMITS ON THE NONRESONANT COSMIC-RAY CURRENT-DRIVEN INSTABILITY

We investigate the so-called nonresonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the nonresonant mode in SN1006 and RX J1713.7-3946 are presented.

On the maximum energy of shock-accelerated cosmic rays at ultra-relativistic shocks

The maximum energy to which cosmic rays can be accelerated at weakly magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropized in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganized field is required on larger scales. The growth of cosmic ray-induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly magnetized ultra-relativistic shocks are disfavoured as high-energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.

Irradiation of water ice by C+ ions in the cosmic environment

We present a first principles molecular dynamics (FPMD) study of the interaction of low energy, positively charged, carbon (C+) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 eV and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the neutral isoformyl radical, COH, and carbon monoxide, CO, are the dominant products of these reactions. All these reactions are accompanied by the transfer of a proton from the reacting water molecule to the ice, where it forms a hydronium ion. We find that COH is formed either via a direct, "knock-out", mechanism following the impact of the C+ projectile upon a water molecule or by creation of a COH_2^+ intermediate. The direct mechanism is more prominent at higher energies. CO is generally produced following the dissociation of COH. More frequent production of the formyl radical, HCO, is observed here than in gas phase calculations. A less commonly occurring product is the dihydroxymethyl, CH(OH)_2, radical. Although a minor result, its existence gives an indication of the increasing chemical complexity which is possible in such heterogeneous environments.

Astrophysics: Cosmic explosions in the young Universe

The discovery of two superluminous supernovae at large distances from Earth pushes the frontier of supernova studies to just 1.5 billion years after the Big Bang, and suggests that they may be common in the young Universe. See Letter p.228

Recent Developments in the Analysis of the Black Mat Layer and Cosmic Impact at 12.8ka

Recent analyses of sediment samples from "black mat" sites in South America and Europe support previous interpretations of an ET impact event that reversed the Late Glacial demise of LGM ice during the Bølling Allerød warming, resulting in a resurgence of ice termed the Younger Dryas (YD) cooling episode. The breakup or impact of a cosmic vehicle at the YD boundary coincides with the onset of a 1-kyr long interval of glacial resurgence, one of the most studied events of the Late Pleistocene. New analytical databases reveal a corpus of data indicating that the cosmic impact was a real event, most possibly a cosmic airburst from Earth's encounter with the Taurid Complex comet or unknown asteroid, an event that led to cosmic fragments exploding interhemispherically over widely dispersed areas, including the northern Andes of Venezuela and the Alps on the Italian/French frontier. While the databases in the two areas differ somewhat, the overall interpretation is that microtextural evidence in weathering rinds and in sands of associated paleosols and glaciofluvial deposits carry undeniable attributes of melted glassy carbon and Fe spherules, planar deformation features, shock-melted and contorted quartz, occasional transition and platinum metals, and brecciated and impacted minerals of diverse lithologies. In concert with other black mat localities in the Western USA, the Netherlands, coastal France, Syria, Central Asia, Peru, Argentina and Mexico, it appears that a widespread cosmic impact by an asteroid or comet is responsible for deposition of the black mat at the onset of the YD glacial event. Whether or not the impact caused a 1-kyr interval of glacial climate depends upon whether or not the Earth had multiple centuries-long episodic encounters with the Taurid Complex or asteroid remnants; impact-related changes in microclimates sustained climatic forcing sufficient to maintain positive mass balances in the reformed ice; and/or inertia in the Atlantic thermohaline circulation system persisted for 1kyr. 

On the electron impact ionization of multi-charged ions: a theoretical investigation of relativistic effects

The electron impact ionization of highly charged ions is investigated. Using a relativistic distorted wave Born approximation, we explore the possible effects that should be observable in a high-energy electron impact ionization (e, 2e) coincidence experiment involving multi-charged ions. We present calculations of triple-differential cross sections. We will focus on relativistic and distortion effects and consider geometries where these may be easily observed.

Detection of a Red Supergiant Progenitor Star of a Type II-Plateau Supernova

We present the discovery of a red supergiant star that exploded as supernova 2003gd in the nearby spiral galaxy M74. The Hubble Space Telescope (HST) and the Gemini Telescope imaged this galaxy 6 to 9 months before the supernova explosion, and subsequent HST images confirm the positional coincidence of the supernova with a single resolved star that is a red supergiant of 8+4-2 solar masses. This confirms both stellar evolution models and supernova theories predicting that cool red supergiants are the immediate progenitor stars of type II-plateau supernovae.

Multiwavelength observations of the M15 intermediate-velocity cloud

We present Westerbork Synthesis Radio Telescope HI images, Lovell telescope multibeam H I wide-field mapping, William Herschel Telescope long-slit echelle Ca II observations, Wisconsin Halpha Mapper (WHAM) facility images, and IRAS ISSA 60- and 100-mum co-added images towards the intermediate- velocity cloud (IVC) at + 70 km s(-1), located in the general direction of the M15 globular cluster. When combined with previously published Arecibo data, the H I gas in the IVC is found to be clumpy, with a peak H I column density of similar to1.5 x 10(20) cm(-2), inferred volume density (assuming spherical symmetry) of similar to24 cm(-3)/D (kpc) and a maximum brightness temperature at a resolution of 81 x 14 arcsec(2) of 14 K. The major axis of this part of the IVC lies approximately parallel to the Galactic plane, as does the low- velocity H I gas and IRAS emission. The H I gas in the cloud is warm, with a minimum value of the full width at half-maximum velocity width of 5 km s(-1) corresponding to a kinetic temperature, in the absence of turbulence, of similar to540 K. From the H I data, there are indications of two-component velocity structure. Similarly, the Ca II spectra, of resolution 7 km s(-1), also show tentative evidence of velocity structure, perhaps indicative of cloudlets. Assuming that there are no unresolved narrow-velocity components, the mean values of log(10)[N(Ca II K) cm(2)] similar to 12.0 and Ca II/H I similar to2 5 x 10(-8) are typical of observations of high Galactic latitude clouds. This compares with a value of Ca II/H I>10(-6) for IVC absorption towards HD 203664, a halo star of distance 3 kpc, some 3.degrees1 from the main M15 IVC condensation. The main IVC condensation is detected by WHAM in Halpha with central local-standard-of-rest velocities of similar to60-70 km s(-1), and intensities uncorrected for Galactic extinction of up to 1.3 R, indicating that the gas is partially ionized. The FWHM values of the Halpha IVC component, at a resolution of 1degrees, exceed 30 km s(-1). This is some 10 km s(-1) larger than the corresponding H I value at a similar resolution, and indicates that the two components may not be mixed. However, the spatial and velocity coincidence of the Halpha and H I peaks in emission towards the main IVC component is qualitatively good. If the Halpha emission is caused solely by photoionization, the Lyman continuum flux towards the main IVC condensation is similar to2.7 x 10(6) photon cm(-2) s(-1). There is not a corresponding IVC Halpha detection towards the halo star HD 203664 at velocities exceeding similar to60 km s(- 1). Finally, both the 60- and 100-mum IRAS images show spatial coincidence, over a 0.675 x 0 625 deg(2) field, with both low- and intermediate-velocity H I gas (previously observed with the Arecibo telescope), indicating that the IVC may contain dust. Both the Halpha and tentative IRAS detections discriminate this IVC from high-velocity clouds, although the H I properties do not. When combined with the H I and optical results, these data point to a Galactic origin for at least parts of this IVC.

Spectroscopy and metastability of CO22+ molecular ions

The spectroscopy and metastability of the carbon dioxide doubly charged ion, the CO22+ dication, have been studied with photoionization experiments: time-of-flight photoelectron photoelectron coincidence (TOF-PEPECO), threshold photoelectrons coincidence (TPEsCO), and threshold photoelectrons and ion coincidence (TPEsCO ion coincidence) spectroscopies. Vibrational structure is observed in TOF-PEPECO and TPEsCO spectra of the ground and first two excited states. The vibrational structure is dominated by the symmetric stretch except in the TPEsCO spectrum of the ground state where an antisymmetric stretch progression is observed. All three vibrational frequencies are deduced for the ground state and symmetric stretch and bending frequencies are deduced for the first two excited states. Some vibrational structure of higher electronic states is also observed. The threshold for double ionization of carbon dioxide is reported as 37.340+/-0.010 eV. The fragmentation of energy selected CO22+ ions has been investigated with TPEsCO ion coincidence spectroscopy. A band of metastable states from similar to38.7 to similar to41 eV above the ground state of neutral CO2 has been observed in the experimental time window of similar to0.1-2.3 mus with a tendency towards shorter lifetimes at higher energies. It is proposed that the metastability is due to slow spin forbidden conversion from bound excited singlet states to unbound continuum states of the triplet ground state. Another result of this investigation is the observation of CO++O+ formation in indirect dissociative double photoionization below the threshold for formation of CO22+. The threshold for CO++O+ formation is found to be 35.56+/-0.10 eV or lower, which is more than 2 eV lower than previous measurements. (C) 2005 American Institute of Physics.

Dissociative charge exchange and ionization of O-2 by fast H+ and O+ ions: Energetic ion interactions in Europa's oxygen atmosphere and neutral torus

Measurements of electron capture and ionization of O-2 molecules in collisions with H+ and O+ ions have been made over an energy range 10 - 100 keV. Cross sections for dissociative and nondissociative interactions have been separately determined using coincidence techniques. Nondissociative channels leading to O-2(+) product formation are shown to be dominant for both the H+ and the O+ projectiles in the capture collisions and only for the H+ projectiles in the ionization collisions. Dissociative channels are dominant for ionizing collisions involving O+ projectiles. The energy distributions of the O+ fragment products from collisions involving H+ and O+ have also been measured for the first time using time-of-flight methods, and the results are compared with those from other related studies. These measurements have been used to describe the interaction of the energetic ions trapped in Jupiter's magnetosphere with the very thin oxygen atmosphere of the icy satellite Europa. It is shown that the ionization of oxygen molecules is dominated by charge exchange plus ion impact ionization processes rather than photoionization. In addition, dissociation is predominately induced through excitation of electrons into high-lying repulsive energy states ( electronically) rather than arising from momentum transfer from knock-on collisions between colliding nuclei, which are the only processes included in current models. Future modeling will need to include both these processes.

Galactic Edge Clouds I. Molecular Line Observations and Chemical Modelling of Edge Cloud 2

Edge Cloud 2 (EC2) is a molecular cloud, about 35 pc in size, with one of the largest galactocentric distances known to exist in the Milky Way. We present observations of a peak CO emission region in the cloud and use these to determine its physical characteristics. We calculate a gas temperature of 20 K and a density of n(H2)~10^4 cm-3. Based on our CO maps, we estimate the mass of EC2 at around 10^4 Msolar and continuum observations suggest a dust-to-gas mass ratio as low as 0.001. Chemical models have been developed to reproduce the abundances in EC2, and they indicate that heavy element abundances may be reduced by a factor of 5 relative to the solar neighborhood (similar to dwarf irregular galaxies and damped Lya systems), very low extinction (A_V <4 mag) due to a very low dust-to-gas mass ratio, an enhanced cosmic-ray ionization rate, and a higher UV field compared to local interstellar values. The reduced abundances may be attributed to the low level of star formation in this region and are probably also related to the continuing infall of primordial (or low-metallicity) halo gas since the Milky Way formed. Finally, we note that shocks from the old supernova remnant GSH 138-01-94 may have determined the morphology and dynamics of EC2.

The 90-degree problem of scattering theory revisited: dynamical turbulence effects versus nonlinear effects

The 90° problem of cosmic-ray transport theory is revisited in this paper. By using standard forms of the wave spectrum in the solar wind, the pitch-angle Fokker–Planck coefficient and the parallel mean free path are computed for different resonance functions. A critical comparison is made of the strength of 90° scattering due to plasmawave effects, dynamical turbulence effects and nonlinear effects. It is demonstrated that, only for low-energy cosmic particles, dynamical effects are usually dominant. The novel results presented here are essential for an effective comparison of heliospheric observations for the parallel mean free path with the theoretical model results.