Particle Acceleration by an Electromagnetic Wave

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.


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   =  0 cos ( − ) ,   =  0 sin ( − ) ,   = 0,   = − 0 sin ( − ) ,   =  0 cos ( − ) ,   = 0 (1) 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.

Solution of the Equation of Motion
The relativistic equation of motion Erschienen in: ISRN High Energy Physics ; (2013).

-Article ID 760983
Konstanzer Online-Publikations-System (KOPS) URL: http://nbn-resolving.de/urn:nbn:de:bsz:352-257676 for a charged point particle (rest mass  0 , 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 ( = 0 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 3 . (17 The maximum velocity ż will be reached for  −  =  and amounts to 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 zdirection.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  max 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 ż = 0 according to (15) for  = 2 corresponding to  = 2(1 + ( 2  2 0 / 2 0  2  2 )) and  = 2( 2  2 0 / 2 0  3 ).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, ẋ = 0 is valid for  = , where ẋ changes its sign and  goes back to zero ( = 0) for  = 2.In contrast to this, ẏ reaches its maximum for  =  with the value ẏ = ( 0 / 0 )(2/(1 + 2( 2  2 0 / 2 0  2  2 ))), and ẏ = 0 is valid for  = 2, where  = 2( 0 / 0  2 ).For ( 0 / 0 ) ≫ 1, the maximum velocity is small compared with  and (2)/(2) ≫ 1 is valid so that the acceleration happens mainly in the direction of the wave (-direction).

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,  = /2 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 barmagnet 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)) ⃗  ⊥ = | ⃗  ⊥ |(− sin , cos , 0).(Vector in the ; -plane perpendicular to the rotation axis (-axis).)Herewith and with | ⃗ | = (/4) ⃗  2 , 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 9, 64 ⋅ 10 38 sin 2  erg/sec) as follows: for protons as particles and the star data (29).For the acceleration time (24), one obtains in the case of protons Otherwise, the negative half-wave will be important and the particle will be retarded, cf.(17).However, because of  max ≫ /2, 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, sin 2  = 10 −7 result in reasonable values for  max ,  max , and  max for protons: The comparison with  max , (38), shows once more that the accelerated motion results mostly in the -direction.
In Table 1, we show the values for  and  in the range 0 <  < 10 −2 for protons and the neutron star data (29) (see also (35)).It is interesting that up to distances of 10 15 km  ≪ 1 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.

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.
2   0 3  sin  = 3, 3 ⋅ 10 9 sin , 64 ⋅ 10 16 sin 2  sec Obviously, the path length of the light  max is larger than  max only by /2 of the wave as follows: