We present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with the circumbinary environment. We focus on studying the impact of the high wind magnetization on the geometrical structure of the shocked flow. This work is an extension of our earlier studies that focused on a purely hydrodynamic interaction and weak wind magnetization. We consider the collision in the axisymmetric approximation, that is, the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star. The increase of the magnetization results in the expansion of the opening cone in which the shocked pulsar wind propagates. This effect is explained in the frameworks of the conventional theory of collimation of magnetized winds. This finding has a direct implication for scenarios that involve Doppler boosting as the primary mechanism behind the GeV flares detected with the Fermi/LAT from PSR B1259-63/LS2883. The maximum enhancement of the apparent emission is determined by the ratio of $4\pi$ to the solid in which the shocked pulsar wind propagates. Our simulations suggest that this enhancement factor is decreased by the impact of the magnetic field.

The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations imposed by the process of conversion of VHE gamma rays into $e^{\pm}$ pairs in the magnetic field of the pulsar magnetosphere. Photons with energies exceeding 1~TeV are effectively absorbed even at large distances from the surface of the neutron star. Our calculations of magnetic absorption in the force-free magnetosphere show that the twisting of the magnetic field due to the pulsar rotation makes the magnetosphere more transparent compared to the dipole magnetosphere. The gamma-ray absorption appears stronger for photons emitted in the direction of rotation than in the opposite direction. There is a small angular cone inside which the magnetosphere is relatively transparent and photons with energy $1.5$~TeV can escape from distances beyond $0.1$~light cylinder radius ($R_{\rm{lc}}$). The emission surface from where photons can be emitted in the observer's direction further restricts the sites of VHE gamma-ray production. For the observation angle $57^{\circ}$ relative to the Crab pulsar axis of rotation and the orthogonal rotation, the emission surface in the open field line region is located as close as $0.4\,R_{\rm{lc}}$ from the stellar surface for a dipole magnetic field, and $0.1\,R_{\rm{lc}}$ for a force-free magnetic field.

Disk accretion onto black holes is accompanied by collimated outflows (jets). In active galactic nuclei (AGN), the kinetic energy flux of the jet may exceed the bolometric luminosity of the disk a few orders of magnitude. This phenomena can be explained in frameworks of so called "cold" disk accretion when the only source of energy of AGN is the energy released at the accretion. The radiation from the disk is suppressed because the wind from the disk carries out almost all the angular momentum and the gravitational energy of the accreted material. In this paper, we calculate an unavoidable radiation from the "cold" disk and the ratio of the kinetic energy luminosity of the outflow to the bolometric luminosity of the accretion disk around a super massive black hole in framework of the paradigm of the optically thick alpha-disk of Shakura \& Sunyaev. The ratio of the luminosities is defined predominantly by the ratio of the magnetic field pressure inside the disk to the magnetic field pressure at the base of the wind. The obtained equations applied to the jet of M87 demonstrate good agreement with observations. In the case of Sgr A*, these equations allow us to predict the kinetic energy flux from the disk around Galactic SMBH.

The very high energy (VHE) gamma-ray emission reported from a number of pulsar wind nebulae (PWNe) is naturally explained by the inverse Compton scattering of multi-TeV electrons. However, the physical dimensions of some gamma-ray-emitting PWNe significantly exceed the scales anticipated by the standard hydrodynamical paradigm of PWN formation. The most "disturbing" case in this regard is HESS J1825-137, which extends to distances $r\approx70\rm\,pc$ from the central pulsar PSR J1826-1334. If the gamma-ray emission is indeed produced inside the PWN, but not by electrons that escaped the nebula and diffuse in the interstellar medium (ISM), the formation of such an anomalously extended plerion could be realized, in a diluted environment with the hydrogen number density $n_{ism}\le10^{-2}\rm\,cm^{-3}$. In this paper, we explore an alternative scenario assuming that the pulsar responsible for the formation of the nebula initially had a very short rotation period. In this case, the sizes of both the PWN and the surrounding supernova remnant depend on the initial pulsar period, the braking index, and the ISM density. To check the feasibility of this scenario, we study the parameter space that would reproduce the size of HESS J1825-137. We show that this demand can be achieved if the braking index is small, $n\leq2$ and the pulsar birth period is short, $P_{\rm b}\simeq1\rm\,ms$. This scenario can reproduce the wind termination position, which is expected at $R_{ts}\simeq0.03\rm\,pc$, only in a dense environment with $n_{ism}\geq\rm1\,cm^{-3}$. The requirement of the dense surrounding gas is supported by the presence of molecular clouds found in the source vicinity.

