Relativistic Effects of Rotation in γ-ray Pulsars—Invited Review

In this paper, we consider the relativistic effects of rotation in the magnetospheres of γ-ray pulsars. The paper reviews the progress achieved in this field during the last three decades. For this purpose, we examine the direct centrifugal acceleration of particles and the corresponding limiting factors: the constraints due to the curvature radiation and the inverse Compton scattering of electrons against soft photons. Based on the obtained results, the generation of parametrically excited Langmuir waves and the corresponding Landau–Langmuir centrifugal drive are studied.


Introduction
One of the fundamental problems of γ-ray pulsar astrophysics is the generation of VHE emission. In this context, it is worth noting that the last decade has been very fruitful in detecting the VHE γ-rays of pulsars [1][2][3][4][5]. On the other hand, it is clear that the generation of VHE emission requires the existence of relativistic particles, and the process of their acceleration is still under debate.
Generally speaking, it is observationally evident that pulsars are rotating neutron stars with an increasing period of rotation, P, characterized by the slow down rate, P ≡ dP/dt > 0. The corresponding slow down powerẆ = IΩ|Ω| then becomes: where the slow down rate is normalized by the typical value for normal period pulsars (with a period of the order of one second) , and we have taken into account that the moment of the inertia of a sphere is expressed as I = 2M ns R 2 /5, where M ns denotes the neutron star's mass, R 10 6 cm is its radius, and M 2 × 10 33 g is the solar mass. It is clear from the aforementioned expression that the slowing down of pulsar's rotation provides enormous power, which might be even higher for millisecond pulsars. Therefore, there is no question about a source of energy, but the problem concerns a particular mechanism that might be responsible for the acceleration of particles up to relativistic energies.
In the framework of the so-called polar-cap models by means of the strong electric field parallel to magnetic field lines (parallel electric field), the particles are uprooted from the surface of the neutron star [6], and the generation of electric fields, in turn, is a direct result of the neutron star's rotation. In the vicinity of the neutron stars, the electrons might undergo further acceleration provided either by magnetic field line curvature [7], or space-charge-limited flow [8], or the effects of the inertial frame [9]. On the other hand, the generation of a parallel electric field might be strongly suppressed either by means of the field curvature [10] or the IC scattering of electrons against photons [11]. Therefore, to explain the origin of VHE particles, outer gap models became actual. In the framework of this approach, the electrons gain energy in the outer magnetosphere in the so-called vacuum gap [12,13], but as it turned out, in these scenarios, the gap size is not enough to explain the observed VHE emission from pulsars. Some attempts to enlarge the gap size [9,14,15] although resulted in a more efficient acceleration process, but still could not solve the problem. The mechanism of energization along the magnetic field lines in the force-free pulsar magnetospheres was studied in detail by Contopulous et al. (see [16]), where the authors pointed out that acceleration is significantly suppressed by the force-free topology of the magnetic field. The acceleration process of a two component electronpositron plasma in a monopole configuration near the LC (a hypothetical zone, where the linear velocity of rotation exactly equals the speed of light) area was examined in [17], and it has been shown that the maximum Lorentz factors of particles are of the order of σ ≡ B 2 /(8πnm e c 2 ), where B is the magnetic field, n is the particle number density, m e is the electron's mass, and c is the speed of light.
A rather different mechanism of acceleration that might take place in the LC zone was proposed by Gold [18,19], who suggested that an extremely strong magnetic field forces particles to slide along the field lines, which are the trajectories of particles in the rotating frame of reference. This is an effective picture: the electrons gyrate in the magnetic field, and a trajectory averaged by time for long time scales is a certain trajectory, which in a very strong magnetic field approximation is a line, along which the particles slide. In particular, the magnetic field in the pulsar's magnetosphere is assumed to be dipolar √ PṖ Gauss is the magnetic induction on the neutron star's surface [20] and r s 10 km is the pulsar's radius. Then, for electrons with relativistic factors of the order of 10 5 , the gyro-radius in the LC area (a hypothetical zone where the linear velocity of rigid rotation exactly equals the speed of light) for normal period pulsars is of the order of 0.0016 cm, which is by many orders of magnitude smaller than the corresponding length scale of the LC radius: 10 9 cm. Therefore, the electrons are frozen into the magnetic field and follow the co-rotating field lines. The corresponding condition of the frozen-in state E + 1 c v × B = 0 explains the drift of particles in the direction of rotation as follows: by means of the radial magnetic field and the generated electric field perpendicular to the equatorial plane, the E × B/c 2 drift occurs, which actually is the co-rotation velocity, v, of the particle.
On the other hand, it is clear that the total velocity cannot be greater than the speed of light, and therefore, because of the co-rotation of the pulsar's magnetosphere [21], by reaching the LC surface (if the rigid rotation is preserved), the longitudinal velocity will tend to zero, and electrons will have only the azimuthal component, leading to the emission along the tangential direction ( Figure 1). Here, Ω denotes the angular velocity of the neutron star's rotation, r L is the radius of the LC, and α is the angle between the dipole magnetic moment and the axis of rotation.
In this review paper, we consider the centrifugally generated relativistic electrons for normal period pulsars, as well as for millisecond ones and present the progress achieved in this field during almost the last three decades. The paper is organized in the following way: Section 2 gives a detailed description of different phenomena that take place due to the centrifugal effects. We apply the model of centrifugal acceleration to the wellknown Crab and Vela pulsars and normal period pulsars, and in Section 3, we outline the obtained results. Figure 1. Sketch of the magnetospheric pulsar model from [22]. The centrifugally accelerated co-rotating particles in the nearby zone of the LC area emit in a relatively narrow cone along the tangential direction.

