• Abstract and Figures
• Public Full-text

Generation of RF by Atmospheric Filaments Travis Garrett,1 Jennifer Elle,1 Michael White,1 Remington Reid,1 Alexander Englesbe,2 Ryan Phillips,1 Peter Mardahl,1 Erin Thornton,1 James Wymer,1 Anna Janicek,3 Oliver Sale,3 and Andreas Schmitt-Sody1 1Air Force Research Laboratory, Directed Energy Directorate, Albuquerque, NM 87123, USA 2Naval Research Laboratory, Plasma Physics Division, Washington, DC 20375, USA 3Leidos Innovations Center, Albuquerque, NM 87106, USA Recent experiments have shown that femtosecond filamentation plasmas generate ultra-broadband RF radiation. We show that a novel physical process is responsible for the RF: a plasma wake field develops behind the laser pulse, and this wake excites (and copropagates with) a surface wave on the plasma column. The surface wave proceeds to detach from the end of the plasma and propagates forward as the RF pulse. We have developed a four stage model of these plasma wake surface waves and find that it accurately predicts the RF from a wide range of experiments, including both 800 nm and 3.9 µm laser systems. The development of chirped pulse amplification [1] has tion. They revealed that a hot shell of electrons expands enabled an array of exotic physics [2{4], including atmo- off the plasma surface into the surrounding atmosphere spheric filamentation [5]. Filamentation in turn produces over roughly 50ps (as anticipated by [26]), and that the interesting and useful effects, from supercontinuum gen- corresponding current density Jr depends strongly on the eration [6,7], broadband THz radiation [8,9], and RF electron-neutral collision frequency. Subsequent simula- pulses, which have recently been explored in great detail tions demonstrated that this Plasma Wake field (PW) [10{15]. We have discovered that a novel combination excites a Surface Wave (SW) on the outer boundary of of a plasma wake field and a surface wave is responsible the plasma (namely, a long wavelength Surface Plasmon for this filament RF. Femtosecond lasers are becoming Polariton (SPP)[27, 28], with a speed approaching c). important spectroscopic tools (both in the lab [16] and Finally the SW was found to effectively detach from the for remote sensing [17]), and the added production of end of the plasma, with most of the energy being then pulsed RF promises to expand the utility of these sys- converted into a broadband RF pulse. tems [18{20], and to provide insight into the physics of The PWSW numerical model is split into 4 stages, as filamentation. length scales spanning 7 orders of magnitude need to be During typical λ = 800 nm filamentation Kerr self- resolved. First the initial electron distribution function focusing and plasma defocusing are in rough balance, re- Ne is determined. In a typical λ = 800 nm experiment 22 −3 sulting a plasma column with a density of ne ∼ 10 m a linearly polarized, single color, 40 mJ pulse is focused [21, 22]. This electron density corresponds to a plasma into a quartz cell with a f=60 mirror, which sets the 12 frequency of !pl ' 6 × 10 rad/s, which is closely linked beam waist !0 to 30µm (10 mJ and f=15 focusing are to the generation of THz radiation. In the single color common for λ = 3:9µm). The resulting Keldysh parame- p THz theory of [8] the freshly ionized electrons receive a ter γ = !laser 2meIpot=(qeE) is close to 1, so we use the ponderomotive push from the laser pulse, thus exciting general ionization rate W from [29], which is supported by a coherent longitudinal current Iz in the plasma (which experimental data [30] near atmospheric pressure (other is quickly damped by the high collision frequency). This ionization mechanisms will be necessary for 10 µm lasers short current pulse translates at c behind the laser pulse, [31]). A laser intensity of I ∼ 5 × 1017 W/m2 generates thereby producing a conical shell of radially polarized enough plasma at 1 atm for intensity clamping [21] to 12 −1 THz radiation. occur (with W ∼ 10 s ), which leads to an estimated The RF resembles the THz radiation in some ways, plasma length Lpl of 25 cm, and a radius rpl of ∼ 1 mm. and differs critically in others. It has a similar conical At 1 Torr I increases to ∼ 2 × 1018 W/m2, which ionizes shell spatial profile and radial polarization, which indi- much of the O2, and the plasma is reduced to Lpl ' 13 arXiv:2102.12657v2 [physics.plasm-ph] 9 Sep 2021 cates that a GHz scale longitudinal current pulse is also cm and rpl ' 0:5 mm. In turn, the lower energy 3.9 µm generated within the plasma. However it also has a well laser plasmas are shorter, with Lpl ranging from about 1 defined broadband peak in the 5-15 GHz range (and is to 2 cm. not the tail of the THz radiation [11, 13, 14]) and the These estimates match well (within 25%) with pho- RF amplitude grows strongly with decreasing pressure tographic measurements of the plasma taken at differ- [11, 13], in contrast to the THz which falls off sharply ent pressures (see the supplementary materials). Simu- [23]. A distinct physical mechanism thus must be re- lations of filamentation with these parameters have also sponsible for the RF. been performed using a Unidirectional Pulse Propaga- Exploratory Particle In Cell (PIC) simulations [24, 25] tion Equations (UPPE) solver [32]. These agree with the provided the key insights that explain the RF genera- overall dimensions of the plasma column, but also indi- 2 cate that there is rich fine scale structure throughout the 100 Torr 10 Torr 1 Torr 0.70 1020 plasma (as expected as we are in the multifilamentary 3 0.65 19 m 10 regime [33, 34]). For this letter we simplify and do not , y t i include fine scale plasma perturbations. The plasma col- 0.60 1018 s n e D umn is approximated as having a constant density central 0.55 1017 n o r Radius, mm core due to intensity clamping, and then a 100µm radial t 16 0.50 10 c e l ramp down to zero density at the outer radius rpl. E 15 The velocity distribution is determined by our PIC 0.45 10 0 25 50 75 0 25 50 75 0 25 50 75 100 code: given E(~x;t) the air is ionized at the W rate [29] Time, ps Time, ps Time, ps and the new electrons are accelerated through the re- mainder of the pulse [35]. In general the velocity distri- FIG. 1. Log plots of the electron density ne(r; t) showing the bution is influenced by both the details of strong field plasma wake variation at pressures of 100, 10, and 1 Torr (for ionization near γ = 1 (e.g. [36]), and by the deformation a λ = 800 nm). The plasma outer radius has been set to of the laser pulse during filamentation. For this work we rpl = 0:5 mm for all 3 pressures for comparison purposes. further simplify and assume that the electrons are ion- 19 3 ized with zero initial velocity and are then accelerated Radial current density, ne = 10 m 0.70 by the remainder of a Gaussian pulse (withx ^ polariza- 0.65 106 tion and propagating in +z). As a whole the initial Ne 0.60 105 is highly non-Maxwellian, and at 100 Torr has a peak 0.55 2 4 m 10 / A kinetic energy K ' 5eV, and an average kinetic en- 0.50 tail 3 , Radius, mm 10 y t 0.45 i ergy Kavg ' 0:6 eV, while these increase at 1 Torr to 2 s 10 n 0.40 e K ' 16eV and K ' 2 eV. For 3.9 µm lasers the d tail avg 20 3 0 t Radial current density, ne = 10 m n kinetic energy is ∼ 25 times larger, as the laser intensities 0.70 e 2 r 10 r 2 u 0.65 are comparable and the energy scales as λ . 3 c 10 l a 0.60 i We next consider the evolution of the plasma column. 104 d a 0.55 R We construct thin transverse slices of plasma given Ne, 105 0.50 6 with periodic boundary conditions used in the longitu- Radius, mm 10 dinal directionz ^ (which is valid to leading order as the 0.45 0.40 electron velocities are a small fraction of c) and simu- 0 20 40 60 80 100 late the radial evolution with our PIC code. The Debye Time, ps length is quite small: λDebye ' 10 nm, so we use an energy conserving method [37] to calculate the Lorentz FIG. 2. Log plots of the radial current density Jr(r; t) at 1 Torr (800 nm source), as it evolves in a transverse slice of force. The electron-neutral elastic collision frequency νeN plasma. The initial density at the edge of the plasma has depends on the cross sections for O2 and N2, which for 19 2 been given a step function profile, with a value of ne = 10 our energies is roughly 10 A˚ . In turn the electron-ion −3 20 −3 m in the top plot and ne = 10 m on the bottom. The momentum-transfer collision frequency is given by νei = profiles of the plasma wakes outside the column (r > rpl = 0:5 −12 3=2 3 7:7×10 ne ln(ΛC )=KeV , where ΛC = 6πneλDebye [38]. mm) are quite similar, while at the surface SPPs are excited The resulting radial current density Jr and electron den- at roughly ∼20 and ∼60 GHz respectively, along with higher frequency volume plasmons in the interior.

Load more

  6

Copy Link