We discuss the physics of stochastic particle acceleration in relativistic magnetohydrodynamic (MHD) turbulence, combining numerical simulations of test-particle acceleration in synthetic wave turbulence spectra with detailed analytical estimates. In particular, we study particle acceleration in wavelike isotropic fast mode turbulence, in Alfv\'en and slow Goldreich-Sridhar type wave turbulence (properly accounting for anisotropy effects), including resonance broadening due to wave decay and pitch-angle randomization. At high particle rigidities, the contributions of those three modes to acceleration are comparable to within an order of magnitude, as a combination of several effects (partial disappearance of transit-time damping for fast modes, increased scattering rate for Alfv\'en and slow modes due to resonance broadening). Additionally, we provide analytical arguments regarding acceleration beyond the regime of MHD wave turbulence, addressing the issue of nonresonant acceleration in a turbulence comprised of structures rather than waves, as well as the issue of acceleration in small-scale parallel electric fields. Finally, we compare our results to the existing literature and provide ready-to-use formulas for applications to high-energy astrophysical phenomenology.

Weakly magnetized, relativistic collisionless shock waves are not only the natural offsprings of relativistic jets in high-energy astrophysical sources, they are also associated with some of the most outstanding displays of energy dissipation through particle acceleration and radiation. Perhaps their most peculiar and exciting feature is that the magnetized turbulence that sustains the acceleration process, and (possibly) the secondary radiation itself, is self-excited by the accelerated particles themselves, so that the phenomenology of these shock waves hinges strongly on the microphysics of the shock. In this review, we draw a status report of this microphysics, benchmarking analytical arguments with particle-in-cell simulations, and extract consequences of direct interest to the phenomenology, regarding in particular the so-called microphysical parameters used in phenomenological studies.

The most powerful persistent accelerators in the Universe are jetted active galaxies. Galaxies whose jets are directed towards Earth, so called blazars, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress is expected in the decade to come, particularly with the deployment of the Cherenkov Telescope Array.

In this first paper of a series dedicated to the microphysics of unmagnetized, relativistic collisionless pair shocks, we discuss the physics of the Weibel-type transverse current filamentation instability (CFI) that develops in the shock precursor, through the interaction of an ultrarelativistic suprathermal particle beam with the background plasma. We introduce in particular the notion of "Weibel frame", or scattering center frame, in which the microturbulence is of mostly magnetic nature. We calculate the properties of this frame, using first a kinetic formulation of the linear phase of the instability, relying on Maxwell-J\"uttner distribution functions, then using a quasistatic model of the nonlinear stage of the instability. Both methods show that: (i) the Weibel frame moves at subrelativistic velocities relative to the background plasma, therefore at relativistic velocities relative to the shock front; (ii) the velocity of the Weibel frame relative to the background plasma scales with $\xi_{\rm b}$, i.e., the pressure of the suprathermal particle beam in units of the momentum flux density incoming into the shock; and (iii), the Weibel frame moves slightly less fast than the background plasma relative to the shock front. Our theoretical results are found to be in satisfactory agreement with the measurements carried out in dedicated large-scale 2D3V PIC simulations.

We develop a comprehensive theoretical model of relativistic collisionless pair shocks mediated by the current filamentation instability. We notably characterize the noninertial frame in which this instability is of a mostly magnetic nature, and describe at a microscopic level the deceleration and heating of the incoming background plasma through its collisionless interaction with the electromagnetic turbulence. Our model compares well to large-scale 2D3V PIC simulations, and provides an important touchstone for the phenomenology of such plasma systems.

In this second paper of a series, we discuss the dynamics of a plasma entering the precursor of an unmagnetized, relativistic collisionless pair shock. We discuss how this background plasma is decelerated and heated through its interaction with a microturbulence that results from the growth of a current filamentation instability (CFI) in the shock precursor. We make use, in particular, of the reference frame $\mathcal R_{\rm w}$ in which the turbulence is mostly magnetic. This frame moves at relativistic velocities towards the shock front at rest, decelerating gradually from the far to the near precursor. In a first part, we construct a fluid model to derive the deceleration law of the background plasma expected from the scattering of suprathermal particles off the microturbulence. This law leads to the relationship $\gamma_{\rm p}\,\sim\,\xi_{\rm b}^{-1/2}$ between the background plasma Lorentz factor $\gamma_{\rm p}$ and the normalized pressure of the beam $\xi_{\rm b}$; it is found to match nicely the spatial profiles observed in large-scale 2D3V particle-in-cell simulations. In a second part, we model the dynamics of the background plasma at the kinetic level, incorporating the inertial effects associated with the deceleration of $\mathcal R_{\rm w}$ into a Vlasov-Fokker-Planck equation for pitch-angle diffusion. We show how the effective gravity in $\mathcal R_{\rm w}$ drives the background plasma particles through friction on the microturbulence, leading to efficient plasma heating. Finally, we compare a Monte Carlo simulation of our model with dedicated PIC simulations and conclude that it can satisfactorily reproduce both the heating and the deceleration of the background plasma in the shock precursor, thereby providing a successful 1D description of the shock transition at the microscopic level.

In this third paper of a series, we discuss the physics of the population of accelerated particles in the precursor of an unmagnetized, relativistic collisionless pair shock. In particular, we provide a theoretical estimate of their scattering length $l_{scatt}(p)$ in the self-generated electromagnetic turbulence, as well as an estimate of their distribution function. We obtain $l_{scatt}(p) \simeq (\gamma_p /\epsilon_B)(p/\gamma_{\infty} mc)^2 (c/\omega_p)$, with p the particle momentum in the rest frame of the shock front, $\epsilon_B$ the strength parameter of the microturbulence, $\gamma_p$ the Lorentz factor of the background plasma relative to the shock front and $\gamma_{\infty}$ its asymptotic value outside the precursor. We compare this scattering length to large-scale PIC simulations and find good agreement for the various dependencies.

In highly conducting astrophysical plasmas, charged particles are generically accelerated through Fermi-type processes involving repeated interactions with moving magnetized scattering centers. The present paper proposes a generalized description of these acceleration processes, by following the momentum of the particle through a continuous sequence of accelerated frames, defined in such a way that the electric field vanishes at each point along the particle trajectory. In each locally inertial frame, the Lorentz force affects the direction of motion of the particle, but the energy changes solely as a result of inertial corrections. This unified description of Fermi acceleration applies equally well in sub- and ultrarelativistic settings, in Cartesian or non-Cartesian geometries, flat or nonflat space-time. Known results are recovered in a variety of regimes -- shock, turbulent and shear acceleration -- and new results are derived in lieu of applications, e.g. nonresonant acceleration in relativistic turbulence, stochastic unipolar inductive acceleration and centrifugo-shear acceleration close to the horizon of a black hole.

The nonlinear evolution of current filaments generated by the Weibel-type filamentation instability is a topic of prime interest in space and laboratory plasma physics. In this paper, we investigate the stability of a stationary periodic chain of nonlinear current filaments in counterstreaming pair plasmas. We make use of a relativistic four-fluid model and apply the Floquet theory to compute the two-dimensional unstable eigenmodes of the spatially periodic system. We examine three different cases, characterized by various levels of nonlinearity and asymmetry between the plasma streams: a weakly nonlinear symmetric system, prone to purely transverse merging modes; a strongly nonlinear symmetric system, dominated by coherent drift-kink modes whose transverse periodicity is equal to, or an integer fraction of the unperturbed filaments; a moderately nonlinear asymmetric system, subject to a mix of kink and bunching-type perturbations. The growth rates and profiles of the numerically computed eigenmodes agree with particle-in-cell simulation results. In addition, we derive an analytic criterion for the transition between dominant filament-merging and drift-kink instabilites in symmetric two-beam systems.

The generation of turbulence at magnetized shocks and its subsequent interaction with the latter is a key question of plasma- and high-energy astrophysics. This paper presents two-dimensional magnetohydrodynamic simulations of a fast shock front interacting with incoming upstream perturbations, described as harmonic entropy or fast magnetosonic waves, both in the relativistic and the sub-relativistic regimes. We discuss how the disturbances are transmitted into downstream turbulence and we compare the observed response for small amplitude waves to a recent linear calculation. In particular, we demonstrate the existence of a resonant response of the corrugation amplitude when the group velocity of the outgoing downstream fast mode matches the velocity of the shock front. We also present simulations of large amplitude waves to probe the non-linear regime.

Relativistic astrophysical collisionless shocks represent outstanding dissipation agents of the huge power of relativistic outflows produced by accreting black holes, core collapsed supernovae and other objects into multi-messenger radiation (cosmic rays, neutrinos, electromagnetic radiation). This article provides a theoretical discussion of the fundamental physical ingredients of these extreme phenomena. In the context of weakly magnetized shocks, in particular, it is shown how the filamentation type instabilities, which develop in the precursor of pair dominated or electron-ion shocks, provide the seeds for the scattering of high energy particles as well as the agent which preheats and slows down the incoming precursor plasma. This analytical discussion is completed with a mesoscopic, non-linear model of particle acceleration in relativistic shocks based on Monte Carlo techniques. This Monte Carlo model uses a semi-phenomenological description of particle scattering which allows it to calculate the back-reaction of accelerated particles on the shock structure on length and momentum scales which are currently beyond the range of microscopic particle-in-cell (PIC) simulations.