Results
In this section, we consider the effects of the rotation of the pulsar's magnetosphere in the LC area and study several problems: centrifugal acceleration of particles; rotationally driven Langmuir waves; the Landau-Langmuir centrifugal drive and the reconstruction of pulsars' magnetospheres by means of the centrifugal effects.

Centrifugal Acceleration
As we already mentioned in the Introduction, the longitudinal velocity component should decrease in the nearby zone of the LC. Deceleration of the radial velocity and thus the centrifugal force reversal has been investigated in the Schwarzschild black-hole metrics [23]. The corresponding problem, but in a special relativistic case, was studied in [24], where the authors considered a gedanken experiment: a bead freely sliding on a rotating wire, and studied the particle's dynamics. It has been demonstrated that the effect of "relativistic mass" leads to a similar phenomenon: if the rigid rotation is preserved, the radial velocity gradually decelerates and vanishes on the LC zone. For studying the problem of acceleration, one can perform calculations in the co-rotating reference frame. then, if one introduces an interval for a straight field line inclined by the angle α with respect to the rotating axis [24,25]: where Ω e f ≡ Ω sin α is the effective angular velocity of rotation for a field line inclined by the angle α with respect to the rotation axis. The particles moving along the co-rotating trajectories are characterized by the following Lagrangian: and the corresponding equation of motion: Considering the radial component of the equation of motion, one can show that the corresponding equation in the dimensionless form is given by: where x ≡ r/r L , r is the radial coordinate, r L = c/Ω denotes the LC radius, τ ≡ tΩ is the dimensionless time, dx/dτ = υ r /c ≡ u, and υ r denotes the radial velocity. In Figure 2, we show the behavior of radial velocity, u(x), and total velocity, υ(x) ≡ x + u(x) 2 , in the dimensionless forms for the following initial values x 0 = 0.001 and u 0 = 0 (note that the linear velocity of rotation in the mentioned dimensionless form coincides with the radial coordinate). As is clear from the plots, initially, the radial velocity increases, but in due time, by reaching the LC (x = 1), it starts decelerating, and the total velocity of the particle tends to the speed of light (υ → 1). By considering a plasma flow in the co-rotating pulsar's magnetosphere, it has been shown that in the limit of a strong magnetic field, the system of fluid equations reduces to a single particle approach described by Equation (5). Equation (3) for µ = 0 leads to energy-constant of motion: where m is the particle's mass, υ ϕ is the azimuthal velocity component, and γ = is the Lorentz factor of the particle. If one imposes the condition of constancy, the radial behavior of the relativistic factor is written as [22,26]: where γ 0 and r 0 denote the initial values of the Lorentz factor and radial coordinate, respectively. Equation (7) confirms that the particle's energy asymptotically increases close to the LC. From the observations, it is clear that particles originating in the pulsar magnetospheres leave these regions and dynamically become force-free. This means that the trajectory of the particle in the laboratory frame of reference is a straight line. Since the particles are in the frozen-in condition and follow the magnetic field lines, the shape of the latter in the co-rotating frame of reference should be given by Archimedes' spiral, lagging behind the rotation [27]. On the other hand, some of the γ-ray pulsars exhibit jet-like structures [28]; therefore, one should introduce a 3D configuration of magnetic field lines co-rotating with the neutron star. It is worth noting that the field lines will become more and more twisted in a narrow region of the LC. This particular problem has been considered in detail in [29]. It was shown that in the nearby zone of the LC, the curvature drift of plasma particles driving perpendicular to the equatorial plane leads to the generation of a toroidal component. This, in turn, changes the shape of the field lines into the configuration of Archimedes' spiral. For this purpose, in the co-rotating frame of reference, we examine the field line configuration: Imposing the aforementioned prescribed trajectory on the the space-time metric: will reduce it to the form: with: where ϕ = dϕ/dr, f = d f /dr, c = 1. Then, in the same manner as before, one can derive an expression for the energy for the non-diagonal metrics: where v ≡ dr/dt is the radial velocity and: is an expression for the relativistic factor in the general case-3D curved trajectories. One can explicitly solve Equations (12) and (13) for the radial velocity: where different signs correspond to different initial conditions. From Figure 2, it is clear that a particle sliding along rectilinear co-rotating field lines never crosses the LC zone. The reason follows from the nature of the rigid rotation of straight field lines. On the other hand, as we have already mentioned, the particles originating in the pulsars' magnetospheres might leave this region, becoming dynamically free. This, in turn, leads to a specific shape of the field lines in the local frame of reference-Archimede's spiral [28]. u(x) Figure 2. Here, we plot u and υ versus the dimensionless radial distance, x. The set of parameters is x 0 = 0.001 and u 0 = 0.
As a first example, we consider curved field lines having the shape of the Archimedean spiral located in the equatorial plane: For a spiral that lags behind the rotation (a < 0), one can find a critical velocity when there is no "reaction force" acting from the curved trajectory. In this case, the angular velocity of the particle in the laboratory frame, Ω l = Ω + aυ, equals zero, and therefore, the corresponding critical radial velocity is defined as υ cr = −Ω/a.
In Figure 3, we plot the behavior of the Lorentz factor versus the radial coordinate for different initial relativistic factors γ 0 . The calculations are performed for an Archimedean spiral with the critical Lorentz factor γ cr = 7. All particles are launched from r 0 = 0. It is clear from the figure that if the initial value of γ equals seven, the particle does not experience the reaction force, and therefore, it stays constant. Unlike this case, if the initial values are different from the critical value of the relativistic factor, it asymptotically tends to γ cr . From the plots, it is evident that the particles go beyond the LC zone. The reason is the following: since in due time, the Lorentz factor tends to its critical value, the effective angular velocity, Ω l , tends to zero, and consequently, the corresponding LC radius "inflates", while the particle always stays inside c/Ω l , but outside the distance, c/Ω.  All particles are launched from r 0 = 0. As is evident from the plot, the curve γ cr attracts all other curves. The figure is taken from [22].
A second example we would like to consider is the field lines located in a common plane with the rotation axis: After imposing this condition, from Equation (14), one can derive a longitudinal velocity (along a field line) as follows: where: and γ 0 and r 0 are initial values quantities. The centrifugal mechanism of acceleration can explain the Crab pulsar's jet-like wind's velocity. From observations, it is evident that jets are collimated flows. Therefore, one can model them as the particles moving along the field lines with the following asymptotic structure: r → R and z → ∞ (R < r L ). The aforementioned velocity in units of the speed of light has its maximum value, which can be straightforwardly obtained from Equation (17) (see for details [28]), v max = 0.5. On the other hand, the velocity averaged across the jet should be less than the possible maximum value, which is known from observations: υ Jet 0.4 [28]. As a last example, it is natural to consider a field configuration given by: (19) which, on the one hand, is characterized by a 3D geometry and, on the other hand, has the property of the Archimedean spiral. After straightforward mathematical manipulations, from Equation (14), one can derive an asymptotic behavior of the radial velocity: It is evident from this expression that the effective angular velocity Ω l ∝ 1/r 2 is a continuously decreasing function of the distance and asymptotically tends to zero, implying that the particle becomes dynamically force-free.
In Figure 4, we show the trajectories of particles in (a) the rotational frame (RF) and (b) the laboratory frame (LF), respectively. we also show (c) the behavior of the curvature (normalized by the initial value) of the trajectory in the LF and (d) the jet velocity, thus the velocity along the z axis. The set of parameters is: a = −28, f = 11, v 0 = 0.01, and r 0 = 0. It is clear that despite the fact that the particle's trajectory in the RF is given by a regularly shaped 3D-Archimedean spiral, in the LF, the trajectory is more complex (see (b)). Initially, the trajectory is more curved (see (c)), but in due time, the curvature vanishes, and as a result, the trajectory in the LF becomes more rectilinear (see (b)). This means that the dynamics tends to the force-free regime. Indeed, from Figure 4d, one can see that the z component of the velocity initially increases, but asymptotically saturates, indicating the transition to the zero force regime.
Since the magnetospheres of pulsars are full of the photon field, similar problems have been studied for a situation when a co-rotating system is imbedded in the isotropic photon field, leading to a certain friction force influencing the dynamics of particle acceleration. Like the previous cases, it is shown that when the trajectories are presented by the Archimedean spiral with arms lagging behind the rotation, the particles reach the force-free regime whatever the value of the friction force is. An interesting feature has been found for the Archimedean spiral with arms oriented in the opposite direction. In this case, the particles do not tend to the force-free regime, but instead reach a certain equilibrium location. Similar results have been obtained for 3D configurations of the field lines (for details, the reader is referred to [30,31]).  [28]. In the top panel, we show the particle trajectories in (a) the rotational frame (RF) of reference and (b) the laboratory frame (LF) of reference, respectively. In Graphs (c,d), the behavior of the curvature (ρ ≡ 1/R c , where R c is the curvature radius of the magnetic field lines) and jet velocity is shown, respectively. The set of parameters is: a = −28, f = 11, v 0 = 0.01, and r 0 = 0.

