Thomson and collisional regimes of in-phase coherent microwave scattering off gaseous microplasmas (original) (raw)

Abstract

The total number of electrons in a classical microplasma can be non-intrusively measured through elastic in-phase coherent microwave scattering (CMS). Here, we establish a theoretical basis for the CMS diagnostic technique with an emphasis on Thomson and collisional scattering in short, thin unmagnetized plasma media. Experimental validation of the diagnostic is subsequently performed via linearly polarized, variable frequency (10.5–12 GHz) microwave scattering off laser induced 1–760 Torr air-based microplasmas (287.5 nm O2 resonant photoionization by ~ 5 ns, < 3 mJ pulses) with diverse ionization and collisional features. Namely, conducted studies include a verification of short-dipole-like radiation behavior, plasma volume imaging via ICCD photography, and measurements of relative phases, total scattering cross-sections, and total number of electrons Ne in the generated plasma filaments following absolute calibration using a dielectric scattering sample. Findings of the paper suggest an ideality of CMS in the Thomson “free-electron” regime—where a detailed knowledge of plasma and collisional properties (which are often difficult to accurately characterize due to the potential influence of inhomogeneities, local temperatures and densities, present species, and so on) is unnecessary to extract Ne from the scattered signal. The Thomson scattering regime of microwaves is further experimentally verified via measurements of the relative phase between the incident electric field and electron displacement.

Subject terms: Plasma physics, Laser-produced plasmas, Techniques and instrumentation

Introduction

The scattering of electromagnetic radiation by an ionized gas is coupled with a multitude of plasma parameters—establishing a fundamental pathway for non-intrusive diagnostics with high spatial and temporal resolution1. In a manifestation devised by Shneider and Miles2, the total number of electrons (Ne) present in a collisional, weakly-ionized, and unmagnetized small plasma object can be measured through the constructive coherent elastic scattering of microwaves. By irradiating the object with a linearly polarized, uniformly distributed long-wavelength microwave field (for a plasma column of length L and diameter D, the skin depth δ≫D and microwave wavelength λ≫D,L), the plasma is periodically polarized and consequently emits short-dipole-like (short relative to the wavelength) radiation. The scattering is thus deemed in-phase coherent for the considered case of minimal phase variation amongst wavelets (scattered electric field contributions from individual electrons) at a far-field detecting location. Emitted radiation can then be measured/isolated using temporally resolved frequency mixer technologies and attributed to an absolute electron count after system calibration3. Both calibrated and uncalibrated variants of in-phase coherent microwave scattering (CMS) have become valuable non-intrusive diagnostic techniques in applications ranging from photoionization and electron-loss rate measurements310 to trace species detection11, gaseous mixture and reaction characterization1214, molecular spectroscopy15, and standoff measurement of local vector magnetic fields in gases through magnetically-induced depolarization16. Notable advantages of the technique include a high sensitivity, good temporal resolution, low shot noise11, non-intrusive probing, and the capability of time gating due to continuous scanning. The advantages of CMS over other available diagnostic techniques (microwave interferometry, laser interferometry, time-of-flight mass spectrometry, laser Thomson scattering, Langmuir probes, hairpin resonators) refer to limited sensitivities, lack of absolute calibration, necessity of signal accumulation, and plasma perturbations of the latter as previously reported elsewhere5, 17.

A careful understanding of electromagnetic-plasma and scattering interactions is essential for the proper extraction of information on total electron numbers and subsequent derivable quantities (e.g., the electron number density ne) from CMS. Vast historical art exists for weakly-ionized, strongly collisional microplasmas in which electron motion is restricted by electron-neutral collisions4, 5, 9, 10. However, application of the diagnostic to collisionless plasmas has remained relatively unexplored (no Ne correspondence is given by11, 13, 18)—despite garnering some attention in recent years19 due to a prevalence in studies on electric propulsion devices20, photoionization in low-pressure conditions21, 22, etc. Such plasmas are particularly advantageous for CMS in that a detailed knowledge of plasma and collisional properties (often susceptible to inhomogeneities, local temperatures and densities, present species, and so on) may not be required to characterize scattering.

Here, we establish a theoretical basis for coherent microwave scattering in the Thomson “free-electron” regime at reduced gas pressures and its transition to the conventional collisional scattering regime at atmospheric pressures—encompassing classical small plasma objects with diverse ionization and collisional properties. We subsequently present an experimental verification of this basis via variable frequency microwave scattering off laser-ionized air-based filaments for a range of pressures 1–760 Torr.

Results and discussion

Theory: coherent plasma scattering of microwaves

To establish a physical intuition for CMS, we first consider scattering contributions at a distance R≫V1/3 (V1/3 is a characteristic scale of the scattering volume V) from N nonrelativistic electron scatterers assuming each electron (distinguished by indices j,l) emits an electric field amplitude ES,0 at the observation location with frequency ω and relative phase Φj. Scattering of microwave radiation from relatively massive ions will be assumed negligible in the context of this work. The time-averaged scattered power per unit area at the detector location dPS/dA in SI units is subsequently given in Eq. 1, in which c is the speed of light in vacuum and ε0 is the vacuum dielectric permittivity1:

dPSdA=cε0∑j=1NEjS·∑l=1NElS=cε02NES,02+cε0∑j,l=1l≠jNEjS·ElS 1

