We outline recent progress in understanding the accretion of plasma to rotating magnetized stars obtained from global axisymmetric (2D) and 3D magnetohydrodynamic (MHD) simulations in three main areas: (1.) Formation of jets from disk accretion onto rotating magnetized stars: From simulations where the viscosity and magnetic diffusivity within the disk are described by alpha models, we find long-lasting conical outflows/jets from the disk/magnetosphere boundary in both the case where the star is slowly rotating and where it is rapidly rotating (the “propeller regime”). Most of the mass flux in the outflows is in a hollow cone but inside this cone there is a low-density high-velocity magnetically dominated flow along the open polar field lines of the star. The outflows occur under conditions where the poloidal magnetic flux of the star is bunched up by the accretion disk near the disk/magnetosphere boundary. Recent simulations show that the conical outflows become well-collimated for axial distances of ≲ 20 times the inner disk radius. Exploratory 3D simulations show that conical winds are axisymmetric about the rotational axis (of the star and the disk), even when the dipole field of the star is significantly misaligned. (2.) Formation of intrinsically one-sided jets from disk accretion to rotating magnetized stars: There is strong observational evidence for an asymmetry between the approaching and receding jets from a number of young stars. We discuss the first MHD simulations of the formation asymmetric or one-sided jets arising from disk accretion to a rotating star with an asymmetric (dipole plus quadrupole) magnetic field. (3.) Global axisymmetric and 3D simulations of the magnetorotational instability (MRI) in disk accretion onto magnetized stars: In the axisymmetric simulations we observe cases where there is episodic or quasi-periodic burst of accretion similar to that observed in one X ray source. In 3D MHD simulations of accretion onto stars with tilted dipole fields using our Godunov-type code based on the “cubed sphere” grid we find that the density distribution is much less smooth than in the case of the laminar accretion flow described by α–viscosity. Instead, large turbulent cells dominate the flows and are strongly elongated in the azimuthal direction.

Within the framework of laboratory astrophysics, we form a qualified multidisciplinary group in radiative hydrodynamics. Since 10 years, we have developed laboratory experiments as radiative shocks and plasma jets in connection to astrophysics. Such laboratory experiments provide a unique opportunity to validate models and numerical schemes introduced in radiative hydrodynamics codes. Here we summarize our experimental researches about plasma jets. Laboratory astrophysical experiments have been performed using LULI2000 (France), VULCAN (UK) and GEKKO XII (Japan) intense lasers. The goal of these experiments is to investigate some of the complex features of jets from Young Stellar Objects (YSO), and in particular its interaction with the interstellar medium (ISM).

Magnetic helicity has received considerable attention in the area of fluid dynamics. Recently, this quantity is attracting the interest of solar physicists and much research has been carried out related to magnetic helicity generation and transport through different solar layers, starting from the interior and the convection zone, towards the photosphere, the corona and finally into the heliosphere. Taking into account the global importance of supergranular cells in convection theories, we study the motion of magnetic features into such a geometrical element simplified as hexagonal cell and we analyse the results in terms of the accumulated magnetic helicity. We compute the emergence of a bipole inside the hexagonal cell and its motion from the centre of the cell towards its sides and its vertices, where the magnetic elements are considered to be sinking down. Multiple bipoles are also considered and phenomena such as cancellation, coalescence and fragmentation are also investigated. We find that the most important process for the accumulation of magnetic helicity is the shear motion between the polarities. The magnetic helicity accumulation changes its trend when one polarity reaches the side of the hexagon, and later the vertex. It has zero value when there is no shear motion inside the hexagonal cell, and it is constant when there is no shear between the two polarities during their motion along the cell sides.

