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 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. sity ne are then recorded as functions of radius and time. These are well resolved with a mesh resolution of ∆x = ∆y = 2µm, and a macro-particle weight of ∼10 sity at a small distance r off the surface of the plasma on the outer edge of the plasma. ∆ at radius r . Using Gauss’s law and approximating The electron number density for λ = 800 nm simula- pl ln((r + r )/r ) as r /r for r  r we find the that tions of the PW at 100, 10, and 1 Torr can be seen in pl ∆ pl ∆ pl ∆ pl the number of electrons per unit length n that escape Fig.1. The outer edge of the plasma at t = 0 has a L 20 −3 to r∆ scales as: simplified step function profile, with ne = 10 m at radius rpl = 0.5 mm. Thus, in addition to the PW be- 2π0KeV rpl ing launched off of the edge of this plasma, a coherent nL(r∆) ' , (1) qe r∆ radial plasma frequency oscillation of ∼90GHz is stimu- lated in the interior and a SPP at ∼63GHz is excited on where KeV ' Ktail is the kinetic energy of the hot the surface. The PW that expands out into the neutral electrons that lead the escape. The radial electrostatic atmosphere (r > rpl) is insensitive to the density of the field Er,stat that corresponds to nL has a simple form: outer edge of the plasma (see Fig.2), as opposed to the Er,stat(r∆) = KeV /r∆. At λ = 800 nm and 1 Torr 5 excited surface and volume plasmons. Ktail ' 16 eV and r∆ ' 60µm, giving Er,stat ' 2.7 × 10 Analytic approximations provide a useful complement V/m. This allows for an estimate of the PW evolution to the simulations. Consider the late time electron den- timescale: an electron with energy Ktail = 16 eV in an 3

