Optimizing electron trajectories in combined helical wiggler and solenoidal magnetic fields for effective particle acceleration
محورهای موضوعی : Journal of Theoretical and Applied Physics
1 - Plasma Science and Technology Laboratory, Department of Physics, Indian Institute of Technology Delhi, New Delhi, India
کلید واژه: Solenoidal magnetic field, Pitch, Helical Wiggler magnetic field, Electron trajectories, Larmor radius,
چکیده مقاله :
Charged particle acceleration is the subject of great interest because of its applications in various fields such as thermonuclear fusion, nuclear physics, radiation generation, coherent harmonic generation, probing materials, medical science, food preservation, etc. The particle acceleration is usually done by interacting the particles with strong electric fields. Since the magnetic field diverts the particle motion, this also plays a vital role in the particle interaction with the electromagnetic fields. In the present work, a combined configuration of helical wiggler and solenoidal magnetic fields has been used to optimize the trajectory of the electron for effective particle acceleration. In this concept, the solenoidal field controls the transverse components of the electron velocities and wiggler field confines the helical motion of the electron. The optimized values of magnitudes of solenoidal field and wiggler magnetic field and its period / wavelength make this configuration useful for particle acceleration in waveguide.
Optimizing Electron Trajectories in Combined Helical Wiggler and Solenoidal Magnetic Fields for Effective Particle Acceleration
Hitendra K. Malik*
Plasma Science and Technology Laboratory, Department of Physics, Indian Institute of Technology Delhi, New Delhi – 110 016, India
*Email: hkmalik@physics.iitd.ac.in and hkmalik@hotmail.com
Abstract
Charged particle acceleration is the subject of great interest because of its applications in various fields such as thermonuclear fusion, nuclear physics, radiation generation, coherent harmonic generation, probing materials, medical science, food preservation, etc. The particle acceleration is usually done by interacting the particles with strong electric fields. Since the magnetic field diverts the particle motion, this also plays a vital role in the particle interaction with the electromagnetic fields. In the present work, a combined configuration of helical wiggler and solenoidal magnetic fields has been used to optimize the trajectory of the electron for effective particle acceleration. In this concept, the solenoidal field controls the transverse components of the electron velocities and wiggler field confines the helical motion of the electron. The optimized values of magnitudes of solenoidal field and wiggler magnetic field and its period / wavelength make this configuration useful for particle acceleration in waveguide.
Keywords
Solenoidal magnetic field, Helical Wiggler magnetic field, electron trajectories, Larmor radius, Pitch
1. Introduction
Several schemes have been developed in laboratories for charged particle acceleration and to achieve a high-quality beam of such particles [1 - 9]. The particles have been found to accelerate by the shocks [10]; even the acceleration of galactic cosmic rays is accounted by the shocks emanating from supernova explosions and propagating through the interstellar medium [11]. In laboratory / computer simulations also, the use of shocks has been made to accelerate the particles [10]. The recent acceleration techniques involve photonic chip-based particle acceleration [12], relativistic reconnection [13], bubble acceleration [14], etc. In most of the schemes, this can be seen that the laser beams have been largely used for the purpose of particle acceleration and very sophisticated equipment are required. However, the microwave pulses with special shapes have proved their effectiveness while being used in the plasma-filled waveguides [8, 15, 16]. In such mechanisms, the magnetic field has been employed for the extraction of accelerated particles [17]. Not only this, but the magnetic field is also found to control the radiation emission and to enhance the efficiency of the schemes [18 - 20].
In the mechanisms of particle acceleration and the radiation generation, the electron motion or its trajectory needs to be controlled as it plays a vital role for the effective interaction of the electron with the electromagnetic field. The calculations for the electron trajectories in free electron laser (FEL) with helical wiggler and axial guide fields has been done by several researchers [21 - 24]. For example, Aamodt [21] has calculated the expressions for the electron orbits analytically in the presence of a longitudinal wiggler magnetic field. Freund and Ganguly [22] have developed a theory for non-steady trajectories in three-dimensional helical wiggler magnetic field. Fajans et al. [23] have analyzed off axis electron orbits in FEL beam in the presence of combination of helical wiggler and axial guide field. El-Bahi et al. [24] have described the electron trajectories in helical free electron laser. On the other hand, Punia and Malik [25] have used the combined configuration of helical wiggler and solenoidal magnetic fields to generate polarization-tunable terahertz radiation. Recently, three-coil set ups with rectangular cross section coils [26] or tapered coils [27] have been proposed to control the magnetic field topology in a magnetic nozzle, which is a device for space propulsion. Since the separation of electrons from the ions in these devices take place due to the magnetic field distribution, the electron motion here played an important role in creating the ambipolar field and hence, in controlling the thrust. In the present letter, considering the combined configuration of helical wiggler and solenoidal magnetic fields, the expressions for the electron velocity components have been derived analytically and the relevant equations are solved numerically in MATLAB using fourth-order Runge-Kutta method (ODE45) for investigating the electron trajectories under different injection positions of the electron. The optimization of the same is expected to produce effective electron acceleration when it is subjected to the electromagnetic field of the microwaves.
2. Mathematical equations
A combined configuration of solenoidal and wiggler magnetic fields is considered in the present work, with the following expression
(1)
Here is the wave number corresponding to the wavelength of the wiggler magnetic field, given as . is the peak value of the helical wiggler field and is the magnitude of the solenoidal field. The solenoidal field has been applied in the z direction, so it will provide a circular motion to the electron and control the transverse (x and y) components of its velocity. The wiggler field is circularly polarised in nature and it will rotate around z axis as the advancement in z results in a helical magnetic field (Figure 1). So, this field is expected to provide helical motion to the electron with higher acceleration if the electron is subjected to external electric field. It can be seen from Figure 1 that the diameter of the circle is enhanced when the peak value of the wiggler field is increased from 0.4 T to 0.5 T. On the other hand, the periodicity or the wavelength of the field decreases with the increasing values of (please see Figure 2, where the graph marked with red stars corresponds to and the graph with blue dots corresponds to . The addition of the solenoidal field increases the peak value of the combined magnetic field, as shown in Figure 3.
Figure 1: Profile of the helical wiggler magnetic field for different values of , when .
Figure 2: Profile of the helical wiggler magnetic field for different values of ( for the graph marked with red stars and for the graph marked with blue dots), when .
Figure 3: Profile of the combined magnetic field having solenoidal () and helical wiggler () magnetic fields.
The equation of motion for the electron under the said configuration of the magnetic field is written as
(2)
Here is the velocity of the electron, is the electronic charge and is the mass of the electron. Equations (1) and (2) result in the following set of three equations:
(3)
(4)
(5)
These equations have been solved by taking three different sets of initial conditions to get the velocity components in each case. Assuming that the electron enters in the magnetic field (at t = 0) with the initial velocity in x-y plane from the origin and , we get
(7)
(8)
For the case the electron enters in the magnetic field (at t = 0) with initial velocity in y-z plane from the origin and , we get
(9)
(10)
(11)
When the electron enters in the magnetic field (at t = 0) with initial velocity in z-x plane from the origin and, we get
(12)
(13)
(14)
3. Results and discussion
The effect of magnitudes of the solenoidal and wiggler magnetic fields is analysed in greater detail in this section by taking their different values and then comparing the electron trajectories in both the cases. This comparison enables us to get the controlling parameters for the trajectory. This effect has been analysed by considering three sets of different initial conditions.
3.1. Electron injection in x-y plane
The particle enters the interaction region initially having velocity 0.2c in the x- and y- directions (x-y plane). The electron trajectory obtained numerically [28 - 31] based on the equations (6) – (8) is shown in Figure 4. Since both the velocity of electron and helical wiggler field are in the x-y plane in this case but the solenoidal field is in the z-direction, the solenoidal field causes the electron to move in a circular path only. On the other hand, the wiggler field’s y-component alters the x-component of the velocity and produces sinusoidal oscillations in the x-z plane. The x-component of the field alters the y-component of the velocity and produces sinusoidal oscillations in the y-z plane. The sinusoidal oscillations are shown in the inset of the figure. Here the helical nature of the field also causes the electron to move in the helical path.
Figure 4: Trajectory of electron when and .
Figure 5: Effect of solenoidal field on trajectory whenand .
Figure 6: Effect of wiggler magnetic field on trajectory when and .
If we look at the role of the solenoidal field (), the radius of the trajectory is found to decrease with the increase in the strength of this field (Figure 5). This can be understood based on the relation of the Larmor radius and the magnetic field, i.e. , where is the perpendicular component of the velocity of the electron. However, the enhancement of the magnitude of the wiggler magnetic field () causes no change in the radius of the trajectory, as the radius in this case is independent of . Hence, both the graphs in Figure 6 overlap each other for the two values of as 0.3T and 0.4T.
Figure 7: Effect of solenoidal field on the trajectory when and .
3.2. Electron injection in y-z plane
Now we discuss the case where the electron enters the interaction region initially having velocity 0.2c in the y- and z-directions. For this, we need to solve the set of equations (9) – (11) to find the x-, y-, and z-coordinates for evaluating the electron trajectory. The results are shown in Figure 7. Since the electron is injected in the y-z plane having velocity components in the x- and z-directions, the solenoidal field (due to Lorentz force) causes the helical motion. As the wiggler field and velocity components are lying in different directions, the wiggler field also contributes to the radius of the electron motion trajectory in this case as well as in the helical motion. The Larmor radius is found to decrease with the increasing strength of the solenoidal field (Figure 7). While this time radius is dependent on the wiggler magnetic field too, so on changing the radius of the motion also changes. Please see this result in Figure 8. An increase in the magnitude of wiggler field increases the radius of the helix of helical wiggler field due to which the radius of trajectory also increases. This is due to the combined effect of the two fields.
Figure 8: Effect of wiggler field on trajectory when and .
3.3. Electron injection in x-z plane
Finally, we consider the situation where the electron enters the interaction region initially having velocity 0.2c in the x- and z-directions, i.e., the x-z plane. For this case, we need to solve the set of equations (12) – (14). The resultant electron trajectory is shown in Figure 9. Due to and components, the solenoidal field (applied in the z-direction) contributes to the helical motion of the electron. Similar to the previous case, the wiggler field and velocity are in different directions (planes), so the wiggler field contributes to the radius of the motion as well as to the helical motion (circular polarisation of the wiggler field). Under this situation, larger the solenoidal field is, smaller is the radius of the trajectory helix (please see Figure 9). On the other hand, the radius of the helical path is also found to change with the change in the wiggler field ; this is evident from Figure 10.
Figure 9: Effect of solenoidal field on trajectory when and .
Figure 10: Effect of wiggler field on trajectory when and .
Finally, we discuss the situation of relativistic speed of the electron when it is launched in the said configuration of magnetic field. The present calculations were carried out by neglecting the relativistic effects. However, there will be a variation in the electron mass when it propagates at the relativistic speed that will modify the cyclotron frequency. Also, the relativistic factor is expected to appear in the calculations and it will become more complicated to separate out all the velocity components in the methodology we developed under the effect of combined solenoidal and helical-wiggler magnetic fields. With regard to achieving the said magnetic field, we mention that a transverse periodic helical magnetic field can be produced on the axis of a double-helix-wound bifilar magnet with equal and opposite currents in each helix [32]. Helical magnetic field has been produced during preheating at the electroweak scale [33]. When the transverse and constant magnetic field rotates in helical fashion the field is called helical-wiggler magnetic field. Wiggler magnetic field with exponential entrance has been produced using the permanent magnets and aluminium holders [34].
4. Conclusions
Different cases were discussed for the electron trajectory under the impact of combined solenoidal magnetic field applied along the z-direction and wiggler magnetic field applied in the x-y plane. A comparison of the electron trajectories under various initial conditions and the magnitudes of both kinds of the magnetic fields enables one to optimize the parameters to control the motion of the electron. The radius of trajectory decreases with increasing strength of the solenoidal field, controlling the transverse components of the electron velocity. For its magnitude of 0.3 T and the electron initial velocity as 6×106 m/s, the radius of motion is of the order of 10-4 m, while the pitch is of the order of 10-3 m. At the entrance of the wiggler, some fraction of axial velocity (z-component) is found to be converted into perpendicular velocity by the wiggler, resulting in the helical trajectory of the electron. The interaction time between the electron and the field is increased by using the helical wiggler magnetic field, which is expected to help providing the electron more intense acceleration in the presence of external electric field. The configuration proposed in this article can be used in a rectangular waveguide for electron acceleration using microwaves, where the radius of circular motion is controlled by the solenoidal field and the helical motion is confined by the wiggler field. This is in view of the fact that the magnetic fields with periodicity / wavelength in centimeters or fraction in centimeters can be easily produced and the wavelength of the microwaves also lies in centimeters. Also by filling this waveguide with plasma, one can additionally control the cutoff frequency and interaction length for achieving the resonance and effective particle acceleration.
References
1. Malik, H.K. Laser-Matter Interaction for Radiation and Energy. CRC Press, 1st edition, 2021.
2. Malik, H.K., Kumar, S. and Singh, K.P., 2008. Electron acceleration in a rectangular waveguide filled with unmagnetized inhomogeneous cold plasma. Laser and Particle Beams, 26(2), pp.197-205.
3. Kumar, S. and Malik, H.K., 2006. Electron acceleration in a plasma filled rectangular waveguide under obliquely applied magnetic field. Journal of plasma physics, 72(6), pp.983-987.
4. Malik, H.K., Kumar, S., Dhaka, V. and Singh, D.K., 2021. Analytical approach for the use of different gauges in bubble wakefield acceleration. Journal of Theoretical and Applied Physics, 15(1), pp.1-12.
5. Malik, H.K., 2004. Effect of plasma density on proton acceleration in a rectangular waveguide. Plasma Science and Technology, 6(5), p.2456.
6. Kumar, S., Gupta, D.N., Malik, H.K., Singh, D., Lee, J. and Yoon, M., 2020. Effect of a tightly focused chirped Gaussian laser pulse on electron acceleration in helical undulator. Physics of Plasmas, 27(4), p.043105.
7. Yadav, M., Gupta, D.N. and Sharma, S.C., 2020. Electron plasma wave excitation by a q-Gaussian laser beam and subsequent electron acceleration. Physics of Plasmas, 27(9), p.093106.
8. Malik, H.K., 2008. Analytical calculations of wake field generated by microwave pulses in a plasma filled waveguide for electron acceleration. Journal of applied physics, 104(5), p.053308.
9. Comisso, L. and Sironi, L., 2021. Pitch-Angle Anisotropy Controls Particle Acceleration and Cooling in Radiative Relativistic Plasma Turbulence. Physical Review Letters, 127(25), p.255102.
10. Decker, R.B., 1988. Computer modeling of test particle acceleration at oblique shocks. Space Science Reviews, 48(3), pp.195-262.
11. Blandford, R.D. and Ostriker, J.P., 1978. Particle acceleration by astrophysical shocks. The Astrophysical Journal, 221, pp.L29-L32.
12. Shiloh, R., Illmer, J., Chlouba, T., Yousefi, P., Schönenberger, N., Niedermayer, U., Mittelbach, A. and Hommelhoff, P., 2021. Electron phase-space control in photonic chip-based particle acceleration. Nature, 597(7877), pp.498-502.
13. Zhang, H., Sironi, L. and Giannios, D., 2021. Fast particle acceleration in three-dimensional relativistic reconnection. The Astrophysical Journal, 922(2), p.261.
14. Ren, Z., Zuo, Z., Wu, S. and Liu, S., 2022. Particulate Projectiles Driven by Cavitation Bubbles. Physical Review Letters, 128(4), p.044501.
15. Malik, H.K., 2003. Energy gain by an electron in the fundamental mode of a rectangular waveguide by microwave radiation. Journal of plasma physics, 69(1), pp.59-67.
16. Aria, A.K., Malik, H.K. and Singh, K.P., 2009. Excitation of wakefield in a rectangular waveguide: Comparative study with different microwave pulses. Laser and Particle Beams, 27(1), pp.41-47.
17. Malik, H.K., 2002. Energy gain by an electron in the fundamental mode of a rectangular waveguide by microwave radiation. Journal of plasma physics, 68(3), pp.211-219.
18. Malik, H.K. and Punia, S., 2019. Dark hollow beams originating terahertz radiation in corrugated plasma under magnetic field. Physics of Plasmas, 26(6), p.063102.
19. Malik, H.K., 2015. Terahertz radiation generation by lasers with remarkable efficiency in electron–positron plasma. Physics Letters A, 379(43-44), pp.2826-2829.
20. Bernáth, B., Gogoi, P., Marchese, A., Kamenskyi, D., Engelkamp, H., Arslanov, D., Redlich, B., Christianen, P.C. and Maan, J.C., 2022. Nonlinear terahertz transmission spectroscopy on Ga-doped germanium in high magnetic fields. Physical Review B, 105(20), p.205204.
21. Aamodt, R.E., 1983. Electron orbits in a free-electron laser with a longitudinal magnetic wiggler. Physical Review A, 28(5), p.2895.
22. Freund, H. and Ganguly, A., 1985. Electron orbits in free-electron lasers with helical wiggler and axial guide magnetic fields. IEEE journal of quantum electronics, 21(7), pp.1073-1079.
23. Fajans, J., Kirkpatrick, D.A. and Bekefi, G., 1985. Off-axis electron orbits in realistic helical wigglers for free-electron-laser applications. Physical Review A, 32(6), p.3448.
24. El-Bahi, R., Rhimi, M.N. and Cheikhrouhou, A.W., 2002. Electron trajectories in a helical free-electron laser. Brazilian Journal of Physics, 32, pp.790-797.
25. Punia, S. and Malik, H.K., 2019. Controlling polarization of THz radiation in pair plasma. Plasma Sources Science and Technology, 28(11), p.115018.
26. Malik, L., Kumar, M. and Singh, I.V., 2021. A three-coil setup for controlled divergence in magnetic nozzle. IEEE Transactions on Plasma Science, 49(7), pp.2227-2237.
27. Malik, L., 2022. Tapered coils system for space propulsion with enhanced thrust: a concept of plasma detachment. Propulsion and Power Research, 11(2), pp.171-180.
28. Malik, L. and Tevatia, A., 2021. Comparative Analysis of Aerodynamic Characteristics of F16 and F22 Combat Aircraft using Computational Fluid Dynamics. Defence Science Journal, 71(2).
29. Malik, L., Rawat, S., Kumar, M. and Tevatia, A., 2021. Simulation studies on aerodynamic features of Eurofighter Typhoon and Dassault Rafale combat aircraft. Materials Today: Proceedings, 38, pp.191-197.
30. Malik, L., Escarguel, A., Kumar, M., Tevatia, A. and Sirohi, R.S., 2021. Uncovering the remarkable contribution of lasers peak intensity region in holography. Laser Physics Letters, 18(8), p.086003.
31. Malik, L., 2020. Dark hollow lasers may be better candidates for holography. Optics & Laser Technology, 132, p.106485.
32. Kincaid, B.M., 1977. A short‐period helical wiggler as an improved source of synchrotron radiation. Journal of Applied Physics, 48(7), pp.2684-2691.
33. Diaz-Gil, A., Garcia-Bellido, J., Perez, M.G. and Gonzalez-Arroyo, A., 2008. Magnetic field production during preheating at the electroweak scale. Physical Review Letters, 100(24), p.241301.
34. Al-Shamma'a, A., Stuart, R.A. and Lucas, J., 1996. A wiggler magnet for a CW-FEM. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 375(1-3), pp.424-426.