The centrifugal acceleration is due to the rotating poloidal magnetic field in the magnetosphere creates the electric field which is orthogonal to the magnetic field. Charged particles with finite cyclotron radii can move along the electric field and receive energy. Centrifugal acceleration pushes particles to the periphery, where their azimuthal velocity reaches the light speed. We have calculated particle trajectories by numerical and analytical methods. The maximum obtained energies depend on the parameter of the particle magnetization $ \kappa $, which is the ratio of rotation frequency of magnetic field lines in the magnetosphere $ \Omega_F $ to non-relativistic cyclotron frequency of particles $ \omega_c $, $ \kappa = \Omega_F /\omega_c << 1 $, and from the parameter $ \alpha $ which is the ratio of toroidal magnetic field $ B_T $ to the poloidal one $ B_P $, $ \alpha = B_T / B_P $. It is shown that for small toroidal fields, $ \alpha <\kappa^{1/4} $, the maximum Lorentz factor $ \gamma_m $ is only the square root of magnetization, $ \gamma_m = \kappa^{-1/2} $, while for large toroidal fields, $ \alpha >\kappa^{1/4} $, the energy increases significantly, $ \gamma_m = \kappa^{-2/3} $. However, the maximum possible acceleration, $ \gamma_m = \kappa^{-1} $, is not achieved in the magnetosphere. For a number of active galactic nuclei, such as M87, maximum values of Lorentz factor for accelerated protons are found. Also for special case of Sgr. A* estimations of the maximum proton energy and its energy flux are obtained. They are in agreement with experimental data obtained by HESS Cherenkov telescope.

Thermal effects on the properties of four electromagnetic waves propagating in the pulsar magnetosphere are analyzed. It is shown that thermal effects change only quantitatively the dispersion properties of superluminal ordinary O-mode freely escaping the pulsar magnetosphere; properties of the extraordinary X-mode remain unchanged. As for two subluminal waves propagating along magnetic field lines, for them thermal effects result in essential absorption. However, this attenuation occurs at considerable distances from the neutron star, so there is no doubt in their existence.

Acceleration of protons in the active galactic nuclei is considered. The largest energy is achieved by protons during centrifugal acceleration in the magnetosphere of the central machine. When the proton accelerated in the magnetosphere of a black hole approaches light cylinder surface, acceleration occurs mainly in the azimuthal direction, i.e. the acceleration is centrifugal. In this paper the acceleration of a proton having smaller synchrotron losses compared to the electron is considered. As a proton experiences the highest energy increase while accelerating near the light surface, a partial solution for the maximum Lorentz factor can be obtained there. In the analysis the obtained dependence of the maximum energy on the parameter of particle magnetization $ \kappa $ and parameter $ \alpha $ which reflects the relation of toroidal $ B_\phi $ and poloidal $ B_T $ magnetic fields , has led to the conclusion that the achievement of theoretical maximum limit of Lorentz factor value $ \gamma_m=\kappa^{-1}$ is not possible for an accelerated particle in the magnetosphere of a black hole due to restrictions of the topology of toroidal and poloidal magnetic fields imposed. The analysis of special cases of the relation of toroidal and poloidal magnetic field has shown that in the presence of magnetic field that is significantly more toroidal the maximum Lorentz factor value reaches $\gamma_m = \kappa^ {-2/3} $, in case when toroidal field becomes smaller in comparison to poloidal field the maximum Lorentz factor value does not exceed $\gamma_m = \kappa^ {-1/2} $. For a number of objects, such as M87 and Sgr. A *, maximum Lorentz factor values for accelerated protons for scenarios of existence or lack of toroidal magnetic field have been derived. The obtained results for magnetosphere of Sgr. A * has confirmed by the experimental data obtained on the massive HESS of Cherenkov telescopes.