Acceleration of bulk motion of plasma due to magneto centrifugal mechanism is investigated for different shapes of the field lines of 3-D magnetic field. It is shown that this mechanism can be efficient provided that the field line is twisted into direction of rotation. According to the last results of the numerical modelling a fraction of the field lines with such geometrical shape apparently exist in the pulsar magnetosphere close to the last closed field line. In this case the Lorentz factor $\gamma$ increases with the radius $r$ along these field lines as $\gamma\sim (1-(r\Omega/c)^2)^{-1}$, where $\Omega$ - angular velocity of the pulsar rotation. The magneto centrifugal mechanism provides acceleration of the electrons to the energy which is sufficient to explain the observed VHE pulsed $\gamma$-rays from Crab pulsar.

We argue that the bright flare of the binary pulsar \object{PSR B1259$-$63/LS2883} detected by the {\it Fermi} Large Area Telescope (LAT), is due to the inverse Compton (IC) scattering of the unshocked electron-positron pulsar wind with a Lorentz factor $\Gamma_0 \approx 10^4$. The combination of two effects both linked to the circumstellar disk (CD), is a key element in the proposed model. The first effect is related to the impact of the surrounding medium on the termination of the pulsar wind. Inside the disk, the "early" termination of the wind results in suppression of its gamma-ray luminosity. When the pulsar escapes the disk, the conditions for termination of the wind undergo significant changes. This would lead to a dramatic increase of the pulsar wind zone, and thus to the proportional increase of the gamma-ray flux. On the other hand, if the parts of the CD disturbed by the pulsar can supply infrared photons of density high enough for efficient Comptonization of the wind, almost the entire kinetic energy of the pulsar wind would be converted to radiation, thus the gamma-ray luminosity of the wind could approach to the level of the pulsar's spin-down luminosity as reported by the {\it Fermi} collaboration.

We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation mechanisms. We show that the fragmented envelope of the star can be accelerated to Lorentz factors up to 100 and radiate effectively the available energy in gamma-rays predominantly through proton synchrotron radiation or external inverse Compton scattering of electrons. The model can readily explain the minute-scale TeV flares on top of longer (typical time-scales of days) gamma-ray variability as observed from the blazar PKS 2155-304. In the framework of the proposed scenario, the key parameters of the source are robustly constrained. In the case of proton synchrotron origin of the emission a mass of the central black hole of $M_{\rm BH}\approx 10^8 M_{\odot}$, a total jet power of $L_{\rm j} \approx 2\times 10^{47} \, \rm erg\,s^{-1}$ and a Doppler factor, of the gamma-ray emitting blobs, of $\delta\geq 40$ are required. Whilst for the external inverse Compton model, parameters of $M_{\rm BH}\approx 10^8 M_{\odot}$, $L_{\rm j} \approx 10^{46} \, \rm erg\,s^{-1}$ and the $\delta\geq 150$ are required.

In this paper, we present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with its environment in binary system. We compare the impact of both the magnetic field and the wind anisotropy to the benchmark case of a purely hydrodynamical (HD) interaction of isotropic winds, which has been studied in detail by Bogovalov et al. (2008). We consider the interaction in axisymmetric approximation, i.e. the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star and the effects related to the pulsar orbiting are neglected. The impact of the magnetic field is studied for the case of weak magnetization (with magnetization parameter $\sigma<0.1$), which is consistent with conventional models of pulsar winds. The effects related to anisotropy in pulsar winds are modeled assuming that the kinetic energy flux in a non-magnetized pulsar wind is strongly anisotropic, with the minimum at the pulsar rotation axis and the maximum in the perpendicular direction. We show that, although both considered effects change the shape of the region occupied by the terminated pulsar wind, their impact appears to be small. In particular, for the magnetization of the pulsar wind below 0.1, the magnetic field pressure remains well below the plasma pressure in the post-shock region. Thus, in the case of interaction of a pulsar with the stellar wind environment (opposite to the case of plerions, i.e. the pulsar interaction with the interstellar medium, when the magnetic field becomes dynamically important independently on the wind magnetization) the HD approach represents a feasible approximation for numerical modelling.

