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

The magnetorotational instability is thought to play an important role in enabling accretion in sufficiently ionized astrophysical disks. The rate at which MRI-driven turbulence transports angular momentum is related to both the strength of the amplitudes of the fluctuations on various scales and the degree of anisotropy of the underlying turbulence. This has motivated several studies of the distribution of turbulent power in spectral space. In this paper, we investigate the anisotropic nature of MRI-driven turbulence using a pseudo-spectral code and introduce novel ways to robustly characterize the underlying turbulence. We show that the general flow properties vary in a quasi-periodic way on timescales comparable to 10 inverse angular frequencies motivating the temporal analysis of its anisotropy. We introduce a 3D tensor invariant analysis to quantify and classify the evolution of the anisotropic turbulent flow. This analysis shows a continuous high level of anisotropy, with brief sporadic transitions towards two- and three-component isotropic turbulent flow. This temporal-dependent anisotropy renders standard shell-average, especially when used simultaneously with long temporal averages, inadequate for characterizing MRI-driven turbulence. We propose an alternative way to extract spectral information from the turbulent magnetized flow, whose anisotropic character depends strongly on time. This consists of stacking 1D Fourier spectra along three orthogonal directions that exhibit maximum anisotropy in Fourier space. The resulting averaged spectra show that the power along each of the three independent directions differs by several orders of magnitude over most scales, except the largest ones. Our results suggest that a first-principles theory to describe fully developed MRI-driven turbulence will likely have to consider the anisotropic nature of the flow at a fundamental level.

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here if the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here for that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.

Gamma ray bursts (GRBs) are thought to originate from highly relativistic jets. The fireball model predicts internal shocks in the jets, causing magnetic field to be amplified and particles to be accelerated. We model the effects of an asymmetric density configuration for an internal plasma collision in a quasi-parallel magnetic field. We measured electron acceleration and found that a tenuous population of electrons is accelerated to Lorentz factors of \sim 300 - close to energy equipartition with ions. We found that the filaments did not remain static, but were deflected by the Lorentz force and rolled up into small vortices, which themselves merge to form a larger current vortex.

The energetic electromagnetic eruptions observed during the prompt phase of gamma-ray bursts are attributed to synchrotron emissions. The internal shocks moving through the ultrarelativistic jet, which is ejected by an imploding supermassive star, are the likely source of this radiation. Synchrotron emissions at the observed strength require the simultaneous presence of powerful magnetic fields and highly relativistic electrons. We explore with one and three-dimensional relativistic particle-in-cell simulations the transition layer of a shock, that evolves out of the collision of two plasma clouds at a speed 0.9c and in the presence of a quasi-parallel magnetic field. The cloud densities vary by a factor of 10. The number densities of ions and electrons in each cloud, which have the mass ratio 250, are equal. The peak Lorentz factor of the electrons is determined in the 1D simulation, as well as the orientation and the strength of the magnetic field at the boundary of the two colliding clouds. The relativistic masses of the electrons and ions close to the shock transition layer are comparable as in previous work. The 3D simulation shows rapid and strong plasma filamentation behind the transient precursor. The magnetic field component orthogonal to the initial field direction is amplified in both simulations to values that exceed those expected from the shock compression by over an order of magnitude. The forming shock is quasi-perpendicular due to this amplification. The simultaneous presence of highly relativistic electrons and strong magnetic fields will give rise to significant synchrotron emissions.

The prompt emissions of gamma-ray bursts are seeded by radiating ultrarelativistic electrons. Internal shocks propagating through a jet launched by a stellar implosion, are expected to amplify the magnetic field & accelerate electrons. We explore the effects of density asymmetry & a quasi-parallel magnetic field on the collision of plasma clouds. A 2D relativistic PIC simulation models the collision of two plasma clouds, in the presence of a quasi-parallel magnetic field. The cloud density ratio is 10. The densities of ions & electrons & the temperature of 131 keV are equal in each cloud. The mass ratio is 250. The peak Lorentz factor of the electrons is determined, along with the orientation & strength of the magnetic field at the cloud collision boundary. The magnetic field component orthogonal to the initial plasma flow direction is amplified to values that exceed those expected from shock compression by over an order of magnitude. The forming shock is quasi-perpendicular due to this amplification, caused by a current sheet which develops in response to the differing deflection of the incoming upstream electrons & ions. The electron deflection implies a charge separation of the upstream electrons & ions; the resulting electric field drags the electrons through the magnetic field, whereupon they acquire a relativistic mass comparable to the ions. We demonstrate how a magnetic field structure resembling the cross section of a flux tube grows in the current sheet of the shock transition layer. Plasma filamentation develops, as well as signatures of orthogonal magnetic field striping. Localized magnetic bubbles form. Energy equipartition between the ion, electron & magnetic energy is obtained at the shock transition layer. The electronic radiation can provide a seed photon population that can be energized by secondary processes (e.g. inverse Compton).

