We present a new formulation to compute numerically stationary and axisymmetric equilibria of magnetized and self-gravitating astrophysical fluids. Under the assumption of ideal MHD, the stream function for the flow can be chosen as a basic variable with which the Euler-Maxwell equations are cast into a set of basic equations, i.e. a generalized Bernoulli equation and a Grad-Shafranov-like equation by employing various integral conditions. A novel feature of this formulation is that systems with stars, disks and winds are treated in a simple unified picture and the magnetic field structures can contain both poloidal and toroidal components.

We have found evidence for interaction between a standing and a traveling shock in the jet of the blazar CTA 102. Our result is based in the study of the spectral evolution of the turnover frequency-turnover flux density (νm, Sm) plane. The radio/mm light curves were taken during a major radio outburst in April 2006.

Understanding the physical mechanisms that play a role in the saturation of the magnetorotational instability (MRI) has been an outstanding problem in accretion physics since the early 90's. Here, we present the summary of a study of the parasitic modes that feed off exact viscous, resistive MRI modes. We focus on the situation in which the amplitude of the magnetic field produced by the MRI is such that the instantaneous growth rate of the fastest parasitic mode matches that of the fastest MRI mode. We argue that this "saturation" amplitude provides an estimate of the magnetic field that can be generated by the MRI before the secondary instabilities suppress its growth significantly. We show that there exist two regimes, delimited by a critical Elsasser number of order unity, in which saturation is achieved via secondary instabilities that correspond to either Kelvin-Helmholtz or tearing modes.

A spherical hydrodynamical expansion flow can be described as the gradient of a potential. In that case no vorticity should be produced, but several additional mechanisms can drive its production. Here we analyze the effects of baroclinicity, rotation and shear in the case of a viscous fluid. Those flows resemble what happens in the interstellar medium. In fact in this astrophysical environment supernovae explosion are the dominant flows and, in a first approximation, they can be seen as spherical. One of the main difference is that in our numerical study we examine only weakly supersonic flows, while supernovae explosions are strongly supersonic.

Jets are produced by young stellar objects (YSOs), by black hole binary star system “microquasars” (μQSOs), by active galactic nuclei (AGN), are associated with neutron stars and pulsar wind nebulae (PWNe), and are thought responsible for the gamma-ray bursts (GRBs). An understanding of these outflows must include how they are launched and collimated into jets, and how they propagate to large distances. Jets be they Poynting flux and/or kinetic flux dominated are current driven (CD) and/or Kelvin-Helmholtz (KH) velocity shear driven unstable. Here I present some of the work that is leading to a better understanding of the properties required for the observed relative stability of astrophysical jets.

The Solar System formation PFO–CFO hypothesis is developed in the direction of creation of a phenomenological model focused on solution of a number of paradoxes and answering to a number of mysterious questions under the same cover. For explanation of the events and processes that occurred over the period from the middle ages of the pre-solar star to the Solar System formation, original approaches are applied.

Using a new formulation to compute structures of stationary and axisymmetric magnetized barotropic stars in Newtonian gravity, we have succeeded in obtaining numerically exact models of stars with extremely high interior magnetic fields. In this formulation, there appear four arbitrary functions of the magnetic flux function from the integrability conditions among the basic equations. Since in our new formulation these arbitrary functions appear in the expression of the current density, configurations with different current distributions can be specified by choosing the forms of the arbitrary functions.

By choosing appropriate forms for the four arbitrary functions, we have solved many kinds of equilibrium configurations both with poloidal and toroidal magnetic fields. Among them, by choosing special form for the toroidal current density, we have been able to obtain magnetized stars which have extremely strong poloidal magnetic fields deep inside the core region near the symmetric axis. By adopting the appropriate model parameters for the neutron stars, the magnetic fields could be 1014 ~ 1015 G on the surfaces and be about 1017 G in the deep interior regions. For other model parameters appropriate for white dwarfs, the magnetic fields could be around 107 ~ 108 G (surface regions) and 109 ~ 1010 G (core regions). It is remarkable that the regions with very strong interior magnetic fields are confined to a very narrow region around the symmetric axis in the central part of the stars. The issues of stability of these configurations and of evolutionary paths to reach such configurations need to be investigated in the future work.

3D-MHD numerical simulations of bipolar, hypersonic, weakly magnetized jets and synthetic synchrotron observations are presented to study the structure and evolution of magnetic fields in FR II radio sources. The magnetic field setup in the jet is initially random. The power of the jets as well as the observational viewing angle are investigated. We find that synthetic polarization maps agree with observations and show that magnetic fields inside the sources are shaped by the jets' backflow. Polarimetry statistics correlates with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than for ones with fatter cocoons. Jets increase the magnetic energy in cocoons, in proportion to the jet velocity. Both, filaments in synthetic emission maps and 3D magnetic power spectra suggest that turbulence develops in evolved sources.