Binary pulsar systems emit potentially detectable components of gamma ray emission due to Comptonization of the optical radiation of the companion star by relativistic electrons of the pulsar wind, both before and after termination of the wind. The recent optical observations of binary pulsar system PSR B1259-63/LS 2883 revealed radiation properties of the companion star which differ significantly from previous measurements. In this paper we study the implications of these observations for the interaction rate of the unshocked pulsar wind with the stellar photons and the related consequences for fluxes of high energy (HE) and very high energy (VHE) gamma rays. We show that the signal should be strong enough to be detected with Fermi close to the periastron passage, unless the pulsar wind is strongly anisotropic or the Lorentz factor of the wind is smaller than $10^3$ or larger that $10^5$. The higher luminosity of the optical star also has two important implications: (i) attenuation of gamma rays due to photon-photon pair production, and (ii) Compton drag of the unshocked wind. While the first effect has an impact on the lightcurve of VHE gamma rays, the second effect may significantly decrease the energy available for particle acceleration after termination of the wind.

We present self-consistent global, steady-state MHD models and synthetic optically thin synchrotron emission maps for the jet of M87. The model consist of two distinct zones: an inner relativistic outflow, which we identify with the observed jet, and an outer cold disk-wind. While the former does not self-collimate efficiently due to its high effective inertia, the latter fulfills all the conditions for efficient collimation by the magneto-centrifugal mechanism. Given the right balance between the effective inertia of the inner flow and the collimation efficiency of the outer disk wind, the relativistic flow is magnetically confined into a well collimated beam and matches the measurements of the opening angle of M87 over several orders of magnitude in spatial extent. The synthetic synchrotron maps reproduce the morphological structure of the jet of M87, i.e. center-bright profiles near the core and limb-bright profiles away from the core. At the same time, they also show a local increase of brightness at some distance along the axis associated to a recollimation shock in the MHD model. Its location coincides with the position of the optical knot HST-1. In addition our best fitting model is consistent with a number of observational constraints such as the magnetic field in the knot HST-1, and the jet-to-counterjet brightness ratio.

The hydrodynamics of the interaction of pulsar and stellar winds in binary systems harboring a pulsar and its impact on the nonthermal radiation of the binary pulsar PSR B1259-63/SS2883 is discussed. The collision of an ultrarelativistic pulsar wind with a nonrelativistic stellar outflow results in significant bulk acceleration of the shocked material from the pulsar wind. Already at distances comparable to the size of the binary system, the Lorentz factor of the shocked flow can be as large as $\gamma$~4. This results in significant anisotropy of the inverse Compton radiation of accelerated electrons. Because of the Doppler boosting of the produced radiation, one should expect a variable gamma-ray signal from the system. In particular, this effect may naturally explain the reported tendency of a decrease of TeV gamma-ray flux close to the periastron. The modeling of the interaction of pulsar and stellar winds allows self-consistent calculations of adiabatic losses. Our results show that adiabatic losses dominate over the radiative losses. These results have direct impact on the orbital variability of radio, X-ray and gamma-ray signals detected from the binary pulsar PSR 1259-63/SS2883.

We consider an extreme case of disc accretion onto a gravitating centre when the viscosity in the disc is negligible. The angular momentum and the rotational energy of the accreted matter is carried out by a magnetized wind outflowing from the disc. The outflow of matter from the disc occurs due to the Blandford & Payne(1982) centrifugal mechanism. The disc is assumed to be cold. Accretion and outflow are connected by the conservation of the energy, mass and the angular momentum. The basic properties of the outflow, angular momentum flux and energy flux per particle in the wind, do not depend on the details of the structure of the accretion disc. In the case of selfsimilar accretion/outflow, the dependence of the rate of accretion $\dot M$ in the disc depends on the disc radius $r$ on the law $\dot M \sim r^{{1\over2(\alpha^2-1)}}$, where $\alpha$ is a dimensionless Alfvenic radius. In the case of $\alpha \gg 1$, the accretion in the disc is provided by very weak matter outflow from the disc and the outflow predominantly occurs from the very central part of the disc. The solution obtained in the work provides mechanism which transforms the gravitational energy of the accreted matter into the energy of the outflowing wind with efficiency close to 100%. The final velocity can essentially exceed Kepler velocity at the site of the wind launch. This mechanism allows us to understand the nature of the astrophysical objects with low luminosity discs and energetic jet-like outflows.