where the first and second terms in Eq. 1 reflect noncollective and collective scattering phenomena, respectively. The collective scattering term is alternatively phrased as the “coherent” expression and can yield various contributions to the total scattered power as outlined below. First, in the case of full constructive interference between wavelets (deemed in-phase coherent), the collective term yields cε0∑j,l=1l≠jNEjS·ElS=cε02NN-1ES,02 and thus a maximum in the scattered power per unit area is reached dPSdAmax=cε02N2ES,02. Second, for phase-locked electrons exhibiting certain phase shifts at the detector location, the contribution of the collective term to the total power would be smaller dPSdA<dPSdAmax. The summation term is consequently zero ∑j,l=1l≠jNEjS·ElS=0 for phase-locked electrons exhibiting uniform and continuous relative phase distributions of Φj on the interval 0,2π. Third, if we are to consider significant phase randomization over the temporal resolution of the detecting equipment (i.e., from thermal motion, Doppler broadening, etc.)23, each term in the summation yields zero EjS·ElSl≠j=0. This case refers to the incoherent scattering regime and, from Eq. 1, the time-averaged scattered power per unit area is subsequently evaluated as dPSdA=cε02NES,02.

The position, velocity, and acceleration of individual scatterers in an ensemble (with respect to the detector) are central factors in determining coherency1. For Thomson scattering of electromagnetic waves on free plasma electrons at optical frequencies (electrons are collisionless on the period of optical field oscillations), the following simplified criteria is commonly utilized. If the parameter α=1kλD≪1 (where k=kI-kS and kI,kS are wavevectors of the incident and scattered radiation and λD is the Debye radius), then the wave resolves the motion of individual uncorrelated free-electrons within the Debye radius and the scattering is incoherent. Otherwise, for α≫1, the waves “see” the collective motion of different plasma elements and the scattering is coherent1, 24, 25.

In this work, we consider irradiation with microwave frequencies—ensuring that scattering occurs in the coherent regime. This conclusion generally follows from the fact that the microwave wavelength greatly exceeds the optical wavelength by about 5 orders of magnitude. Therefore, the signals generated at the detector from each individual electron in the plasma volume are coherent (e.g., phase randomization due to thermal motion and associated electron displacements comparable to λ occurs at a very slow time scale, Doppler broadening is unable to randomize phases for reasonable detector distances). More specifically, this work deals with a weakly-ionized plasma ellipsoid of length L≤1cm and diameter D≤0.5mm—where the linearly polarized microwave field (D≪L<λ∼3cm) is oriented along the ellipsoid length and the skin depth δ≫D. Therefore, all plasma electrons “see” the same intensity and phase of the incident microwave field. Phase differences generated by electrons from different portions of the plasma volume are then negligible as the considered experimental detector distance R = 10 cm greatly exceeds L2λ∼0.3cm. In this situation, full constructive interference between wavelets occurs and indicates that CMS scattering refers to the aforementioned in-phase coherent case.

Theory: model for microwave scattering from a short, thin plasma ellipsoid

In this section, we consider scattering of linearly polarized electromagnetic radiation off a plasma ellipsoid with the major axis (of length L) oriented parallel to the incident field. The considered ellipsoid will be regarded as short (L < λ), homogeneous, and thin (D≪λ,δ). The corresponding radiation characteristics can be sufficiently described by short dipole theory26, 27 and an accompanying plasma fluid model for the motion of plasma electrons. Generally, the equation of motion for a nonrelativistic electron in the swarm is:

where s is the electron displacement vector, e is the absolute electron charge, m is the electron mass, E and B are local electric and magnetic fields (respectively), and the effective collisional frequency for momentum transfer νm=∑tνm,t reflects the sum total of electron collisions with specie index t. The collisional frequency (νm=νm,eg+νm,ei=ngσtrvTe+neσCoulvTe, for which vTe is the electron thermal velocity and σ are cross-sections) is governed by electron–gas collisions νm,e-g for ionization degrees on the order of <10-2-10-3, while electron–ion collisions νm,e-i may become appreciable at higher ionization degrees28. νm resultingly exhibits complicated dependencies on temperatures, densities, and present species.

Equation 2 is subsequently reduced to the following form noting a unidirectional electric field and assuming negligible B-field interactions (similar to the approach used in the Lorentz oscillator model29):

where the total electric field inside the plasma volume (E) is composed of an incident irradiating field (EI) and a depolarizing field (Edep) due to an accumulation of charges at the ellipsoid ends.

The problem can be restated in complex formalism for an eiωt-process considering the incident electric field in the form EI=EI,0cosωt=ReE~I,0eiωt and the electron displacement in the form s=s0cosωt+Φ=Res~0eiωt to account for the phase-lag (Φ) between the electron displacement and the incident electric field (E~I,0 and s~0 refer to complex amplitudes). Then, since we examine scattering off a thin (D≪λ,δ), homogeneous plasma ellipsoid with the incident field oriented along its length, the established complex amplitude of the total electric field inside the plasma can be written as follows: E~0=E~I,01+ξε~-1 where ε~=ε′-iε″=1-ωp2ω2+νm2-iωp2ω2+νm2νmω is the complex dielectric permittivity for the plasma, ωp=nee2mε0 is the plasma frequency, and ξ is the depolarization factor for an ellipsoid30, 31. The established solution for the complex electron displacement amplitude is correspondingly written as s~0=-eE~I,0m1+ξε′-1-iξε″1-ω2+iωνm=-emE~I,0ξωp2-ω2+iνmω. Finally, the actual displacement amplitude s0 and phase with respect to the incident field Φ is reduced to the form:

