Emission of electromagnetic radiation is a natural classical phenomenon appearing because of the acceleration of any charged particle  . This radiation implies emission of energy of the charge and, consequently, a damping effect on its motion, which is modeled through a force called “radiation reaction force”. An alternative formulation to the usual Abraham-Lorentz-Dirac radiation reaction force    has appeared recently  , where the radiation reaction force is given in terms of the same external force which brings about the acceleration of the charge. This approach intents to solve some difficulties of the previous formulation  - , where the radiation reaction is given in terms of the acceleration which brings about pre and post acceleration of a charged particle after the external force appears and disappears, having no experimental or consistent theoretical justification. The novel approach to radiation reaction has not been experimentally tested yet, although some theoretical studies have been done already  . Because of this experimental requirement, we make in this paper a study of the 1-D motion of an electron under a periodic external force. The modified relativistic equation of 1-D motion of a charged particle is given by 
where with being the normalized velocity of the charge, and F is the external force. This is the equation of motion that we will use in our study. We will study first the dynamics under the simple sinusoidal force
where , and are the amplitude, the angular frequency and the phases of the electric field, and then we will extend our study to the position and time depending force
2. Dynamical Equations
Using the external force (2) and (3) in (1), and after some rearrangements, we get the following dynamical systems
where the term in (4a) and (5a) appears from the differentiation of with respect the time. In the first dynamical system the motion in the coordinates “x” and “ ” are disconnected, but on the second dynamical system these coordinates are coupled. These dynamical systems are solved using Runge-Kutta method at fourth order and are solved for the cases without radiation reaction force ( , the solution is denoted as ) and with radiation reaction force ( , the solution is denoted as ). Then, we calculate the difference
3. Numerical Analysis
We solved the first dynamical system, Equations (4a) and (4b), for different but fixed (electric field intensities) to see the change in the electron’s dynamics behavior. This is done for , , and for , , simultaneously at each time “t”, and the difference (6) is calculated. The inital conditions are the same, position and normalized velocity . Figure 1 shows the expression , solution of first dynamical system with Equations (4a) and (4b), as a function of time for several values. As one can see, this figure has some peaks and discontinuity, the one that is in the middle can be associated with a change of force’s direction, this will happen when (for this case and ). The peaks at the very beginning are due to the fact that and are very close together in magnitude, that is and as is known . In addition, we can see that the difference increases with the intensity of electric field. From Equation (4a) we notice that when the electron’s normalized velocity is closed to one ( ), the external force approaches to zero leading to a almost constant, and this is observed for . As stronger the electric field
Figure 1. Dynamical system (4): Logarithm of as a function of time for different electric field intensities in units of Dynes/ues with a fixed frequency of . Each curve is associated to the shown intensity of the electric field .
is, the difference between non-radiation reaction dynamics and with radiation reaction dynamics become more separated.
Figure 2 shows the same function as before but with the solutions for the second dynamical system (5a) and (5b). For this case, the electron’s behavior is a little more complicated, as we can see there are more peaks and discontinuities. As increases the difference increases for times above where almost a constant value is reached. The peaks appear because the difference between normalized velocities starts to decrease at those points, and this occurs when the is very close to the speed of light. Thus from (5a), the external force starts to decrease leading to .
Finally we studied the case for different frequencies leaving the electric field intensity fixed at . The results for the dynamical system (4a) and (4b) are shown in Figure 3. We can see that for higher frequencies it appears more peaks and discontinuities at the beginning of the motion of the electron, and the reason is due to the same as the previous case, electron speed leads to a difference . Also we have three cases of interest here, first one for frequencies of order 1 GHz to 10 GHz, the second one for 100 GHz to 1000 GHz and the last one when . For the first case, we observe that is almost constant, this is because electron’s speed reaches out to the relativistic case very soon and the external force disappears. For the second one, force’s frequency leads to a change of force’s direction in a way that the electron doesn’t
Figure 2. Dynamical system (5): Logarithm of as a function of time for fixed frequency . Each curve is associated to the shown intensity of the electric field .
Figure 3. Dynamical system (4): Logarithm of as a function of time for different frequencies (in Hz) and fixed .
Figure 4. Dynamical system (5): Logarithm of as a function of time for fixed and several frequencies (in Hz).
reach the relativistic case for the simulation time considered. And for the last case, the frequency makes the force changes direction so rapidly that we start to have a lot of discontinuities until electron speed is stable.
For the dynamical system (5) in Figure 4, we have fixed the electric field intensity at and we change the frequencies from . One can see that electron’s behavior is very complicated, we have more discontinuities while the frequency is higher until the electron reaches a relativistic speed at frequency and the external force disappears. Those discontinuities can be explained similarly as we have done with all the previous cases.
We studied 1-D electron motion under two periodic forces and including the novel radiation reaction force. When the periodic external force has a simple time dependence (dynamical system (4)), we found that the difference increases with . For the periodic force depending on position and time (dynamical system (5)), the relation between and is more complicated but the same dependence is observed after some time ( ) when the difference becomes more stable. On the other hand and for both dynamical systems, the relation between the frequency of the electric field and is different since several peaks and discontinuities appear. The observed value for for both cases is of the order of 10−31 to 10−14 which could be a concern to detect experimentally, but we think it is worthy to try to make an experiment to see whether or not this approach for radiation reaction force points to the right experimental direction.
 Dirac, P.A.M. (1938) Classical Theory of Radiating Electrons. Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 167, 148-169.
 Jackson, J.D. (2007) Comment on “Preacceleration without Radiation: The Nonexistence of Preradiation Phenomenon,” by J.A. Heras. American Journal of Physics, 75, 844-845.
 Hnizdo, V. (2007) Comment on “Preacceleration without Radiation: The Nonexistence of Preradiation Phenomenon,” by J.A. Heras [Am. J. Phys. 74(11), 1025-1030 (2006)]. American Journal of Physics, 75, 845-846.