We study diffusion of charged particles in stationary stochastic magnetic field ${\bf B}$ with zero mean, $\langle {\bf B} \rangle = 0 $. In the case when electric current is carried by electrons, the field is force-free, $\mathrm{curl} \,{\bf B} = \alpha{\bf B} $, where $\alpha({\bf r})$ is an arbitrary scalar function. In a small region where the function $\alpha $ and the field magnitude $|{\bf B}|$ are approximately constant, the equations of motion of charged particles are integrated and reduced to the equation of mathematical pendulum. The transition from trapped to untrapped particles is continuously traced. Averaging over the magnetic field spectrum gives the spatial diffusion coefficient $D$ of particles as a function of the Larmor radius $r_L$ in the large-scale magnetic fields ($B_{LS}$) and magnetic field correlation length $L_0$. The diffusion coefficient turns out to be proportional to the Larmor radius, $D\propto r_L $, for $r_L <L_0 / 2\pi $, and to the Larmor radius squared, $ D \propto r_L^2 $, for $ r_L> L_0 /2\pi $. We apply obtained results to the diffusion of cosmic rays in the Galaxy, which contains a large number of independent regions with parameters $L_0$ and $B_{LS}$ varying in wide range. We average over $B_{LS}$ with the Kolmogorov spectrum and over $L_0$ with the distribution function $f(L_0)\propto L_0^{- 1+ \sigma}$. For the practically flat spectrum $\sigma = 1/15$, we have $ D\propto r_m^{0.7}$, which is consistent with observations.

The flow of a matter, accreting onto a magnetized neutron star, is accompanied by an electric current. The closing of the electric current occurs in the crust of a neutron stars in the polar region across the magnetic field. But the conductivity of the crust along the magnetic field greatly exceeds the conductivity across the field, so the current penetrates deep into the crust down up to the super conducting core. The magnetic field, generated by the accretion current, increases greatly with the depth of penetration due to the Hall conductivity of the crust is also much larger than the transverse conductivity. As a result, the current begins to flow mainly in the toroidal direction, creating a strong longitudinal magnetic field, far exceeding an initial dipole field. This field exists only in the narrow polar tube of $r$ width, narrowing with the depth, i.e. with increasing of the crust density $\rho$, $r\propto \rho^{-1/4}$. Accordingly, the magnetic field $B$ in the tube increases with the depth, $B\propto \rho^{1/2}$, and reaches the value of about $10^{17}$ Gauss in the core. It destroys super conducting vortices in the core of a star in the narrow region of the size of the order of ten centimeters. Because of generated density gradient of vortices they constantly flow into this dead zone and the number of vortices decreases, the magnetic field of a star decreases as well. The attenuation of the magnetic field is exponential, $B=B_0(1+t/\tau)^{-1}$. The characteristic time of decreasing of the magnetic field $\tau$ is equal to $\tau\simeq 10^3$ years. Thus, the magnetic field of accreted neutron stars decreases to values of $10^8 - 10^9$ Gauss during $10^7-10^6$ years.

From the analysis of the flux of high energy particles, $E>3\cdot 10^{18}eV$, it is shown that the distribution of the power density of extragalactic rays over energy is of the power law, ${\bar q}(E)\propto E^{-2.7}$, with the same index of $2.7$ that has the distribution of Galactic cosmic rays before so called 'knee', $E<3\cdot 10^{15}eV$. However, the average power of extragalactic sources, which is of ${\cal E}\simeq 10^{43}erg \,s^{-1}$, at least two orders exceeds the power emitted by the Galaxy in cosmic rays, assuming that the density of galaxies is estimated as $N_g\simeq 1 Mpc^{-3}$. Considering that such power can be provided by relativistic jets from active galactic nuclei with the power ${\cal E}\simeq 10^{45} - 10^{46} erg \, s^{-1}$, we estimate the density of extragalactic sources of cosmic rays as $N_g\simeq 10^{-2}-10^{-3}\, Mpc^{-3}$. Assuming the same nature of Galactic and extragalactic rays, we conclude that the Galactic rays were produced by a relativistic jet emitted from the Galactic center during the period of its activity in the past. The remnants of a bipolar jet are now observed in the form of bubbles of relativistic gas above and below the Galactic plane. The break, observed in the spectrum of Galactic rays ('knee'), is explained by fast escape of energetic particle, $E>3\cdot 10^{15}eV$, from the Galaxy because of the dependence of the coefficient of diffusion of cosmic rays on energy, $D\propto E^{0.7}$. The obtained index of the density distribution of particles over energy, $N(E)\propto E^{-2.7-0.7/2}=E^{-3.05}$, for $E>3\cdot 10^{15}eV$ agrees well with the observed one, $N(E)\propto E^{-3.1}$. Estimated time of termination of the jet in the Galaxy is $4.2\cdot 10^{4}$ years ago.