s0=s~0=emEI,0ξωp2-ω22+νmω2 4

Note that an identical result can be obtained by evaluating the depolarizing field as Edep=ξmωp2es (or by including a restoring force term ξωp2s in the left side of Eq. 3)16.

The total dipole moment of the oscillating plasma channel can then be characterized by the amplitude d0=es0∫ner,z2πrdrdz=es0Ne. In the far-field (short dipole criteria R = 2_λ_26, 32), the total time-averaged scattered power can be treated as a point dipole and related to the dipole moment amplitude as PS=d¨26πε0c3→PS=ω4d0212πε0c3 and, using Eq. 4 and the incident microwave field intensity II=12ε0cEI,02, finally reduced to the form:

PS=e46πm2ε02c4IIω4ξωp2-ω22+νmω2Ne2 6

with associated far-field total and differential cross-sections for a point dipole:

σTot=PSII=e46πm2ε02c4Ne2ω4ξωp2-ω22+νmω2=σThω4ξωp2-ω22+νmω2Ne2=σeNe2 7
dσTotdΩ=38πσTotsin2θ;dPSdΩ=dσTotdΩII 8

where σTh=e46πm2ε02c4 is the Thomson cross-section33, 34 and σe=σThω4ξωp2-ω22+νmω2 is the total cross-section of an individual electron.

Three distinct coherent scattering regimes can be interpreted from Eq. 7 corresponding to the dominant denominator term—as detailed below and summarized in Table 1. First, the Thomson scattering regime can be identified which is associated with free plasma electrons oscillating in-phase with the incident microwave field (the ω4 term is dominant in the denominator of Eq. 7). The total cross-section of an individual electron then coincides with the classical Thomson cross-section σe=σTh and is thus independent of the microwave wavelength λ. Second, the collisional scattering regime (often referred to as the quasi-Rayleigh regime) refers to collision dominated electron motion (νm2ω2 term is dominant) with displacement oscillations shifted 90∘ with respect to the EI-field oscillations. The individual electron cross-section σe=σThωνm2 correspondingly scales as ∝1λ2. Previous experiments on microwave scattering off microplasmas have been conducted in this regime as discharges were operated at elevated pressures35, 9, 10, 35. A detailed knowledge of νm is consequently required for an accurate characterization of scattering. Third, the Rayleigh scattering regime can be observed which is associated with restoring-force-dominated electron motion (ξ2ωp4 term is dominant) with s oscillating in antiphase with EI. The total individual electron cross-section then evaluates as σe=σThω2ξωp22 and classically scales with the wavelength as ∝1λ4. Rayleigh scattering is not comprehensively explored in this work due to relatively low plasma densities and depolarization factors for the encountered experimental conditions (this regime can be somewhat exotic for practical implementation). However, the current work reports the first demonstration of in-phase coherent Thomson scattering of microwaves off miniature plasma objects.

Table 1.

Coherent microwave scattering regimes.

Regime Condition Coherent scattering cross-section, σe Phase shift, Φ
Thomson ω2≫νmω,ξωp2 σTh
Collisional νmω≫ξωp2-ω2 σThωνm2∝1λ2 90°
Rayleigh ξωp2≫ω2,νmω σThω2ξωp22∝1λ4 180°

A visualization of scattering regimes is illustrated in Fig. 1 for 11 GHz microwaves (as utilized in the experiment) using contour plots of individual electron scattering cross-sections (σe) and phase-lags between s and EI (Φ). An air-based, singly ionized thin/short plasma ellipsoid with depolarization factor ξ=0.01 and electron temperature Te=0.9 eV is considered for Fig. 1—where the collisional frequency was subsequently calculated as νm=nO2σeO2vTe+nN2σeN2vTe+neσCoulvTe (for values and assumptions on electron-neutral collisions, see the “Methodology” section—Electron number measurements). The Coulomb cross-section for electron–ion collisions was evaluated as σCoul=4π9e4lnΛkTe2,Λ=32πkTe3/2e3ne1/328.

Figure 1.

Figure 1

Scattering regimes. Contour plots for f = 11 GHz (a) total individual electron scattering cross-sections σe—calculated from Eq. 7 and (b) phase-lag between electron displacement and the incident microwave field—calculated from Eq. 5. The corresponding scattering regimes (Thomson, collisional, and Rayleigh) are additionally depicted. An air-based, singly ionized short/thin plasma ellipsoid with electron temperature T e = 0.9 eV (relevant to the considered experimental plasmas) and depolarization factor ξ=0.01 is considered.

Experimental results and discussions