Most of the plasma microphysics which shapes the acceleration process of particles at collisionless shock waves takes place in the cosmic-ray precursor, through the interaction of accelerated particles with the unshocked plasma. Detecting directly or indirectly the synchrotron radiation of accelerated particles in this precursor would open a new window on the microphysics of acceleration and of collisionless shock waves. We provide analytical estimates of the spectrum and of the polarization fraction of the synchrotron precursor for both relativistic and non-relativistic collisionless shock fronts, accounting for the self-generation or amplification of magnetic turbulence. In relativistic sources, the spectrum of the precursor is harder than that of the shocked plasma, because the upstream residence time increases with particle energy, leading to an effectively hard spectrum of accelerated particles in the precursor volume. At high frequencies, typically in the optical to X-ray range, the contribution of the precursor becomes sizeable, but we find that in most cases studied, it remains dominated by the synchrotron or inverse Compton contribution of the shocked plasma; its contribution might be detectable only in trans-relativistic shock waves. Non-relativistic sources offer the possibility of spectral imaging of the precursor by viewing the shock front edge-on. We calculate this spectro-morphological contribution for various parameters. The synchrotron contribution is also sizeable at the highest frequencies (X-ray range). If the turbulence is tangled in the plane transverse to the shock front, the resulting synchrotron radiation should be nearly maximally linearly polarized; polarimetry thus arises as an interesting tool to reveal this precursor.[Abridged]

Successful phenomenological models of pulsar wind nebulae assume efficient dissipation of the Poynting flux of the magnetized electron-positron wind as well as efficient acceleration of the pairs in the vicinity of the termination shock, but how this is realized is not yet well understood. The present paper suggests that the corrugation of the termination shock, at the onset of non-linearity, may lead towards the desired phenomenology. Non-linear corrugation of the termination shock would convert a fraction of order unity of the incoming ordered magnetic field into downstream turbulence, slowing down the flow to sub-relativistic velocities. The dissipation of turbulence would further preheat the pair population on short length scales, close to equipartition with the magnetic field, thereby reducing the initial high magnetization to values of order unity. Furthermore, it is speculated that the turbulence generated by the corrugation pattern may sustain a relativistic Fermi process, accelerating particles close to the radiation reaction limit, as observed in the Crab nebula. The required corrugation could be induced by the fast magnetosonic modes of downstream nebular turbulence; but it could also be produced by upstream turbulence, either carried by the wind or seeded in the precursor by the accelerated particles themselves.

As a shock front interacts with turbulence, it develops corrugation which induces outgoing wave modes in the downstream plasma. For a fast shock wave, the incoming wave modes can either be fast magnetosonic waves originating from downstream, outrunning the shock, or eigenmodes of the upstream plasma drifting through the shock. Using linear perturbation theory in relativistic MHD, this paper provides a general analysis of the corrugation of relativistic magnetized fast shock waves resulting from their interaction with small amplitude disturbances. Transfer functions characterizing the linear response for each of the outgoing modes are calculated as a function of the magnetization of the upstream medium and as a function of the nature of the incoming wave. Interestingly, if the latter is an eigenmode of the upstream plasma, we find that there exists a resonance at which the (linear) response of the shock becomes large or even diverges. This result may have profound consequences on the phenomenology of astrophysical relativistic magnetized shock waves.

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebul\ae, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in-situ observations, analytical and numerical developments. A particular emphasize is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics

Recent analyses of multiwavelength light curves of gamma-ray bursts afterglows point to values of the magnetic turbulence well below the canonical $\sim1\,$\% of equipartition, in agreement with theoretical expectations of a micro-turbulence generated in the shock precursor, which then decays downstream of the shock front through collisionless damping. As a direct consequence, the Compton parameter $Y$ can take large values in the blast. In the presence of decaying micro-turbulence and/or as a result of the Klein-Nishina suppression of inverse Compton cooling, the $Y$ parameter carries a non-trivial dependence on the electron Lorentz factor, which modifies the spectral shape of the synchrotron and inverse Compton components. This paper provides detailed calculations of this synchrotron-self-Compton spectrum in this large $Y$ regime, accounting for the possibility of decaying micro-turbulence. It calculates the expected temporal and spectral indices $\alpha$ and $\beta$ customarily defined by $F_\nu\,\propto\,t_{\rm obs}^{-\alpha}\nu^{-\beta}$ in various spectral domains. This paper also makes predictions for the very high energy photon flux; in particular, it shows that the large $Y$ regime would imply a detection rate of gamma-ray bursts at $>10\,$GeV several times larger than currently anticipated.

Pulsar wind nebulae (PWNe) are outstanding accelerators in Nature, in the sense that they accelerate electrons up to the radiation reaction limit. Motivated by this observation, this paper examines the possibility that young pulsar wind nebulae can accelerate ions to ultra-high energies at the termination shock of the pulsar wind. We consider here powerful PWNe, fed by pulsars born with $\sim\,$millisecond periods. Assuming that such pulsars exist, at least during a few years after the birth of the neutron star, and that they inject ions into the wind, we find that protons could be accelerated up to energies of the order of the Greisen-Zatsepin-Kuzmin cut-off, for a fiducial rotation period $P\,\sim\,1\,$msec and a pulsar magnetic field $B_{\star}\,\sim\,10^{13}\,$G, implying a fiducial wind luminosity $L_{\rm p}\,\sim\,10^{45}\,$erg/s and a spin-down time $t_{\rm sd}\,\sim 3\times 10^7\,$s. The main limiting factor is set by synchrotron losses in the nebula and by the size of the termination shock; ions with $Z\geq 1$ may therefore be accelerated to even higher energies. We derive an associated neutrino flux produced by interactions in the source region. For a proton-dominated composition, our maximum flux lies slightly below the 5-year sensitivity of IceCube-86 and above the 3-year sensitivity of the projected Askaryan Radio Array. It might thus become detectable in the next decade, depending on the exact level of contribution of these millisecond pulsar wind nebulae to the ultra-high energy cosmic ray flux.

We review the physics of relativistic shocks, which are often invoked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possible sources of ultra-high energy cosmic-rays. We focus on particle acceleration and magnetic field generation, and describe the recent progress in the field driven by theory advances and by the rapid development of particle-in-cell (PIC) simulations. In weakly magnetized or quasi parallel-shocks (where the magnetic field is nearly aligned with the flow), particle acceleration is efficient. The accelerated particles stream ahead of the shock, where they generate strong magnetic waves which in turn scatter the particles back and forth across the shock, mediating their acceleration. In contrast, in strongly magnetized quasi-perpendicular shocks, the efficiencies of both particle acceleration and magnetic field generation are suppressed. Particle acceleration, when efficient, modifies the turbulence around the shock on a long time scale, and the accelerated particles have a characteristic energy spectral index of ~ 2.2 in the ultra-relativistic limit. We discuss how this novel understanding of particle acceleration and magnetic field generation in relativistic shocks can be applied to high-energy astrophysical phenomena, with an emphasis on PWNe and GRB afterglows.

The Weibel/filamentation instability is known to play a key role in the physics of weakly magnetized collisionless shock waves. From the point of view of high energy astrophysics, this instability also plays a crucial role because its development in the shock precursor populates the downstream with a small-scale magneto-static turbulence which shapes the acceleration and radiative processes of suprathermal particles. The present work discusses the physics of the dissipation of this Weibel-generated turbulence downstream of relativistic collisionless shock waves. It calculates explicitly the first-order non-linear terms associated to the diffusive nature of the particle trajectories. These corrections are found to systematically increase the damping rate, assuming that the scattering length remains larger than the coherence length of the magnetic fluctuations. The relevance of such corrections is discussed in a broader astrophysical perspective, in particular regarding the physics of the external relativistic shock wave of a gamma-ray burst.

We report here on a fast current-driven instability at relativistic collisionless shocks, triggered by the perpendicular current carried by the supra-thermal particles as they gyrate around the background magnetic field in the shock precursor. We show that this instability grows faster than any other instability studied so far in this context, and we argue that it is likely to shape the physics of the shock and of particle acceleration in a broad parameter range.

This paper summarizes recent progresses in our theoretical understanding of particle acceleration at relativistic shock waves and it discusses two salient consequences: (1) the maximal energy of accelerated particles; (2) the impact of the shock-generated micro-turbulence on the multi-wavelength light curves of gamma-ray burst afterglows.

The physics of collisionless relativistic shocks with a moderate magnetization is presented. Micro-physics is relevant to explain the most energetic radiative phenomena of Nature, namely that of the termination shock of Gamma Ray Bursts. A transition towards Fermi process occurs for decreasing magnetization around a critical value which turns out to be the condition for the scattering to break the mean field inhibition. Scattering is produced by magnetic micro-turbulence driven by the current carried by returning particles, which had not been considered till now, but turns out to be more intense than Weibel's one around the transition. The current is also responsible for a buffer effect on the motion of the incoming flow, on which the threshold for the onset of turbulence depends.

The physics of instabilities in the precursor of relativistic collisionless shocks is of broad importance in high energy astrophysics, because these instabilities build up the shock, control the particle acceleration process and generate the magnetic fields in which the accelerated particles radiate. Two crucial parameters control the micro-physics of these shocks: the magnetization of the ambient medium and the Lorentz factor of the shock front; as of today, much of this parameter space remains to be explored. In the present paper, we report on a new instability upstream of electron-positron relativistic shocks and we argue that this instability shapes the micro-physics at moderate magnetization levels and/or large Lorentz factors. This instability is seeded by the electric current carried by the accelerated particles in the shock precursor as they gyrate around the background magnetic field. The compensation current induced in the background plasma leads to an unstable configuration, with the appearance of charge neutral filaments carrying a current of the same polarity, oriented along the perpendicular current. This ``current-driven filamentation'' instability grows faster than any other instability studied so far upstream of relativistic shocks, with a growth rate comparable to the plasma frequency. Furthermore, the compensation of the current is associated with a slow-down of the ambient plasma as it penetrates the shock precursor (as viewed in the shock rest frame). This slow-down of the plasma implies that the ``current driven filamentation'' instability can grow for any value of the shock Lorentz factor, provided the magnetization \sigma <~ 10^{-2}. We argue that this instability explains the results of recent particle-in-cell simulations in the mildly magnetized regime.