We investigate the magnetic field which is generated by turbulent motions of a weakly ionized gas. Galactic molecular clouds give us an example of such a medium. As in the Kazantsev-Kraichnan model we assume a medium to be homogeneous and a neutral gas velocity field to be isotropic and delta-correlated in time. We take into consideration the presence of a mean magnetic field, which defines a preferred direction in space and eliminates isotropy of magnetic field correlators. Evolution equations for the anisotropic correlation function are derived. Isotropic cases with zero mean magnetic field as well as with small mean magnetic field are investigated. It is shown that stationary bounded solutions exist only in the presence of the mean magnetic field for the Kolmogorov neutral gas turbulence. The dependence of the magnetic field fluctuations amplitude on the mean field is calculated. The stationary anisotropic solution for the magnetic turbulence is also obtained for large values of the mean magnetic field.

It was as early as the 1980s that A V Gurevich and his group proposed a theory to explain the magnetosphere of radio pulsars and the mechanism by which they produce coherent radio emission. The theory has been sharply criticized and is currently rarely mentioned when discussing the observational properties of radio pulsars, even though all the criticisms were in their time disproven in a most thorough and detailed manner. Recent results show even more conclusively that the theory has no internal inconsistencies. New observational data also demonstrate the validity of the basic conclusions of the theory. Based on the latest results on the effects of wave propagation in the magnetosphere of a neuron star, we show that the developed theory does indeed allow quantitative predictions of the evolution of neutron stars and the properties of the observed radio emission.

The problem of interaction of the rotating magnetic field, frozen to a star, with a thin well conducting accretion disk is solved exactly. It is shown that a disk pushes the magnetic field lines towards a star, compressing the stellar dipole magnetic field. At the point of corotation, where the Keplerian rotation frequency coincides with the frequency of the stellar rotation, the loop of the electric current appears. The electric currents flow in the magnetosphere only along two particular magnetic surfaces, which connect the corotation region and the inner edge of a disk with the stellar surface. It is shown that the closed current surface encloses the magnetosphere. Rotation of a disk is stopped at some distance from the stellar surface, which is 0.55 of the corotation radius. Accretion from a disk spins up the stellar rotation. The angular momentum transferred to the star is determined.

A conducting disk significantly changes the generation of the electromagnetic radiation excited by the rotation of the magnetic field frozen to a star. Due to the reflection of waves from a disk there appear waves propagating toward a star, not only outward a star as it takes place for the magneto-dipole radiation. Because that the angular momentum can be transformed from a disk to a star when the inner edge of a disk approaches the light surface of a rotating star. This is purely electromagnetic effect. At some distance of a disk from a star, $r_d=r^*\simeq c/\omega_s$, the stellar angular momentum losses due to the electromagnetic radiation become zero. It results the stable stellar rotation.

The paper deals with the one possible mechanism of the pulsar radio emission, i.e., with the collective curvature radiation of the relativistic particle stream moving along the curved magnetospheric magnetic field lines. It is shown that the electromagnetic wave containing one cylindrical harmonic exp{is{\phi}} can not be radiated by the curvature radiation mechanism, that corresponds to radiation of a charged particle moving along curved magnetic field lines. The point is that the particle in vacuum radiates the triplex of harmonics (s, s \pm 1), so for the collective curvature radiation the wave polarization is very important and cannot be fixed a priori. For this reason the polarization of real unstable waves must be determined directly from the solution of wave equations for the media. Its electromagnetic properties should be described by the dielectric permittivity tensor \^{\epsilon}({\omega},k,r), that contains the information on the reaction on all possible types of radiation.

The production of electron-positron pairs in a vacuum neutron star magnetosphere is investigated for both low (compared to the Schwinger one) and high magnetic fields. The case of a strong longitudinal electric field where the produced electrons and positrons acquire a stationary Lorentz factor in a short time is considered. The source of electron-positron pairs has been calculated with allowance made for the pair production by curvature and synchrotron photons. Synchrotron photons are shown to make a major contribution to the total pair production rate in a weak magnetic field. At the same time, the contribution from bremsstrahlung photons may be neglected. The existence of a time delay due to the finiteness of the electron and positron acceleration time leads to a great reduction in the electron-positron plasma generation rate compared to the case of a zero time delay. The effective local source of electron-positron pairs has been constructed. It can be used in the hydrodynamic equations that describe the development of a cascade after the absorption of a photon from the cosmic gamma-ray background in a neutron star magnetosphere.