We will now experimentally demonstrate the Thomson CMS regime and utilize the one-dimensional model described above for measurements of electron numbers and electron number densities in the Thomson scattering regime. To this end, we first consider unmagnetized, laser induced ≤ 760 Torr weakly ionized air-based microplasmas predominantly realized via the O2 resonance-enhanced multiphoton ionization (REMPI) mechanism (see “Methodology” section for details). More information on parameters of REMPI-induced plasma filaments at 1 atm can be found elsewhere7, 9, 36. An Abel inverted image (gate width 100 ns, 250 accumulations) of the produced non-equilibrium plasma for a background pressure of 760 Torr is illustrated in Fig. 2a. One can see that the plasma can be regarded as an oblate spheroid with upper bound (due to imaging uncertainties and diffusion) semi-axis estimates of a = b = 280 μm and c = 4.5 mm. For the probing microwave frequencies used, the filament geometrically corresponds to a one-third wavelength dipole antenna with a longitudinal depolarization factor of ξ≈0.0131. Two conclusions can be drawn in the context of this work (based on plasma imaging and ne-measurements described below). First, the depolarization factor appearing in the denominators of Eqs. 4, 5, 6, 7 and 12 can be neglected restricting scattering to the collisional and Thomson regimes (see “Methodology” section for details). Second, the approximation of a short (L<λ) and thin (D≪λ,δ) plasma ellipsoid considered in the theoretical model above is satisfied in experiments (the 11 GHz skin depth condition allows for ne measurements up to 1020 m−3 and 4×1021 m−3 for background pressures of 1 and 760 Torr, respectively, while actual measured ne-values were lower as shown below).

Figure 2.

Figure 2

Coherent microwave scattering from short, thin weakly-ionized plasmas. (a) Relative emission function (representative of an axisymmetric electron number density) of the laser-induced plasma at a background pressure of 760 Torr, calculated from an Abel-inversion of the corresponding ICCD image (gate width 100 ns, 250 accumulations). (b) Illustration of in-phase coherent microwave scattering and short-dipole-like far-field radiation behavior. (c) Short dipole (SD, line) and experimental f = 10.5 GHz plasma scattered (8 and 760 Torr, scatter) radiation patterns at t = 5 ns for a Tx/Rx distance of 10 cm. The microwave scattering signal (normalized VS) is proportional to the scattered E-field amplitude, and angular variation is based on default positions ϕ = 270°, θ = 90°, and β = 180°.

We then scatter 10.5–12 GHz microwaves off the REMPI-induced microplasma—as schematically shown in Fig. 2b, where the microwave electric field is oriented parallel to the plasma channel (see the “Methodology” section for details). The radiation pattern measured by an I/Q mixer-based microwave system at t = 5 ns (after initiation of the laser pulse) is depicted in Fig. 2c. Both the theoretical short-dipole (SD, line) and experimental (8 and 760 Torr, scatter) f = 10.5 GHz E-field amplitude ∝VS radiation patterns at a far-field Tx/Rx distance of 10 cm are illustrated. It was observed that the dependence of the measured scattered signal on θ and ϕ closely resembles the short dipole for both the collisionless (8 Torr) and collisional (760 Torr) plasmas over the early timescales presented in this work. Specifically, the radiation patterns in Fig. 2c demonstrate a sinusoidal dependence with respect to θ and isotropic with respect to ϕ. Further, the experimental pattern is indeed linearly polarized, as indicated by the β plot.

Temporal evolution of total electron numbers Ne in the plasma filament measured by 11 GHz microwave scattering at 1, 10, 100, and 760 Torr are shown in Fig. 3a—along with the normalized laser photodiode signal. The associated line-of-sight ICCD images at 100 and 760 Torr are additionally illustrated in Fig. 3b. One can see that Ne rapidly increases during the ionizing laser pulse and then decays shortly after (governed by dissociative recombination, electron attachment, and diffusion28, 3739). Observed decay times (> 5–7 ns, after accounting for the logarithmic nature of the plot) are consistent with others’ computational and experimental work for air plasmas4, 7, 9, 36, 38. Note that plasmas may persist well after the conclusion of the laser pulse—particularly for lower pressures where the rates of dissociative recombination and electron attachment are reduced.

Figure 3.

Figure 3

Coherent MS measurements of total electron numbers in variable pressure, short-dipole microplasmas assuming negligible depolarization and skin depth effects. (a) Temporal evolution of the total electron number measured via 11 GHz microwave scattering with the dielectric scatterer calibration technique at 1, 10, 100, and 760 Torr (with error bars)—alongside (b) line-of-sight intensified charge-coupled device images (gate width 100 ns) of the corresponding laser-induced plasmas. (c) Total number of electrons produced by the laser (approximated by the peak Ne value, which holds particularly well in the solid black-line region P≤ 100 Torr) as a function of background pressure for 10.5, 11, 11.5, and 12 GHz microwave scattering (with error bars).

The total number of electrons produced by the laser pulse Ne,laser (approximated by the peak Ne-value) is displayed on a logarithmic scale for four discrete microwave frequencies in Fig. 3c. It can be seen that Ne,laser scales nearly linearly with pressure in the range 1–100 Torr and begins to saturate (less-steep slope) for 300–760 Torr. This saturation onset at higher pressures can be explained by electron losses during the laser pulse and, as a result, underestimation of the actual Ne,laser-value for P ≥ 300 Torr. For example, an electron decay time of ~ 7 ns was observed at 760 Torr which is comparable with the laser pulse duration of ~ 5 ns. In contrast, observed electron decay times substantially exceeded the laser pulse duration for P ≤ 100 Torr and, therefore, Ne,laser can be reliably approximated by the peak Ne-value in that range.