Jet formation MHD simulations are presented considering a variety of model setups. The first approach investigates the interrelation between the disk magnetisation profile and jet collimation. Our results suggest (and quantify) that outflows launched from a very concentrated region at the inner disk tend to be weakly collimated. In the second approach, jet formation is investigated from a magnetic field configuration consisting of a stellar dipole superposed by a strong disk field. We find that the central dipole considerably de-collimates the disk wind. In addition, reconnection flares are launched in the interaction region of disk and stellar magnetic field, subsequently changing the outflow mass flux by factors of two. The time interval between flare ejection is about 1000 Keplerian periods - surprisingly similar to the observed time lag between jet knots. The third approach considers radiative pressure effects on jet collimation - an environment which is interesting mainly for outflows from massive young stars (but also for relativistic jets). Finally we present relativistic MHD simulations of jet formation from accretion disks extenting the previous non-relativistic approaches.

In this work, we present the analysis results using UMRAO preliminary data base. We used the light curves 1) to get the shortest timescales and then to get the brightness temperature so that we can estimate the Doppler factors; 2) to investigate the periodicity and discuss the variability index. We also used the data base to discuss the polarization properties of blazars. We found that the periodicity distribution in BL Lacs and that in the flat spectrum radio quasars should be from the same distribution. The Doppler factor in FSRQs is higher than that in BL. The polarization in BLs are higher than that in the flat spectrum radio quasars

We report our recent results from multiwavelengths studies of microquasars, focusing on X-ray data of GX 339–4 and GRS 1915+105 obtained with Suzaku and other observatories. The broad band coverage and high energy resolution achieved with Suzaku (or a combination of Chandra/HETGS and RXTE) enable us to perform the most reliable spectral analysis both on iron-K features and continuum, and thus to best constrain the accretion disk structure of microquasars and its relation to the jet formation at various mass accretion rates.

Simultaneous multi-wavelength observations are crucial for understanding the physics of microquasars, especially the accretion disk/jet connection. The enigmatic microquasar Cygnus X-3 exhibits strong, relativistic jet ejection events producing radio flares up to 20 Jy. These events are preceded by a very soft X-ray state with quenched emission in the radio and hard X-ray bands. Recently, GeV flux was observed by the AGILE and Fermi γ-ray observatories during the newly-identified hypersoft state. By using an extensive database of simultaneous multi-wavelength observations gathered from Cygnus X-3 we can form a more unified picture of the nature of the source and show how the recent γ-ray observations fit into it.

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneities, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in our two-dimensional simulations. The magnetic energy spectrum is flatter than Kolmogorov and indicates that a so-called small-scale dynamo is operating in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field.

In this published note I attempt to sketch my understanding of the universal working scheme of all the astrophysical jet sources, or ‘bipolar flows’, on both stellar and galactic scales, also called ‘microquasars’, and ‘quasars’. A crucial building block will be their medium: extremely relativistic e±-pair plasma performing quasi loss-free E × B-drifts through self-rammed channels, whose guiding equi-partition E- and B-fields convect the electric potential necessary for eventual single-step post-acceleration, at their terminating ‘knots’, or ‘hotspots’. The indispensible pair plasma is generated in magnetospheric reconnections of the central rotator. Already for this reason, black holes cannot serve as jet engines.

GRS 1915+105 is a very peculiar black hole binary that exhibits accretion-related states that are not observed in any other stellar-mass black hole system. One of these states, however – referred to as the plateau state – may be related to the canonical hard state of black hole X-ray binaries. Both the plateau and hard state are associated with steady, relatively lower X-ray emission and flat/inverted radio emission, that is sometimes resolved into compact, self-absorbed jets. To investigate the relationship between the plateau and the hard state, we fit two multi-wavelength observations using a steady-state outflow-dominated model, developed for hard state black hole binaries. The data sets consist of quasi-simultaneous observations in radio, near-infrared and X-ray bands. Interestingly, we find both significant differences between the two plateau states, as well as between the best-fit model parameters and those representative of the hard state. We discuss our interpretation of these results, and the possible implications for GRS 1915+105's relationship to canonical black hole candidates.