Cold steady-state disk wind theory from near Keplerian accretion disks requires a large scale magnetic field at near equipartition strength. However the minimum magnetization has never been tested. We investigate the time evolution of an accretion disk threaded by a weak vertical magnetic field. The strength of the field is such that the disk magnetization falls off rapidly with radius. Four 2.5D numerical simulations of viscous resistive accretion disk are performed using the magnetohydrodynamic code PLUTO. In these simulations, a mean field approach is used and turbulence is assumed to give rise to anomalous transport coefficients (alpha prescription). The large scale magnetic field introduces only a small perturbation to the disk structure, with accretion driven by the dominant viscous torque. A super fast magnetosonic jet is observed to be launched from the innermost regions and remains stationary over more than 953 Keplerian orbits. The self-confined jet is launched from a finite radial zone in the disk which remains constant over time. Ejection is made possible because the magnetization reaches unity at the disk surface, due to the steep density decrease. However, no ejection is reported when the midplane magnetization becomes too small. The asymptotic jet velocity remains nevertheless too low to explain observed jets due to the negligible power carried away by the jet. Astrophysical disks with superheated surface layers could drive analogous outflows even if their midplane magnetization is low. Sufficient angular momentum would be extracted by the turbulent viscosity to allow the accretion process to continue. The magnetized outflows would be no more than byproducts, rather than a fundamental driver of accretion. However, if the midplane magnetization increases towards the center, a natural transition to an inner jet dominated disk could be achieved.

A 2D particle simulation models the collision of two electron-ion plasma clouds along a quasi-parallel magnetic field. The collision speed is 0.9c and the density ratio 10. A current sheet forms at the front of the dense cloud, in which the electrons and the magnetic field reach energy equipartition with the ions. A structure composed of a solenoidal and a toroidal magnetic field grows in this sheet. It resembles that in the cross-section of the torus of a force-free spheromak, which may provide the coherent magnetic fields in gamma-ray burst (GRB) jets needed for their prompt emissions.

Plasma processes close to SNR shocks result in the amplification of magnetic fields and in the acceleration of electrons, injecting them into the diffusive acceleration mechanism. The acceleration of electrons and the B field amplification by the collision of two plasma clouds, each consisting of electrons and ions, at a speed of 0.5c is investigated. A quasi-parallel guiding magnetic field, a cloud density ratio of 10 and a plasma temperature of 25 keV are considered. A quasi-planar shock forms at the front of the dense plasma cloud. It is mediated by a circularly left-hand polarized electromagnetic wave with an electric field component along the guiding magnetic field. Its propagation direction is close to that of the guiding field and orthogonal to the collision boundary. It has a low frequency and a wavelength that equals several times the ion inertial length, which would be indicative of a dispersive Alfven wave close to the ion cyclotron resonance frequency of the left-handed mode (ion whistler), provided that the frequency is appropriate. However, it moves with the super-alfvenic plasma collision speed, suggesting that it is an Alfven precursor or a nonlinear MHD wave such as a Short Large-Amplitude Magnetic Structure (SLAMS). The growth of the magnetic amplitude of this wave to values well in excess of those of the quasi-parallel guiding field and of the filamentation modes results in a quasi-perpendicular shock. We present evidence for the instability of this mode to a four wave interaction. The waves developing upstream of the dense cloud give rise to electron acceleration ahead of the collision boundary. Energy equipartition between the ions and the electrons is established at the shock and the electrons are accelerated to relativistic speeds.

Stellar jets are normally constituted by chains of knots with some periodicity in their spatial distribution, corresponding to a variability of order of several years in the ejection from the protostar/disk system. A widely accepted theory for the presence of knots is related to the generation of internal working surfaces due to variations in the jet ejection velocity. In this paper we study the effect of variations in the inner disk-wind radius on the jet ejection velocity. We show that a small variation in the inner disk-wind radius produce a variation in the jet velocity large enough to generate the observed knots. We also show that the variation in the inner radius may be related to a variation of the stellar magnetic field.

Magnetic reconnection plays a critical role in many astrophysical processes where high energy emission is observed, e.g. particle acceleration, relativistic accretion powered outflows, pulsar winds and probably in dissipation of Poynting flux in GRBs. The magnetic field acts as a reservoir of energy and can dissipate its energy to thermal and kinetic energy via the tearing mode instability. We have performed 3d nonlinear MHD simulations of the tearing mode instability in a current sheet. Results from a temporal stability analysis in both the linear regime and weakly nonlinear (Rutherford) regime are compared to the numerical simulations. We observe magnetic island formation, island merging and oscillation once the instability has saturated. The growth in the linear regime is exponential in agreement with linear theory. In the second, Rutherford regime the island width grows linearly with time. We find that thermal energy produced in the current sheet strongly dominates the kinetic energy. Finally preliminary analysis indicates a P(k) 4.8 power law for the power spectral density which suggests that the tearing mode vortices play a role in setting up an energy cascade.