Regions of Thomson, mixed, and collisional scattering regimes are indicated on Fig. 3c and correspond to the pressure ranges 0–10 Torr, 10–100 Torr, and 100–760 Torr, respectively. The boundaries of these regions can be estimated using the theoretical model above corresponding to the value of the effective collisional frequency for momentum transfer νm (νm∝ pressure) in relation to the microwave angular frequency ω. The plasma dynamics measured at 760 Torr corresponds to the collisional scattering regime and is consistent with previous works. A maximum total electron number of Ne = 3×1011 is correspondingly observed and, based on the previous Abel-inverted upper bound on the volume, the peak spatially averaged electron number density can be estimated on the order 10201m3. This density agrees well with similar studies in literature7.

In this work, in-phase coherent scattering of microwaves off a microplasma at reduced gas pressures in the Thomson and mixed collisional-Thomson regimes is demonstrated for the first time. The Thomson regime is particularly advantageous in that a detailed knowledge of collisional properties is unnecessary to characterize scattering. That is, an accurate description of νm (often difficult to determine due to unknown local temperatures and species densities, inhomogeneities, and so on) is not required to accurately extract Ne from the measured scattered signal VS since the following simple relation holds in the Thomson regime: Ne=mωAe2VS (see Eq. 14). Additionally, the Thomson regime has higher total electron number sensitivities than the collisional regime due to larger individual electron scattering cross-sections.

The Thomson regime is particularly advantageous for the evaluation of photoionization rates and cross-sections in comparison with collisional scattering at 1 atm used previously4, 5, 9, 10. In addition to νm-independent measurement results in the Thomson regime (as discussed in previous paragraph), it also enables measurements for an extended range of laser intensities by delaying the onset of nonlinear optical phenomena to even higher laser intensities. This refers to the fact that nonlinear optical phenomena including Kerr and plasma nonlinear terms in the refractive index reduce proportionally with a decrease in pressure (n2∝ng∝P,ωp2∝ne∝ng∝P).

Finally, we will experimentally verify the theoretical model of in-phase CMS by direct measurements of total scattering cross-sections of the dielectric calibration scatterer and measurements of the relative phase ∆Φ (see “Methodology” section for details). First, the total cross-section of the dielectric calibration sample was experimentally determined using direct measurements of the total power scattered by the sample and the incident microwave intensity at the bullet’s location to calculate σTot,D=PSII. This experimental value was then compared with the one predicted by the theoretical model for the case of a Teflon dielectric scatterer σTot,D=ω26πε02c4ε0εD-1ωVD2, where εD is the sample’s dielectric constant, and VD is the sample’s volume (electrical conductivity and depolarization factor of the bullet are assumed to be 0). Measurements matched the theoretical prediction within 20%. Second, relative measurements of the phase between the electron displacement and the incident electric field were conducted. Specifically, the parameter ΔΦ=Φ-Φ1Torr (representing a relative phase which is set to be 0 at P=1 Torr) was measured by the I/Q mixer-based microwave system in the pressure range 1–760 Torr. This was accomplished via measurement of the scattered radiation phase with respect to the local oscillator ΦS through the I¯ and Q¯ mixer channels. ΦS remained relatively invariant over the timescales presented in this work. The measured dependence of ΔΦ (at t = 3 ns) on background pressure is depicted in Fig. 4 for 11 GHz CMS, alongside the theoretical relative phase of electron displacement from Eq. 5. One can see that the phase shift changes about 90 degrees during the transition from the Thomson to the collisional scattering regime. This is fundamentally rooted in the fact that s and EI oscillate in-phase for free electrons in the Thomson regime (s¨∝-EI), while a 90-degree phase shift is introduced between s and EI in the collisional regime (s˙∝EI). Relatively good agreement between experimental data and theoretical predictions is additionally observed in Fig. 4.

Figure 4.

Figure 4

Relative phases for coherent microwave scattering as a function of background pressure. The I/Q-mixer measured relative phase of scattered radiation ΔΦ and theoretical prediction for 11 GHz microwave scattering at t = 3 ns.

Concluding remarks

We report the first demonstration of in-phase coherent scattering of microwaves off a microplasma at reduced gas pressures in the Thomson and mixed collisional-Thomson regimes. Ideality of the Thomson regime (< 10 Torr) refers to the fact that a precise knowledge of the collisional frequency νm is unnecessary to extract Ne from the scattered signal in contrast with the collisional regime. In-phase coherent microwave scattering offers substantially higher sensitivities in comparison with its incoherent counterpart in the optical frequency domain (laser Thomson scattering). This is due to the fact that the total scattering cross-section of a microplasma (and, correspondingly, total scattered power for a fixed irradiating intensity) is Ne times higher for the case of in-phase coherent scattering than for the case of incoherent scattering. The Thomson regime is particularly advantageous for the evaluation and tabulation of photoionization rates and cross-sections at extended range of laser intensities.

Methodology

A schematic of the general experimental setup is depicted in Fig. 5.

Figure 5.

Figure 5

Experimental set-up. Conceptual illustration of in-phase coherent microwave scattering in the far-field by a short, thin plasma object—coupled with a general schematic of the photoionization (purple), plasma imaging (blue), and microwave scattering (yellow) systems.

Vacuum chamber