3.9 micron, 10 Torr Surface Wave electric field of this magnitude follows a parabolic tra- 3 105 jectory over 100 ps. This compares well with the spread of electrons seen in Fig.1. The magnitude of the radial current density can likewise be approximated as: 2 104 ( 3 qenL vtail Radius, cm 1 10 Jr ' veff (t) ; veff (t) = min (2) 2πrplr∆ vdiff (t) Electric Field Strength, V/m 0 102 where the effective velocity veff (t) is the minimum of the 2 1 0 1 2 3 4 Length, cm tail velocity and the Fickian diffusion speed vdiff (t) = (K /(m νt))1/2. At 1 Torr this gives J ' 1 × 105 tail e r FIG. 3. Log plot of electric field strength shortly after the A/m2. Finally we note that caution is needed for small detachment of a λ = 3.9 µm, 10 Torr surface wave from a 2 r∆ values: Eq.s (1) and (2) only hold up to densities that cm plasma column, as simulated by the axisymmetric FDTD- are comparable to the original plasma edge density. Drude model. We next switch from these transverse PIC simulations to a continuum FDTD-Drude model in a 2D axisymmet- ric coordinate system [39] that spans the length of the typical plasma conductivity σ a complex ξ can be defined: ~ ~ 2 ~ plasma. The Drude model: ∂tJpl + νplJpl = 0ωplE −2γ r ke euler 0ω is integrated into the code via an auxiliary differential ξ = − rpl(1 + i), (4) equation [40]. A mixed resolution of ∆r = 2µm and 2 2σ ∆z = 50µm suffices to resolve the SWs. The PW cur- (where γeuler is the Euler-Mascheroni constant), and h is rent profile Jr is driven across the surface of the plasma then given by: column in the longitudinal direction at the speed of light. In this work only the external radial currents J are s r>rpl 4ξe−2γeuler used to drive the SW and subsequent RF. The higher fre- h = k2 + , (5) r2 W (ξ) quency volume plasmons can also excite waves in neigh- pl 1 boring lower density plasma, but we suspect that includ- where the first branch of W is used. At 1 Torr the plasma ing fine scale plasma structure will act as surface rough- has a radius r ' 0.5 mm and conductivity σ ' 2 × 103 ness [41] and cause these to continuously detach. pl S/m, and thus for a 10 GHz SW we find h ' 210 + 3.2i. As expected the PW current excites a broadband The dispersive SW thus travels at about 0.988c, atten- surface wave with similar frequency content. The fre- uates by a factor of 1/e every 30 cm (if not being ac- quency ω of the wave is considerably lower than the tively excited), and has a 1/r transverse envelope out to limit SPP frequency and it thus travels at velocity r ' 2.4 cm, all of which compares well with simula- c((2ω2−ω2 )/(ω2−ω2 ))1/2 which approaches the speed of outer pl pl tions. light. The resulting phase-matched copropagation leads We next consider the excitation of the surface wave by to a steady growth in the surface wave intensity (see also the plasma wake currents. An isolated Hertzian dipole [42, 43]). with uniform current I along an of length L will The SW is well approximated by the Sommerfeld- radiate electromagnetic waves at a wavelength λ Goubau [44] solution for SPPs on a cylinder of finite con- Dipole with power: P = (πη /3)I2L2/λ2 . In our ductivity. The external radial component of the electric Dipole 0 Dipole case we have a radial current distribution of approximate field Er has the form: length r∆ that is normal to a conducting cylinder of ra-

πrplE0 (1) p dius rpl. When there is a separation of length scales E (r, z, t) = − ei(ωt−hz)H (r k2 − h2), (3) r 1 r∆  rpl  λSW and an overall radial current of Ir then 2router the power that is primarily radiated into SWs is: for waves of frequency ω propagating in the +z direction, πη r2 with an amplitude of E0 at the surface. The complex h 0 2 ∆ PSW ' I , (6) wavenumber describes the SW wavelength and attenua- 60 r 5/3 1/3 rpl λSW tion length scale, k is the free space wavenumber k = ω/c, (1) (where we have simplified with λ = 2π/ e(h)). We H1 is a Hankel√ function of the first kind with order 1, SW R 2 2 2 5/3 1/3 and router = 1/| k − h |. For radii rpl < r < router note that the dimensionless term r∆/(rpl λSW ) (com- the Hankel function with complex argument is approxi- 2 2 pare to L /λDipole for the Hertzian dipole), has different (1) mately H1 ' −2router/(πr), and for larger radii it falls asymptotic scalings in other length scale regimes. Equa- off exponentially (as is typical for SPPs). tion (6) applies both to a stationary radial antenna and Recent work [45, 46] has simplified the calculation of the ±z wavetrains it emits, and to our case, where the h through the use of the Lambert W function. Given a antenna translates at c with its comoving SW pulse. 4

The total effective current Ir(z) is found by integrating = 800 nm, 40mJ, fN = 60 = 3.9 m, 12mJ, fN = 15 3.5 the current density Jr from 0 to z, and absent dissipation 80 3.0 the amplitude E0 of the SW in (3) would grow linearly 60 2.5 (for constant Jr). However there is dissipation due to 2.0 the finite conductivity of the plasma, and the attenuation 40 1.5 1.0 coefficient α = m(h) is similar in magnitude to 1/L . 20 I pl Angle, Degrees The SW thus grows linearly initially and then saturates 0.5 0 0.0 E Field Amplitude, V/m after several 1/α, and the effective total radial current 0 10 20 30 40 0 10 20 30 40 50 after propagating a distance z is: Frequency, GHz Frequency, GHz

1 − e−αz FIG. 4. Plots of simulated electric field strength as a function Ir(z) = 2πrplJr . (7) α of frequency and angle for RF pulses from a 40 mJ, fN = 60, 800 nm laser source [11] and a 12 mJ, fN = 15, 3.9 µm We express the amplitude of the final SW in terms of its source [13] (both run at 10 Torr). The numerical profiles and longitudinal current, as Iz = 2πrplE0/η0. Integrating (3) amplitudes closely match the experimental results. to find the SW peak power and setting this equal to (6) we find: RF voltage, laboratory data RF voltage, simulated and Eqn. (9)

2 −αz  1/6 90° 1 Torr, Lab 90° 1 Torr, Sim. 2π r∆Jr(1 − e ) rpl 10 Torr, Lab 10 Torr, Sim. Iz(z) = . (8) 100 Torr, Lab 100 Torr, Sim. 1/2 60° 60° 1 Torr, Eqn. (9) α(30 ln(router/rpl)) λSW 10 Torr, Eqn. (9) 100 Torr, Eqn. (9) (deg.) (deg.) For the λ = 800 nm laser this gives Iz ' 0.6 A at the end of the plasma at 1 Torr, a current of I ' 0.2 A at z 30° 30° 10 Torr, and Iz ' 0.05 A at 100 Torr, which are close to the measured currents in the corresponding simulations. When the surface wave reaches the end of plasma the majority of the energy detaches and propagates forward as the RF pulse, with the reflected wave decaying as it 0° 0° travels back down the plasma. Fig.3 shows an example 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 V (V) V (V) of this, where the SW has just detached from a 2 cm pp pp plasma column produced by a λ = 3.9 µ laser source at 10 Torr. In this way the system resembles a surface wave FIG. 5. Polar plots of the detector peak to peak Voltage Vpp for the λ = 800 nm laser source as a function of angle end-fire antenna [47, 48]. The recorded Iz(z, t) currents 2 and pressure as measured in the lab (left), and as calculated from the SW simulation are used in a ∼1 m laboratory by both the simulation framework (right, triangles), and by scale axisymmetric FDTD model, with resolution ∆r = the associated analytic approximations (right, × ticks) that ∆z = 0.5 mm. Plots of the resulting Fourier transformed culminate in (9). RF field strengths at 10 Torr for both λ = 800 nm and λ = 3.9 µm simulations can be seen in Fig.4. To directly compare with lab results the simulated elec- bit lower than the lab data for the runs shown in Fig.5, tric fields are converted to peak to peak voltages Vpp via but is within the errors for both the simulations and in a frequency dependent calibration factor CE/V (fGHz) ' the experimental data. Subtle changes in Ti:Sapphire 50 + 10fGHz [11]. We adapt an equation for energy spec- laser alignment generates ±50% amplitude variations at tral density radiated per unit solid angle derived for THz a given pressure (with the spatial and frequency content radiation [8]: being much more stable). It is suspected that the slight laser pulse variations are amplified during filamentation, d2U |I (ω)|2 sin2 θ L ω  = z sin2 pl (1 − cos θ) , (9) leading to fine scale changes in the plasma radius, density, dωdΩ 4π c(1 − cos θ)2 2c 0 and electron velocities. Plausible variations in the same to construct voltage profiles for our analytic approxima- parameters lead to comparable changes in the simulated tions. The average amplitude of Iz (8) over the length RF amplitude as well. of the plasma is applied to a ∼10 GHz pulse, which is We have demonstrated that a plasma wake field that Fourier transformed for use with (9) (along with the cal- excites a surface wave is responsible for the generation of ibration factor), and then converted back to an effective filament RF. The rich internal dynamics of the plasma Vpp. merits further research: we expect that stochastic higher In general simulations with varying pressures, focal ge- frequency RF may inform on the fine scale structure of ometries, and laser pulse energies and wavelengths all the plasma and the nonlinear optics that generate it. The produce RF that nicely matches recent experimental re- ionization model needs to be expanded to include more sults [11–14]. The magnitude of the simulated RF is a physical effects (e.g. [31, 36, 49]), and strong field ioniza- 5 tion can be further explored with noble gas filament RF. [11] A. Englesbe, J. Elle, R. Reid, A. Lucero, H. Pohle, The plasma can also be driven into the strongly coupled M. Domonkos, S. Kalmykov, K. Krushelnick, and regime with harder focusing, and we intend to model re- A. Schmitt-Sody, Gas pressure dependence of sulting RF with the use of [50]. Finally we note that it pulses generated by laser-produced filament plasmas, Op- tics Letters 43, 4953 (2018). may be possible to recreate the core mechanism – an an- [12] A. Janicek, E. Thornton, T. Garrett, A. Englesbe, J. Elle, tenna that copropagates with the wave that it excites – and A. Schmitt-Sody, Length dependence on broadband in other settings and thereby drive strong coherent radi- microwave emission from laser-generated plasmas, IEEE ation. Transactions on Plasma Science 48, 1979 (2020). [13] A. Mitrofanov, D. Sidorov-Biryukov, M. Nazarov, A. Voronin, M. Rozhko, A. Fedotov, and A. Zheltikov, The authors thank Edward Ruden, Serge Kalmykov Coherently enhanced microwave pulses from midinfrared- and Charles Armstrong for useful discussions, and the Air driven laser plasmas, Optics Letters 46, 1081 (2021). Force Office of Scientific Research (AFOSR) for support [14] A. Englesbe, J. Elle, R. Schwartz, T. Garrett, D. Wood- via contracts #FA9550-18RDCOR0009 and #FA9550- bury, D. Jang, K.-Y. Kim, H. Milchberg, R. Reid, 19RDCOR025. Approved for public release; distribution A. Lucero, et al., Ultrabroadband microwave radiation from near-and mid-infrared laser-produced plasmas in is unlimited. Public Affairs release approval #AFRL- air, Physical Review A 104, 013107 (2021). 2020-0489. [15] A. Mitrofanov, A. Voronin, M. Rozhko, D. Sidorov- Biryukov, M. Nazarov, A. Fedotov, and A. Zheltikov, Polarization and spatial mode structure of mid-infrared- driven terahertz-to-microwave radiation, ACS Photonics (2021). [1] D. Strickland and G. Mourou, Compression of amplified [16] T. Kampfrath, K. Tanaka, and K. A. Nelson, Resonant chirped optical pulses, Optics communications 56, 219 and nonresonant control over matter and light by intense (1985). terahertz transients, Nature Photonics 7, 680 (2013). [2] J. Faure, Y. Glinec, A. Pukhov, S. Kiselev, S. Gordienko, [17] T. A. Labutin, V. N. Lednev, A. A. Ilyin, and E. Lefebvre, J.-P. Rousseau, F. Burgy, and V. Malka, A A. M. Popov, Femtosecond laser-induced breakdown laser–plasma accelerator producing monoenergetic elec- spectroscopy, Journal of Analytical Atomic Spectrome- tron beams, Nature 431, 541 (2004). try 31, 90 (2016). [3] X. Liu, D. Du, and G. Mourou, Laser ablation and mi- [18] G. B. Park and R. W. Field, Perspective: The first ten cromachining with ultrashort laser pulses, IEEE journal years of broadband chirped pulse fourier transform mi- of quantum electronics 33, 1706 (1997). crowave spectroscopy, The Journal of chemical physics [4] A. Di Piazza, C. M¨uller,K. Hatsagortsyan, and C. H. 144, 200901 (2016). Keitel, Extremely high-intensity laser interactions with [19] A. Doll, S. Pribitzer, R. Tschaggelar, and G. Jeschke, fundamental quantum systems, Reviews of Modern Adiabatic and fast passage ultra-wideband inversion in Physics 84, 1177 (2012). pulsed epr, Journal of magnetic resonance 230, 27 (2013). [5] A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and [20] J. Capmany and D. Novak, Microwave photonics com- G. Mourou, Self-channeling of high-peak-power femtosec- bines two worlds, Nature photonics 1, 319 (2007). ond laser pulses in air, Optics letters 20, 73 (1995). [21] A. Couairon and A. Mysyrowicz, Femtosecond filamen- [6] A. Brodeur and S. Chin, Ultrafast white-light continuum tation in transparent media, Physics reports 441, 47 generation and self-focusing in transparent condensed (2007). media, JOSA B 16, 637 (1999). [22] S. L. Chin, Femtosecond laser filamentation, Vol. 55 [7] V. Kandidov, O. Kosareva, I. Golubtsov, W. Liu, (Springer, 2010). A. Becker, N. Akozbek, C. Bowden, and S. Chin, Self- [23] G. Rodriguez and G. L. Dakovski, Scaling behavior of ul- transformation of a powerful femtosecond laser pulse into trafast two-color terahertz generation in plasma gas tar- a white-light laser pulse in bulk optical media (or su- gets: energy and pressure dependence, Optics express 18, percontinuum generation), Applied Physics B 77, 149 15130 (2010). (2003). [24] C. K. Birdsall and A. B. Langdon, Plasma physics via [8] C. Amico, A. Houard, S. Akturk, Y. Liu, J. Le Bloas, computer simulation (CRC press, 2004). M. Franco, B. Prade, A. Couairon, V. Tikhonchuk, and [25] R. E. Peterkin and J. W. Luginsland, A virtual proto- A. Mysyrowicz, Forward thz radiation emission by fem- typing environment for directed-energy concepts, Com- tosecond filamentation in gases: theory and experiment, puting in Science & Engineering 4, 42 (2002). New Journal of Physics 10, 013015 (2008). [26] B. Zhou, A. Houard, Y. Liu, B. Prade, A. Mysyrow- [9] V. Andreeva, O. Kosareva, N. Panov, D. Shipilo, icz, A. Couairon, P. Mora, C. Smeenk, L. Arissian, and P. Solyankin, M. Esaulkov, P. G. de Alaiza Mart´ınez, P. Corkum, Measurement and control of plasma oscil- A. Shkurinov, V. Makarov, L. Berg´e, et al., Ultrabroad lations in femtosecond filaments, Physical review letters terahertz spectrum generation from an air-based filament 106, 255002 (2011). plasma, Physical review letters 116, 063902 (2016). [27] J. Pitarke, V. Silkin, E. Chulkov, and P. Echenique, The- [10] B. Forestier, A. Houard, M. Durand, Y.-B. Andr´e, ory of surface plasmons and surface-plasmon polaritons, B. Prade, J.-Y. Dauvignac, F. Perret, C. Pichot, M. Pel- Reports on progress in physics 70, 1 (2006). let, and A. Mysyrowicz, Radiofrequency conical emission [28] S. A. Maier, S. R. Andrews, L. Martin-Moreno, and from femtosecond filaments in air, Applied Physics Let- F. Garcia-Vidal, Terahertz surface plasmon-polariton ters 96, 141111 (2010). propagation and focusing on periodically corrugated 6

metal wires, Physical Review Letters 97, 176805 (2006). [39] A. Taflove and S. C. Hagness, Computational elec- [29] S. Popruzhenko, V. Mur, V. Popov, and D. Bauer, Strong trodynamics: the finite-difference time-domain method field ionization rate for arbitrary laser frequencies, Phys- (Artech house, 2005). ical review letters 101, 193003 (2008). [40] M. Okoniewski, M. Mrozowski, and M. Stuchly, Simple [30] A. Sharma, M. N. Slipchenko, M. N. Shneider, X. Wang, treatment of multi-term dispersion in fdtd, IEEE Mi- K. A. Rahman, and A. Shashurin, Measurements of elec- crowave and Guided Wave Letters 7, 121 (1997). tron numbers in femtosecond laser induced plasmas using [41] E. Kretschmann, The angular dependence and the po- rayleigh microwave scattering, in 2018 AIAA Aerospace larisation of light emitted by surface plasmons on metals Sciences Meeting (2018) p. 0177. due to roughness, Optics Communications 5, 331 (1972). [31] S. Tochitsky, E. Welch, M. Polyanskiy, I. Pogorelsky, [42] S. J. Smith and E. Purcell, Visible light from localized P. Panagiotopoulos, M. Kolesik, E. M. Wright, S. W. surface charges moving across a grating, Physical Review Koch, J. V. Moloney, J. Pigeon, et al., Megafilament in 92, 1069 (1953). air formed by self-guided terawatt long-wavelength in- [43] R. Dahan, S. Nehemia, M. Shentcis, O. Reinhardt, frared laser, Nature Photonics 13, 41 (2019). Y. Adiv, X. Shi, O. Be’er, M. H. Lynch, Y. Kur- [32] M. Kolesik and J. V. Moloney, Nonlinear optical pulse man, K. Wang, et al., Resonant phase-matching between propagation simulation: From maxwell’s to unidirec- a light wave and a free-electron wavefunction, Nature tional equations, Physical Review E 70, 036604 (2004). Physics 16, 1123 (2020). [33] L. Berg´e,S. Skupin, F. Lederer, G. M´ejean, J. Yu, J. Kas- [44] G. Goubau, Surface waves and their application to trans- parian, E. Salmon, J. Wolf, M. Rodriguez, L. W¨oste, mission lines, Journal of Applied Physics 21, 1119 (1950). et al., Multiple filamentation of terawatt laser pulses in [45] D. Jaisson, Simple formula for the wave number of the air, Physical Review Letters 92, 225002 (2004). goubau line, Electromagnetics 34, 85 (2014). [34] W. Ettoumi, J. Kasparian, and J.-P. Wolf, Laser fila- [46] J. R. G. Mendon¸ca,Electromagnetic surface wave prop- mentation as a new phase transition universality class, agation in a metallic wire and the lambert w function, Physical Review Letters 114, 063903 (2015). American Journal of Physics 87, 476 (2019). [35] J. P. Boris, Relativistic plasma simulation-optimization [47] J. Andersen, Radiation from surface-wave antennas, of a hybrid code, in Proc. Fourth Conf. Num. Sim. Plas- Electronics Letters 3, 251 (1967). mas (1970) pp. 3–67. [48] M. Sporer, R. Weigel, and A. Koelpin, A 24 ghz dual- [36] P. B. Corkum, Plasma perspective on strong field mul- polarized and robust rod antenna, IEEE Trans- tiphoton ionization, Physical review letters 71, 1994 actions on Antennas and Propagation 65, 6952 (2017). (1993). [49] K. Amini, J. Biegert, F. Calegari, A. Chac´on,M. F. [37] T. D. Pointon, Second-order, exact charge conservation Ciappina, A. Dauphin, D. K. Efimov, C. F. de Moris- for electromagnetic particle-in-cell simulation in complex son Faria, K. Giergiel, P. Gniewek, et al., Symphony geometry, Computer Physics Communications 179, 535 on strong field approximation, Reports on Progress in (2008). Physics 82, 116001 (2019). [38] K. S. Thorne and R. D. Blandford, Modern Classical [50] S. D. Baalrud and J. Daligault, Effective potential theory Physics: Optics, Fluids, Plasmas, Elasticity, Relativ- for transport coefficients across coupling regimes, Physi- ity, and Statistical Physics (Princeton University Press, cal review letters 110, 235001 (2013). 2017).