Co-Rotation Constraints
It is natural to suppose that for almost straight field lines, rigid rotation cannot be maintained up to the LC surface; therefore, in order not to violate the relativity principle, the pulsar's magnetosphere must reconstruct in a proper way. In particular, the rigid rotation is maintained by means of the strong magnetic field. Therefore, the co-rotation condition is satisfied if the magnetic field energy density, B 2 /8π, exceeds the plasma energy density, γn GJ Mm e c 2 , where B denotes the magnetic induction, n GJ = ΩB/(2πec) is the Goldreich-Julian number density of electrons in the pulsar's magnetosphere [32], e is the electron's charge, and M is the multiplicity factor of particles. In this context, one has to note that in the pulsar models, it is assumed that primary particles, accelerated by the parallel electric field, radiate due to the curvature radiation. The interaction of the emitted photons with the strong magnetic field will inevitably lead to the very efficient pair creation process (γ + B → e + + e − + B), and newly generated secondaries repeat the same cascading mechanism until the plasma particles screen the electric field [33,34].
By taking the magnetic field strength and the pulsar's radius r s 10 km into account, combined with the aforementioned condition and Equation (7), one can show that for r 0 << r, the maximum attainable relativistic factor of electrons moving on the open field lines in the LC zone is written as [22,35]: where the initial Lorentz factor is normalized by its typical value in the Crab-like pulsars [35], P andṖ are normalized by the Crab pulsar's parameters, P Cr 0.0332 s, anḋ P Cr 4.21 × 10 −13 ss −1 . The condition r r L / sin α has been taken into account, and we considered the open field lines with curvature radii exceeding the LC length scale, which means that in the framework of this approximation, the field lines are assumed to be straight. As is clear from the aforementioned expression in the Crab-like pulsar's magnetospheres, the particles might achieve extremely high energies. One also can straightforwardly show that for normal pulsars (P 1 s,Ṗ 10 −15 ss −1 ), the co-rotation condition might guarantee the maximum attainable relativistic factor of the order of 2.6 × 10 5 .