The joint analysis of anisotropy signals and chemical composition of ultra-high energy cosmic rays offers strong potential for shedding light on the sources of these particles. Following up on an earlier idea, this paper studies the anisotropies produced by protons of energy >E/Z, assuming that anisotropies at energy >E have been produced by nuclei of charge Z, which share the same magnetic rigidity. We calculate the number of secondary protons produced through photodisintegration of the primary heavy nuclei. Making the extreme assumption that the source does not inject any proton, we find that the source(s) responsible for anisotropies such as reported by the Pierre Auger Observatory should lie closer than ~20-30, 80-100 and 180-200 Mpc if the anisotropy signal is mainly composed of oxygen, silicon and iron nuclei respectively. A violation of this constraint would otherwise result in the secondary protons forming a more significant anisotropy signal at lower energies. Even if the source were located closer than this distance, it would require an extraordinary metallicity >120, 1600, 1100 times solar metallicity in the acceleration zone of the source, for oxygen, silicon and iron respectively, to ensure that the concomitantly injected protons not to produce a more significant low energy anisotropy. This offers interesting prospects for constraining the nature and the source of ultra-high energy cosmic rays with the increase in statistics expected from next generation detectors.

The origin of magnetic fields that permeate the blast waves of gamma-ray bursts (GRBs) is a long-standing problem. The present paper argues that in four GRBs revealing extended emission at >100 MeV, with follow-up in the radio, optical and X-ray domains at later times, this magnetization can be described as the partial decay of the micro-turbulence that is generated in the shock precursor. Assuming that the bulk of the extended emission >100 MeV can be interpreted as synchrotron emission of shock accelerated electrons, we model the multi-wavelength light curves of GRB 090902B, GRB 090323, GRB 090328 and GRB 110731A, using a simplified then a full synchrotron calculation with power-law-decaying microturbulence \epsilon_B \propto t^{\alpha_t} (t denotes the time since injection through the shock, in the comoving blast frame). We find that these models point to a consistent value of the decay exponent -0.5 < \alpha_t < -0.4.

Fermi/LAT has detected long-lasting high-energy photons (>100 MeV) from gamma-ray bursts (GRBs), with the highest energy photons reaching about 100 GeV. One proposed scenario is that they are produced by high-energy electrons accelerated in GRB forward shocks via synchrotron radiation. We study the maximum synchrotron photon energy in this scenario, considering the properties of the microturbluence magnetic fields behind the shock, as revealed by recent Particle-in-Cell simulations and theoretical analyses of relativistic collisionless shocks. Due to the small-scale nature of the micro-turbulent magnetic field, the Bohm acceleration approximation breaks down at such high energies. This effect leads to a typical maximum synchrotron photon of a few GeV at 100 s after the burst and this maximum synchrotron photon energy decreases quickly with time. We show that the fast decrease of the maximum synchrotron photon energy leads to a fast decay of the synchrotron flux. The 10-100 GeV photons detected after the prompt phase can not be produced by the synchrotron mechanism. They could originate from the synchrotron self-Compton emission of the early afterglow if the circum-burst density is sufficiently large, or from the external inverse-Compton process in the presence of central X-ray emission, such as X-ray flares and prompt high-latitude X-ray emission.

Collisionless relativistic shocks have been the focus of intense theoretical and numerical investigations in recent years. The acceleration of particles, the generation of electromagnetic microturbulence and the building up of a shock front are three interrelated essential ingredients of a relativistic collisionless shock wave. In this paper we investigate two issues of importance in this context: (1) the transport of suprathermal particles in the excited microturbulence upstream of the shock and its consequences regarding particle acceleration; (2) the preheating of incoming background electrons as they cross the shock precursor and experience relativistic oscillations in the microturbulent electric fields. We place emphasis on the importance of the motion of the electromagnetic disturbances relatively to the background plasma and to the shock front. This investigation is carried out for the two major instabilities involved in the precursor of relativistic shocks, the filamentation instability and the oblique two stream instability. Finally, we use our results to discuss the maximal acceleration at the external shock of a gamma-ray burst; we find in particular a maximal synchrotron photon energy of the order of a few GeV.

The origin of the highest energy cosmic rays represents one of the most conspicuous enigmas of modern astrophysics, in spite of gigantic experimental efforts in the past fifty years, and of active theoretical research. The past decade has known exciting experimental results, most particularly the detection of a cut-off at the expected position for the long sought Greisen-Zatsepin-Kuzmin suppression as well as evidence for large scale anisotropies. This paper summarizes and discusses recent achievements in this field.

Microphysics of weakly magnetized relativistic collisionless shock waves, corroborated by recent high performance numerical simulations, indicate the presence of a microturbulent layer of large magnetic field strength behind the shock front, which must decay beyond some hundreds of skin depths. The present paper discusses the dynamics of such microturbulence, borrowing from these same numerical simulations, and calculates the synchrotron signature of a powerlaw of shock accelerated particles. The decaying microturbulent layer is found to leave distinct signatures in the spectro-temporal evolution of the spectrum $F_\nu \propto t^{-\alpha}\nu^{-\beta}$ of a decelerating blast wave, which are potentially visible in early multi-wavelength follow-up observations of gamma-ray bursts. This paper also discusses the influence of the evolving microturbulence on the acceleration process, with particular emphasis on the maximal energy of synchrotron afterglow photons, which falls in the GeV range for standard gamma-ray burst parameters. Finally, this paper argues that the evolving microturbulence plays a key role in shaping the spectra of recently observed gamma-ray bursts with extended GeV emission, such as GRB090510.

In this review we confront the current theoretical understanding of particle acceleration at relativistic outflows with recent observational results on various source classes thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy cosmic rays.

Relativistic sources, e.g. gamma-ray bursts, pulsar wind nebulae and powerful active galactic nuclei produce relativistic outflows that lead to the formation of collisionless shock waves, where particle acceleration is thought to take place. Our understanding of relativistic shock acceleration has improved in the past decade, thanks to the combination of analytical studies and high level numerical simulations. In ultra-relativistic shocks, particle acceleration is made difficult by the generically transverse magnetic field and large advection speed of the shocked plasma. Fast growing microturbulence is thus needed to make the Fermi process operative. It is thought, and numerical simulations support that view, that the penetration of supra-thermal particles in the shock precursor generates a magnetic turbulence which in turn produces the scattering process needed for particle acceleration through the Fermi mechanism. Through the comparison of the growth timescale of the microinstabilities in the shock precursor and the precursor crossing timescale, it is possible to delimit in terms of magnetization and shock Lorentz factor the region in which micro-turbulence may be excited, hence whether and how Fermi acceleration is triggered. These findings are summarized here and astrophysical consequences are drawn.

Recent analytical and numerical work argue that successful relativistic Fermi acceleration requires a weak magnetization of the unshocked plasma, all the more so at high Lorentz factors. The present paper tests this conclusion by computing the afterglow of a gamma-ray burst outflow propagating in a magnetized stellar wind using "ab initio" principles regarding the microphysics of relativistic Fermi acceleration. It is shown that in magnetized environments, one expects a drop-out in the X-ray band on sub-day scales as the synchrotron emission of the shock heated electrons exits the frequency band. At later times, Fermi acceleration becomes operative when the blast Lorentz factor drops below a certain critical value, leading to the recovery of the standard afterglow light curve. Interestingly, the observed drop-out bears resemblance with the fast decay found in gamma-ray bursts early X-ray afterglows.

Fermi acceleration can develop efficiently at relativistic collisionless shock waves provided the upstream (unshocked) plasma is weakly magnetized. At low magnetization, the large size of the shock precursor indeed provides enough time for electromagnetic micro-instabilities to grow and such micro-instabilities generate small scale turbulence that in turn provides the scattering required. The present paper extends our previous analysis on the development of these micro-instabilities to account for the finite angular dispersion of the beam of reflected and accelerated particles and to account for the expected heating of the upstream electrons in the shock precursor. We show that the oblique two stream instability may operate down to values of the shock Lorentz factor gamma_{sh}~10 as long as the electrons of the upstream plasma remain cold, while the filamentation instability is strongly inhibited in this limit; however, as electrons get heated to relativistic temperatures, the situation becomes opposite and the two stream instability becomes inhibited while the filamentation mode becomes efficient, even at moderate values of the shock Lorentz factor. The peak wavelength of these instabilities migrates from the inertial electron scale towards the proton inertial scale as the background electrons get progressively heated during the crossing of the shock precursor. We also discuss the role of current driven instabilities upstream of the shock. In particular, we show that the returning/accelerated particles give rise to a transverse current through their rotation in the background magnetic field. We find that the compensating current in the background plasma can lead to a Buneman instability which provides an efficient source of electron heating. [Abridged]

Various astrophysical studies have motivated the investigation of the transport of high energy particles in magnetic turbulence, either in the source or en route to the observation sites. For strong turbulence and large rigidity, the pitch-angle scattering rate is governed by a simple law involving a mean free path that increases proportionally to the square of the particle energy. In this paper, we show that perpendicular diffusion deviates from this behavior in the presence of a mean field. We propose an exact theoretical derivation of the diffusion coefficients and show that a mean field significantly changes the transverse diffusion even in the presence of a stronger turbulent field. In particular, the transverse diffusion coefficient is shown to reach a finite value at large rigidity instead of increasing proportionally to the square of the particle energy. Our theoretical derivation is corroborated by a dedicated Monte Carlo simulation. We briefly discuss several possible applications in astrophysics.

Detecting and characterizing the anisotropy pattern of the arrival directions of the highest energy cosmic rays are crucial steps towards the identification of their sources. We discuss a possible distortion of the cosmic ray flux induced by the anisotropic and inhomogeneous distribution of extragalactic magnetic fields in cases where sources of ultrahigh energy cosmic rays are rare transient phenomena, such as gamma-ray bursts and/or newly born magnetars. This distortion does not involve an angular deflection but the modulation of the flux related to the probability of seeing the source on an experiment lifetime. To quantify this distortion, we construct sky maps of the arrival directions of these highest energy cosmic rays for various magnetic field configurations and appeal to statistical tests proposed in the literature. We conclude that this distortion cannot affect present experiments but should be considered when performing anisotropy studies with future large-scale experiments that record as many as hundreds of events above 6x10^19 eV.