Microwave scattering experiments were performed in a 0.66 m3 vacuum chamber lined with carbon-coated polyurethane foam (insertion loss ~ 22 dB, reflectivity level < − 18 dB). An MKS 925 Micro Pirani transducer (uncertainty 5%) was utilized for pressure monitoring. To ensure consistent air composition (P ≤ 760 Torr, T = 298 K, specific humidity ≤ 7 g/kg), the chamber was first mechanically pumped, purged and outgassed using a diffusion pump at 10–5 Torr, sealed, then filled with ambient air through a needle valve for experiments. Minimization of outgassing was confirmed after purging through an MKS Microvision 2 residual gas analyzer. UV fused silica and N-BK7 windows were installed for laser delivery and plasma imaging, respectively. Hermetically sealed SMA feedthroughs were inserted for microwave line connections and compression vacuum tube fittings were utilized for external angular positioning of the Rx microwave horn.

Laser-induced plasma generation

4.93 ± 0.25 ns FWHM laser pulses at 287.5 nm (FWHM 0.1 nm) were first produced by a 10 Hz, broadly tunable Ekspla NT342 Q-switched Nd:YAG laser system (1064 nm → SHG → THG → OPO → OPO-SHG). The laser was measured to be nearly linearly-polarized with an axial ratio of > 10—where the major axis is aligned parallel to the resultant of orthogonal unit vectors d^H,Tx and d^k,Tx (ref. Figure 2). The 7 mm diameter top-hat beam with full divergence angle < 2 mrad was then focused through a f = 175 mm plano-convex spherical lens to a peak intensity on the order of IL ~ 109 W/cm2 for plasma generation (maximum pulse energy of 3 ± 0.38 mJ after the lens, measured with a pyroelectric energy meter). A silicon photodiode with a rise/fall time of < 1 ns was used for relative intensity measurements at the position of plasma (assuming a constant beam waist), and an electrical output signal from the laser (jitter 0.3 ns) was utilized for external device triggering. Pulse energy could be varied via Q-switch timing modifications. The laser was finally quenched by a beam dump positioned sufficiently-far from the focal point.

A notable ionization process is observed for diatomic oxygen at 287.5 nm through the (2 + 1) REMPI mechanism. Specifically, via 2-photon excitation O2C3Πgv′=2,J′←O2X3Σg-v″=0,J″ followed by 1-photon absorption to ionization continuum O2+15. The effective ionization cross-section for (2 + 1) REMPI of O2 is significantly greater than that for tunneling and non-resonant multiphoton ionization of air constituents4043. This feature was confirmed experimentally by significantly reduced ICCD and microwave scattering signals off resonance (282 and 290 nm). The resulting photoelectrons naturally exhibit angular and kinetic energy distribution signatures of the underlying REMPI process41, 44. Thermalization, reactions, and supplemental radiation heating quickly modify these distributions. For an order of magnitude estimate, the electron temperature can be evaluated through multiphoton excess Te ≈ 3ћωL − UI,O2=0.868eV8. Supplemental electron heating is minimally expected for our laser conditions and multiphoton ionization can thus be interpreted as the primary contributor to electron temperature.

For the low energy depositions encountered in this work (non-detectable changes in laser pulse energy across the plasma on the precision 0.02 mJ) and nature of the REMPI mechanism, minimal optical perturbations to the neutral gas density and temperature are expected for the timescales presented8, 36. The plasma can be further regarded as predominantly weakly-ionized8. Despite REMPI-enabled selectivity, several species may arise over the laser pulse duration at sufficiently-high pressures due to secondary ionization (e.g. electron collisions with metastable species) and a myriad of reactions13, 3739—a set of processes hastened by high transport frequencies. A nonequilibrium air-based, weakly-ionized plasma is consequently considered for decay rate comparisons.

Plasma imaging

Photographs of the photoionized plasmas were taken using a Princeton Instruments PI-MAX 4: 1024i-RB ICCD camera with a lens of focal length fICCD = 55 mm (f stop fICCD/2.8). The camera was triggered using the laser electrical output pulse and surveilled with a gate monitoring signal. In reality, the signals observed can be highly intricate—exhibiting effects of straylight, laser-induced fluorescence, laser light scattering (elastic and inelastic), plasma emission, etc.45, 46. Straylight contributions were minimized through short gate times and a proper background environment, as validated through near-zero signals far-off the main laser line. The bulk of laser-induced fluorescence and laser light scattering processes were additionally found to be weak—as justified / confirmed through low N-BK7 transmissivity and camera quantum efficiencies in the ultraviolet range, imaging of off-resonant laser wavelengths (282 and 290 nm), and relative intensity comparisons with the corresponding microwave scattering signals for various photoionization conditions and ICCD gate features (delay and width). Some competition to plasma-related emission (see47 for relevant processes) is expected from two-photon-excited fluorescence (TPEF). However, overlap and intensity scaling similarities between TPEF and the (2 + 1) REMPI process roughly preserve the perceived plasma size. It should be noted that emission does not explicitly indicate the presence of electrons but can provide a reliable order-of-magnitude estimate of the plasma size9. An improved method for an estimation would entail the isolation of electron-induced fluorescence (see studies on the second positive system of N248). Further, the filament volume and distribution will be unique for a given pressure and timestamp due to optical effects45, 49, diffusion, and diverse decay/reaction mechanisms28.