We present a self-consistent MHD model for the jet of M87. The model consist of two distinct zones: an inner relativistic outflow, which we identify with the observed jet, and an outer cold disk-wind. While the former does not self-collimate efficiently due to its high effective inertia, the latter fulfills all the conditions for efficient collimation by the magneto-centrifugal mechanism. Given the right balance between the effective inertia of the inner flow and the collimation efficiency of the outer disk wind, the relativistic flow is magnetically confined into a well collimated beam for a wide range of parameters and matches the measurements of the opening angle of M87 over several orders of magnitude in spatial extent. In the second part of this work, we present synthetic synchrotron emission maps for our MHD models. In principle the two-zone model can reproduce the morphological structure seen in radio observations, as central-peaked profiles across the jet close the the source, limb-bright further down the jet, and a bright knot close to the position of HST-1. However it is difficult to reconcile all features into a single set of parameters.

We consider disk accretion resulting purely from the loss of angular momentum due to the outflow of plasma from a magnetized disk. In this limiting case, the dissipation due to the viscosity and finite electrical conductivity of the plasma can be neglected. We have obtained self-consistent, self-similar solutions for dissipationless disk accretion. Such accretion may result in the formation of objects whose bolometric luminosities are lower than the flux of kinetic energy in the ejected material.

In this paper, we present a detailed hydrodynamical study of the properties of the flow produced by the collision of a pulsar wind with the surrounding in a binary system. This work is the first attempt to simulate interaction of the ultrarelativistic flow (pulsar wind) with the nonrelativistic stellar wind. Obtained results show that the wind collision could result in the formation of an "unclosed" (at spatial scales comparable to the binary system size) pulsar wind termination shock even when the stellar wind ram pressure exceeds significantly the pulsar wind kinetical pressure. Moreover, the post-shock flow propagates in a rather narrow region, with very high bulk Lorentz factor ($\gamma\sim100$). This flow acceleration is related to adiabatical losses, which are purely hydrodynamical effects. Interestingly, in this particular case, no magnetic field is required for formation of the ultrarelativistic bulk outflow. The obtained results provide a new interpretation for the orbital variability of radio, X-ray and gamma-ray signals detected from binary pulsar system PSR 1259-63/SS2883.

The inverse Compton scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. While this process can explain the energy spectrum of the observed TeV emission, the detected gamma-ray fluxes do not agree with the published theoretical predictions of the TeV lightcurve. The main objective of this paper is to show that the HESS results can be explained, under certain reasonable assumptions, by inverse Compton scenarios of gamma-ray production in the system. In this paper we study evolution of the energy spectra of relativistic electrons under different assumptions about the acceleration and energy-loss rates of electrons, and the impact of these processes on the lightcurve of IC gamma-rays. We demonstrate that the observed TeV lightcurve can be explained (i) by adiabatic losses which dominate over the entire trajectory of the pulsar with a significant increase towards the periastron, or (ii) by the sub-TeV cutoffs in the energy spectra of electrons due to the enhanced rate of Compton losses close to the periastron. The Compton deceleration of the pulsar wind contributes to the decrease of the nonthermal power released in the accelerated electrons, and thus to the reduction of the IC and synchrotron components of radiation close to the periastron. Although this effect alone cannot explain the observed lightcurves, the Comptonization of the pulsar wind leads to the formation of gamma-radiation with a line-type energy spectrum. While the HESS data already constrain the Lorentz factor of the wind, $\Gamma \le 10^6$, future observations of this object with GLAST should allow a deep probe of the wind Lorentz factor in the range between $10^4$ and $10^6$.

We apply a two-zone MHD model to the jet of M87. The model consists of an inner relativistic outflow, which is surrounded by a non-relativistic outer disk-wind. The outer disk-wind collimates very well through magnetic self-collimation and confines the inner relativistic jet into a narrow region around the rotation axis. Further, we show by example, that such models reproduce very accurately the observed opening angle of the M87 jet over a large range from the kiloparsec scale down to the sub-parsec scale.

The formation of the toroidal and jet-like structures in the central part of the Crab Nebula is explained in the framework of Kennel & Coroniti theory. The only new element introduced by us in this theory is the initial anisotropy of the energy flux in the wind. We estimate the X-ray surface brightness of the Crab Nebula from the region of interaction of this wind with the interstellar medium and compare it with observations.

We discuss the potential of future sub-10 GeV threshold imaging atmospheric Cherenkov telescope arrays for exploring the physics of rotation powered pulsars and their interactions with the ambient medium through relativistic winds and termination shocks.