The injection of ultra-high energy cosmic rays in the intergalactic medium leads to the production of a GeV-TeV gamma-ray halo centered on the source location, through the production of a high electromagnetic component in the interactions of the primary particles with the radiation backgrounds. This paper examines the prospects for the detectability of such gamma ray halos. We explore a broad range of astrophysical parameters, including the inhomogeneous distribution of magnetic fields in the large scale structure as well as various possible chemical compositions and injection spectra; and we consider the case of a source located outside clusters of galaxies. With respect to the gamma-ray flux associated to synchrotron radiation of ultra-high energy secondary pairs, we demonstrate that it does not depend strongly on these parameters and conclude that its magnitude ultimately depends on the energy injected in the primary cosmic rays. Bounding the cosmic ray luminosity with the contribution to the measured cosmic ray spectrum, we then find that the gamma-ray halo produced by equal luminosity sources is well below current or planned instrument sensitivities. Only rare and powerful steady sources, located at distances larger than several hundreds of Mpc and contributing to a fraction > 10% of the flux at 10^19 eV might be detectable. We also discuss the gamma-ray halos that are produced by inverse Compton/pair production cascades seeded by ultra-high energy cosmic rays. This latter signal strongly depends on the configuration of the extragalactic magnetic fields; it is dominated by the synchrotron signal on a degree scale if the filling factor of magnetic fields with B > 10^(-14) G is smaller than a few percents. Finally, we discuss briefly the case of nearby potential sources such as Centaurus A.

In the pre-reheating era, following cosmic inflation and preceding radiation domination, the energy density may be dominated by an oscillating massive scalar condensate, such as is the case for quadratic chaotic inflation. We have found in a previous paper that during this period, a wide range of sub-Hubble scale perturbations are subject to a preheating instability, leading to the growth of density perturbations ultimately collapsing to form non-linear structures. We compute here the gravitational wave signal due to these structures in the linear limit and present estimates for emission in the non-linear limit due to various effects: the collapse of halos, the tidal interactions, the evaporation during the conversion of the inflaton condensate into radiation and finally the ensuing turbulent cascades. The gravitational wave signal could be rather large and potentially testable by future detectors.

After cosmic inflation and before the transition to radiation domination, the cosmic energy density may have been dominated during an extended period by an oscillating massive scalar condensate. We show that during this period, sub-Hubble scale perturbations are subject to a metric preheating instability in the narrow resonance regime. We analyze in detail both, quadratic and quartic potentials. The instability leads to the growth of density perturbations which in many cases become non-linear already before the beginning of a radiation dominated Universe. This is particularly the case when requiring a phenomenologically preferred low reheat temperature. These early structures may lead to the emission of gravitational waves and the production of primordial black holes. Furthermore, it is not clear if they could modify the prediction of linear curvature perturbations on very large scales.

(Abridged) This paper addresses the issue of magnetic field generation in a relativistic shock precursor through micro-instabilities. The level of magnetization of the upstream plasma turns out to be a crucial parameter, notably because the length scale of the shock precursor is limited by the Larmor rotation of the accelerated particles in the background magnetic field and by the speed of the shock wave. We discuss in detail and calculate the growth rates of the following beam plasma instabilities seeded by the accelerated and reflected particle populations: for an unmagnetized shock, the Weibel and filamentation instabilities, as well as the Cerenkov resonant longitudinal and oblique modes; for a magnetized shock, the Weibel instability and the resonant Cerenkov instabilities with the longitudinal electrostatic modes, as well as the Alfven, Whisler and extraordinary modes. All these instabilities are generated upstream, then they are transmitted downstream. The modes excited by Cerenkov resonant instabilities take on particular importance with respect to the magnetisation of the downstream medium since, being plasma eigenmodes, they have a longer lifetime than the Weibel modes. We discuss the main limitation of the wave growth associated with the length of the precursor and the magnetisation of the upstream medium for both oblique and parallel shock waves. We also characterize the proper conditions to obtain Fermi acceleration. We recover some results of most recent particle-in-cell simulations and conclude with some applications to astrophysical cases of interest, pulsar winds and gamma-ray burst external shock waves in particular. (Abridged)

This paper proposes and discusses a test of the chemical composition of ultra-high energy cosmic rays that relies on the anisotropy patterns measured as a function of energy. In particular, we show that if one records an anisotropy signal produced by heavy nuclei of charge Z above an energy E_{thr}, one should record an even stronger (possibly much stronger) anisotropy at energies >E_{thr}/Z due to the proton component that is expected to be associated with the sources of the heavy nuclei. This conclusion remains robust with respect to the parameters characterizing the sources and it does not depend at all on the modelling of astrophysical magnetic fields. As a concrete example, we apply this test to the most recent data of the Pierre Auger Observatory. Assuming that the anisotropy reported above 55EeV is not a statistical accident, and that no significant anisotropy has been observed at energies <10EeV, we show that the apparent clustering toward Cen A cannot be attributed to heavy nuclei. Similar conclusions are drawn regarding the apparent excess correlation with nearby active galactic nuclei. We then discuss a robust lower bound to the magnetic luminosity that a source must possess in order to be able to accelerate particles of charge Z up to 100EeV, L_B>10^{45}Z^{-2}erg/s. Using this bound in conjunction with the above conclusions, we argue that the current PAO data does not support the model of cosmic ray origin in active radio-quiet or even radio-loud galaxies. Finally, we demonstrate that the apparent clustering in the direction of Cen A can be explained by the contribution of the last few gamma-ray bursts or magnetars in the host galaxy thanks to the scattering of the cosmic rays on the magnetized lobes.

Moduli fields generically produce strong dark matter -- radiation and baryon -- radiation isocurvature perturbations through their decay if they remain light during inflation. We show that existing upper bounds on the magnitude of such fluctuations can thus be translated into stringent constraints on the moduli parameter space m_\sigma (modulus mass) -- \sigma_{inf} (modulus vacuum expectation value at the end of inflation). These constraints are complementary to previously existing bounds so that the moduli problem becomes worse at the perturbative level. In particular, if the inflationary scale H_{inf}~10^{13} GeV, particle physics scenarios which predict high moduli masses m_\sigma > 10-100 TeV are plagued by the perturbative moduli problem, even though they evade big-bang nucleosynthesis constraints.

We set constraints on moduli cosmology from the production of dark matter -- radiation and baryon -- radiation isocurvature fluctuations through modulus decay, assuming the modulus remains light during inflation. We find that the moduli problem becomes worse at the perturbative level as a significant part of the parameter space m_\sigma (modulus mass) -- \sigma_{inf} (modulus vev at the end of inflation) is constrained by the non-observation of significant isocurvature fluctuations. We discuss in detail the evolution of the modulus vev and perturbations, in particular the consequences of Hubble scale corrections to the modulus potential and the stochastic motion of the modulus during inflation. We show, in particular, that a high modulus mass scale m_\sigma > 100 TeV, which allows the modulus to evade big-bang nucleosynthesis constraints is strongly constrained at the perturbative level. We find that generically, solving the moduli problem requires the inflationary scale to be much smaller than 10^{13} GeV.

Fermi acceleration can take place at ultra-relativistic shock waves if the upstream or downstream magnetic field has been remodeled so that most of the magnetic power lies on short spatial scales. The relevant conditions under which Fermi acceleration become efficient in the presence of both a coherent and a short scale turbulent magnetic field are addressed. Within the MHD approximation, this paper then studies the amplification of a pre-existing magnetic field through the streaming of cosmic rays upstream of a relativistic shock wave. The magnetic field is assumed to be perpendicular in the shock front frame, as generally expected in the limit of large shock Lorentz factor. In the MHD regime, compressive instabilities seeded by the net cosmic-ray charge in the shock precursor (as seen in the shock front frame) develop on the shortest spatial scales but saturate at a moderate level deltaB/B ~ 1, which is not sufficient for Fermi acceleration. As we argue, it is possible that other instabilities outside the MHD range provide enough amplification to allow successful Fermi acceleration.

After a successful development of theoretical and numerical works on Fermi acceleration at relativistic shocks, some difficulties recently raised with the scattering issue, a crucial aspect of the process. Most pioneering works were developed assuming the scattering off magnetic fluctuations as given. Even in that case, when a mean field is considered, its orientation is mostly perpendicular to the shock normal in the front frame, and this tends to quench the scattering process. Solving this difficulty leads to address the issue of the generation of very intense magnetic fluctuations at short wave lengths. The relativistic motion of the shock front let the cosmic rays to visit upstream during a very short time only, making this generation of magnetic fluctuations very challenging. Anyway there is some hope to solve the problem. Thanks to a recent work by Spitkovsky (2008) \cite{AS}, we know that the process works without any mean field and now we have to investigate up to which intensity the mean field can be amplified for allowing Fermi process with appropriate fast instabilities. In this presentation, the collisionless shock structure in relativistic regime is sketched, the scattering issue is presented, and the instabilities that can provide the expected magnetic field amplification are presented as well. Although there exists observational evidence that particles are accelerated in relativistic flows and are distributed according to a power law suggesting a Fermi process, the drastic conditions for Fermi process to work are not always clearly fulfilled.

This paper calculates the amount of baryon/radiation isocurvature fluctuation produced through the decay of a curvaton field. It is shown in particular that if curvaton decay preserves baryon number and the curvaton dominates the energy density at the time of decay, the initial curvaton/radiation isocurvature mode is entirely transfered into a baryon/radiation isocurvature mode. This situation is opposite to that previously studied in three fluid models of curvaton decay; this difference is related to the conservation of the pre-existing baryon asymmetry and to the efficiency of the annihilation of all baryon/anti-baryon pairs produced in the decay. We study in detail the relevant cases in which the curvaton decay preserves or not baryon number and provide analytical and numerical calculations for each situation.