In the case where emission is linearly proportional to the electron number density, a line-of-sight ICCD image represents the graphical projection of a three-dimensional plasma object rather than ne. For improved insight on the electron number density distribution and plasma volume, an Abel inversion (gr,z=-1π∫yRdhy,zdy·dyy2-r2, where g(r,z) is an emission function approximately proportional to ne and h(y,z) is the corresponding projection) was applied to the photograph in Fig. 2a under a reasonable assumption of axial symmetry. A Fourier-based algorithm was utilized for the inversion.

Microwave scattering (MS) system

The microwave field scattered off the plasma volume was measured by a receiving (Rx) horn coupled with an in-phase/quadrature (I/Q) mixer-based homodyne detection system that provides an output voltage VS∝ES,0 (where ES,0 is the amplitude of the scattered field at the Rx horn location)—as illustrated in Fig. 5 and described in detail below. First, an Anritsu 68369B synthesized signal generator was used to produce a continuous microwave signal at frequencies 10.5–12 GHz (FWHM 1 MHz). Such frequencies weakly interact with the neutral species considered in this work50, 51. The microwaves were then amplified and split. One arm was connected to a MITEQ IRM0812LC2Q I/Q mixer local oscillator (LO) port, with a confirmed operational power in the linear mixer regime (10–13 dBm). The other arm was then further amplified, isolated, and sent to a pyramidal 20 dB transmitting (Tx) horn for plasma irradiation. The linearly polarized microwave intensity at the plasma interface was measured to be ~ 3.5 W/m2 (relatively low incident microwave intensities were utilized to minimize plasma perturbations40 and nonlinear interactions49). Negligible heating by probing microwaves were supplementally confirmed though the invariance of plasma decay behavior on modest increases in supplied microwave power (0.1W→6W). A 15 dB receiving horn was then installed and connected to the I/Q mixer radio frequency (RF) port to detect the scattered electric field. Both the transmitting and receiving horns were placed inside the vacuum chamber at a proper far-field distance of 10 cm from the plasma ellipsoid center—where sealed SMA feedthroughs were installed for external connections to the I/Q mixer. The in-phase (Ῑ) and quadrature (Q¯) channels of the mixer were subsequently connected to a WavePro 735Zi 3.5 GHz oscilloscope with 50 Ohm DC termination. Calibrated amplification of the MW scattering signals (before the RF port) and the I/Q channels were used as needed to improve system sensitivity (at the consequence of uncertainty) and linear, undistorted amplifier operation was ensured throughout experiments.

A schematic of the plasma-microwave scattering experiment is shown in Fig. 2b. The corresponding notation is defined as follows: R—distance from the plasma ellipsoid center to the Rx detector central point, EI—incident electric field at the location of plasma, ES—scattered electric field at the observation location, Tx—transmitting horn, Rx—receiving horn, d^k,Tx/Rx—horn unit vector parallel to the wave vector when operating as a transmitter, d^E,Tx/Rx—E-plane horn flare unit vector, d^H,Tx/Rx—H-plane horn flare unit vector. Three angles are further defined: θ is the polar angle between the unit dipole vector (=d^E,Tx) and -R^ (=d^k,Rx), β is the angle between the Rx horn E-field flare vector d^E,Rx and R^×R^×E^I=T⇀, and the azimuthal angle ϕ (in which ES,0 is ideally invariant under rotations of ϕ for a short dipole) is additionally defined between d^k,Tx and the projection projUd^k,Rx=W⇀—where U is the plane formed by vectors d^k,Tx and d^H,Tx. Default positions ϕ = 270°, θ = 90°, and β = 180° are referenced for the subsequent angular variation of the Rx horn in experiment.

The amplitude of the scattered electric field (ES,0) can be related to the linear I/Q voltage expressions as follows:

VI¯=κBLO2ηBBcosΦB+κBLO2ηBScosΦSt=VI¯,0+ΔVI¯ 9
VQ¯=κBLO2ηBBsinΦB+κBLO2ηBSsinΦSt=VQ¯,0+ΔVQ¯ 10
ΔVI¯2+ΔVQ¯2=VS∝BS∝ES,0 11

where κ is the conversion factor, B is the local oscillator (BLO)/background (BB)/scattered radiation (BS) signal amplitude, η is the MW circuit attenuation factor, and ΦB/ΦS are corresponding phases of background/scattered signals relative to the local oscillator. ΦS can be determined through Eqs. 9 and 10 as ΦS=tan-1ΔVQ¯/ΔVI¯. For an improved signal-to-noise ratio and trial consistency (mitigating variations in laser pulse energies, vibrations, etc.), averaging over 120 events was accomplished via low-jitter laser triggering of the oscilloscope. Microwave scattering data at various probing frequencies were collected for diverse laser-induced air-based plasmas with background pressures in the range of 1–760 Torr.

Electron number measurements

The total number of electrons in the produced small plasma objects were determined from VS measurements coupled with MS system calibration via dielectric scattering samples. For CMS, scaling of VS with Ne can be generally derived based on results obtained in previous sections (see Eqs. 6, 7, 8, 11) as:

VS∝ES,0cosβ=3PS4πε0csinθcosβR=e2ω2EI,04πRmε0c2sinθcosβξωp2-ω22+νmω2Ne 12

