The relativistic equation of motion for an electrically charged massive point particle is integrated exactly for the case of a circularly polarized electromagnetic plane wave. The acceleration of the particle by the wave is calculated in detail, however, under neglection of the radiation damping. In spite of this simplification, the results are applied to a typical pulsar, where the maximum kinetic energy of accelerated protons reaches values of the order of the fastest cosmic ray particles.

1. Introduction

The present model for a pulsar is that of an oblique rotator, where the magnetic moment of the neutron star is inclined against the rotation axis of the star (-axis) by an angle . Then, in consequence of the time variable magnetic moment perpendicular to the rotation axis, a magnetic dipole radiation is emitted, which is circularly polarized (see Figure 1).

For this reason it may be interesting to calculate the acceleration of an electrically charged particle by a circularly polarized electromagnetic wave propagating in the -direction. For simplicity, we assume a plane wave with the field strength satisfying the sourceless Maxwell equations. An additional constant phase parameter is not necessary in (1), because the - and -axis can be chosen in any case in such a way that (1) is valid by turning around the -axis. (This simplification is not possible in the case of a linearly polarized wave.) The application of the ansatz (1) to a neutron star is possible only with reservation, first because the electromagnetic wave field (1) represents only the electromagnetic far field and second because the decrease of the intensity of the wave with increasing distance from the star is not taken into account. However, the maximum of the particle energy attainable by such an acceleration process may be evaluable, which is of interest for explaining the energy observed in cosmic rays. Beyond this, the results may be applicable in the accelerator physics.

Such calculations are already performed for pulsars first by Gunn and Ostriker [1, 2] in 1969 as well as later by Grewing and Heintzmann [3] and for accelerators by Holkundkar et al. [4]. But all these calculations are not performed by an exact analytical integration of the relativistic equation of motion of the charged particle. To do that is the aim of the present paper.

2. Solution of the Equation of Motion

The relativistic equation of motion for a charged point particle (rest mass , electric charge ) reads in the case of the field (1) in cartesian coordinates The radiation damping is neglected for simplicity. Accordingly, a particle at rest will be accelerated at first by the electric field in - and -direction and then by the magnetic field in -direction. Furthermore, from (2) the energy law follows immediately: which takes the following form according to (1):

Insertion of (1) into the first two equations of (3) gives () Because of the dispersion relation of the electromagnetic wave, the integrals of (6) read using the following initial condition ( may be identified later with the surface of the star): which means acceleration out of the rest state. Herewith, the third equation of (3) and (5) take the following form using (1): Comparison of these two equations gives with the integral taking into account the initial values (9) By squaring of (8) and (12) and by addition, we find the following kinetic energy: and herewith from (12) analogously to (8) we have The substitution of the square root in (14) by the use of (13) results finally in and from (8) and (14), one finds the following quotient: For the acceleration in -direction, one obtains from (15) immediately The maximum velocity will be reached for and amounts to For , it is valid: .

Obviously, the maximum kinetic energy , which can be reached by the acceleration of the electromagnetic wave, is given according to (13) also for by Evidently, the energy transfer is very effective for low frequencies and high intensities. However, the intensities are correlated with at least the fourth power of the frequencies so that pulsars with high frequencies will produce the particles with the highest energies; see Section 3.

Finally, the total integration of the equations of motion is possible; however, of especial interest is the motion in z-direction. With the substitution one obtains from (12) and (15) with the integral for being valid for the moving particle using the initial values of (9). Subsequently, it follows for from (19):

The acceleration of the particle in positive -direction happens according to (17), see also (18), within the range Then, the end of the acceleration phase is reached according to (21) for at (see (22)) If is very large, the intensity of a real spherical wave is decreased drastically so that the interaction between particle and wave can be considered as finished. Otherwise the particle will be decelerated and reaches the initial velocity according to (15) for corresponding to and .

For completeness, we integrate also the motion in - and -direction. Equation (20) gives Then, it follows from (8) that with the integrals using the initial condition (9) where is given by (21). It is interesting that the very easy integration is possible because of the relativistic factor (26) in (8) (velocity dependence of the mass). Obviously, is valid for , where changes its sign and goes back to zero () for . In contrast to this, reaches its maximum for with the value , and is valid for , where . For , the maximum -velocity is small compared with and is valid so that the acceleration happens mainly in the direction of the wave (-direction).

3. Application to a Neutron Star

For the comparison with a real situation, we choose a neutron star with the following data ( radius of the star, rotating frequency): For the dipole radiation intensity (Poynting vector), it is valid: where is the angle between and the magnetic moment . For radiation in the direction of the rotation axis (-axis) of the pulsar, holds (maximum intensity). On the other hand, we use for the connection between the magnetic field strength and the magnetic moment the relation for a bar-magnet outside, but on the axis (-axis) of the magnet. At the magnetic pole () we find Furthermore, for the oblique rotator it is valid ( angle between and -axis (see Figure 1) as follows: where (cf. (1) and (32)) . (Vector in the ; -plane perpendicular to the rotation axis (-axis).) Herewith and with , we find from (30) for the wave intensity the following at the rotation pole of the star . Insertion of (32) into (33) and of (33) into (34) results in Herewith, we find from (18) the maximum kinetic energy (the total radiation power of the dipole-radiation of the pulsar amounts to ) as follows: for protons as particles and the star data (29). For the acceleration time (24), one obtains in the case of protons and for the acceleration distance (25) one obtains Obviously, the path length of the light is larger than only by of the wave as follows: Otherwise, the negative half-wave will be important and the particle will be retarded, cf. (17). However, because of , the relation (39) means that the particle runs nearly with the speed of light surfriding on the wave. For the ratio (16) of the velocities in z- and y-direction one gets so that the motion results mostly in the z-direction, that is, in the direction of the wave. Finally, we see that also very small angles, for example, result in reasonable values for , and for protons: At the end of the acceleration phase () it follows: The comparison with , (38), shows once more that the accelerated motion results mostly in the -direction.

Besides these maximum values, it is of especial interest to determine the attainable energy in the distance from the star in dependence of the phase . From (21) to (22), and from (13) to (19) it follows immediately: In Table 1, we show the values for and in the range for protons and the neutron star data (29) (see also (35)). It is interesting that up to distances of  km is valid, which demonstrates once more that the protons surfride on the wave very early, so that changes its values only very slowly. On the other hand, the large -values for higher energies show clearly that the neglection of the decrease of the wave intensity with increasing distance is indeed a very rough approximation.

4. Final Remark

Some critical remarks are necessary concerning the application of the results of Section 2 to neutron stars. The reached maximum energy of protons is indeed of the order of magnitude of cosmic rays. However, we have neglected the enormous gravitational attraction of the particle by the neutron star. Furthermore, we have used the following simplifications:(1)plane wave instead of a spherical one, (2)using only the radiation far field, (3)neglection of the influences on the wave by the surrounding plasma (see e.g., Kegel [5] and Thielheim [6]) and by the strong gravitational field, (4)neglection of the radiation damping. Therefore, the results for the neutron star can be considered only as very rough estimations.