This paper provides an analytical description of the transport of ultrahigh energy cosmic rays in an inhomogeneously magnetized intergalactic medium. This latter is modeled as a collection of magnetized scattering centers such as radio cocoons, magnetized galactic winds, clusters or magnetized filaments of large scale structure, with negligible magnetic fields in between. Magnetic deflection is no longer a continuous process, it is rather dominated by scattering events. We study the interaction between high energy cosmic rays and the scattering agents. We then compute the optical depth of the Universe to cosmic ray scattering and discuss the phenomological consequences for various source scenarios. For typical parameters of the scattering centers, the optical depth is greater than unity at 5x10^{19}eV, but the total angular deflection is smaller than unity. One important consequence of this scenario is the possibility that the last scattering center encountered by a cosmic ray be mistaken with the source of this cosmic ray. In particular, we suggest that part of the correlation recently reported by the Pierre Auger Observatory may be affected by such delusion: this experiment may be observing in part the last scattering surface of ultrahigh energy cosmic rays rather than their source population. Since the optical depth falls rapidly with increasing energy, one should probe the arrival directions of the highest energy events beyond 10^{20}eV on an event by event basis to circumvent this effect.

Various experiments indicate the existence of a second knee around energy E=3.10^{17} eV in the cosmic ray spectrum. This feature could be the signature of the end of the galactic component and of the emergence of the extragalactic one, provided that the latter cuts off at low energies. Recent analytical calculations have shown that this cut-off could be a consequence of the existence of extragalactic magnetic fields: low energy protons diffuse on extragalactic magnetic fields and cannot reach the observer within a given time. We study the influence of inhomogeneous magnetic fields on the magnetic horizon, using a new semi-analytical propagation code. Our results indicate that, at a fixed value of the volume averaged magnetic field <B>, the amplitude of the low energy cut-off is mainly controled by the strength of magnetic fields in the voids of the large scale structure distribution.

We build a realistic model of curvaton cosmology, in which the energy content is described by radiation, WIMP dark matter and a curvaton component. We calculate the curvature and isocurvature perturbations, allowing for arbitrary initial density perturbations in all fluids, following all species and their perturbations from the onset of dark matter freeze-out onto well after curvaton decay. We provide detailed numerical evaluations as well as analytical formulae which agree well with the latter. We find that substantial isocurvature perturbations, as measured relatively to the total curvature perturbation, can be produced even if the curvaton energy density is well underdominant when it decays; high precision measurements of cosmic microwave background anisotropies may thus open a window on underdominant decoupled species in the pre-nucleosynthesis early Universe. We also find that in a large part of parameter space, curvaton decay produces enough dark matter particles to restore WIMP annihilations, leading to the partial erasure of any pre-existing dark matter - radiation isocurvature perturbation.

It is shown that gravitinos with mass m_{3/2} ~ 0.1-1 MeV may provide suitable cold dark matter candidates in scenarios of gauge mediated supersymmetry breaking (GMSB) under SO(10) grand unification coupled to supergravity, which accommodate a messenger sector of mass scale M_X ~ 10^6 GeV. This is due to the combined effects of renormalizable loop-suppressed operators and generic non-renormalizable ones governing the dilution of a pre-existing equilibrium gravitino abundance via messenger decay. The above range of gravitino and messenger masses can be accommodated in indirect GMSB scenarios. The gravitino abundance does not depend on the post-inflationary reheat temperature and it is shown that leptogenesis can generate successfully the baryon asymmetry.

It is shown that Fermi acceleration at an ultra-relativistic shock wave cannot operate on a particle for more than 1 1/2 Fermi cycle (i.e., u -> d -> u -> d) if the particle Larmor radius is much smaller than the coherence length of the magnetic field on both sides of the shock, as is usually assumed. This conclusion is shown to be in excellent agreement with recent numerical simulations. We thus argue that efficient Fermi acceleration at ultra-relativistic shock waves requires significant non-linear processing of the far upstream magnetic field with strong amplification of the small scale magnetic power. The streaming or transverse Weibel instabilities are likely to play a key role in this respect.

Stable particle dark matter may well originate during the decay of long-lived relic particles, as recently extensively examined in the cases of the axino, gravitino, and higher-dimensional Kaluza-Klein (KK) graviton. It is shown that in much of the viable parameter space such dark matter emerges naturally warm/hot or mixed. In particular, decay produced gravitinos (KK-gravitons) may only be considered cold for the mass of the decaying particle in the several TeV range, unless the decaying particle and the dark matter particle are almost degenerate. Such dark matter candidates are thus subject to a host of cosmological constraints on warm and mixed dark matter, such as limits from a proper reionisation of the Universe, the Lyman-alpha forest, and the abundance of clusters of galaxies.. It is shown that constraints from an early reionsation epoch, such as indicated by recent observations, may potentially limit such warm/hot components to contribute only a very small fraction to the dark matter.

This paper investigates the parameter space of theories with gauge mediated supersymmetry breaking leading to gravitino (cold) dark matter with mass m_{3/2}= 1keV - 10MeV. We pay particular attention to the cosmological role of messenger fields. Cosmology requires that these messengers decay to the visible sector if the lightest messenger mass M_X > 30TeV. We then examine the various possible messenger number violating interactions allowed by the symmetries of the theory and by phenomenology. Late messenger decay generally results in entropy production hence in the dilution of pre-existing gravitinos. We find that in SU(5) grand unification only specific messenger-matter couplings allow to produce the required amount of gravitino dark matter particles. Gravitino dark matter with the correct abundance is however expected in larger gauge groups such as SO(10) for generic non-renormalizable messenger-matter interactions and for arbritrarily high post-inflationary reheating temperatures.

Supernovae remnant shock waves could be at the origin of cosmic rays up to energies in excess of the knee ($E\simeq3\cdot 10^{15} $eV) if the magnetic field is efficiently amplified by the streaming of accelerated particles in the shock precursor. This paper follows up on a previous paper \citep{pell05} which derived the properties of the MHD turbulence so generated, in particular its anisotropic character, its amplitude and its spectrum. In the present paper, we calculate the diffusion coefficients, also accounting for compression through the shock, and show that the predicted three-dimensional turbulence spectrum $k_\perp S_{3\rm d}(k_\parallel,k_\perp)\propto k_\parallel^{-1}k_\perp^{-\alpha}$ (with $k_\parallel$ and $k_\perp$ the wavenumber components along and perpendicular to the shock normal) generally leads to Bohm diffusion in the parallel direction. However, if the anisotropy is constrained by a relation of the form $k_\parallel \propto k_\perp^{2/3}$, which arises when the turbulent energy cascade occurs at a constant rate independent of scale, then the diffusion coefficient loses its Bohm scaling and scales as in isotropic Kolmogorov turbulence. We show that these diffusion coefficients allow to account for X-ray observations of supernova remnants. This paper also calculates the modification of the Fermi cycle due to the energy lost by cosmic rays in generating upstream turbulence and the concomittant steepening of the energy spectrum. Finally we confirm that cosmic rays can produced an amplified turbulence in young SNr during their free expansion phase such that the maximal energy is close to the knee and the spectral index is close to 2.3 in the warm phase of the interstellar medium

This paper investigates the nature of the MHD turbulence excited by the streaming of accelerated cosmic rays in a shock wave precursor. The two recognised regimes (non-resonant and resonant) of the streaming instability are taken into account. We show that the non-resonant instability is very efficient and saturates through a balance between its growth and non-linear transfer. The cosmic-ray resonant instability then takes over and is quenched by advection through the shock. The level of turbulence is determined by the non-resonant regime if the shock velocity $V_{\rm sh}$ is larger than a few times $\xi_{\rm CR} c$, where $\xi_{\rm CR}$ is the ratio of the cosmic-ray pressure to the shock kinetic energy. The instability determines the dependence of the spectrum with respect to $k_\parallel$ (wavenumbers along the shock normal). The transverse cascade of Alfv\'en waves simultaneously determines the dependence in $k_{\perp}$. We also study the redistribution of turbulent energy between forward and backward waves, which occurs through the interaction of two Alfv\'en and one slow magneto-sonic wave. Eventually the spectra at the longest wavelengths are found almost proportional to $k_{\parallel}^{-1}$. Downstream, anisotropy is further enhanced through the compression at shock crossing.

This paper presents numerical simulations of test particle Fermi acceleration at relativistic shocks of Lorentz factor Gamma_sh = 2-60, using a realistic downstream magnetic structure obtained from the shock jump conditions. The upstream magnetic field is described as pure Kolmogorov turbulence; the corresponding downstream magnetic field lies predominantly in the plane tangential to the shock surface and the coherence length is smaller along the shock normal than in the tangential plane. Acceleration is nonetheless efficient and leads to powerlaw spectra with index s = 2.6-2.7 at large shock Lorentz factor Gamma_sh >> 1, markedly steeper than for isotropic scattering downstream. The acceleration timescale t_acc in the upstream rest frame becomes a fraction of Larmor time t_L in the ultra-relativistic limit, t_acc ~ 10 t_L/Gamma_sh. Astrophysical applications are discussed, in particular the acceleration in gamma-ray bursts internal and external shocks.

Recent work suggests that the cosmic ray spectrum may be dominated by Galactic sources up to ~10^{17.5} eV, and by an extra-Galactic component beyond, provided this latter cuts off below the transition energy. Here it is shown that this cut-off could be interpreted in this framework as a signature of extra-galactic magnetic fields with equivalent average strength B and coherence length l_c such that B\sqrt{l_c} ~ 2-3.10^{-10} G.Mpc^{1/2}, assuming l_c < r_L (Larmor radius at 10^{17} eV) and continuously emitting sources with density 10^{-5}/Mpc^3. The extra-Galactic flux is suppressed below 10^{17} eV as the diffusive propagation time from the source to the detector becomes larger than the age of the Universe.

This paper summarizes recent developments in the study of the propagation of ultra-high energy cosmic ray protons or nuclei in the Universe, with emphasis on the consequences of energy losses and on the influence of extra-galactic magnetic fields.