We present results of recent observations and theoretical modeling of data from black holes accreting at very low luminosities (L/LEdd ≲ 10−8). We discuss our newly developed time-dependent model for episodic ejection of relativistic plasma within a jet framework, and a successful application of this model to describe the origin of radio flares seen in Sgr A*, the Galactic center black hole. Both the observed time lags and size-frequency relationships are reproduced well by the model. We also discuss results from new Spitzer data of the stellar black hole X-ray binary system A0620–00. Complemented by long term SMARTS monitoring, these observations indicate that once the contribution from the accretion disk and the donor star are properly included, the residual mid-IR spectral energy distribution of A0620–00 is quite flat and consistent with a non-thermal origin. The results above suggest that a significant fraction of the observed spectral energy distribution originating near black holes accreting at low luminosities could result from a mildly relativistic outflow. The fact that these outflows are seen in both stellar-mass black holes as well as in supermassive black holes at the heart of AGNs strengthens our expectation that accretion and jet physics scales with mass.

In recent years, numerous efforts have been devoted to unravel the connection between accretion flow and jets in accreting compact objects. Here we report new constraints on these issues, through the long term study of the radio and X-ray behaviour of the black hole candidate H 1743–322. This source is known to be one of the “outliers” of the universal radio/X-ray correlation, i.e. a group of stellar mass accreting black holes displaying fainter radio emission for a given X-ray luminosity, than expected from the correlation. In this work we find, at high X-ray luminosity in the hard state, a tight radio/X-ray correlation with an unusual steep slope of b = 1.38 ± 0.03. This correlation then breaks below ~5 × 10−3LEdd (M/10M⊙)−1 in X-rays and becomes shallower. When compared with radio/X-ray data from other black hole X-ray binaries, we see that the deviant points of H 1743–322 join the universal correlation and seem to follow it at low luminosity. Based on these results, we investigate several hypotheses that could explain both the b ~ 1.4 slope and the transition toward the universal correlation.

The fractal shape and multi-component nature of the interstellar medium together with its vast range of dynamical scales provides one of the great challenges in theoretical and numerical astrophysics. Here we will review recent progress in the direct modelling of interstellar hydromagnetic turbulence, focusing on the role of energy injection by supernova explosions. The implications for dynamo theory will be discussed in the context of the mean-field approach.

Results obtained with the test field-method are confronted with analytical predictions and estimates from quasilinear theory. The simulation results enforce the classical understanding of a turbulent Galactic dynamo and, more importantly, yield new quantitative insights. The derived scaling relations enable confident global mean-field modelling.

Of all possible black hole sources, the event horizon of the Galactic Center black hole, Sgr A*, subtends the largest angular scale on the sky. It is therefore a prime candidate to study and image plasma processes in strong gravity and it even allows imaging of the shadow cast by the event horizon. Recent mm-wave VLBI and radio timing observations as well as numerical GRMHD simulations now have provided several breakthroughs that put Sgr A* back into the focus. Firstly, VLBI observations have now measured the intrinsic size of Sgr A* at multiple frequencies, where the highest frequency measurements have approached the scale of the black hole shadow. Moreover, measurements of the radio variability show a clear time lag between 22 GHz and 43 GHz. The combination of size and timing measurements, allows one to actually measure the flow speed and direction of magnetized plasma at some tens of Schwarzschild radii. This data strongly support a moderately relativistic outflow, consistent with an accelerating jet model. This is compared to recent GRMHD simulation that show the presence of a moderately relativistic outflow coupled to an accretion flow Sgr A*. Further VLBI and timing observations coupled to simulations have the potential to map out the velocity profile from 5-40 Schwarzschild radii and to provide a first glimpse at the appearance of a jet-disk system near the event horizon. Future submm-VLBI experiments would even be able to directly image those processes in strong gravity and directly confirm the presence of an event horizon.

The strength and structure of cosmic magnetic fields is best studied by observations of radio continuum emission, its polarization and its Faraday rotation. Fields with a well-ordered spiral structure exist in many types of galaxies. Total field strengths in spiral arms and bars are 20–30 μG and dynamically important. Strong fields in central regions can drive gas inflows towards an active nucleus. The strongest regular fields (10–15 μG) are found in interarm regions, sometimes forming “magnetic spiral arms” between the optical arms. The typical degree of polarization is a few % in spiral arms, but high (up to 50%) in interarm regions. The detailed field structures suggest interaction with gas flows. Faraday rotation measures of the polarization vectors reveals large-scale patterns in several spiral galaxies which are regarded as signatures of large-scale (coherent) fields generated by dynamos. – Polarization observations with the forthcoming large radio telescopes will open a new era in the observation of magnetic fields and should help to understand their origin. Low-frequency radio synchrotron emission traces low-energy cosmic ray electrons which can propagate further away from their origin. LOFAR (30–240 MHz) will allow us to map the structure of weak magnetic fields in the outer regions and halos of galaxies, in galaxy clusters and in the Milky Way. Polarization at higher frequencies (1–10 GHz), to be observed with the EVLA, MeerKAT, APERTIF and the SKA, will trace magnetic fields in the disks and central regions of galaxies in unprecedented detail. All-sky surveys of Faraday rotation measures towards a dense grid of polarized background sources with ASKAP and the SKA are dedicated to measure magnetic fields in distant intervening galaxies and clusters, and will be used to model the overall structure and strength of the magnetic field in the Milky Way.