Emission Constraints
The approach presented in the previous section is simplified for several reasons: in realistic astrophysical scenarios, the dynamics of particles is strongly affected either by the IC scattering of relativistic particles against soft photons or curvature radiation. Unlike them, the synchrotron mechanism cannot impose constraints on the motion of electrons. In particular, by taking the single particle synchrotron power into account, P syn 2e 2 ω 2 B γ 2 /3c, where ω B = eB/mc is the cyclotron frequency, and the synchrotron cooling time scale, t syn γmc 2 /P syn is of the order of 6 × 10 −9 s for γ = 10 5 . This value is very small compared to the kinematic time scale of pulsars (rotation period); therefore, immediately after becoming relativistic, the electrons lose their perpendicular momentum, go to the ground Landau level, and continue sliding along the field lines, which, in turn, guarantees the centrifugal mechanism of acceleration.
In the previous subsection, we showed that the efficiency of rotation energy pumping significantly increases in the LC area. If one assumes a straight co-rotating field line inclined by the angle α with respect to the rotation axis, then the acceleration time scale on the LC: reduces to: From Equation (21), it is evident that in the nearby zone of the LC, the relativistic factor increases. This, in turn, leads to the decrease of the acceleration time scale (see Equation (23)), indicating the high efficiency of the acceleration process. On the other hand, the pulsar magnetosphere is full of soft-thermal-photons originating from the neutron star's surface, and the IC scattering seems to be important. In this context, one has to note that very often, the role of the IC mechanism with thermal photons is supposed to be negligible because it has been accepted that the photons and relativistic electrons travel in the same-radial-direction [36]. Contrary to this, following the idea, originally proposed by Gold (see [18,19]), a strong magnetic field (see the previous subsection) guarantees the frozen-in condition of electrons almost up to the LC zone, resulting in the co-rotation of particles, having only the azimuthal velocity component; therefore, the angle of interaction is ∼90 • , and consequently, the IC scattering of electrons against thermal photons should be significant.
In the standard theory, the temperature of a neutron star's surface strongly depends on its age, τ = P/(2Ṗ) [37] T(τ) 5.9 × 10 5 × therefore, for the Crab-like pulsars, the age is of the order of 10 3 y, and the surface temperature is T 1.2 × 10 6 K; for the normal pulsars (P = 1 s,Ṗ = 10 −15 ss −1 ), the age is three orders of magnitude bigger, and consequently, T 2.4 × 10 5 K. For the peak frequencies ν 2.8 kT/h and the corresponding photon energies, one has ν 6.7 × 10 16 K and max 0.3 keV for the Crab-type pulsars and ν 4.1 × 10 15 K and max 0.02 keV for the normal pulsars, respectively (here, h denotes Planck's constant and k is Boltzmann's constant). In general, the IC scattering might occur in two extreme regimes, in the so-called Thomson regime for very low energy particles and in the Klein-Nishina regime for very energetic electrons. In particular, the former takes place if the following condition is satisfied γ ph /(mc 2 ) << 1 [38] and in the latter case γ ph /(mc 2 ) >> 1, where ph is the photon's energy. Then, one can straightforwardly check that for both classes of pulsars, where we have already shown that γ 10 7 (Crab-like pulsars) and γ 10 5 (normal pulsars), the IC process occurs in the KN regime when the photons gain almost the total energy of electrons. Then, after taking into account the spherically symmetric photon field, the KN-Compton power is given by [38]: Logarithmic dependence on the Lorentz factor implies that the IC cooling time scale in the LC area: is a continuously increasing function of γ, which indicates that the IC process does not impose any significant constraints on the maximum attainable energies. Therefore, the maximum energy is limited by the co-rotation constraint, and consequently, after scattering, photons will have energies: , TeV (27) in Crab-type pulsars and 130 MeV in normal pulsars. Despite the fact that the IC mechanism does not affect the maximum achievable energy of electrons, the produced VHE luminosity might be quite high. For calculating this particular quantity, one should emphasize that the acceleration is extremely efficient in a narrow thickness, where the relativistic factor asymptotically increases and the corresponding scaling factor-rate of change, dγ/dr-increases as well. Therefore, acceleration is more efficient for highly nonuniform γ's, and hence, the corresponding thickness length scale can be estimated to be d γ dγ/dr γ 0 r L /(2γ), where we have taken into account Equation (21) for r 0 = 0. The azimuthal length scale involved in the emission process is of the order of δl r L θ, where θ 2πξ/10 (ξ < 10), and consequently, for the active volume one has ∆V (δl) 2 d.
We discussed in the previous subsection that the pair cascading is very efficient. The corresponding number density of electron-positron pairs is significantly increased, and as was shown in [32], the multiplicity factor equals 1/(1 − r 2 /r 2 L ), which combined with Equation (7) leads to M = γ/γ 0 . After taking into account all these quantities, the expression for the bolometric IC luminosity is given by: where we have taken into account that in the generation of the TeV emission, the most energetic photons are involved. The normal pulsars (with T 2.4 × 10 5 K) in the VHE domain (∼100 MeV) might provide luminosity of the order of ∼1.6 × 10 21 ergs s −1 .
Even if the field lines are almost straight, due to the co-rotation condition, the particles follow the field lines, and therefore, in the laboratory frame of reference, their trajectories might be significantly curved, leading to the mechanism of curvature radiation. In the context of pulsars, this mechanism was investigated in [39,40], where the authors discussed in detail the generation of radio emission by means of the curvature radiation mechanism. We will see that this process might account for the generation of much higher energies.
The energy loss rate by means of the curvature radiation of the electron is given by: where R c denotes the radius of the curvature. The characteristic cooling time scale, t c ≡ γmc 2 /P c , then becomes: It is clear from this expression that t c is vanishing faster than the acceleration time scale, t acc ∝ 1/γ 1/2 (see Equation (22)); therefore, the maximum attainable Lorentz factor, γ c , is achieved when t acc t c , leading to the following expression of γ c : where we have taken into account that the curvature radius, R c , of the trajectory in the nearby zone of the LC is of the order of r L , which is clearly seen from Equation (7). In particular, close to the mentioned area, the Lorentz factor asymptotically increases, and the linear velocity of rotation tends to the speed of light; therefore, the radial velocity tends to zero, and consequently, the curvature radius equals the LC radius.
Since the derived value of the maximum relativistic factor is higher than the Lorentz factor limited by the co-rotation constraint (see Equation (21)), the latter is responsible for reaching the maximum energy of electrons, leading to the following energy of curvature photons [35]: GeV.
As we see, the curvature radiation provides the emission in the GeV range for the millisecond pulsars. For the total curvature luminosity, one obtains: Despite the fact that the curvature emission mechanism is not responsible for limiting the maximum attainable energy of particles, most of the energy is emitted by means of it.
For the normal pulsars (P 1 s,Ṗ 10 −15 ss −1 ), one can straightforwardly check that the co-rotation still remains the major mechanism limiting the maximum energy, and the curvature radiation provides emission in the upper UV domain ( c 100 eV) with the luminosity of the order of 2.4 × 10 22 ergs s −1 .