This paper presents a new method of Monte-Carlo simulations of test particle Fermi acceleration at relativistic shocks. The particle trajectories in tangled magnetic fields are integrated out exactly from entry to exit through the shock, and the conditional probability of return as a function of ingress and egress pitch angles is constructed by Monte-Carlo iteration. These upstream and downstream probability laws are then used in conjunction with the energy gain formula at shock crossing to reproduce Fermi acceleration. For pure Kolmogorov magnetic turbulence upstream and downstream, the spectral index is found to evolve smoothly from s=2.09 +/- 0.02 for mildly relativistic shocks with Lorentz factor Gamma=2 to s=2.26 +/- 0.04 in the ultra-relativistic limit Gamma >> 1. The energy gain is ~Gamma^2 at first shock crossing, and ~2 in all subsequent cycles as anticipated by Gallant & Achterberg (1999). The acceleration timescale is found to be as short as a fraction of Larmor time when Gamma >> 1.

It has recently been proposed by Mersini et al. 01, Bastero-Gil and Mersini 02 that the dark energy could be attributed to the cosmological properties of a scalar field with a non-standard dispersion relation that decreases exponentially at wave-numbers larger than Planck scale (k_phys > M_Planck). In this scenario, the energy density stored in the modes of trans-Planckian wave-numbers but sub-Hubble frequencies produced by amplification of the vacuum quantum fluctuations would account naturally for the dark energy. The present article examines this model in detail and shows step by step that it does not work. In particular, we show that this model cannot make definite predictions since there is no well-defined vacuum state in the region of wave-numbers considered, hence the initial data cannot be specified unambiguously. We also show that for most choices of initial data this scenario implies the production of a large amount of energy density (of order M_Planck^4) for modes with momenta of order M_Planck, far in excess of the background energy density. We evaluate the amount of fine-tuning in the initial data necessary to avoid this back-reaction problem and find it is of order H/M_Planck. We also argue that the equation of state of the trans-Planckian modes is not vacuum-like. Therefore this model does not provide a suitable explanation for the dark energy.

We extend existing work on the propagation of ultra-high energy cosmic rays in extragalactic magnetic fields to a possible component of heavy nuclei, taking into account photodisintegration, pion production, and creation of e+e- pairs. We focus on the influence of the magnetic field on the spectrum and chemical composition of observed ultra-high energy cosmic rays. We apply our simulations to the scenarios proposed by Anchordoqui et al, in which Iron nuclei are accelerated in nearby starburst galaxies, and show that it is in marginal agreement with the data. We also show that it is highly unlikely to detect He nuclei from M87 at the highest energies observed $\sim3 10^{20} $eV as required for the scenario of Ahn et al. in which the highest energy cosmic rays originate from M87 and are deflected in a Parker spiral Galactic magnetic field.

It is shown that fusion of neutrons and protons to He-4 nuclei occurs in gamma-ray burst outflows in a process similar to big-bang nucleosynthesis in the early Universe. Only the surviving free neutrons can then decouple kinematically from the charged fluid so that the multi-GeV neutrino signal predicted from inelastic nuclear n-p collisions is significantly reduced. It is also argued that a sizeable fraction of ultra-high energy cosmic rays accelerated in gamma-ray bursts should be He-4 nuclei.

We present an analysis of interstellar absorption along the line of sight to the nearby white dwarf star HZ43A. The distance to this star is 68+/-13 pc, and the line of sight extends toward the north Galactic pole. Column densities of OI, NI, and NII were derived from spectra obtained by the Far Ultraviolet Spectroscopic Explorer (FUSE), the column density of DI was derived from a combination of our FUSE spectra and an archival HST GHRS spectrum, and the column density of HI was derived from a combination of the GHRS spectrum and values derived from EUVE data obtained from the literature. We find the following abundance ratios (with 2-sigma uncertainties): DI/HI = (1.66 +/- 0.28) x 10^-5, OI/HI = (3.63 +/- 0.84) x 10^-4, and NI/HI = (3.80 +/- 0.74) x 10^-5. The NII column density was slightly greater than that of NI, indicating that ionization corrections are important when deriving nitrogen abundances. Other interstellar species detected along the line of sight were CII, CIII, OVI, SiII, ArI, MgII, and FeII; an upper limit was determined for NIII. No elements other than HI were detected in the stellar photosphere.

High resolution far-ultraviolet spectra of the O-type subdwarf BD +28 4211 were obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) to measure the interstellar deuterium, nitrogen, and oxygen abundances in this direction. The interstellar D I transitions are analyzed down to Lyman iota at 920.7 A. The star was observed several times at different target offsets in the direction of spectral dispersion. The aligned and coadded spectra have high signal-to-noise ratios (S/N =50-100). D I, N I, and O I transitions were analyzed with curve-of-growth and profile fitting techniques. A model of interstellar molecular hydrogen on the line of sight was derived from H2 lines in the FUSE spectra and used to help analyze some features where blending with H2 was significant. The H I column density was determined from high resolution HST/STIS spectra of Lyman alpha to be log(N HI) = 19.846+/-0.035 (2 sigma), which is higher than is typical for sight lines in the local ISM studied for D/H. We found that D/H =(1.39+/-0.21) E-5 (2 sigma) and O/H = (2.37+/-0.55) E-4 (2 sigma). O/H toward BD +28 4211 appears to be significantly below the mean O/H ratio for the ISM and the Local Bubble.

Observations obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) have been used to determine the column densities of D I, O I, and N I along seven sight lines that probe the local interstellar medium (LISM) at distances from 37 pc to 179 pc. Five of the sight lines are within the Local Bubble and two penetrate the surrounding H I wall. Reliable values of N(H I) were determined for five of the sight lines from HST data, IUE data, and published EUVE measurements. The weighted mean of D I/H I for these five sight lines is (1.52 +/- 0.08) x10-5 (1 sigma uncertainty in the mean). It is likely that the D I/H I ratio in the Local Bubble has a single value. The D I/O I ratio for the five sight lines within the Local Bubble is (3.76 +/- 0.20) x10-2. It is likely that the O I column densities can serve as a proxy for H I in the Local Bubble. The weighted mean for O I/H I for the seven FUSE sight lines is (3.03 +/-0.21) x10-4, comparable to the weighted mean (3.43 +/- 0.15) x10-4 reported for 13 sight lines probing larger distances and higher column densities (Meyer et al. 1998, Meyer 2001). The FUSE weighted mean of N I/H I for the five sight lines is half that reported by Meyer et al. (1997) for seven sight lines with larger distances and higher column densities. This result combined with the variability of O I/N I (six sight lines) indicates that at the low column densities found in the LISM, nitrogen ionization balance is important. Thus, unlike O I, N I cannot be used as a proxy for H I or as a metallicity indicator in the LISM. Subject Headings: cosmology: observations- ISM: abundances- ISM: evolution - Galaxy:abundances-Ultraviolet:ISM

We use Far Ultraviolet Spectrocopic Explorer (FUSE) observations to study interstellar absorption along the line of sight to the white dwarf WD1634-573 (d=37.1+/-2.6 pc). Combining our measurement of D I with a measurement of H I from Extreme Ultraviolet Explorer data, we find a D/H ratio toward WD1634-573 of D/H=(1.6+/-0.5)e-5. In contrast, multiplying our measurements of D I/O I=0.035+/-0.006 and D I/N I=0.27+/-0.05 with published mean Galactic ISM gas phase O/H and N/H ratios yields D/H(O)=(1.2+/-0.2)e-5 and D/H(N)=(2.0+/-0.4)e-5, respectively. Note that all uncertainties quoted above are 2 sigma. The inconsistency between D/H(O) and D/H(N) suggests that either the O I/H I and/or the N I/H I ratio toward WD1634-573 must be different from the previously measured average ISM O/H and N/H values. The computation of D/H(N) from D I/N I is more suspect, since the relative N and H ionization states could conceivably vary within the LISM, while the O and H ionization states will be more tightly coupled by charge exchange.

We present a deuterium abundance analysis of the line of sight toward the white dwarf WD2211-495 observed with the Far Ultraviolet Spectroscopic Explorer (FUSE). Numerous interstellar lines are detected on the continuum of the stellar spectrum. A thorough analysis was performed through the simultaneous fit of interstellar absorption lines detected in the four FUSE channels of multiple observations with different slits. We excluded all saturated lines in order to reduce possible systematic errors on the column density measurements. We report the determination of the average interstellar D/O and D/N ratios along this line of sight at the 95% confidence level: D/O = 4.0 +/-1.2 *10^-2; D/N = 4.4 +/-1.3 *10^-1. In conjunction with FUSE observations of other nearby sight lines, the results of this study will allow a deeper understanding of the present-day abundance of deuterium in the local interstellar medium and its evolution with time.

High-resolution spectra of the hot white dwarf G191-B2B, covering the wavelength region 905-1187A, were obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE). This data was used in conjunction with existing high-resolution Hubble Space Telescope STIS observations to evaluate the total HI, DI, OI and NI column densities along the line of sight. Previous determinations of N(DI) based upon GHRS and STIS observations were controversial due to the saturated strength of the DI Lyman-alpha line. In the present analysis the column density of DI has been measured using only the unsaturated Lyman-beta and Lyman-gamma lines observed by FUSE. A careful inspection of possible systematic uncertainties tied to the modeling of the stellar continuum or to the uncertainties in the FUSE instrumental characteristics has been performed. The column densities derived are: log N(DI) = 13.40 +/-0.07, log N(OI) = 14.86 +/-0.07, and log N(NI) = 13.87 +/-0.07 quoted with 2-sigma uncertainties. The measurement of the HI column density by profile fitting of the Lyman-alpha line has been found to be unsecure. If additional weak hot interstellar components are added to the three detected clouds along the line of sight, the HI column density can be reduced quite significantly, even though the signal-to-noise ratio and spectral resolution at Lyman-alpha are excellent. The new estimate of N(HI) toward G191-B2B reads: log N(HI) = 18.18 +/-0.18 (2-sigma uncertainty), so that the average (D/H) ratio on the line of sight is: (D/H) = 1.66 (+0.9/-0.6) *10^-5 (2-sigma uncertainty).

