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.

We use particle-in-magnetohydrodynamics-cells to model particle acceleration and magnetic field amplification in a high Mach, parallel shock in three dimensions and compare the result to 2-D models. This allows us to determine whether 2-D simulations can be relied upon to yield accurate results in terms of particle acceleration, magnetic field amplification and the growth rate of instabilities. Our simulations show that the behaviour of the gas and the evolution of the instabilities are qualitatively similar for both the 2-D and 3-D models, with only minor quantitative differences that relate primarily to the growth speed of the instabilities. The main difference between 2-D and 3-D models can be found in the spectral energy distributions (SEDs) of the non-thermal particles. The 2-D simulations prove to be more efficient, accelerating a larger fraction of the particles and achieving higher velocities. We conclude that, while 2-D models are sufficient to investigate the instabilities in the gas, their results have to be treated with some caution when predicting the expected SED of a given shock.

The rather elusive High-Frequency Quasi-Periodic Oscillations(HFQPO) observed in the X-ray lightcurve of black holes have been seen in a wide range of frequencies, even within one source. It is also notable to have been detected in "pairs" of HFQPOs with a close to integer ratio between the frequencies. The aim of this paper is to investigate some of the possible observable that we could obtain from having the Rossby Wave Instability (RWI) active in the accretion disc surrounding the compact object. Using the newly developed GR-AMRVAC code able to follow the evolution of the RWI in a full general relativistic framework, we explore how RWI can reproduce observed HFQPO frequencies ratios and if it is compatible with the observations. In order to model the emission coming from the disc we have linked our general relativistic simulations to the general relativistic ray-tracing GYOTO code and delivered synthetic observables that can be confronted to actual data from binary systems hosting HFQPOs.}{We have demonstrated in our study that some changes in the physical conditions prevailing in the part of the disc where RWI can be triggered leads to various dominant RWI modes whose ratio recovers frequency ratios observed in various X-ray binary systems. In addition, we have also highlighted that when RWI is triggered near the last stable orbit of a spinning black hole, the amplitude of the X-ray modulation increases with the spin of the black hole. Revisiting published data on X-ray binary systems, we show that this type of relationship actually exists in the five systems where an indirect measure of the spin of the black hole is available.

We have performed general relativistic hydrodynamics (GRHD) simulations of 2D discs orbiting around spinning black holes and prone to the Rossby Wave Instability (RWI). We show that the RWI can develop at any location in the disc and for any spin parameter. After recovering the overall patterns of the instability in this general relativistic context, we have analysed its development and identified some modifications induced by the combined effects of the relativistic rotational profile of the disc and local time dilatation that affects the propagation of waves in the disc. We have found in particular that the saturation level of the instability increases significantly when RWI occurs in the very close vicinity of fast-rotating black holes where general relativistic effects are strong. Such finding suggests that even more strongly than in the case of Schwarzschild black-hole, it is necessary to complement such GRHD simulations with a full GR ray-tracing processing in order to provide synthetic observations of the disc in the distant observer frame.

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.

We present simulations of magnetized astrophysical shocks taking into account the interplay between the thermal plasma of the shock and supra-thermal particles. Such interaction is depicted by combining a grid-based magneto-hydrodynamics description of the thermal fluid with particle in cell techniques devoted to the dynamics of supra-thermal particles. This approach, which incorporates the use of adaptive mesh refinement features, is potentially a key to simulate astrophysical systems on spatial scales that are beyond the reach of pure particle-in-cell simulations. We consider in this study non-relativistic shocks with various Alfvenic Mach numbers and magnetic field obliquity. We recover all the features of both magnetic field amplification and particle acceleration from previous studies when the magnetic field is parallel to the normal to the shock. In contrast with previous particle-in-cell-hybrid simulations, we find that particle acceleration and magnetic field amplification also occur when the magnetic field is oblique to the normal to the shock but on larger timescales than in the parallel case. We show that in our simulations, the supra-thermal particles are experiencing acceleration thanks to a pre-heating process of the particle similar to a shock drift acceleration leading to the corrugation of the shock front. Such oscillations of the shock front and the magnetic field locally help the particles to enter the upstream region and to initiate a non-resonant streaming instability and finally to induce diffuse particle acceleration.

In an early work the Rossby Wave Instability (RWI) has been proposed to explain variability thought to originate in the close vicinity of black-holes but this was done in the pseudo-Newtonian approach. Here we present the first general relativistic hydrodynamics simulations of this instability not only proving its theorized existence in a full general relativistic (GR) environment but also studying the effect of the strong gravity on the instability. To that end we performed a set of simulations increasingly closer to the black hole with our new GR version of the MPI-AMRVAC code.This allows us to study the minute changes in the behaviour of the instability.We found that the pseudo-Newtonian approach gives adequate results provided that time shifting induced by the black hole gravity is taken into account. Hence, to view the disc as a distant observer would a full GR ray-tracing post treatment of the simulations is a must.

Classical Supergiant X-ray Binaries host a neutron star orbiting a supergiant OB star and display persistent X-ray luminosities of 10$^{35}$ to 10$^{37}$ erg/s. The stellar wind from the massive companion is believed to be the main source of matter accreted by the compact object. With this first paper, we introduce a ballistic model to characterize the structure of the wind at the orbital scale as it accelerates, from the stellar surface to the vicinity of the accretor. Thanks to the parametrization we retained and the numerical pipeline we designed, we can investigate the supersonic flow and the subsequent observables as a function of a reduced set of characteristic numbers and scales. We show that the shape of the permanent flow is entirely determined by the mass ratio, the filling factor, the Eddington factor and the $\alpha$-force multiplier which drives the stellar wind acceleration. Provided scales such as the orbital period are known, we can trace back the observables to evaluate the mass accretion rates, the accretion mechanism, the shearing of the inflow and the stellar parameters. We confront our model to three persistent Supergiant X-ray Binaries (Vela X-1, IGR J18027-2016, XTE J1855-026) and discuss the likelihood of wind-formed accretion discs around the accretors in each case, further investigated in a following paper.

Bondi-Hoyle accretion configurations occur as soon as a gravitating body is immersed in an ambient medium with a supersonic relative velocity. From wind-accreting X-ray binaries to runaway neutron stars, such a regime has been witnessed many times and is believed to account for shock formation, the properties of which can be only marginally derived analytically. In this paper, we present the first results of the numerical characterization of the stationary flow structure of Bondi-Hoyle accretion onto a compact object, from the large scale accretion radius down to the vicinity of the compact body. For different Mach numbers, we study the associated bow shock. It turns out that those simulations confirm the analytical prediction by Foglizzo & Ruffert (1996) concerning the topology of the inner sonic surface with an adiabatic index of 5/3. They also enable us to derive the related mass accretion rates, the position and the temperature of the bow shock, as function of the flow parameters, along with the transverse density and temperature profiles in the wake.

The black hole at the Galactic centre exhibits regularly flares of radiation, the origin of which is still not understood. In this article, we study the ability of the near-future GRAVITY infrared instrument to constrain the nature of these events. We develop realistic simulations of GRAVITY astrometric data sets for various flare models. We show that the instrument will be able to distinguish an ejected blob from alternative flare models, provided the blob inclination is >= 45deg, the flare brightest magnitude is 14 <= mK <= 15 and the flare duration is >= 1h30.

We present in this paper both a linear study and numerical relativistic MHD simulations of the non-resonant streaming instability occurring in the precursor of relativistic shocks. In the shock front restframe, we perform a linear analysis of this instability in a likely configuration for ultra-relativistic shock precursors. This considers magneto-acoustic waves having a wave vector perpendicular to the shock front and the large scale magnetic field. Our linear analysis is achieved without any assumption on the shock velocity and is thus valid for all velocity regimes. In order to check our calculation, we also perform relativistic MHD simulations describing the propagation of the aforementioned magneto-acoustic waves through the shock precursor. The numerical calculations confirm our linear analysis, which predicts that the growth rate of the instability is maximal for ultra-relativistic shocks and exhibits a wavenumber dependence $\propto k_x^{1/2}$. Our numerical simulations also depict the saturation regime of the instability where we show that the magnetic amplification is moderate but nevertheless significant ($\delta B/B\leq 10$). This latter fact may explain the presence of strong turbulence in the vicinity of relativistic magnetized shocks. Our numerical approach also introduces a convenient means to handle isothermal (ultra-)relativistic MHD conditions.

We investigate under which conditions cold, fan-shaped winds can be steadily launched from thin (Keplerian) accretion discs. Such winds are magneto-centrifugal winds launched from a thin annulus in the disc, along open magnetic field lines that fan out above the disc. In principle, such winds could be found in two situations: (1) at the interface between an inner Jet Emitting Disc, which is itself powering magneto-centrifugally driven winds, and an outer standard accretion disc; (2) at the interface between an inner closed stellar magnetosphere and the outer standard accretion disc. We refer to Terminal or T-winds to the former kind and to Magnetospheric or M-winds to the latter. The full set of resistive and viscous steady state MHD equations are analyzed for the disc (the annulus), which allow us to derive general expressions valid for both configurations. We find that, under the framework of our analysis, the only source of energy able to power any kind of fan-shaped winds is the viscous transport of rotational energy coming below the inner radii. Using standard local $\alpha$ prescriptions for the anomalous (turbulent) transport of angular momentum and magnetic fields in the disc, we derive the strength of the transport coefficients that are needed to steadily sustain the global configuration. It turns out that, in order for these winds to be dynamically relevant and explain observed jets, the disc coefficients must be far much larger than values expected from current knowledge of turbulence occurring inside proto-stellar discs. Either the current view on MHD turbulence must be deeply reconsidered or steady-state fan-shaped winds are never realized in Nature. The latter hypothesis seems to be consistent with current numerical simulations.

In the context of planet formation, anticyclonic vortices have recently received lots of attention for the role they can play in planetesimals formation. Radial migration of intermediate size solids toward the central star may prevent their growth to larger solid grains. On the other hand, vortices can trap the dust and accelerate this growth, counteracting fast radial transport. Multiple effects have been shown to affect this scenario, such as vortex migration or decay. The aim of this paper is to study the formation of vortices by the Rossby wave instability and their long term evolution in a full three dimensional protoplanetary disc. We use a robust numerical scheme combined with adaptive mesh refinement in cylindrical coordinates, allowing to affordably compute long term 3D evolutions. We consider a full disc stratified both radially and vertically that is prone to formation of vortices by the Rossby wave instability. We show that the 3D Rossby vortices grow and survive over hundreds of years without migration. The localized overdensity which initiated the instability and vortex formation survives the growth of the Rossby wave instability for very long times. When the vortices are no longer sustained by the Rossby wave instability, their shape changes toward more elliptical vortices. This allows them to survive shear-driven destruction, but they may be prone to elliptical instability and slow decay. When the conditions for growing Rossby wave-related instabilities are maintained in the disc, large-scale vortices can survive over very long timescales and may be able to concentrate solids.

Molecular counterparts to atomic jets have been detected within 1000 AU of young stars. Reproducing them is a challenge for proposed ejection models. We explore whether molecules may survive in an MHD disk wind invoked to reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. The coupled ionization, chemical and thermal evolution along dusty flow streamlines is computed for a prescribed MHD disk wind solution, using a method developed for magnetized shocks in the interstellar medium. Irradiation by wind-attenuated coronal X-rays and FUV photons from accretion hot spots is included, with self-shielding of H2 and CO. Disk accretion rates of 5e-6, 1e-6 and 1e-7 solar masses per year are considered, representative of low-mass young protostars (Class 0), evolved protostars (Class I) and very active T Tauri stars (Class II). The disk wind has an onion-like thermo-chemical structure, with streamlines launched from larger radii having lower temperature and ionisation, and higher H2 abundance. The coupling between charged and neutral fluids is sufficient to eject molecules from the disk out to 9 AU. The launch radius beyond which most H2 survives moves outward with evolutionary stage. CO survives in the Class 0 but is significantly photodissociated in the Class I/II. Balance between ambipolar heating and molecular cooling establishes an asymptotic temperature 700-3000 K, with cooler jets at earlier protostellar stages. Endothermic formation of H2O is efficient with abundances up to 1e-4, while CH+ and SH+ can exceed 1e-6 in the Class I/II winds. A centrifugal MHD disk wind launched from beyond 0.2-1 AU can produce molecular jets/winds up to speeds 100 km/s in young low-mass stars. The model predicts a high ratio H2/CO and an increase of molecular launch radius, temperature, and flow width as the source evolves, in agreement with current observed trends.

Planet formation and migration in accretion discs is a very active topic. Among the many aspects related to that question, dead zones are of particular importance as they can influence both the formation and the migration of planetary embryos. The ionisation level in the disc is the key element in determining the existence and the location of the dead zone. This has been studied either within the Standard Accretion Disc (SAD) framework or using parametric discs. In this paper, we extend this study to the case of Jet Emitting Discs (JED), the structure of which strongly differ from SADs because of the new energy balance and angular momentum extraction imposed by the jets. We make use of the (r,z) density distributions provided by self-similar accretion-ejection models, along with the JED thermal structure derived in a previous paper, to create maps of the ionisation structure of JEDs. We compare the ionisation rates we obtain to the critical value required to trigger the magneto-rotational instability. It is found that JEDs have a much higher ionisation degree than SADs which renders very unlikely the presence of a dead zone in these discs. As JEDs are believed to occupy the inner regions of accretion discs, the extension of the dead zones published in the literature should be re-considered for systems in which a jet is present. Moreover, since JEDs require large scale magnetic fields close to equipartition, our findings raise again the question of magnetic field advection in circumstellar accretion discs.

Context. The formation of vortices in accretion disks is of high interest in various astrophysical contexts, in particular for planet formation or in the disks of compact objects. But despite numerous attempts it has thus far not been possible to produce strong vortices in fully three-dimensional simulations of disks. Aims. The aim of this paper is to present the first 3D simulation of a strong vortex, established across the vertically stratified structure of a disk by the Rossby Wave Instability. Methods. Using the Versatile Advection Code (VAC), we set up a fully 3D cylindrical stratified disk potentially prone to the Rossby Wave Instability. Results. The simulation confirms the basic expectations obtained from previous 2D analytic and numerical works. The simulation exhibits a strong vortex that grows rapidly and saturates at a finite amplitude. On the other hand the third dimension shows unexpected additional behaviours that could be of strong importance in the astrophysical roles that such vortices can play.

The present article investigates magnetic amplification in the upstream medium of SNR blast wave through both resonant and non-resonant regimes of the streaming instability. It aims at a better understanding of the diffusive shock acceleration (DSA) efficiency considering various relaxation processes of the magnetic fluctuations in the downstream medium. Multi-wavelength radiative signatures coming from the SNR shock wave are used in order to put to the test the different downstream turbulence relaxation models. We confirm the result of Parizot et al (2006) that the maximum CR energies should not go well beyond PeV energies in young SNRs where X-ray filaments are observed. In order to match observational data, we derive an upper limit on the magnetic field amplitude insuring that stochastic particle reacceleration remain inefficient. Considering then, various magnetic relaxation processes, we present two necessary conditions to achieve efficient acceleration and X-ray filaments in SNRs: 1/the turbulence must fulfil the inequality $2-\beta-\delta_{\rm d} \ge 0$ where $\beta$ is the turbulence spectral index while $\delta_d$ is the relaxation length energy power-law index; 2/the typical relaxation length has to be of the order the X-ray rim size. We identify that Alv\'enic/fast magnetosonic mode damping does fulfil all conditions while non-linear Kolmogorov damping does not. Confronting previous relaxation processes to observational data, we deduct that among our SNR sample, the older ones (SN1006 & G347.3-0.5) fail to verify all conditions which means that their X-ray filaments are likely controlled by radiative losses. The younger SNRs, Cas A, Tycho and Kepler, do pass all tests and we infer that the downstream magnetic field amplitude is lying in the range of 200-300 $\mu$ Gauss.

Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind the Mach disk, for various jet Lorentz factors and magnetic field helicities. We investigate the influence of helical magnetic fields on jet beam propagation in overdense external medium. We use the AMRVAC code, employing a novel hybrid block-based AMR strategy, to compute ideal plasma dynamics in special relativity. The helicity of the beam magnetic field is effectively transported down the beam, with compression zones in between diagonal internal cross-shocks showing stronger toroidal field regions. In comparison with equivalent low-relativistic jets which get surrounded by cocoons with vortical backflows filled by mainly toroidal field, the high speed jets demonstrate only localized, strong toroidal field zones within the backflow vortical structures. We find evidence for a more poloidal, straight field layer, compressed between jet beam and backflows. This layer decreases the destabilizing influence of the backflow on the jet beam. In all cases, the jet beam contains rich cross-shock patterns, across which part of the kinetic energy gets transferred. For the high speed reference jet considered here, significant jet deceleration only occurs beyond distances exceeding ${\cal O}(100 R_j)$, as the axial flow can reaccelerate downstream to the internal cross-shocks. This reacceleration is magnetically aided, due to field compression across the internal shocks which pinch the flow.

We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. Our simulations show that the inner outflow is accelerated from the central object hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force. We also show that when a dense inner stellar wind occurs, the resulting disc-driven jet have a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass ejection rate in the disc-driven jet and a larger radial extension which is in better agreement with the observations besides being more consistent.

We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $\gamma=100$ up to the almost Newtonian $\gamma\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.

We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar magnetospheric winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. These innovative MHD simulations are devoted to the study of the interplay between a stellar wind (having different ejection mass rates) and an MHD disc-driven jet embedding the stellar wind. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. We first investigate the various angular momentum transports occurring in the magneto-viscous accretion disc. We then analyze the modifications induced by the interaction between the two components of the outflow. Our simulations show that the inner outflow is accelerated from the central object's hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force.}{The simulations show that the MHD disc-driven outflow extracts angular momentum more efficiently than do viscous effects in near-equipartition, thin-magnetized discs where turbulence is fully developed. We also show that, when a dense inner stellar wind occurs, the resulting disc-driven jet has a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass-ejection rate in the disc-driven jet and a larger radial extension that is in better agreement with the observations, besides being more consistent.

Gamma-ray bursts (GRB), at least those with a duration longer than a few seconds are the most energetic events in the Universe and occur at cosmological distances. The ECLAIRs micro-satellite, to be launched in 2009, will provide multi-wavelength observations of GRB, to study their astrophysics and to use them as cosmological probes. Furthermore in 2009 ECLAIRs is expected to be the only space borne instrument capable of providing a GRB trigger in near real-time with sufficient localization accuracy for GRB follow-up observations with the powerful ground based spectroscopic telescopes available by then. A "Phase A study" of the ECLAIRs project has recently been launched by the French Space Agency CNES, aiming at a detailed mission design and selection for flight in 2006. The ECLAIRs mission is based on a CNES micro-satellite of the "Myriade" family and dedicated ground-based optical telescopes. The satellite payload combines a 2 sr field-of-view coded aperture mask gamma-camera using 6400 CdTe pixels for GRB detection and localization with 10 arcmin precision in the 4 to 50 keV energy band, together with a soft X-ray camera for onboard position refinement to 1 arcmin. The ground-based optical robotic telescopes will detect the GRB prompt/early afterglow emission and localize the event to arcsec accuracy, for spectroscopic follow-up observations.

The present paper deals with the yield and transport of high-energy particle within extragalactic jet terminal shocks, also known as hotspots. We investigate in some details the cosmic ray, neutrinos and high-energy photons yield in hotspots of powerful FRII radio-galaxies by scanning all known spatial transport regimes, adiabatic and radiative losses as well as Fermi acceleration process. Since both electrons and cosmic rays are prone to the same type of acceleration, we derive analytical estimates of the maximal cosmic ray energy attainable in both toroidal and poloidal magnetic field dominated shock structures by using observational data on synchrotron emission coming from various hot-spots. One of our main conclusions is that the best hot-spot candidates for high energy astroparticle production is the extended ($L_{HS}\geq 1kpc$), strongly magnetized ($B> 0.1mG$) terminal shock displaying synchrotron emission cut-off lying at least in the optical band. We found only one object (3C273A) over the six objects in our sample being capable to produce cosmic rays up to $10^{20}$ eV. Secondly, we investigate the astroparticle spectra produced by two characteric hot-spots (Cygnus A and 3C273 A) by applying a multi-scale MHD-kinetic scheme, coupling MHD simulations to kinetic computations using stochastic differential equations. We show that 3C273 A, matching the previous properties, may produce protons up to $10^{20}$ eV in a Kolmogorov type turbulence by both computing electron and cosmic ray acceleration. We also calculate the high-energy neutrino and gamma-ray fluxes on Earth produced through p-$\gamma$ and p-p processes and compare them to the most sensitive astroparticle experiments.

Motivated by the first detections of rotation signatures in the DG Tau jet (Bacciotti et al. 2002), we examine possible biases affecting the relation between detected rotation signatures and true azimuthal velocity for self-similar MHD disc winds, taking into account projection, convolution as well as excitation gradients effects. We find that computed velocity shifts are systematically smaller than the true underlying rotation curve. When outer slower streamlines dominate the emission, we predict observed shifts increasing with transverse distance to the jet axis, opposite to the true rotation profile. Determination of the full transverse rotation profile thus requires high angular resolution observations (< 5 AU) on an object with dominant inner faster streamlines. Comparison of our predictions with HST/STIS observations of DG Tau clearly shows that self-similar, warm MHD disc wind models with lambda = 13 and an outer radius of the disc ~ 3 AU are able to reproduce detected velocity shifts, while cold disc wind models (lambda > 50) are ruled out for the medium-velocity component in the DG Tau jet.

We present self-similar models of resistive viscous Keplerian disks driving non-relativistic magnetohydrodynamics (MHD) jets becoming super fast-magnetosonic. We show that in order to obtain such solutions, the thermal pressure must be a sizeable fraction of the poloidal magnetic pressure at the Alfven surface. These steady solutions which undergo a recollimation shock causally disconnected from the driving engine, account for structures with a high temperature plasma in the sub-Alfvenic region. We suggest that only unsteady outflows with typical time-scales of several disk dynamical time-scales can be produced if the suitable pressure conditions are not fulfilled.

We present magnetohydrodynamic simulations of a resistive accretion disk continuously launching transmagnetosonic, collimated jets. We time-evolve the full set of magnetohydrodynamic equations, but neglect radiative losses in the energetics (radiatively inefficient). Our calculations demonstrate that a jet is self-consistently produced by the interaction of an accretion disk with an open, initially bent large-scale magnetic field. A constant fraction of heated disk material is launched in the inner equipartition disk regions, leading to the formation of a hot corona and a bright collimated, super-fastmagnetosonic jet. We illustrate the complete dynamics of the ``hot'' near steady-state outflow (where thermal pressure $\simeq$ magnetic pressure) by showing force balance, energy budget and current circuits. The evolution to this near stationary state is analyzed in terms of the temporal variation of energy fluxes controlling the energetics of the accretion disk. We find that unlike advection-dominated accretion flow, the energy released by accretion is mainly sent into the jet rather than transformed into disk enthalpy. These magnetized, radiatively inefficient accretion-ejection structures can account for under-luminous thin disks supporting bright fast collimated jets as seen in many systems displaying jets (for instance M87).

Multidimensional magneto-hydrodynamical (MHD) simulations coupled with stochastic differential equations (SDEs) adapted to test particle acceleration and transport in complex astrophysical flows are presented. The numerical scheme allows the investigation of shock acceleration, adiabatic and radiative losses as well as diffusive spatial transport in various diffusion regimes. The applicability of SDEs to astrophysics is first discussed in regards to the different regimes and the MHD code spatial resolution. The procedure is then applied to 2.5D MHD-SDE simulations of kilo-parsec scale extragalactic jets. The ability of SDE to reproduce analytical solutions of the diffusion-convection equation for electrons is tested through the incorporation of an increasing number of effects: shock acceleration, spatially dependent diffusion coefficients and synchrotron losses. The SDEs prove to be efficient in various shock configuration occurring in the inner jet during the development of the Kelvin-Helmholtz instability. The particle acceleration in snapshots of strong single and multiple shock acceleration including realistic spatial transport is treated. In chaotic magnetic diffusion regime, turbulence levels $\eta_T=<\delta B^2>/(B^2+<\delta B^2>)$ around $0.2-0.3$ are found to be the most efficient to enable particles to reach the highest energies. The spectrum, extending from 100 MeV to few TeV (or even 100 TeV for fast flows), does not exhibit a power-law shape due to transverse momentum dependent escapes. Out of this range, the confinement is not so efficient and the spectrum cut-off above few hundreds of GeV, questioning the Chandra observations of X-ray knots as being synchrotron radiation. The extension to full time dependent simulations to X-ray extragalactic jets is discussed.

We present numerical magnetohydrodynamic (MHD) simulations of a magnetized accretion disk launching trans-Alfvenic jets. These simulations, performed in a 2.5 dimensional time-dependent polytropic resistive MHD framework, model a resistive accretion disk threaded by an initial vertical magnetic field. The resistivity is only important inside the disk, and is prescribed as eta = alpha_m V_AH exp(-2Z^2/H^2), where V_A stands for Alfven speed, H is the disk scale height and the coefficient alpha_m is smaller than unity. By performing the simulations over several tens of dynamical disk timescales, we show that the launching of a collimated outflow occurs self-consistently and the ejection of matter is continuous and quasi-stationary. These are the first ever simulations of resistive accretion disks launching non-transient ideal MHD jets. Roughly 15% of accreted mass is persistently ejected. This outflow is safely characterized as a jet since the flow becomes super-fastmagnetosonic, well-collimated and reaches a quasi-stationary state. We present a complete illustration and explanation of the `accretion-ejection' mechanism that leads to jet formation from a magnetized accretion disk. In particular, the magnetic torque inside the disk brakes the matter azimuthally and allows for accretion, while it is responsible for an effective magneto-centrifugal acceleration in the jet. As such, the magnetic field channels the disk angular momentum and powers the jet acceleration and collimation. The jet originates from the inner disk region where equipartition between thermal and magnetic forces is achieved. A hollow, super-fastmagnetosonic shell of dense material is the natural outcome of the inwards advection of a primordial field.

A complete analytical and numerical treatment of all magnetohydrodynamic waves and instabilities for radially stratified, magnetized accretion disks is presented. The instabilities are a possible source of anomalous transport. While recovering results on known hydrodynamicand both weak- and strong-field magnetohydrodynamic perturbations, the full magnetohydrodynamic spectra for a realistic accretion disk model demonstrates a much richer variety of instabilities accessible to the plasma than previously realized. We show that both weakly and strongly magnetized accretion disks are prone to strong non-axisymmetric instabilities.The ability to characterize all waves arising in accretion disks holds great promise for magnetohydrodynamic spectroscopic analysis.

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 summarize results on the thermal and ionization structure of self-similar, magnetically-driven, atomic disk winds heated by ambipolar diffusion. We improve upon earlier work by Safier by considering (1) new MHD solutions consistent with underlying cold keplerian disk equilibrium, (2) a more accurate treatment of the micro-physics, and (3) predictions for spatially resolved forbidden line emission (maps, long-slit spectra). The temperature plateau $\simeq 10^4$ K found earlier is recovered, but ionization fractions are revised downward by a factor of 10, due to previous omission of thermal speeds in ion-neutral momentum-exchange rates. The physical origin of the temperature plateau is outlined. Predictions are then compared with T~Tauri star observations, with emphasis on the necessity of suitable beam convolution. Jet widths and variations in line profiles with distance and line tracer are well reproduced. However, predicted maximum velocities are too high, total densities too low, and the low-velocity [OI] component is too weak. Denser, slower MHD winds from warm disks might resolve these discrepancies.

We present a self-similar, steady-state model describing both a magnetized accretion disc and a magnetohydrodynamic jet. We focus on the role of a hot corona in such a structure. This corona enables the disc to launch various types of jets. By considering the energy conservation, we also present a diagnostic of the luminosity of the magnetized disc, which could explain some observational signatures of galactic objects.

In this paper, steady-state MHD calculations of non-relativistic magnetized accretion discs driving jets are presented. For the first time, an energy equation describing the effects of entropy generation along streamlines is included. Using a simplified approach, we showed that thermal effects have a tremendous influence on the disc capability to feed jets with mass. The disc ejection efficiency is measured by the parameter $\xi= d \ln \dot{M}_{a}/d \ln{r}$, where $ \dot{M}_{a}(r)$ is the local disc accretion rate. While previous self-similar solutions were only able to produce jets with $\xi \sim 0.01$, solutions with a coronal heating display huge efficiencies up to $\xi \sim 0.5$. Moreover, taking thermal effects into account allows to obtain both fast and slow magnetic rotators. Since most of the jet properties (like asymptotic velocity or degree of collimation) depend on the mass load, it arises from this study that any quantitative result requires a detailled analysis of the disc energetics.

We present steady-state calculations of self-similar magnetized accretion discs driving cold, adiabatic, non-relativistic jets. For the first time, both the magnetic torque due to the jets and a turbulent "viscous" torque are taken into account. This latter torque allows a dissipation of the accretion power as radiation at the disc surfaces, while the former predominantly provides jets with power. The parameter space of these structures has been explored. It is characterized by four free parameters, namely the disc aspect ratio and three MHD turbulence parameters, related to the anomalous magnetic diffusivities and viscosity. It turns out that launching cold jets from thin, dissipative discs implies anisotropic turbulent dissipation. Jets that asymptotically reach a high Alfvenic Mach number are only produced by weakly dissipative discs. We obtained general analytical relations between disc and jet quantities that must be fulfilled by any steady-state model of cold jets, launched from a large radial extension of thin discs. We also show that such discs cannot have a dominant viscous torque. This is because of the chosen geometry, imposing the locus of the Alfven surface. Some observational consequences of these cold magnetized accretion-ejection structures are also briefly discussed.