The connection between the radio emission from "lightnings" produced by the absorption of high-energy photons from the cosmic gamma-ray background in a neutron star magnetosphere and radio bursts from rotating radio transients (RRATs) is investigated. The lightning length reaches 1000 km; the lightning radius is 100 m and is comparable to the polar cap radius. If a closed magnetosphere is filled with a dense plasma, then lightnings are efficiently formed only in the region of open magnetic field lines. For the radio emission from a separate lightning to be observed, the polar cap of the neutron star must be directed toward the observer and, at the same time, the lightning must be formed. The maximum burst rate is related to the time of the plasma outflow from the polar cap region. The typical interval between two consecutive bursts is ~100 s. The width of a single radio burst can be determined both by the width of the emission cone formed by the lightning emitting regions at some height above the neutron star surface and by a finite lightning lifetime. The width of the phase distribution for radio bursts from RRATs, along with the integrated pulse width, is determined by the width of the bundle of open magnetic field lines at the formation height of the radio emission. The results obtained are consistent with the currently available data and are indicative of a close connection between RRATs, intermittent pulsars, and extreme nullers.

The absorption of a high-energy photon from the external cosmic gamma-ray background in the inner neutron star magnetosphere triggers the generation of a secondary electron-positron plasma and gives rise to a lightning - a lengthening and simultaneously expanding plasma tube. It propagates along magnetic fields lines with a velocity close to the speed of light. The high electron-positron plasma generation rate leads to dynamical screening of the longitudinal electric field that is provided not by charge separation but by electric current growth in the lightning. The lightning radius is comparable to the polar cap radius of a radio pulsar. The number of electron-positron pairs produced in the lightning in its lifetime reaches 10^28. The density of the forming plasma is comparable to or even higher than that in the polar cap regions of ordinary pulsars. This suggests that the radio emission from individual lightnings can be observed. Since the formation time of the radio emission is limited by the lightning lifetime, the possible single short radio bursts may be associated with rotating radio transients (RRATs).

It is shown that the relativistic jet, emitted from the center of the Galaxy during its activity, possessed power and energy spectrum of accelerated protons sufficient to explain the current cosmic rays distribution in the Galaxy. Proton acceleration takes place on the light cylinder surface formed by the rotation of a massive black hole carring into rotation the radial magnetic field and the magnetosphere. Observed in gamma, x-ray and radio bands bubbles above and below the galactic plane can be remnants of this bipolar get. The size of the bubble defines the time of the jet's start, $\simeq 2.4\cdot 10^7$ years ago. The jet worked more than $10^7$ years, but less than $2.4\cdot10^7$ years.

We consider the electron-positron plasma generation processes in the magnetospheres of magnetars - neutron stars with strong surface magnetic fields, B = 10^(14) - 10^(15) G. We show that the photon splitting in a magnetic field, which is effective at large field strengths, does not lead to the suppression of plasma multiplication, but manifests itself in a high polarization of gamma-ray photons. A high magnetic field strength does not give rise to the second generation of particles produced by synchrotron photons. However, the density of the first-generation particles produced by curvature photons in the magnetospheres of magnetars can exceed the density of the same particles in the magnetospheres of ordinary radio pulsars. The plasma generation inefficiency can be attributed only to slow magnetar rotation, which causes the energy range of the produced particles to narrow. We have found a boundary in the P - Pdot diagram that defines the plasma generation threshold in a magnetar magnetosphere.

The motion of electrons and positrons in the vacuum magnetosphere of a neutron star with a surface magnetic field of B~10^12 G is considered. Particles created in the magnetosphere or falling into it from outside are virtually instantaneously accelerated to Lorentz factors gamma~10^8. After crossing the force-free surface, where the projection of the electric field onto the magnetic field vanishes, a particle begins to undergo ultra-relativistic oscillations. The particle experiences a regular drift along the force-free surface simultaneous with this oscillatory motion.

The accumulation of electrons and positrons in the vacuum magnetosphere of a neutron star with a surface magnetic field of B~10^12 G is considered. It is shown that particles created in the magnetosphere or falling into the magnetosphere from outside undergo ultra-relativistic oscillations with a frequency of 10-100 MHz. These oscillations decay due to energy losses to curvature radiation and bremsstrahlung, with their frequencies reaching 1-10 GHz. Simultaneously, the particles undergo regular motion along the force-free surface along closed trajectories. This leads to the gradual accumulation of particles at the force-free surface and the formation of a fully charge-separated plasma layer with a density of the order of the Goldreich-Julian density. The presence of a constant source of electron-positron pairs in the magnetosphere due to the absorption of energetic cosmic gamma-rays leads to the growth of this layer, bringing about a rapid filling of the pulsar magnetosphere with electron-positron plasma if the pair-creation multiplication coefficient is sufficiently high.