This article presents the derivation of the stress-energy tensor of a free scalar field with a general non-linear dispersion relation in curved spacetime. This dispersion relation is used as a phenomelogical description of the short distance structure of spacetime following the conventional approach of trans-Planckian modes in black hole physics and in cosmology. This stress-energy tensor is then used to discuss both the equation of state of trans-Planckian modes in cosmology and the magnitude of their backreaction during inflation. It is shown that gravitational waves of trans-Planckian momenta but subhorizon frequencies cannot account for the form of cosmic vacuum energy density observed at present, contrary to a recent claim. The backreaction effects during inflation are confirmed to be important and generic for those dispersion relations that are liable to induce changes in the power spectrum of metric fluctuations. Finally, it is shown that in pure de Sitter inflation there is no modification of the power spectrum except for a possible magnification of its overall amplitude independently of the dispersion relation.

The transport of charged particles in disorganised magnetic fields is an important issue which concerns the propagation of cosmic rays of all energies in a variety of astrophysical environments, such as the interplanetary, interstellar and even extra-galactic media, as well as the efficiency of Fermi acceleration processes. We have performed detailed numerical experiments using Monte-Carlo simulations of particle propagation in stochastic magnetic fields in order to measure the parallel and transverse spatial diffusion coefficients and the pitch angle scattering time as a function of rigidity and strength of the turbulent magnetic component. We confirm the extrapolation to high turbulence levels of the scaling predicted by the quasi-linear approximation for the scattering frequency and parallel diffusion coefficient at low rigidity. We show that the widely used Bohm diffusion coefficient does not provide a satisfactory approximation to diffusion even in the extreme case where the mean field vanishes. We find that diffusion also takes place for particles with Larmor radii larger than the coherence length of the turbulence. We argue that transverse diffusion is much more effective than predicted by the quasi-linear approximation, and appears compatible with chaotic magnetic diffusion of the field lines. We provide numerical estimates of the Kolmogorov length and magnetic line diffusion coefficient as a function of the level of turbulence. Finally we comment on applications of our results to astrophysical turbulence and the acceleration of high energy cosmic rays in supernovae remnants, in super-bubbles, and in jets and hot spots of powerful radio-galaxies.

We present numerical simulations for energy spectra and angular distributions of nucleons above 10^{19} eV injected by the radio-galaxy Centaurus A at a distance 3.4 Mpc and propagating in extra-galactic magnetic fields in the sub-micro Gauss range. We show that field strengths B~0.3 micro Gauss, as proposed by Farrar and Piran, cannot provide sufficient angular deflection to explain the observational data. A magnetic field of intensity ~1 micro Gauss could reproduce the observed large-scale isotropy and could marginally explain the observed energy spectrum. However, it would not readily account for the E=320 plusminus 93 EeV Fly's Eye event that was detected at an angle 136 degrees away from Cen-A. Such a strong magnetic field also saturates observational upper limits from Faraday rotation observations and X-ray bremsstrahlung emission from the ambient gas (assuming equipartition of energy). This scenario may already be tested by improving magnetic field limits with existing instruments. We also show that high energy cosmic ray experiments now under construction will be able to detect the level of anisotropy predicted by this scenario. We conclude that for magnetic fields B~0.1-0.5 micro Gauss, considered as more reasonable for the local Supercluster environment, in all likelihood at least a few sources within ~10 Mpc from the Earth should contribute to the observed ultra high energy cosmic ray flux.

We discuss the interstellar absorption lines found in FUSE spectra of the Wolf-Rayet binary Sk 108, located in the northeastern part of the main ``bar'' of the Small Magellanic Cloud. The spectra cover the wavelength range 988-1187 Angstroms, at a resolution of about 12,000 and S/N of 20--40. We use detailed component information from higher resolution near-UV and optical spectra to model the far-UV lines of similarly distributed species. Both the Galactic and SMC gas toward Sk 108 seem to be predominantly neutral, though a significant fraction of the SMC gas is ionized. The column densities of P II, S II, and Ar I are consistent with essentially solar ratios, relative to N(Zn II), in both the Galactic and SMC gas; the column density of N I remains somewhat uncertain. Molecular hydrogen is present in the Galactic gas, with properties similar to those found in low mean density Galactic lines of sight and in the Galactic gas toward several other LMC and SMC stars. We report a tentative detection of H_2 in the SMC gas for J = 1 and 3, with rotational level populations consistent with an excitation temperature of order 1000 K -- similar to the H_2 found in diffuse Galactic gas toward zeta Pup. Strong absorption from N III, S III, and Fe III has revealed a significant ionized component, particularly in the SMC; O VI is present, but relatively weak, especially in the Galactic gas. The N(C IV)/N(O VI) ratio varies somewhat within the SMC --- suggesting that several processes may contribute to the observed high ion abundances.

The metal deficient starburst galaxy IZw18 has been observed with FUSE in a search for H2 molecules. The spectrum obtained with an aperture covering the full galaxy shows no absorption lines of diffuse H2 at the radial velocity of the galaxy. The upper limit for the diffuse H2 column density is found to be very low: N(H2) <~ 10^15 cm-2 (10 sigma), unlike our Galaxy where H2 is generally present for even low HI column densities. Although the HI column density is here as high as N(HI) ~ 2.10^21 cm-2, we observe 2N(H2)/N(HI)<< 10^-6. We cannot exclude the possibility that some H2 could be in very dense, small and discrete clumps which cannot be detected with the present observation. However, the remarkable absence of diffuse H2 in this metal-poor galaxy can be explained by the low abundance of dust grains (needed to form this molecule from H-atoms), the high ultraviolet flux and the low density of the HI cloud surrounding the star-forming regions. Thus having eliminated diffuse H2 as a significant contributor to the total mass, it appears that the gas of the galaxy is dominated by HI, and that the high dynamical mass is not composed of cold and diffuse baryonic dark matter.

The Lyman and Werner band systems of deuterated molecular hydrogen (HD) occur in the far UV range below 1200 A. The high sensitivity of the FUSE mission can give access, at moderate resolution, to hot stars shining through translucent clouds, in the hope of observing molecular cores in which deuterium is essentially in the form of HD. Thus, the measurement of the HD/H2 ratio may become a new powerful tool to evaluate the deuterium abundance, D/H, in the interstellar medium. We report here on the detection of HD toward the high extinction star HD 73882 [E(B-V)=0.72]. A preliminary analysis is presented.

We report the results of initial FUSE observations of molecular hydrogen (H2) in translucent clouds. These clouds have greater optical depth than any of the diffuse clouds previously observed for far-UV H2 absorption, and provide new insights into the physics and chemistry of such regions. Our initial results involve observations of HD 73882, a well-studied southern hemisphere star lying behind substantial interstellar material (E(B-V) = 0.72; A_V = 2.44). We find a total H2 column density, N(H2) = 1.2 x 10^{21} cm^{-2}, about three times larger than the values for diffuse clouds previously measured in the far-UV. The gas kinetic temperature indicated by the ratio N(J=1)/N(J=0) is 58 +/- 10 K. With the aid of ground-based data to calculate an appropriate multi-component curve of growth, we have determined column densities for all rotational levels up to J = 7. The J >= 2 states can be reasonably fitted with a rotational excitation temperature of 307 +/- 23 K. Both the kinetic and rotational temperatures are similar to those found in previous investigations of diffuse clouds. The ratios of carbonaceous molecules to hydrogen molecules are also similar to ratios in diffuse clouds, suggesting a similar chemistry for this line of sight.

The amount of late decaying massive particles (e.g., gravitinos, moduli) produced in the evaporation of primordial black holes (PBHs) of mass $\Mbh\la10^9 $g is calculated. Limits imposed by big-bang nucleosynthesis on the abundance of these particles are used to constrain the initial PBH mass fraction $\beta$ (ratio of PBH energy density to critical energy density at formation), as: $\beta\la 5\times10^{-19} (\xp/6 10^{-3})^{-1} (\Mbh/10^9 {\rm g})^{-1/2} (\bar{\Yp}/10^{-14})$; $\xp$ is the fraction of PBH luminosity going into gravitinos or moduli, $\bar{\Yp}$ is the upper bound imposed by nucleosynthesis on the number density to entropy density ratio of gravitinos or moduli. This notably implies that such PBHs should never come to dominate the cosmic energy density.

We present new observations of the binary Sirius A / Sirius B performed with HST-GHRS. Two interstellar clouds are detected on this sightline, one of them being identified as the Local Interstellar Cloud (LIC). The interstellar structure of this sightline is constrained by high spectral resolution data of the species OI, NI, SiII, CII, FeII and MgII. Lyman alpha interstellar lines are also observed toward the two stars. But whereas the deuterium Lyman alpha line is well detected in the LIC with an abundance in agreement with that obtained by Linsky et al. (1993 & 1995), no significant DI line is detected in the other cloud. However, the Lyman alpha lines toward Sirius A and Sirius B are not trivial. An excess of absorption is seen in the blue wing of the Sirius A Lyman alpha line and interpreted as the wind from Sirius A. In its white dwarf companion, an excess in absorption is seen in the red wing and interpreted as the core of the Sirius B photospheric Lyman alpha line. A composite Lyman alpha profile can nonetheless be constructed, and allows one to measure the deuterium abundance in the second cloud 0 < D/H_ISM < 1.6 10^-5, which is marginally in agreement with the Linsky et al. (1993 & 1995) value. This sightline appears consequently as a good candidate for a low D/H_ISM.

We calculate the number density of helicity $\pm3/2$ gravitinos produced out of the vacuum by the non-static gravitational field in a generic inflation scenario. We compare it to the number density of gravitinos produced in particle interactions during reheating.

We predict energy spectra and angular distributions of nucleons above 10**(19) eV that originate from sources distributed in the Local Supercluster, which is also supposed to contain a large scale magnetic field of strength between 0.05 and 0.5 micro Gauss. We show that this model can explain all present-day features of ultra-high energy cosmic rays, at least for field strengths close to 0.5 micro Gauss. The large-scale anisotropy and the clustering predicted by this scenario will allow strong discrimination against other models with next generation experiments.