We interpret the observed X-ray morphology of the central part of the Crab Nebula (torus + jets) in terms of the standard theory by Kennel and Coroniti (1984). The only new element is the inclusion of anisotropy in the energy flux from the pulsar in the theory. In the standard theory of relativistic winds, the Lorentz factor of the particles in front of the shock that terminates the pulsar relativistic wind depends on the polar angle as $\gamma=\gamma_0+\gamma_m\sin^2\theta$, where $\gamma_0 \sim 200$ and $\gamma_m \sim 4.5\times 10^6$. The plasma flow in the wind is isotropic. After the passage of the pulsar wind through the shock, the flow becomes subsonic with a roughly constant (over the plerion volume) pressure $P={1\over 3}n\epsilon$, where $n$ is the plasma particle density and $\epsilon$ is the mean particle energy. Since $\epsilon \sim \gamma mc^2$, a low-density region filled with the most energetic electrons is formed near the equator. A bright torus of synchrotron radiation develops here. Jet-like regions are formed along the pulsar rotation axis, where the particle density is almost four orders of magnitude higher than that in the equatorial plane, because the particle energy there is four orders of magnitude lower. The energy of these particles is too low to produce detectable synchrotron radiation. However, these quasi-jets become comparable in brightness to the torus if additional particle acceleration takes place in the plerion. We also present the results of our study of the hydrodynamic interaction between an anisotropic wind and the interstellar medium. We compare the calculated and observed distributions of the volume intensity of X-ray radiation.

A stationary self-consistent outflow of a magnetised relativistic plasma from a rotating object with an initially monopole-like magnetic field is investigated in the ideal MHD approximation under the condition $\sigma/U_0^2 >1$, where $\sigma$ is the ratio of the Poynting flux over the mass energy flux at the equator and the surface of the star, with $U_0=\gamma_0v_0/c$ and $\gamma_0$ the initial four-velocity and Lorentz factor of the plasma. The mechanism of the magnetocentrifugal acceleration and self-collimation of the relativistic plasma is investigated. A jet-like relativistic flow along the axis of rotation is found in the steady-state solution under the condition $\sigma/U_0^2 >1$ with properties predicted analytically. The amount of the collimated matter in the jet is rather small in comparison to the total mass flux in the wind. An explanation for the weak self-collimation of relativistic winds is given.

Recent numerical simulations and analytical models of magnetically collimated plasma outflows from a uniformly rotating central gravitating object and/or a Keplerian accretion disk have shown that relatively low mass and magnetic fluxes reside in the produced jet. Observations however indicate that in some cases, as in jets of YSO's, the collimated outflow carries higher fluxes than these simulations predict. A solution to this problem is proposed here by assuming that jets with high mass flux originate in a central source which produces a noncollimated outflow provided that this source is surrounded by a rapidly rotating accretion disk. The relatively faster rotating disk produces a collimated wind which then forces all the enclosed outflow from the central source to be collimated too. This conclusion is confirmed by self-consistent numerical solutions of the full set of the MHD equations.

Violation of the axial synnetry of a magnetic field essentially modifies the physics of the plasma outflow in the magnetosphere of rotating objects. In comparison to the axisymmetric outflow, two new affects appear: more efficient magnetocentrifugal acceleration of the plasma along particular field lines and generation of MHD waves. Here, we use an ideal MHD approximation to study the dynamics of a cold wind in the nonaxisymmetric magnetosphere. We obtain a self-consistent analytical solution of the problem of cold plasma outflow from a slowly rotating star with a slightly nonaxisymmetric magnetic field using perturbation theory. In the axisymmetric (monopole-like) magnetic field, the first term in the expansion of the terminating energy of the plasama in powers of $\Omega$ is proportional to $\Omega^4$, where $\Omega$ is the angular velocity of the central source. Violation of the axial symmetry of the magnetic field crucially changes this dependence. The first correction to the energy of the plasma becomes proportional to $\Omega$. Efficient magnetocentrifugal acceleration occurs along the field lines curved initially in the direction of the rotation. I argue that all necessary conditions for the efficient magnetocentrifugal acceleration of the plasma exist in the radio pulsar megnetosphere. We calculated the first correction of the rotational losses due to the generation of the MHD waves and analysed the plasma acceleration by these waves.