where the total electron number measurements were conducted in the default configuration ϕ = 270°, θ = 90°, β = 180°, R=10 cm as shown in Fig. 2b. Highly prolate plasma ellipsoids (ξ<0.01) with moderate-to-low electron number densities were considered such that the ξωp2 term in the denominator of Eq. 12 is negligible and only the collisional or Thomson scattering regimes can be observed. In this case, scaling of VS with Ne can be presented as VS∝e2EI,04πRmε0c2ωω2+νm2Ne. One particularly convenient form of the VS-Ne relation (used in this work) is shown in Eq. 13, where A denotes an overall calibration factor of the MS system.

The calibration factor A was determined for each microwave frequency using a cylindrical PTFE sample with known relative permittivity εD = 2.1 and volume VD (length—12.7 mm, diameter—3.175 mm). The sample was mounted at the position of plasma, where the cylinder’s axis was oriented parallel to the laser wave vector. Assuming the PTFE bullet to be sufficiently short and thin with minimal depolarization effects (ED,0≈EI,0), the dependence of measured VS on dielectric scatterer properties can be written as the lower expression in Eq. 133, 35. This relation was utilized to determine the calibration factors A.

VS=Ae2mω2+νm2Ne-plasmascattererAVDε0εD-1ω-dielectricscatterer 13

The total number of electrons Ne in the produced plasma objects were finally evaluated using the following procedure. First, the calibration factor A was determined for each microwave frequency from VS measurements of the dielectric scattering sample using the lower expression in Eq. 13. For example, the correspondence constant for 11 GHz scattering was measured to be 48.4VΩcm2. Then, the derived calibration factors A and VS measurements of the plasma scatterers were used via the upper expression to determine Ne values. The effective collisional frequency for momentum transfer νm was estimated as a function of pressure under the assumption of a weakly-ionized air plasma with minimal laser-perturbations to the neutral background, a constant electron temperature of Te=0.9eV=10444K, and isotropic, elastic electron-neutral scattering. That is, νm=nO2σeO2vTe+nN2σeN2vTe≈1.88·109·PTorr, where σeg represents the electron–gas momentum transfer cross-section (σeO2≈7.4·10-20m2,σeN2≈9.74·10-20m2 for Te=0.9eV53), vTe=8kbTeπm is the thermal electron velocity, and ng is the neutral gas density. Despite often difficult to characterize dependencies on temperature, present species/reactions, specie densities, electron–ion collisions, etc.28, 35, this estimate remains relatively reliable for the considered plasma and experimental conditions5, 9, 54. A further discussion on the effect of electron temperature, specie densities, etc. on the scattered signal can be seen in our previous work55.

Note that Eq. 13 reduces to the following form when extreme cases of collisional (νm≫ω) and Thomson (νm≪ω) scattering are considered:

VS=Ae2mωNe-plasmascattererThomsonregimeAe2mνmNe-plasmascattererCollisionalregimeAVDε0εD-1ω-dielectricscatterer. 14

The Thomson regime is especially convenient for practical use since Ne (determined from the top expression in Eq. 14) is independent of collisional frequency ν m(and, correspondingly, independent of ng, Te, and σeg). Uncertainty in Ne measurements are attributed to skin depth and collective scattering effects, nonidealities of the MW system, and nonlinear laser phenomena.

Scattered radiation pattern measurements

The microplasma-scattered radiation pattern in the far field was assessed at 8 and 760 Torr for 10.5 GHz microwaves. In these experiments, the MW signal was not amplified to prevent mixer saturation. Further, data was not collected in inaccessible locations or positions where an excessive RF-port power level was experienced. Angular variation is based on the default configuration ϕ = 270°, θ = 90°, β = 180°, R=10 cm, as referenced in Fig. 2, where compression vacuum tube fittings were utilized for external angular positioning of the Rx horn. Data for the plots in Fig. 2c were generated by varying one angle at a time (namely, ϕ, θ, β) starting every time from the default configuration.

Scattering cross-section measurements

The total scattering cross-section σTot=PSII was determined using radiation pattern data and measurements with the same MS detection system at two locations. First, VS,def∝IS,def was measured in the default position (ϕ = 270°, θ = 90°, β = 180°, R=10 cm). Using the experimentally confirmed scattered radiation pattern ISθ,ϕ=IS,defsin2θ, the total time-averaged scattered power PS can be derived as PS=∫∫◯ISθ,ϕdS=∫∫◯IS,defsin2θdS=8π3R2IS,def. Next, the incident microwave intensity VI∝II was measured at the position of plasma (the Rx horn facing the Tx horn at R = 0, β = 180°) when the plasma-generating laser was disabled. The total scattering cross-section σTot was then calculated from measured quantities VI and VS,def (note that VS,def∝1/R):

σTot=PSII=∫∫◯ISdSII=IS,defII8π3R2=VS,defVI28π3R2. 15

Acknowledgements

This work was supported by the National Science Foundation (Grant No. 1903415) and the Directorate for Mathematical and Physical Sciences (No. 1903360).

Author contributions

Experimental design, data collection, post-processing, and figure preparation were performed by A.R.P., A.R., and X.W. Theoretical models and the interpretation of results were composed by A.R.P., A.S., M.N.Sh., and M.N.Sl. All authors contributed to the manuscript. This collaborative effort was led and organized by A.S.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Footnotes

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.