We consider the motion of charged particles in the vacuum magnetospheres of rotating neutron stars with a strong surface magnetic field, B>10^(12) G. The electrons and positrons falling into the magnetosphere or produced in it are shown to be captured by the force-free surface EB=0. Using the Dirac-Lorentz equation, we investigate the dynamics of particle capture and subsequent motion near the force-free surface. The particle energy far from the force-free surface has been found to be determined by the balance between the power of the forces of an accelerating electric field and the intensity of curvature radiation. When captured, the particles perform adiabatic oscillations along the magnetic field lines and simultaneously move along the force-free surface. We have found the oscillation parameters and trajectories of the captured particles. We have calculated the characteristic capture times and energy losses of the particles through the emission of both bremsstrahlung and curvature photons by them. The capture of particles is shown to lead to a monotonic increase in the thickness of the layer of charged plasma accumulating near the force-free surface. The time it takes for a vacuum magnetosphere to be filled with plasma has been estimated.

We show numerically that vector antenna arrays can generate radio beams which exhibit spin and orbital angular momentum characteristics similar to those of helical Laguerre-Gauss laser beams in paraxial optics. For low frequencies (< 1 GHz), digital techniques can be used to coherently measure the instantaneous, local field vectors and to manipulate them in software. This opens up for new types of experiments that go beyond those currently possible to perform in optics, for information-rich radio physics applications such as radio astronomy, and for novel wireless communication concepts.

Interaction of a fast shock wave generated during a supernova explosion with a magnetized star-companion of the supernova precursor produces a current sheet. We consider the evolution of this current sheet and show that a singularity (shock) is formed in finite time within the ideal MHD framework. Charged particles (electrons) are accelerated in the vicinity of the singularity, and their distribution function has a plateau up to the energies of the order of $10^4 mc^2$. These fast particles radiate in the $\gamma$-range in the strong magnetic field of the current sheet ($B\simeq 10^6 G$). Radiation is concentrated within a narrow angle around the current sheet, $\Delta\theta \simeq 3\cdot 10^{-4}$, and its spectrum has the maximum at several hundreds of $keV$. Presented calculations confirm the model of cosmological GRBs proposed by Istomin & Komberg (2002).

It is shown that the slow glitches in the spin rate of the pulsar B1822-09 can be explained by the reconstruction of the neutron star shape, which is not matched with the star rotation axis. Owing to the evolution of the inclination angle, i.e. the angle between the rotation axis and the axis of the magnetic dipole, under the action of the braking torque, there appears the disagreement between the rotation axis and the symmetry axis. After the angle between the axis of symmetry and the axis of the rotation achieves the maximum value of alpha ~ 2x10^-4 the shape of the neutron star becomes matching with the rotation axis. Such reconstruction is observed as the slow glitch.

We consider the propagation and stability of the helical electromagnetic perturbations in a relativistic force-free Pointing flux dominated jet (magnetic helix). Perturbations are found to dissipate energy in the vicinity of the resonant surfaces inside the jet, where the resonance of the perturbation with the Alfven waves occurs. The energy of the perturbations is converted into the energy of accelerated particles, which emit synchrotron radiation. The intensity and polarization of this radiation is modulated with the perturbation. Perturbations with vanishing group velocity are the most important and could be associated with the bright knots in parsec scale jets.

We calculate the optically thin synchrotron emission of fast electrons and positrons in a spiral stationary magnetic field and a radial electric field of a rotating relativistic strongly magnetized force-free jet consisting of electron-positron pair plasma. The magnetic field has a helical structure with a uniform axial component and a toroidal component that is maximal inside the jet and decreasing to zero towards the boundary of the jet. Doppler boosting and swing of the polarization angle of synchrotron emission due to the relativistic motion of the emitting volume are calculated. The distribution of the plasma velocity in the jet is consistent with the electromagnetic field structure. Two spatial distributions of fast particles are considered: uniform, and concentrated in the vicinity of the Alfven resonance surface. The latter distribution corresponds to the regular acceleration by an electromagnetic wave in the vicinity of its Alfven resonance surface inside the jet. The polarization properties of the radiation have been obtained and compared with the existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give a natural explanation of the observed bimodality in the alignment between the electric field vector of the polarized radiation and the projection of the jet axis on the plane of the sky. We interpret the motion of bright knots as a phase velocity of standing spiral eigenmodes of electromagnetic perturbations in a cylindrical jet. The degree of polarization and the velocity of the observed proper motion of bright knots depend upon the angular rotational velocity of the jet. The observed polarizations and velocities of knots indicate that the magnetic field lines are bent in the direction opposite to the direction of the jet rotation.