We resolve the paradox that although magnetic collimation of an isotropic solar wind results in an enhancement of its proton flux along the polar directions, several observations indicate a wind proton flux peaked at the equator. To that goal, we solve the full set of the time-dependent MHD equations describing the axisymmetric outflow of plasma from the magnetized and rotating Sun, either in its present form of the solar wind, or, in its earlier form of a protosolar wind. Special attention is directed towards the collimation properties of the solar outflow at large heliocentric distances. For the present day solar wind it is found that the poloidal streamlines and fieldlines are only slightly focused toward the solar poles. However, even such a modest compression of the flow by the azimuthal magnetic field would lead to an increase of the mass flux at the polar axis by about 20% at 1 AU, relatively to its value at the equator, for an initially isotropic at the base wind, contrary to older and recent (Prognoz, Ulysses, SOHO) observations. For the anisotropic in heliolatitude wind with parameters at the base inferred from in situ observations by ULYSSES/SWOOPS and SOHO/CDS the effect of collimation is almost totally compensated by the initial velocity and density anisotropy of the wind. This effect should be taken into account in the interpretation of the recent SOHO observations by the SWAN instrument. Similar simulations have been performed for a five- and ten-fold increase of the solar angular velocity corresponding presumably to the wind of an earlier phase of our Sun. For such conditions it is found that for initially radial streamlines, the azimuthal magnetic field created by the fast rotation focus them toward the rotation axis and forms a tightly collimated jet.

We show that the relativistic wind in the Crab pulsar, which is commonly thought to be invisible in the region upstream of the termination shock at R < 0.1 pc, in fact could be directly observed through its inverse Compton gamm-ray emission. The search for such specific component of radiation in the gamma-ray spectrum of the Crab can provide unique information about the unshocked pulsar wind that is not accessible at other wavelengths.

I show that the self-consistent solution of the problem of MHD plasma flow in magnetosphere of an oblique rotator with an initially split-monopole magnetic field is reduced to the solution of the similar problem for the axisymmetric rotator. All properties of the MHD cold plasma flows from the axisymmetric rotators with the initially split-monopole magnetic field are valid for the oblique rotators as well. Rotational losses of the oblique rotator do not depend on the inclination angle and there is no temporal evolution of this angle. Self-consistent analytical and numerical solutions for the axisymmetric plasma flows obtained earlier show that the rotators can be divided on fast rotators ($\sigma_0/U_0^2 > 1$) and slow rotators ($\sigma_0/U_0^2 < 1$), where $\sigma_0$ is the ratio of the Poynting flux to the matter energy flux in the flow at the equator on the surface of the star, $U_0 = \gamma_0v_0/c$, $v_0$ and $\gamma_0$ are the initial velocity and Lorentz-factor of the plasma. The self-consistent approximate analytical solution for the plasma flow from the oblique rotator is obtained under the condition $\sigma_0/U_0^2 \ll 1$. Implications of these results for radio pulsars are discussed. In particular, I argue that all radio pulsars are apparently the slow rotators ejecting the Poynting dominated relativistic wind.

The axisymmetric 3-D MHD outflow of a cold plasma from a magnetized and rotating astrophysical object is numerically simulated with the purpose of investigating the outflow's magnetocentrifugal acceleration and eventual collimation. Gravity and thermal pressure are neglected while a split-monopole is used to describe the initial magnetic field configuration. It is found that the stationary final state depends critically on a single parameter alpha expressing the ratio of the corotating speed at the Alfven distance to the initial flow speed along the initial monopole-like magnetic fieldlines. Several angular velocity laws have been used for relativistic and nonrelativistic outflows. The acceleration of the flow is most effective at the equatorial plane and the terminal flow speed depends linearly on alpha. Significant flow collimation is found in nonrelativistic efficient magnetic rotators corresponding to relatively larger than 1 values of alpha while very weak collimation occurs in inefficient magnetic rotators with values of alpha smaller than about 1. Part of the flow around the rotation and magnetic axis is cylindrically collimated while the remaining part obtains radial asymptotics. The transverse radius of the jet is inversely proportional to alpha while the density in the jet grows linearly with alpha. For alpha greater than about 5 the magnitude of the flow in the jet remains below the fast MHD wave speed everywhere. In relativistic outflows, no collimation is found in the supersonic region for parameters typical for radio pulsars. All above results verify the main conclusions of general theoretical studies on the magnetic acceleration and collimation of outflows from magnetic rotators and extend previous numerical simulations to large stellar distances.