Parametric Excitation of Electrostatic Waves and Langmuir-Landau Centrifugal Drive
In the previous subsection, we showed that the centrifugal mechanism of particle acceleration might be very efficient in millisecond, as well as normal pulsars' magnetospheres. It is worth noting that different species of centrifugally accelerated particles experience different forces, leading to charge separation and a consequent generation of the electrostatic field. On the other hand, as was found in [25,41], the relativistic centrifugal force is time dependent, and therefore, the excitation process is parametrically amplified.
For studying the dynamics of a plasma flow following the magnetic field lines, it is assumed that the field lines are almost rectilinear, which is satisfied for distances less than the curvature radius. For simplicity, one can consider the dynamics of plasma particles in the local frame of reference of a rotating pulsar. If one considers the 1 + 1 formalism [42] in the framework of a single particle approximation, for the chosen reference frame, the zero angular momentum observers (ZAMOs) will measure the following proper time dτ ≡ ζdt, where dt is the universal time and ζ ≡ 1 − Ω 2 e f r 2 /c 2 is the so-called lapse function. The equation of motion then is given by the following equation [25]: where E is the electric field, p → p/mc is the dimensionless momentum; γg is the effective centrifugal force [43] with g = −∇ζ/ζ, and in the relativistic factor γ = 1 − V 2 /c 2 1/2 , the velocity is defined according to the 1 + 1 approach V = dr/dτ. In order to rewrite the aforementioned equation for a magnetospheric fluid, one should emphasize a complex composition of plasmas. In Section 2.1.1, we briefly discussed that initially, accelerated particles radiate by means of the curvature process. Newly produced photons under certain conditions might generate secondary electron positron pairs. Finally, the magnetospheric e + , e − population may be divided into two major components: ( for several species of particles denoted by i = b, e, p, where b corresponds to bulk particles (mostly electrons), e-to electrons and p-to positrons. One should note that the ZAMO's momentum is coincident with the momentum measured in the inertial frame of reference. Indeed, by combining the definition of momentum, p = γV, where γ = ζγ and V = dr/dτ (prime denotes the quantities in the inertial frame of reference), one can straightforwardly show p = p . Then, for Equation (35) rewritten in the inertial frame, one has (omitting primes): As we already discussed, the centrifugal force, acting on different species, leads to charge separation, which automatically creates the electrostatic field. Therefore, for describing the physical system, the momentum equation should be complemented by the continuity equation and the Poisson equation: where n i denotes the number density of a corresponding component. It is worth noting that a strong magnetic field guarantees the frozen-in condition of plasmas, E 0 + 1 c v 0i × B 0 = 0, which is supposed to be a leading state of the system. By taking into account this condition, the momentum equation reduces to Equation (5), which in a dimensional form is given by: where υ = dr/dt denotes the radial velocity. In [24], it has been shown that for ultrarelativistic particles (γ >> 1), Equation (39) has the following solution: where φ denotes the phase. After taking into account these solutions, one can straightforwardly show that the analog of the centrifugal force, −ζ i γ i ∇ζ i , is time dependent and drives the parametric excitation of the Langmuir waves. In this paper, we intend to study the linear generation of the electrostatic waves. Therefore, all physical quantities are expanded around the leading state Ψ 0 : where Ψ = (n, v, p, E, B) and Ψ 1 denotes a small perturbation of the corresponding quantity. It is straightforward to show that the ansatz: leads to the following linearized versions of Equations (36)- (38): By introducing the so-called plasma quantities for perturbations n 1 pl = n 1 e − n 1 p , υ 1 pl = υ 1 e − υ 1 p (the plasma component) with the unperturbed quantities of electrons and positrons n 0 e = n 0 p and υ 0 e = υ 0 p , and considering the bulk component and the plasma component, respectively, after straightforward, but tedious mathematical manipulations, one arrives at the coupled non-autonomous differential equations governing the centrifugal generation of the electrostatic instability [44]: where ω 1,2 = 8πe 2 n 1,2 /mγ 3 1,2 and N 1 and N 2 relate to the first-order perturbations of the number densities of two species in the following way: where φ β denote the corresponding phases and Ξ is expressed in the following way: and with 2φ ± = φ 1 ± φ 2 . We have assumed that energy is almost uniformly distributed among the species, n 1 γ 1 n 2 γ 2 n GJ γ b . After Fourier transforming Equations (47) and (48) by applying the tools developed in [25], one can arrive at the following dispersion relation: From the obtained expression, it is clear that the resonance frequency is given by ω r = ω 2 1 + ω 2 2 J 2 0 (b) 1/2 . Following the method described in [25], expanding the frequency, ω = ω r + ∆, and emphasizing that near the resonance, the basic contribution to the r.h.s. of the equation comes form the term ω r = −µ 0 Ω, one obtains the following equation: Solving this cubic equation, one can show that the imaginary part of the solution, thus the growth rate of the instability, is written as [45]: Already amplified electrostatic waves will inevitably Landau damp with the following increment [46]: It is worth noting that the most optimal regime is when both growth rates are of the same order of magnitude. For example, if one considers the Crab pulsar with the possible maximum relativistic factor provided by the direct CA, 1.5 × 10 7 (see Equation (21)), one can see that two streams with γ 1 10 4 and γ 2 7 × 10 4 drive the electrostatic instability with the corresponding time scale, 1/Γ, of the order of 0.01 s. For the aforementioned stream parameters (considering the Stream 1 component as to be the plasma component), both growth rates are equal, and since the time scale is less than the kinematic time scale (period of rotation) of the Crab pulsar, the process is efficient. Therefore, the energy pumped by the primary beam by means of the Landau process should be enormous.
In particular, by considering the force responsible for the particle acceleration, F r 2mcΩ ζ −3 [45], the energy accumulated in the primary beam after Landau damping is given by [45]: n p F r δr n GJ and for the aforementioned two streams is of the order of 3.7 PeV. Here, δr c/Γ is the corresponding length scale where the process of re-acceleration takes place. After applying the same mechanism to normal pulsars (P 1 s,Ṗ 10 −15 ss −1 ) with the beam Lorentz factor of the order of 2.6 × 10 5 , one can straightforwardly check that two streams with γ 1 10 2 and γ 2 5 × 10 2 lead to the parametric instability of the generation of Langmuir waves with a high growth rate Γ ∼ Γ L 46 s −1 . Consequently, by Equation (56), one can check the re-acceleration process in normal pulsars; although, it cannot guarantee such high energies as in the Crab-like pulsars, still, energies will be in the VHE domain, 360 GeV. In Section 2.1.2, we showed that the IC scattering and synchrotron mechanism cannot affect the dynamics of particles. The Langmuir waves propagate along the magnetic field lines, and the resulting re-acceleration will not cause the pitch angle significantly suppressing the synchrotron process. We already discussed that from observations, it is evident that cosmic particles leave the pulsar's magnetospheres, which indicates that they are in the force-free regime. Therefore, as was studied in [28,29,47], the field lines lag behind the rotation, having the shape of the Archimedean spiral characterized by the rectilinear trajectories in the laboratory frame of reference. This means that the curvature emission is also strongly suppressed and does not impose any significant constraints on the maximum achievable energies of electrons.

Conclusions
Centrifugal processes occurring in the pulsar's magnetospheres might lead to new and very interesting phenomena. This may allow new insights into the physics of these extreme environments. The detection of the deceleration of pulsars indicates that a significant fraction of neutron stars' rotation energy is extracted from the central object. A certain fraction of it can be pumped into the magnetospheric electrons and electromagnetic radiation.
In particular, we have demonstrated that direct CA is a very efficient mechanism providing the observed jet velocity of the Crab pulsar. This process is responsible also for the generation of relativistic particles with the Lorentz factors of the order of 10 5 (for normal pulsars) and 10 7 (for millisecond Crab-like pulsars).
It has been found that IC scattering might provide photon energies of the order from 100 MeV up to several TeV depending on a pulsar's type. The total luminosity for Crab-like pulsars might be at least 10 29 ergs s −1 and for relatively slowly rotating pulsars −10 21 ergs s −1 .
Although the curvature mechanism provides electromagnetic radiation in the GeV (Crablike pulsars) and UV (normal pulsars) domains, the corresponding luminosities are higher than in the IC scattering: 10 34 ergs s −1 (Crab-like pulsars) and 10 22 ergs s −1 (normal pulsars).
We also examined the effects of CA on the excitation of Langmuir waves. It was found that by means of the strong centrifugal force in the LC zone, the electrostatic waves are parametrically excited. The corresponding instability time scale is much less than the rotation period of a neutron star, indicating the high efficiency of the driving process.
The amplification of the Langmuir waves will be terminated by means of the Landau damping, resulting in the extremely high energies of particles. We found that the electrons in Crab-like pulsars after re-acceleration via the mechanism of the Langmuir-Landau centrifugal drive might reach energies in the PeV domain. For the normal pulsars, the achieved energies of electrons are of the order of several hundreds of GeV.
The centrifugal parametric excitation of Langmuir waves has a very interesting consequence. In particular, as we showed, the growth rate of the instability is so large that the electrostatic waves are efficiently amplified. During this process, the field can reach the Schwinger limit, and consequently, an efficient pair creation will take place. This, in turn, might significantly influence physical processes in the γ-ray pulsar's magnetospheres, and sooner or later, we are going to study this particular problem.
Funding: This research received no external funding.