We have searched for quasi-periodic oscillations in the hard X-ray emission of solar flares. We have selected 14 flare events which were divided into two groups: a) the events with the X-ray sources located at the flare loop footpoints and b) the events with the X-ray source above the solar limb, i.e. with the loop-top X-ray source. We found that while in the case with the footpoints X-ray sources the quasi-periods of the recorded oscillations were in the interval 2–380 s, in the events with loop-top sources only the quasi-periods longer than 50 s were recognized. These results are probably connected with the MHD oscillation modes of the flaring loop. While the long periods, which are dominant in loop-top sources, are produced by acoustic oscillations along the whole long loop, in the layers close to the loop footpoints also the MHD wave modes in shorter structures with shorter periods are generated.

Instabilities of nonuniform flows is a fundamental problem in dynamics of fluids and plasmas. This presentation outlines atypical dynamics of instabilities for unmagnetized and magnetized astrophysical differentially rotating flows, including, our efforts in the development of general theory of magneto rotation instability (MRI) that takes into account plasma compressibility, pressure anisotropy, dissipative and kinetic effects. Presented analysis of instability (transient growth) processes in unmagnetized/hydrodynamic astrophysical disks is based on the breakthrough of the hydrodynamic community in the 1990s in the understanding of shear flow non-normality induced dynamics. This analysis strongly suggests that the so-called bypass concept of turbulence, which has been developed by the hydrodynamic community for spectrally stable shear flows, can also be applied to Keplerian disks. It is also concluded that the vertical stratification of the disks is an important ingredient of dynamical processes resulting onset of turbulence.

The correlations between the rest frame peak of the νFν spectrum of GRBs (Epeak) and their isotropic energy (Eiso) or luminosity (Liso) could have several implications for the understanding of the GRB prompt emission. These correlations are presently founded on the time–averaged spectral properties of a sample of 95 bursts, with measured redshifts, collected by different instruments in the last 13 years (pre–Fermi). One still open issue is wether these correlations have a physical origin or are due to instrumental selection effects. By studying 10 long and 14 short GRBs detected by Fermi we find that a strong time–resolved correlation between Epeak and the luminosity Liso is present within individual GRBs and that it is consistent with the time–integrated correlation. This result is a direct proof of the existence in both short and long GRBs of a similar physical link between the hardness and the luminosity which is not due to instrumental selection effects. The origin of the Epeak – Liso correlation should be searched in the radiation mechanism of the prompt emission.

Collimated outflows (jets) are ubiquitous in the universe, appearing around sources as diverse as protostars and extragalactic supermassive black holes. Jets are thought to be magnetically collimated, and launched from a magnetized accretion disk surrounding a compact gravitating object. We have developed the first laboratory experiment to address time-dependent, episodic phenomena relevant to the poorly understood jet acceleration and collimation region (Ciardi et al., 2009). The experiments were performed on the MAGPIE pulsed power facility (1.5 MA, 250 ns) at Imperial College. The experimental results show the periodic ejections of magnetic bubbles naturally evolving into a heterogeneous jet propagating inside a channel made of self-collimated magnetic cavities. The results provide a unique view of the possible transition from a relatively steady-state jet launching to the observed highly structured outflows.

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs in the shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and for particle acceleration. These magnetic fields contribute to the electron's transverse deflection behind the shock. The “jitter” radiation from deflected electrons in turbulent magnetic fields has different properties from synchrotron radiation calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure of gamma-ray bursts, relativistic jets in general, and supernova remnants. In order to calculate radiation from first principles and go beyond the standard synchrotron model, we have used PIC simulations. We will present detailed spectra for conditions relevant to various astrophysical sites of collisionless shock formation. In particular we will discuss application to GRBs and SNRs.

K 3-35 is a very young planetary nebula (PN) with a characteristic S-shaped radio emission morphology. It is the first PN where water vapor maser was detected: the emission is located in a torus-like structure with a radius of 100 AU and also at the surprisingly large distance of 5000 AU from the star, in the tips of the bipolar lobes. Several mechanism have been proposed to explain the bipolar morphology of PNe, and in the case of K 3-35 we believe we may be observing several of them at the same time: i) a disk-like structure traced by the H2O masers, ii) a precessing bipolar jet probably due to the presence of a binary companion and iii) circular polarization in the OH 1665 MHz masers, which suggests the presence of a magnetic field. Additional observations and modeling are needed to establish what mechanisms are shaping K 3-35.