We investigate potential models that could explain why multiple proto-stellar systems predominantly show single jets. During their formation, stars most frequently produce energetic outflows and jets. However, binary jets have only been observed in a very small number of systems. We model numerically 3D binary jets for various outflow parameters. We also model the propagation of jets from a specific source, namely L1551 IRS 5, known to have two jets, using recent observations as constraints for simulations with a new MHD code. We examine their morphology and dynamics, and produce synthetic emission maps. We find that the two jets interfere up to the stage where one of them is almost destroyed or engulfed into the second one. We are able to reproduce some of the observational features of L1551 such as the bending of the secondary jet. While the effects of orbital motion are negligible over the jets dynamical timeline, their interaction has significant impact on their morphology. If the jets are not strictly parallel, as in most observed cases, we show that the magnetic field can help the collimation and refocusing of both of the two jets.

We present a 127deg x 63deg extinction map of the Anticenter of the Galaxy, based on <J-H> and <H-K> colour excess maps from 2MASS. This 8001 square degree map with a resolution of 4 arcminutes is provided as online material. The colour excess ratio <J-H> / <H-K> is used to determine the power law index of the reddening law (\beta) for individual regions contained in the area (e.g. Orion, Perseus, Taurus, Auriga, Monoceros, Camelopardalis, Cassiopeia). On average we find a dominant value of \beta=1.8+-0.2 for the individual clouds, in agreement with the canonical value for the interstellar medium. We also show that there is an internal scatter of \beta values in these regions, and that in some areas more than one dominant \beta value is present. This indicates large scale variations in the dust properties. The analysis of the A_V values within individual regions shows a change in the slope of the column density distribution with distance. This can either be attributed to a change in the governing physical processes in molecular clouds on spatial scales of about 1pc or an A_V dilution with distance in our map.

In this work, the gas infall and the formation of outflows around low and high mass protostars are investigated. A radial self-similar approach to model the transit of the molecular gas around the central object is employed. We include gravitational and radiative fields to produce heated pressure-driven outflows with magneto-centrifugal acceleration and collimation. Outflow solutions with negligible or vanishing magnetic field are reported. They indicate that thermodynamics is a sufficient engine to generate an outflow. The magnetized solutions show dynamically significant differences in the axial region, precisely where the radial velocity and collimation are the largest. They compare quantitatively well with observations. The influence of the opacity on the transit solutions is also studied. It is found that, when dust is not the dominant coolant, such as in the primordial universe, mass infall rates have substantial larger values in the equatorial region. This suggests that star forming in a dust-free environment should be able to accrete much more mass and become more massive than present day protostars.It is also suggested that molecular outflows may be dominated by the global transit of material around the protostar during the very early stages of star formation, especially in the case of massive or dust-free star formation.

We present three 14400 square degree relative extinction maps of the Galactic Plane (|b|<20degrees) obtained from 2MASS using accumulative star counts (Wolf diagrams). This method is independent of the colour of the stars and the variation of extinction with wavelength. Stars were counted in 3.5'x3.5' boxes, every 20". 1x1degree surrounding fields were chosen for reference, hence the maps represent local extinction enhancements and ignore any contribution from the ISM or very large clouds. Data reduction was performed on a Beowulf-type cluster (in approximately 120 hours). Such a cluster is ideal for this type of work as areas of the sky can be independently processed in parallel. We studied how extinction depends on wavelength in all of the high extinction regions detected and within selected dark clouds. On average a power law opacity index (\beta) of 1.0 to 1.8 in the NIR was deduced. The index however differed significantly from region to region and even within individual dark clouds. That said, generally it was found to be constant, or to increase, with wavelength within a particular region.

We present a large-scale study of the IC1396 region using new deep NIR and optical images, complemented by 2MASS data. For ten globules in IC1396 we determine (H-K, J-H) colour-colour diagrams and identify the young stellar population. Five of these globules contain a rich population of reddened objects, most of them probably young stellar objects. Two new HH objects (HH865 and HH864) could be identified by means of [SII] emission, one of them a parsec-scale flow. Using star counts based on 2MASS data we create an extinction map of the whole region. This map is used to identify 25 globules and to estimate their mass. The globule masses show a significant increase with the distance from the exciting O6.5V star HD206267. We explain this correlation by the enhanced radiation pressure close to this star, leading to evaporation of the nearby clouds and hence smaller globule masses. We see evidence that the radiation from HD206267 has a major impact on the star formation activity in these globules.