Dense wind of a massive star can be partially captured by a companion neutron star (NS) creating a very turbulent and magnetized transition region at some distance from the NS surface. We consider the consequences of electron and hadron acceleration at such a transition region. Electrons lose energy on the synchrotron process and the inverse Compton (IC) scattering of thermal radiation from the NS surface and/or the massive star. We calculate the synchrotron spectra (from X-rays to soft γ-rays) and IC spectra in the case of sources accreting the matter under the accretor and propeller scenarios. It is argued that a population of accreting massive binaries, recently discovered by the INTEGRAL observatory, can be detectable by the Fermi LAT telescope. On the other hand, TeV γ-ray emission from other class of massive binaries can be interpreted in terms of a magnetar accreting matter in the propeller scenario. We also calculate the expected neutrino event rates in a km2 detector produced by relativistic hadrons accelerated in such scenario.

We present the first applications of a new time-dependent multi-zone jet radiation transfer code to the study the multiwavelength emission of the TeV Blazar Mrk 421. The code couples Fokker-Planck and Monte Carlo methods. For the first time all light travel time effects are fully considered as well as proper self-consistent treatment of Compton cooling, which depends on them. The first tests focus on the March 2001 observations of Mrk 421, still one of the best datasets available for phenomenology and X-ray/TeV data coverage. We summarize the results of scenarios of variability induced by injection of relativistic electrons in a blob encountering a shock, and with different combinations with a second component, either co-spatial or independent from the active region.

In this paper we investigate the quasi periodic oscillation (QPO) behavior of the black hole candidate GX 339-4 during its 2010 outburst using RXTE/PCA data. We perform spectral and timing analysis of the observations, where the QPOs are observed. We analyze the relationship between the centroid frequency of QPO and the spectral parameters. The correlation of spectral and timing properties can be used to estimate the mass of black hole with the scaling method. Using this method we estimate a mass of 7.5 ± 0.8 M⊙ of GX 339-4.

It is analytically shown that passages of comets near the Sun's surface with velocities more than 600 km/s is accompanied by aerodynamic crushing of their nuclei within the solar chromosphere and transversal expansion of the crushed matter. The deceleration of the flattened hypervelocity body within the solar photosphere has sharply impulsive and strongly explosive character. The specific energy release in the explosion zone near the solar surface 10-100 thousand times exceeds the evaporation heat of the nucleus material, so that the process is accompanied by generation of high-temperature plasma and non-stationary explosive phenomena around the photosphere. Spectral observations of these phenomena by SOHO and SDO type space observatories with high spatial and temporal resolutions are of interest for the plasma astrophysics as well as the physics of solar flares.

Large regions of protoplanetary discs are believed to be too weakly ionised to support magnetorotational instabilities, because abundant tiny dust grains soak up free electrons and reduce the conductivity of the gas. At the outer edge of this “dead zone”, the ionisation fraction increases gradually and the resistivity drops until the magnetorotational instability can develop turbulence. We identify a new viscous instability which operates in the semi-turbulent transition region between “dead” and “alive” zones. The strength of the saturated turbulence depends strongly on the local resistivity in this transition region. A slight increase (decrease) in dust density leads to a slight increase (decrease) in resistivity and a slight decrease (increase) in turbulent viscosity. Such spatial variation in the turbulence strength causes a mass pile-up where the turbulence is weak, leading to a run-away process where turbulence is weakened and mass continues to pile up. The final result is the appearance of high-amplitude pressure bumps and deep pressure valleys. Here we present a local linear stability analysis of weakly ionised accretion discs and identify the linear instability responsible for the pressure bumps. A paper in preparation concerns numerical results which confirm and expand the existence of the linear instability.

The key characteristics of molecular and atomic ejection from young stars are summarized, with emphasis on similarities across evolutionary stages, and the need for efficient magnetic collimation and ejection extracting a large fraction of the accretion power. The jet kinematics, and its dust and molecular content, are confronted to steady MHD jet models, and the probable contribution of non-steady processes is pointed out.

We used near-IR integral field spectroscopy, obtained with Gemini NIFS and GNIRS integral field units (IFUs), to map the ionized and molecular flux distributions and kinematics in the central few hundreds of parsecs of Seyfert galaxies. We conclude that the molecular gas emission can be considered a tracer of the feeding of the AGN, while the emission of the ionized gas a tracer of its feedback.