Magnetic reconnection (Parker, 1957; Sweet, 1958; Petschek, 1964; Yamada et al., 2010; Biskamp, 2000; Tsuneta, 1996; Kivelson and Russell, 1995; Yamada, 2007; Birn et al., 2001; Drake et al., 2003) is considered important to many astrophysical phenomena including stellar flares, magnetospheric disruptions of magnetars, and dynamics of galactic lobes. Research on magnetic reconnection started with observations in solar coronae and in the Earths magnetosphere, and a classical theory was developed based on MHD. Recent progress has been made by understanding the two-fluid physics of reconnection, through space and astrophysical observations (Tsuneta, 1996; Kivelson and Russell, 1995), laboratory experiments (Yamada, 2007), and theory and numerical simulations (Birn et al., 2001; Daughton et al., 2006; Uzdensky and Kulsrud, 2006). Laboratory experiments dedicated to the study of the fundamental reconnection physics have tested the physics mechanisms and their required conditions, and have provided a much needed bridge between observations and theory. For example, the Magnetic Reconnection Experiment (MRX) experiment (http://mrx.pppl.gov) has rigorously cross-checked the leading theories though quantitative comparisons of the numerical simulations and space astrophysical observations (Mozer et al., 2002). Extensive data have been accumulated in a wide plasma parameter regime with Lundquist numbers of S = 100 − 3000, where S is a ratio of the magnetic diffusion time to the Alfven transit time.

Several galaxy clusters are known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths. It suggests that the AGN precessing jets play a role in the formation of the misaligned bubbles. Also, X-ray spectra reveal that typically these systems are also able to supress cooling flows, predicted theoretically. The absence of cooling flows in galaxy clusters has been a mistery for many years since numerical simulations and analytical studies suggest that AGN jets are highly energetic, but are unable to redistribute it at all directions. We performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, in which dynamics is coupled to a NFW dark matter gravitational potential. Radiative cooling has been taken into account and the cooling flow problem was studied. We found that precession is responsible for multiple pairs of bubbles, as observed. The misaligned bubbles rise up to scales of tens of kiloparsecs, where the thermal energy released by the jets are redistributed. After ~150 Myrs, the temperature of the gas within the cavities is kept of order of ~107 K, while the denser plasma of the intracluster medium at the central regions reaches T ~ 105 K. The existence of multiple bubbles, at diferent directions, results in an integrated temperature along the line of sight much larger than the simulations of non-precessing jets. This result is in agreement with the observations. The simulations reveal that the cooling flows cessed ~50–70 Myr after the AGN jets are started.

Strong and variable radiation detected over all accessible energy bands in blazar arises from a relativistic, Doppler-boosted jet pointing close to our line of sight. Flat Spectrum Radio Quasar 3C 279 was one of the brightest γ-ray blazars in the sky at the time of the discovery with EGRET. Since the successful launch of the Fermi Gamma-ray Space telescope in 2008, we have organized extensive multi-band observational campaign of 3C 279 from radio to γ-ray bands, also including optical polarimetric observations. The uninterrupted monitoring in the γ-ray band by Fermi-LAT together with the multi-band data provide us with new insights of the relativistic jet of blazar. Here, we present the results of the first-year multi-band campaign of 3C 279 including the discovery of a γ-ray flare event associated with a dramatic change of the optical polarization - as well as a discovery of an “orphan” X-ray flare, unassociated with prominent outbursts in other bands.

We present results of a broadband spectroscopy of the galactic microquasars and black hole candidates XTE J1550-564 and GRO J1655-40, performed with the INTEGRAL and RXTE observatories during strong outbursts in 2003 and 2005, respectively. The spectral parameters evolution was traced during brightening and fading phases of each outburst to search a possible hysteresis and transitions from state to state. We estimated a size and optical depth of different regions around XTE J1550-564, like a hot plasma zone and optically thick accretion disk. Upper limits to the annihilation 511 keV line emission were obtained for both sources using data of the SPI spectrometer onboard the INTEGRAL observatory.