We develop a theory of particle acceleration inside relativistic rotating electron-positron force-free jets with spiral magnetic fields. We considered perturbation of the stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled, i.e. where the local Alfven speed is equal to the phase speed of an eigenmode. The acceleration mechanism is regular acceleration by the electric field of the electromagnetic wave, which is the eigenmode of the force-free cylindrical jet: particles drift out of the region of the large wave amplitude near the Alfven resonant surface and gain energy. Acceleration in the strong electric field near the Alfven resonance and synchrotron losses combine to form a power-law energy spectrum of ultra-relativistic electrons and positrons with index between 2 and 3 depending upon the initial energy of the injected particles. The power law distribution ranges from 10 MeV to 1000 MeV.

The SN explosion in the closed binary can give the magnetospheric flare possessing the properties of GRB. The SN shock, flowing around the magnetosphere of a magnetized neutron star or a white dwarf, produces a narrow magnetic tail 10^9 cm long, 10^8 cm wide and a magnetic field of 10^6 Gauss. Fast particles (Lorentz factor of 10^4), generated in the tail by reconnection processes, radaite gamma rays of the 100 KeV - 1 MeV energies. The duration of radiation T<1 sec corresponds to a short GRB. Apart, the powerful shock can tear and accelerate part of the tail. That is the relativistic, strongly magnetized jet, producing gamma radiation and also X-rays and optic afterglow. That is long (T> 10 sec) GRB. The duration of the afterglow is inversly proportional to the photon energy and is several months for optic.

It is suggested that the accretion disk excites the eigen modes of Alfven oscillations of the magnetic field tubes the ends of which are frozen to the neutron star surface. The resonance takes place when the eigen Alfven frequency coincides with the difference of the Keplerian disk rotation frequency and the frequency of rotation of a neutron star. This model explains the two narrow lines of QPO oscillations in kHz band and permits to calculate the masses and radii of the neutron stars in X-ray binaries.

We developed theory of particle acceleration inside the relativistic rotating electron-positron force-free jet with spiral magnetic field. We considered perturbation of stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled. Acceleration process and synchrotron losses combined together form power law energy spectrum of ultrarelativistic electrons and positrons with index between 2 and 3 depending upon initial energy of injected particles. Synchrotron emission of these electrons and positrons in spiral magnetic field of rotating force-free jet has been calculated. Polarization properties of the radiation has been obtained and compared with existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give explanation of observed bimodality in alignment between electric field vector and jet axis. Degree of polarization and velocity of observed proper motion of bright knots depend upon angular rotational velocity of the jet. Thus, comparing them to each other, we can estimate angular rotational velocity in jets. We obtain that the fact, that generally in BL Lac objects electric field vector is oriented parallel to the jet axis while in quasars perpendicular to the jet axis, may be due to intrinsically larger angular rotational velocity and large winding of magnetic field in BL Lac jets than in quasar jets.

We investigate the linear stability of a hydrodynamic relativistic flow of magnetized plasma in the force--free approximation. We considered the case of light cylindrical jet in cold and dense environment, so the jet boundary remains at rest. Continuous and discrete spectra of frequencies are investigated analytically. An infinite sequence of eigenfrequencies is found near the edge of Alfv\'en continuum. Numerical calculations showed that modes are stable and have attenuation increment $\gamma$ small. The dispersion curves $\omega =\omega (k_\parallel )$ have a minimum for $k_{{\parallel}_0}\simeq 1/R$ ($R$ is the jet radius ). This results in accumulation of perturbations inside the jet with wavelength of the order of the jet radius. The wave crests of the perturbation pattern formed in such a way move along the jet with the velocity exceeding light speed. If one has relativistic electrons emitting synchrotron radiation inside the jet, than this pattern will be visible. This provide us with the new type of superluminal source. If the jet is oriented close to the line of sight, than the observer will see knots moving backward to the core.