We discuss the measurements of deuterium abundances in high redshift quasar absorbers, in the solar system and in the interstellar medium. We present new results that indicate spatial variations of the deuterium abundance in the interstellar medium at the level of 50% over scales possibly as small as 10 pc, and discuss plausible causes for the origin of these variations.

We present detailed numerical simulations and analytical approximations of the propagation of nucleons above 10**(19) eV in the Local Supercluster, assuming that the ambient magnetic field is turbulent, and its strength 0.01 < B_rms < 1 micro-Gauss. In such strong magnetic fields, protons in the low energy part of the spectrum, 10**(19) eV < E < E_C diffuse, while the higher energy particles, with E > E_C propagate along nearly straight lines. The magnitude of the transition energy E_C depends mainly on the strength of the magnetic field, the coherence length, and the distance to the source; for B_rms=0.1 micro-Gauss, a largest eddy of length 10 Mpc, and a distance to the source of 10 Mpc, E_C=100 EeV. Our numerical treatment substantially improves on previous analytical approximations, as it allows to treat carefully the transition between the two propagation regimes, as well as the effects due to inhomogeneities expected on scales of a few Mpc. We show that a turbulent magnetic field B_rms=0.1 micro-Gauss, close to equipartition, would allow to reproduce exactly the observed spectrum of ultra high energy cosmic rays, up to the highest energy observed, for a distance to the source below 10 Mpc, for the geometry of the Local Supercluster, i.e. a sheet of thickness 10 Mpc. Diffusion, in this case, allows to reproduce the high flux beyond the Greisen Zatsepin Kuzmin cut-off, with a soft injection spectrum proportional to E**(-2.4). Moreover, the large deflection angles at the highest energies observed, typically 10 degrees for the above values, would explain why no close-by astrophysical counterpart could be associated with these events.

We present high spectral resolution HST-GHRS observations of the white dwarf G191-B2B, and derive the interstellar D/H ratio on the line of sight. We have observed and analysed simultaneously the interstellar lines of HI, DI, NI, OI, SiII, and SiIII. We detect three absorbing clouds, and we measure an average (D/H) ratio: N(DI)_Total / N(HI)_Total = 1.12 +/- 0.08 x 10**-5, in disagreement with the value of the D/H ratio as measured toward Capella by Linsky et al. (1993, 1995). We re-analyze the GHRS data of the Capella line of sight, and confirm their estimate, as we find: (D/H)_Capella = 1.56 +/- 0.1 x 10**-5, in the Local Interstellar Cloud in which the Sun is embedded. This shows that the D/H ratio varies by at least 30% within the local interstellar medium (not to be confused with the Local Interstellar Cloud). Furthermore, the Local Interstellar Cloud is also detected toward our target G191-B2B, and we show that the D/H ratio in this component, toward G191-B2B, can be made compatible with that derived toward Capella. However, this comes at the expense of a much smaller value for the D/H ratio as averaged over the other two components, of order 0.9 x 10**-5, and in such a way that the D/H ratio, as averaged over all three components, remains at the above value, i.e. 1.12 x 10**-5. We thus conclude that, either the D/H ratio varies from cloud to cloud, and/or the D/H ratio varies within the Local Interstellar Cloud in which the Sun is embedded, although we note that our observations neither prove nor disprove this latter possibility.

We present a detailed Monte Carlo study coupled to a likelihood analysis of the potential of next generation ultra-high energy cosmic ray experiments to reconstruct properties of the sources and the extra-galactic magnetic field. Such characteristics are encoded in the distributions of arrival time, direction, and energy of clusters of charged cosmic rays above a few 10**(19) eV. The parameters we consider for reconstruction are the emission timescale, total fluence (or power), injection spectrum, and distance of the source, as well as the r.m.s. field strength, power spectrum, and coherence length of the magnetic field. We discuss five generic situations which can be identified relatively easily and allow a reasonable reconstruction of at least part of these parameters. Our numerical code is set up such that it can easily be applied to the data from future experiments.

Recent HST/GHRS observations of the z = 0.7010 absorber toward the QSO 1718+4807 (selected because of apparently ideal characteristics for measuring D/H) yield a D/H significantly higher D/H than some recent high-redshift measurements. Our analysis indicates D/H}= 1.8 - 2.5 x 10^{-4}. This may indicate a cosmological inhomogeneity in the deuterium abundance of at least a factor of ten.

Of the light elements, the primordial abundance of deuterium, (D/H)_p, provides the most sensitive diagnostic for the cosmological mass density parameter Omega_B. Recent high redshift (D/H) measurements are highly discrepant, although this may reflect observational uncertainties. The larger (D/H) values, which imply a low Omega_B and require the Universe to be dominated by non-baryonic matter (dynamical studies indicate a higher total density parameter), cause problems for galactic chemical evolution models since they have difficulty in reproducing the large decline down to the lower present-day (D/H). Conversely, low (D/H) values imply an Omega_B greater than derived from ^7Li and ^4He abundance measurements, and may require a deuterium abundance evolution that is too low to easily explain. Here we report the first measurement at intermediate redshift, where the observational difficulties are smaller, of a gas cloud with ideal characteristics for this experiment. Our analysis of the z = 0.7010 absorber toward 1718+4807 indicates (D/H) = 2.0 +/- 0.5 x 10^{-4} which is in the high range. This and other independent observations suggests there may be a cosmological inhomogeneity in (D/H)_p of at least a factor of ten.

Protons of energies up to 10**20 eV can be subject to significant deflection and energy dependent time delay in lage scale extragalactic or halo magnetic fields of strengths comparable to current upper limits. By performing 3-dimensional Monte Carlo simulations of nucleon propagation, we show how observations of arrival direction and time distributions can be used to measure the structure and strength of large-scale magnetic fields, and constrain the nature of the source of ultra-high energy cosmic rays.

We present a numerical code designed to conduct a likelihood analysis for clusters of nucleons above 10**19 eV originating from discrete astrophysical sources such as powerful radio galaxies, gamma-ray bursts or topological defects. The code simulates the propagation of nucleons in a large-scale magnetic field and constructs the likelihood of a given observed event cluster as a function of the average time delay due to deflection in the magnetic field, the source activity time scale, the total fluence of the source, and the power law index of the particle injection spectrum. Other parameters such as the coherence length and the power spectrum of the magnetic field are also considered. We apply it to the three pairs of events above 4X10**19 eV recently reported by the Akeno Giant Air Shower Array (AGASA) experiment, assuming that these pairs were caused by nucleon primaries which originated from a common source. Although current data are too sparse to fully constrain each of the parameters considered, and/or to discriminate models of the origin of ultra-high energy cosmic rays, several tendencies are indicated. If the clustering suggested by AGASA is real, next generation experiments with their increased exposure should detect more than 10 particles per source over a few years and our method will put strong constraints on both the large-scale magnetic field parameters and the nature of these sources.

We present the results of observations of the $\lambda 3130.4$ \AA interstellar absorption line of Be II in the direction of zeta Per. The data were obtained at the Canada-France-Hawaii 3.6m Telescope using the Coud\'e f/4 Gecko spectrograph at a resolving power $\simeq 1.1 \times 10^5$, and a signal-to-noise ratio S/N $\simeq$ 2000. The Be II line is not detected, and we obtain an upper limit on the equivalent width $W_{3130.4}\leq30$ $\mu$\AA. This upper limit is 7 times below the lowest upper limit ever reported hitherto. The derived interstellar abundance is ($^9$Be/H) $\leq 7 \times 10^{-13}$, not corrected for the depletion of Be onto interstellar grains; it corresponds to an upper limit $\delta_{Be} \leq -1.5$ dex on the depletion factor of Be. As such, it argues in favour of models of formation of dust grains in stellar atmospheres.

We examine the possibility that Big Bang Nucleosynthesis (BBN) may produce non-trivial amounts of $^6$Li. If a primordial component of this isotope could be observed, it would provide a new fundamental test of Big-Bang cosmology, as well as new constraints on the baryon density of the universe. At present, however, theoretical predictions of the primordial $^6$Li abundance are extremely uncertain due to difficulties in both theoretical estimates and experimental determinations of the deuterium-plus-alpha radiative capture reaction cross-section. We also argue that present observational capabilities do not yet allow the detection of primeval $^6$Li in very metal-poor stars of the galactic halo. However, if the critical cross section is towards the upper end of its plausible range, then improvements in $^6$Li detection capabilities may allow the establishment of $^6$Li as another product of BBN. It is also noted that a primordial $^6$Li detection could help resolve current concerns about the extragalactic D/H determination.

We review the status of the measurements of the deuterium abundance from the local interstellar medium to the solar system and high redshifts absorbers toward quasars. We present preliminary results toward a white dwarf and a QSO. We conclude that the deuterium evolution from the Big-Bang to now is still not properly understood.

We explore possible depletion factors for Li6 in the three stars where it has been detected. To this end, we assume that 6Li/H scales as 16O/H throughout the galactic evolution, an assumption motivated by the observations of a similar scaling for 9Be. We examine the uncertainties attached to this modeling of 6Li galactic evolution; we also incorporate a recent evaluation of the primordial production of 6Li. The depletion factor D6, for 6Li, in HD84937, is constrained to be less than 4 at 95% c.l.. This depletion factor is an extreme upper limit on the depletion factor for 7Li. We conclude that the primordial abundance of 7Li cannot be higher that twice the observed value, since no current stellar model is able to reproduce the observed trends of 7Li with temperature and metallicity, and yet deplete 7Li by more than a factor 2, while satisfying D6<4. All three 6Li observations are in excellent agreement with all standard expectations: Big-Bang nucleosynthesis with 2<$\eta_{10}$<6.5, and standard isochrones of 6Li survival. We discuss possible deviations from our assumption of 6Li/H scaling with 16O/H, due to alpha-alpha creation of 6Li in various sites. We notably reject stellar flares as a significant source of 6Li in HD84937.

We discuss the observational status of the interstellar lithium isotopic ratio and its significance with respect to the galactic evolution of lithium.