Streamer-induced kinetics of excited states in pure N2: I. Propagation velocity, E/N and vibrational distributions of N2(C 3Πu ) and N 2+ (B 2Σu+ ) states

The emission spectra of a streamer discharge in pure nitrogen provide an important tool for investigating the fundamental kinetics of excited electronic states of N2 and benchmark data for validating advanced kinetic schemes for numerical models. In this work, we characterize a streamer monofilament developed in a dielectric barrier discharge configuration, including electrical characteristics, time-resolved images and N2/N 2+ emission spectra, all acquired with nanosecond temporal resolution. Time-resolved images and emission characteristics provide clear evidence of the formation of a cathode-directed streamer and allow determining the streamer propagation velocity and the typical E/N values using the intensity ratio of nitrogen spectral bands in the center of the discharge gap. We also measure the vibrational distributions of the N2(C 3Πu , v=  0–4) and N 2+ (B 2Σu+ , v=  0–2) states. The population of N2(C 3Πu , v=  0–4) state, initially formed by energetic electrons in the streamer head, changes later significantly due to the decrease in the mean energy and concentration of the streamer channel electrons. After a few tens of nanoseconds, the electron-impact excitation rate of N2(C 3Πu ) becomes negligible compared to its population by the N2( A3Σu+ ) + N2( A3Σu+ ) pooling. The experimental findings are supported and consistent with the 0D state-to-state kinetic model results and reveal the participation of high vibrational levels of N2( A3Σu+ ) in the pooling reactions.

The emission spectra of a streamer discharge in pure nitrogen provide an important tool for investigating the fundamental kinetics of excited electronic states of N 2 and benchmark data for validating advanced kinetic schemes for numerical models.In this work, we characterize a streamer monofilament developed in a dielectric barrier discharge configuration, including electrical characteristics, time-resolved images and N 2 /N + 2 emission spectra, all acquired with nanosecond temporal resolution.Time-resolved images and emission characteristics provide clear evidence of the formation of a cathode-directed streamer and allow determining the streamer propagation velocity and the typical E/N values using the intensity ratio of nitrogen spectral bands in the center of the discharge gap.We also measure the vibrational distributions of the N 2 (C 3 Π u , v = 0-4) and N + 2 (B 2 Σ + u , v = 0-2) states.The population of N 2 (C 3 Π u , v = 0-4) state, initially formed by energetic electrons in the streamer head, changes later significantly due to the decrease in the mean energy and concentration of the streamer channel electrons.After a few tens of nanoseconds, the electron-impact excitation rate of N 2 (C 3 Π u ) becomes negligible compared to its population by the N 2 (A 3 Σ + u ) + N 2 (A 3 Σ + u ) pooling.The experimental findings are supported and consistent with the 0D state-to-state kinetic model results and reveal the participation of high vibrational levels of N 2 (A 3 Σ + u ) in the pooling reactions.
Keywords: dielectric barrier discharge, molecular nitrogen, second positive system of N 2 , streamer discharge, vibrational distributions

Introduction
Emission spectroscopy, although a simple and non-invasive diagnostic method is an essential tool for studying the basic characteristics of non-equilibrium discharges, including the processes of ionization and energy transfers driven by energetic electrons [1][2][3][4].Under streamer discharge conditions, most of the excitation and ionization processes are localized in the streamer-head region, which is characterized by an amplified electric field (E/N > 100 Td) and propagates at a remarkable speed (v ∼ 10 5 -10 6 m s −1 ) in the interelectrode gap.In molecular gases, the kinetics of streamers is very complex.A single-component gas like N 2 may involve hundreds of species in specific quantum states and thousands of elementary processes [5][6][7][8][9][10].A better understanding of the physics and chemistry of streamers then requires an approach based on a systematic comparison of experimental results with numerical modeling [11,12].Simplified approaches based on reduced data sets and models that consider electronic states of molecular species without their vibrational degrees of freedom cannot provide conclusive results regarding, for example, gas heating or dissociation kinetics.Indeed, the production of atomic nitrogen under non-equilibrium discharge conditions can result from direct electron-impact dissociation, or a sequence of electron-to-vibration (e-V), vibration-tovibration (V-V), vibration-to-translation (V-T) and vibrationto-dissociation (V-D) processes [13].While several studies address the importance of the vibrational distribution of ground state N 2 (X 1 Σ + g ) molecules in various phenomena, the role of the vibrational distributions of N 2 electronic states remains largely unexplored.Therefore, their investigation is of utmost interest to improve the predictive power and consistency of numerical models, allowing a deeper understanding and optimization of nitrogen/air-based discharges [14,15].
In previous works [16][17][18], we investigated post-discharge kinetics of metastable N 2 (A 3 Σ + u , v = 0-10) and ground-state N( 4 S) species in the dielectric barrier discharge (DBD) geometry employing laser-induced fluorescence (LIF) and twophoton absorption LIF (TALIF) techniques, respectively.A triggered mono-filamentary streamer was periodically produced in pure nitrogen at 50 and 200 Torr pressures.The investigated species were monitored during the streamer channel decay in the center of the discharge gap.At 200 Torr, densities of the metastable levels N 2 (A 3 Σ + u , v = 6-10) were determined to peak with the densities of approximately 10 14 cm −3 about 100 ns after the discharge onset [16].The population maxima of N 2 (A 3 Σ + u , v < 6) levels occur with a delay of 0.5-10 µs after the streamer onset, and this delay increases with the decreasing vibrational number.This is due to the radiative/quenching cascade N 2 (C 3 Π u ) → N 2 (B' 3 Σ u )/N 2 (B 3 Π g )/N 2 (W 3 ∆ u ) → N 2 (A 3 Σ + u ), populating the N 2 (A 3 Σ + u , v) states with µs delay after the active discharge phase.The N( 4 S) density reaches the value of approximately 1.4× 10 15 cm −3 and peaks 10 µs after the discharge initiation.
The N( 4 S) production is, for low gas temperatures (300-400 K) as in this case, due to the predissociation of molecular nitrogen singlet states and stored mainly in the three lowest states of the atomic nitrogen (N( 4 S), N( 2 D), N( 2 P)), which later interact with the other dominant species e, N 2 (X 1 Σ + g ) and N 2 (A 3 Σ + u ), therefore, affecting also the post-discharge phase.The understanding of the temporal evolution of N( 4 S), N 2 (A 3 Σ + u ) and underlying kinetics requires knowledge of the temporal profile of E/N since the production of all significant species springs from the electron-impact processes in the active discharge phase.However, it is extremely challenging to determine E/N in streamer discharges due to its sub-nanosecond time evolution.The measurement of E/N can be addressed with the non-invasive intensity ratio method of N 2 emission bands [19], which relates E/N with the ratio of intensities of the first negative system (FNS) + and the second positive system (SPS) C 3 Π u → B 3 Π g of N 2 .The intensity ratio is very sensitive to the change of E/N due to different cross-section profiles and energy thresholds for the electron-impact population of the N + 2 (B 2 Σ + u ) and N 2 (C 3 Π u ) states, which are 11.0 and 18.8 eV, respectively.Moreover, the study of the FNS and SPS emissions and their vibrational distribution functions (VDFs) also bring information about the dominant processes populating the excited states and the confirmation of the applicability of the intensity ratio method.The present work consists of two parts.Part I corresponds to this paper, in which we analyze the electrical characteristics and two-dimensional dynamics of the UV emission produced by an N 2 streamer monofilament developing in a point-plane DBD electrode geometry at a pressure of 200 Torr.We propose a method to determine the streamer velocity at near-to-atmospheric pressure using a parallel combination of intensified charge-coupled device (ICCD) and photomultiplier tube (PMT) measurements.Based on electro-optical measurements, we reveal the fundamental processes and obtain the temporal evolution of vibrational distributions of the N 2 (C 3 Π u ) and N + 2 (B 2 Σ + u ) states.Finally, we compare the observed emission characteristics with the results of a detailed 0D model, which uses as an input parameter the electric field determined using the FNS/SPS method.The paper consists of two sections: Experimental setup and methods section 2 includes details of the used discharge reactor, a description of emission spectroscopy and imaging techniques and information about data processing and analysis.The results and discussion section 3 incorporates measured electrical characteristics, streamer imaging, the determination of the propagation velocity, measured ICCD emission spectra, VDFs of N + 2 (B 2 Σ + u ) and N 2 (C 3 Π u ) states, determination of E/N and implications for kinetic modeling.
Forthcoming Part II will focus on the vis-NIR emission characteristics, including the complete vibrational distribution of the N 2 (B 3 Π g ) state and will be presented in a companion paper.The UV-vis-NIR emission characteristics will then be analyzed using an advanced kinetic model.

Discharge reactor
We carried out the present experiments using a diagnostic discharge reactor described in detail in preceding works [16][17][18] and illustrated in figures 1(a) and (c).The DBD consists of a point-plane electrode system placed in a stainlesssteel chamber based on the KF 6-way cross.The copper tips (1 and 16 mm for the high-voltage (HV) and grounded electrode, respectively) were embedded in MACOR glassceramic disks (0.5 mm thick dielectric layers, relative permittivity is approximately 6 at 25 • C and 1 kHz).The discharge was powered by employing an HV waveform made of two sine-waves (f AC = 1 kHz, peak-to-peak the amplitude of 19.4 kV) superposed with a fast positive pulse (with the amplitude of 9 kV) acting on the second positive sinusoidal half-cycle, see panel (b) of figure 1.This complex waveform was applied at a fixed repetition frequency of f M = 10 Hz.Within each waveform, the fast positive pulse produces a triggered streamer mono-filament that is locked with respect to the onset of the pulse (HV rise time of a few nanoseconds).Untriggered discharges regularly occur during the preceding AC phases (typically one streamer mono-filament per one AC half-cycle), leaving residual electrons that serve as initial charges for subsequent triggered events.
Figure 2 shows the block diagram of the setup used to drive the discharge periodically by a specific HV waveform.The AC high-voltage power supply is composed of a TG1010A Function Generator (TTi), a Powertron model 250A RF amplifier, and a high-voltage step-up transformer.The pulsed power supply generates the positive HV pulse (u PPS , peak amplitude 9 kV, 100 ns duration, the typical rising slope of 200 Vns −1 , repetition rate f P = f M ) by discharging a 10 nF capacitor through a Behlke model HTS-81 solid state switch (8 kVdc/30 A).We used two channels of the digital pulse/delay generator (BNC Model 575) as a master trigger to synchronize the AC burst with the positive HV pulse at a suitable phase.Other BNC channels were used to synchronise the diagnostic equipment (ICCD detectors and multichannel photon counter).The discharge was fed with pure nitrogen (99.9999 % N 2 + 0.5 ppmv H 2 O + 0.3 ppmv O 2 + 0.1 ppmv total hydrocarbon) through Bronkhorst HI-TEC mass flow controller (total flow of 85.3 sccm).The pressure of 200 Torr in the chamber was stabilised by using a needle valve backed by an oil-free diaphragm vacuum pump, and measured with a TPG 260 piezo gauge (Pfeiffer Vacuum).The present experimental conditions are dictated by conditions used to perform LIF/TALIF experiments [16][17][18]; the pressure was set as a suitable compromise, which ensures high stability (low onset jitter) of the discharge and a sufficiently large radius of the streamer filament for LIF diagnostics.A Tektronix DPO 5204B oscilloscope was used to record the discharge characteristics (the HV, current, PMT waveforms) and to control the timing of the discharge with respect to the detectors.The discharge voltage waveforms were sampled through a Tektronix P6015 high-voltage probe (1000:1@100 MΩ).The current pulses produced by individual micro-discharges were monitored through the Pearson probe right before the high voltage electrode.

Time-resolved emission spectroscopy and imaging
We used the digital pulse/delay generator (BNC Model 575) as a trigger and delay source to synchronize the HV waveforms with diagnostic equipment (ICCDs, oscilloscope).Two BNC 575 channels were used to generate the triggered discharge, while the other channels were used for synchronous electrical characteristics and spectrometric data acquisition.We simultaneously used two ICCD-based imaging spectrometers to collect spectrometric data in the 200-1100 nm spectral range.To analyze plasma-induced emission (PIE) in 200-850 nm inter- val, we used Jobin-Yvon iHR-320 spectrometer equipped with an Andor DH740 ICCD detector.On the other hand, the Andor Shamrock 303i spectrometer with DH340 ICCD was used to collect spectra in 600-1100 nm range.All emission spectra were calibrated only relatively since the absolute calibration has no effect on VDF analysis and E/N determination.The iHR-320 was used in an imaging mode (1:1 image through a pair of UV-grade lenses).The Shamrock 303i was used in full vertical binning mode because the PIE was delivered through the 10 m long quartz fiber bundle.By directly imaging the DBD gap onto the entrance slit of the iHR-320, we could also take time-resolved images of the gap using the spectrometer's grating in the 0th-order of diffraction.Furthermore, to better follow the evolution of the discharge morphology during a single event, we used a 4-channel XX RapidFrame camera (Stanford Computer Optics, Inc), essentially a four-channel ICCD system based on UV-extended mirror-based image splitting optics and four 4-Picos ICCD cameras.The image splitter splits the arriving flux of photons equally among the four exit channels (Ch1-Ch4), each registered by one ICCD.Each ICCD detector is equipped with an 18 mm dual-stage multichannel plate (MCP) image intensifier (S20 photocathode), allowing the acquisition of time-resolved images taken with a pre-selected MCP gating time (a minimum gate of 1 ns was used in this work).This enables the registration of four successive time-resolved images taken with a required frame rate.The four ICCD detectors' relative positions of the MCP gates were controlled through the internal delay generators of the respective 4-Picos cameras.Complementary PIE characteristics measurements were performed using fast PMTs (Hamamatsu H10721, Photek PMT210) combined with appropriate bandpass filters.

Data processing and analysis
The acquired datasets were analyzed by using synthetic models for N 2 and N + 2 emission spectra reported in detail in [3,20,21].Based on the fact that a typical N 2 /N + 2 emission spectrum consists of a group of vibronic bands arranged in so-called sequences (a group of bands characterized by a constant difference between the vibrational numbers of the upper/lower states participating in the radiative transition), we use an approach based on the calculation of a set of band profiles by considering the actual experimental instrumental function of the ICCD spectrometer and assuming an equilibrium population distribution in the upper rotovibronic states (therefore characterized by a certain rotational temperature T rot ).Then the reproduction/interpretation of the experimental spectra is reduced to finding a suitable multiplicative constant for each band profile such that the superposition of all calculated band profiles (after being multiplied by the corresponding constants) minimizes the residual spectrum obtained as the difference between the experimental and calculated spectra.The set of multiplicative constants then contains information on the vibrational distribution of the upper (emitting) vibronic state.

Electrical characteristics
In order to understand the discharge initiation and the charge transfer in discharge chamber, we consider an equivalent circuit which, as best as possible, represents the circuit.The equivalent electric circuit consists of lumped elements as shown in figure 3(a).The voltage source V 1 represents a device which superposes the low-frequency harmonic waveform and short (100 ns) voltage pulse (∼9 kV) with approximately exponential decay.The internal impedance of V 1 is represented by the resistor R i and the inductor L i .The elements C p and R p represent loading by the Tektronix P6015 high-voltage probe.The discharge cell's inner electrodes are accessible via two short transmission lines T 1 and T 2 .The anode is fed from V 1 via T 1 , whereas the cathode is connected to the ground via T 2 and an inductor L 2 .The capacitor C 1 represents dielectrics covering the electrodes, and C 2 defines the capacity of the gap between the dielectrics.These capacities were estimated based on the geometrical size of the dielectrics and the gap (C = ε 0 ε r S d , where the S denotes the area of the electrodes, d gap or dielectrics thickness and ε 0 , ε r the permittivity of vacuum or the dielectrics, respectively).The streamer monofilament in the gap, and an unavoidable conductive leakage between the HV electrode and the groundings, are defined by the two variable resistors denoted filament and leakage, respectively.The meaning of the leakage resistor is that the filament current is not restricted to flow to the grounded inner walls of the discharge chamber.
Figure 3(b) shows the experimentally measured voltage (red) and current (blue) waveforms at the V 1 , when there is no discharge between the electrodes due to elevated gas pressure.The simulation on the equivalent circuit is performed in the LTspice software [22], where the measured voltage waveform is a simulation input, whereas its output is the simulated current (green).The simulated current's noisy character originates in the weak noise present in the voltage waveform as the load impedance is predominantly capacitive.The noise is not present in the measured current waveform, probably because the noise bandwidth is larger than 1 GHz, whereas the bandwidth of the Pearson probe is limited to 200 MHz. Figure 3(c) shows voltage and current waveforms at the V 1 with the discharge.The measured current waveform (blue) is obviously elevated compared to the previous case in figure 3(b).Although the discharge itself is modelled by a variable filament resistor, to achieve a similar simulated current at V 1 , it is necessary to add the additional leakage resistor to the model.Without this resistor, the agreement between the measured and simulated currents is not good.Individual components of simulated current are plotted in figure 3(d), where they are labeled as follows: filament current (blue), leakage current (magenta), and the current on the dielectrics (I C1 , black).The steep increase of I C1 is due to streamer arrival to the cathode, connected with the charging of the dielectric cathode.This increase begins approximately 14 ns after the high voltage onset.The initial charging of the dielectric cathode, defined by a relatively small capacity C 1 , is after approximately 11 ns taken over by the leakage resistor.The probable explanation is that: although the streamer filament initially gaps the dielectrics, the surface ionization wave later propagates along the dielectric electrodes into remote areas of the chamber, where no insulation exists.Therefore, the leakage current flowing between conductive parts in the chamber cannot be excluded.Figure 3(e) plots the time dependence of estimated filament (blue) and leakage (magenta) resistances, which are found by the simulation to be appropriate to reproduce the measured current.

Streamer imaging: 4-channel ICCD
The fast time-resolved imaging using 4-channel ICCD provides information about the streamer morphology and dynamics.The most important results obtained through the 4-channel XX RapidFrame are summarized in figures 4-6. Figure 4 shows Ch1-Ch4 images acquired during a single event.Channels Ch1-Ch2 integrate the very weak emission produced during the rising slope of the second positive AC half-cycle before the onset of the superimposed HV pulse.These two ICCD images capture the very slow phase of the regeneration of the charges in the anode region after the polarity reversal of the HV electrode (temporary cathode to temporary anode).Channels Ch3 and Ch4 integrate the emission during the first 20 ns of the HV pulse and several microseconds after the pulse.Ch3 provides a characteristic image of a fully developed streamer monofilament in the gap.The image from Ch4 shows the glow disappearance after the streamer decay.Note that in order to visualize the luminosity of the formative pre-streamer phase, we had to use MCP gates with microsecond duration (5 and 2 µs for Ch1 and Ch2, respectively) compared to nanosecond gates during the streamer phase (Ch3, Ch4). Figure 5 shows a series of four consecutive frames captured during the onset of a streamer triggered by a HV pulse.Four images taken during one event in four consecutive time windows allow monitoring the dynamics of monofilament formation.They were obtained using the same MCP gate widths (1 ns) for all channels and were appropriately delayed with respect to the onset of the HV pulse and with respect to each other (∆t = 3 ns).Events 1-3 were selected from all captured single discharge events to illustrate the basic PIE characteristics during the first tens of nanoseconds.
Finally, figure 6 shows a kinetic series of ICCD frames created from Ch2 images and selected from a set of 500 registered events similar to those shown in figure 5 (i.e.all acquired with a 1 ns MCP gate, but with a variable delay relative to the HV pulse).Each frame in figure 6 was created by averaging five equivalent Ch2 images (taken at the same MCP delay and showing similar evolution on all four Ch1-Ch4 channels) and corrected for background, resulting in significant noise reduction.The image in the upper right corner captures the moment immediately after the impact of the streamer on the dielectric surface covering the temporary cathode, and the corresponding MCP delay is denoted as T 0 .The three preceding images in the top row capture the evolution before streamer formation.The images in the middle row demonstrate the formation of the cathode spot over the next ten nanoseconds.The streamer velocity on the cathode surface can be determined by investigating the evolution of the surface streamer front in time, based on the middle row.However, this method is not very precise since figure 6 was obtained by averaging.Nevertheless, the streamer velocity on the cathode surface reaches approximately (2±1)×10 5 m s −1 .The images in the last row show the discharge's upward expansion along the temporary anode's dielectric surface.Based on the observed features of the discharge evolution, the following mechanism of streamer initiation between the dielectric barriers can be concluded, similarly as described in detail in [23]: (i) Townsend avalanche pre-phase (charge accumulation, t < T 0 ).(ii) Avalanche to streamer transition, propagation of the cathode-directed streamer and its impact on the cathode (t ≈ T 0 ).(iii) Propagation of surface streamer on the cathode dielectric (t > T 0 ).(iv) Charge redistribution: electrons, generated during the streamer impact on the cathode, are transferred to the dielectric surface on the anode (t > T 0 + 16 ns).The charge redistribution in the gap causes a decrease in the electric field.

Streamer ICCD imaging: determination of the propagation velocity
The determination of the streamer propagation velocity in a 4 mm DBD gap is a challenging task.On the one hand, the statistical dispersion of the streamer onset time, in this case of ∼ 6 ns, makes it difficult to average the streamer propagation over several events.On the other hand, the minimum imaging gate of 2 ns of the DH740 ICCD, when it is necessary to obtain a sub-nanosecond evolution, turns the task even harder.Herein, we propose a method to determine the streamer propagation velocity that overcomes both aforementioned problems and is described by the following steps.
(i) For different opening times of the 2 ns gate (t gate on ), chosen to catch the streamer propagation phase, 500 ICCD images   ing the first vertical position, starting from the bottom, whose horizontally-averaged intensity is higher than a certain threshold.In this case, a threshold of 1.5 × 10 4 is optimal to determine the SFP.However, note that using thresholds of 10 4 or 2 × 10 4 leads to the same trend of SFP as a function of time, with the exception of a constant time shift, which is not relevant for the velocity determination.(v) Finally, the streamer propagation velocity is calculated by fitting the SFP(t gate corr ) with linear functions.
From the observation of figure 7(b), we find that the streamer takes approximately 3 ns to propagate from the top (anode) to the bottom (cathode) electrode.Additionally, there are clearly two different linear trends in the streamer propagation: initially, the streamer velocity is ∼ 9.8 × 10 5 m s −1 and then, in the lower half of the gap, the velocity increases to ∼ 2.4 × 10 6 m s −1 .After the streamer arrives at the cathode, it propagates around it, leading to the apparent deceleration pointed out by the third linear trend.Light reflections at the cathode surface may also play a role in this third stage.Note that, simultaneously with the streamer propagation over the gap, there is also propagation around the anode, as evidenced by the increased luminosity after 1 ns at positions higher than the end of the top electrode.
We should remark that the study is done with only 500 ICCD images, which corresponds to an effective measurement time of ∼1 min (considering the 10 Hz repetition frequency).A similar study done by varying the opening time of the gate t gate on would require a very large number of averages at each point in order to cancel the effect of the jitter on the streamer onset.For this reason, such a study would take at least several hours.This shows the potential of the method presented in this work.

ICCD emission spectra
The UV part of the spectra contains a blend of the SPS of N 2 and the γ-system of nitric oxide (NO-γ), see figure 8(a).Even in high-purity nitrogen, the dissociation of O 2 impurities subsequently forms nitric oxide in the ground state (NO(X 2 Π)).Electronic excitation of NO(X 2 Π) then occurs either through electron-impact excitation [24]: or resonant energy transfer from the N 2 (A 3 Σ + u ) metastables [25][26][27][28]: Therefore, the NO-γ emission can serve as an important diagnostic tool about the presence of N 2 (A 3 Σ + u ) metastables.The NO-γ emission begins to occur at the time of approximately 60 ns; at these times, the vibrational levels w > 5 of N 2 (A 3 Σ + u ) peak [16].Therefore, we can expect that the onset emission is connected mainly with the energy transfer from the higher levels, which is later replaced by the energy transfer from the lower N 2 (A 3 Σ + u ) vibrational levels.This is also in agreement with the higher rate constants of the process (2) for the higher w-levels; see [27] and references therein.
The observation of the Goldstein-Kaplan and NO-γ bands confirms the presence of N 2 (A 3 Σ + u ) states at the early postdischarge phase, t ∈ (80, 200) ns.

VDFs of N +
2 (B 2 Σ + u ) and N 2 (C 3 Πu) states The investigation of the VDFs allows revealing the specific precursors of the excited N 2 species.The Franck-Condon factors (FCF) are defined using vibrational overlap integrals and serve as a benchmark when evaluating the efficiency of vertical electron-impact excitation processes [32].Therefore, it is usually assumed that when exciting multiple vibronic levels (e.g.N 2 (C 3 Π u , v = 0-4)) from a single ground state level (N 2 (X 1 Σ + g , v = 0)) the population of excited states is proportional to the corresponding FCFs between N 2 (X 1 Σ + g , v = 0) and N 2 (C 3 Π u , v = 0-4).Under present discharge conditions, it provides important evidence justifying the use of the intensity ratio method for the determination of the reduced electric field (confirmation of N 2 (C 3 Π u ) and N + 2 (B 2 Σ + u ) state production due to direct electron-impact ionization from the N 2 (X 1 Σ + g , v = 0) ground state).The analysis of N 2 (C 3 Π u ) VDF is relatively simple due to the high intensity of SPS radiation and many bands/sequences that can be analyzed.Conversely, determination of the VDF of N + 2 (B 2 Σ + u ) is significantly more difficult due to the FNS lower intensity and the omnipresent overlap with SPS bands, especially under conditions of a low-to-medium range of E/N (100-200 Td) [12,33].
The SPS bands belonging to the ∆v = 0, −2, and −3 sequences were analyzed using the technique detailed in [3,21].When calculating synthetic SPS band profiles, we considered the experimental instrumental function given by dispersion grating and entrance slit width.Figure 9(a) shows the UV spectra from the wavelength range 385-400 nm, revealing FNS(0,0) and the SPS bands belonging to the ∆v = −3.(2,4) and FNS(0,0) PMT profiles as well as the band intensities of SPS(0,3), SPS (1,4) and SPS(2,5) obtained by the spectrometer using the technique of kinetic series.These measurements can be converted to the time evolution of N 2 (C 3 Π u , v = 0-4) profiles as shown in figure 9(c).Obviously, the onset of the ICCD intensities (1-8 ns) is more gradual than the onset measured using the PMT, which is due to the 2 ns gate width of ICCD used to collect spectra.Complementary to the preceding, figure 10 shows the time evolution of the VDF within 950 ns after the discharge onset.
(i) Interval 0-8 ns, VDF shows a typical FCF-like shape, which reveals the initiation of the active discharge phase (see the population comparison in table 1 and discussion further below).During this phase, the N 2 (C 3 Π u , v = 0) and N 2 (C 3 Π u , v = 1) densities can be described using the following formulas.We use the shortened notation where the k C0 im and k C1 im denote rate constants of electron impact excitation of the N 2 (C 3 Π u , v = 0) and N 2 (C 3 Π u , v = 1) states, respectively.(ii) Interval 12-20 ns, the N 2 (C 3 Π u , v = 0-4) densities start to reflect quenching by molecular nitrogen and relax from the higher to the lower vibrational levels due to vibrational relaxation within the N 2 (C 3 Π u ) state: The evolution of the lowest two states N 2 (C 3 Π u ) can be described by the following equations: where the k C10 rel is the rate constant of the equation ( 7) for v = 1 and w = 0. τ C0 eff and τ C1 eff denote effective lifetimes of N 2 (C 3 Π u , v = 0) and N 2 (C 3 Π u , v = 1) states, respectively.Note that τ C1 eff includes the influence of the process (7).In the case of N time interval implies that the gain of the excited states population is already balanced by the losses due to quenching by N 2 , confirming that the electric field does not further increase, since k C0 im and k C1 im increase with increasing electric field.(iii) Interval 40-60 ns, the VDF reaches its minimum since the electric field is already very low.Higher vibrational levels v ⩾ 1 are relaxed to the N 2 (C 3 Π u , v = 0).We observe that [N2(C 3 Πu,v=0)] < 10%, corresponding to characteristic N 2 (C 3 Π u ) state vibrational temperature of approximately 1300 K.The equations describing the N 2 (C 3 Π u , v = 0) and N 2 (C 3 Π u , v = 1) densities are the following: (iv) Interval 60-1000 ns, the switching to N 2 (A 3 Σ + u )-poolingcontrolled VDFs.The N 2 (A 3 Σ + u )-pooling reactions affect observed populations of the N 2 (C 3 Π u , v = 3-4) states at first (60-80 ns) since their population by the electron impact was low.In the later times (80-1000 ns) also, the lower N 2 (C 3 Π u , v = 0-2) levels are populated by the pooling reactions [35,36].The VDF still changes when comparing the (140-350 ns) and (500-950 ns) cases.The change is due to different N 2 (A 3 Σ + u , v) vibrational levels populating the N 2 (C 3 Π u , w) states by the pooling reactions.These results will be further discussed in the section 3.7.The equations describing the N 2 (C 3 Π u , v = 0) and N 2 (C 3 Π u , v = 1) densities are the following: The k pool denotes the pooling reaction rate constant from metastable states and The equations can be written in the steady state, since the nanosecond timescale of quenching dominates over the microsecond timescale of pooling reactions.
The conditions for the use of the intensity ratio method are fulfilled in the interval 0-20 ns, where the population of the N 2 (C 3 Π u , v = 0-4) states is due to the electron-impact excitation also with the contribution of the quenching and relaxation processes.
Table 1.Comparison of the experimental and FCF-like VDFs.The SPS errors were estimated using three different distributions obtained from three SPS sequences (∆v = 0, −2 and −3).Analogous to SPS, we analyze the VDF of the N + 2 (B 2 Σ + u ) state.We have chosen 380-480 nm interval containing ∆v = 0, −1 and −2 sequences of the FNS. Figure 11 shows synthetic bands of the FNS (upper trace) and SPS (lower trace), both calculated by assuming FCF-like VDFs (using triangular instrumental function and rotational temperature of 300 K).The ∆v = 0 sequence of FNS shows relatively low intensities of bands originating from v >0 vibrational levels; furthermore, this region is affected by several SPS bands.The ∆v = −1 sequence of FNS shows reasonable intensities of (0,1) and (1,2) bands while ∆v = −2 also predicts weak (2,4) band at 460 nm.Upon inspection of the experimental spectra obtained during the first nanoseconds of the discharge, it is clear that the FNS system cannot be analyzed without separating the FNS bands from the dominant SPS emission as illustrated in figure 12. Upper panels (a),(b) in figure 12 show experimental spectra in 405-480 nm interval.All labeled SPS bands were subtracted from (a),(b) and the resulting residuals are plotted in the lower panels (c), (d).Residuals (c) in the 420-430 nm interval clearly show the shapes of the two FNS(0,1) and FNS(1,2) bands.Residuals (d) in the interval 455-475 nm reveal the shape of the FNS(0,2) band and a small peak at the position of the FNS(2,4) band.The residual peak occurring at 465 nm is due to overlapping almost coinciding SPS(4,10) and FNS (1,3) bands.In order to estimate the VDF, we used intensities of FNS(0,1) and FNS(1,2) from (c) and FNS(0,2) and FNS(2,4) from (d).Experimental VDF and FCF-like VDF are compared in table 1.The value of VDF for N + 2 (B 2 Σ + u , v = 1) level is notably below the FCF-like value.The possible explanation could be that the ratio of rate coefficients for N ) levels is not simply given by a ratio of FCFs; as already reported in the case of the N 2 (C 3 Π u ) state under low E/N conditions [3].The primary reason for such behavior is the non-FCF scaling of N 2 (X 4) cross-section amplitudes in the region of the excitation threshold, which results in an increasing contribution of threshold electrons to the total rate constant with decreasing E/N.Similar trends should be observed in the case of FNS at comparably higher E/N due to the notable difference in the N 2 (X ) ionization thresholds.However, SPS and FNS spectra obtained over a wide range of E/N would be required to provide explicit experimental confirmation of this statement, which is beyond the scope of this work.On the other hand, based on the same reasoning, close to FCF-like scaling should be expected in the case of the electronic excitation process e + state (1:0.44:0.065).Therefore, based on the previous arguments, we believe that the obtained VDF can be considered a signature of a direct e process at present conditions and justifies the use of the FNS/SPS method for E/N evaluation.

E/N estimation
Determining the reduced electric field using the FNS/SPS intensity ratio technique includes many technical details described in the following paragraph.Figure 13(a) presents band intensities obtained by the spectrometer using the technique of kinetic series with 2 nanosecond MCP gate and step.This measurement informs us about the ratio of the FNS(0,0) and SPS (2,5) intensities, also showing that the FNS signal maximum precedes the SPS maximum, which is a common phenomenon caused by the different cross sections for electron impact excitation of the FNS and SPS upper states and the much shorter effective lifetime of the FNS compared to the SPS [37].The photomultiplier (PMT) measurements in  (2,4) intensity and a blend of FNS(0,0) and SPS (2,5) intensities, which are in the range of 388-395 nm.The short response time of PMT enables measurements with a time resolution of down to 100 ps, which is important for capturing small delays between emission curves from different excited states [37].The signals were measured using the slits of 250 µm.The instant of the FNS signal maximum in figure 13(b) was shifted to correspond with the FNS signal maximum in figure 13(a).Then, the onset of the PMT signals begin at 7 ns, and from this instant, the measured signals start to grow continuously.The maximum of the FNS/SPS blend is reached approximately at 12 ns and precedes the maximum of the SPS(2,4) signal, which is reached approximately at 13.5 ns.Both signals are scaled to 1.The FNS-narrow filter is defined by its transmission function: T = T 0 exp(−σ(λ − λ 0 ) 2 ) with parameters T 0 , λ 0 and σ being 33.6%, 391.5 nm and 1.05, respectively.The FNS(0,0) and SPS(2,5) blend needs to be resolved into individual components.Due to the much lower effective lifetime of the FNS(0,0) compared to SPS(2,5), we can consider that the tail of the FNS/SPS blend is created only due to the pure SPS(2,5) signal.The SPS (2,5) profile scaled by a factor 0.68 to fit the tail of the FNS/SPS blend is labeled using blue dashed line.The difference of the FNS/SPS blend and the SPS(2,4) signal scaled to fit the tail of the FNS/SPS blend, then provides the pure FNS(0,0) signal labeled using red solid line.
The dependence R(E/N) relating the intensity ratio and reduced electric field [19]: was used to obtain the reduced electric field.The I FNS,SPS are the measured intensities of respective spectral bands and τ FNS,SPS eff their effective lifetimes.At the investigated pressure of 200 Torr, the values of the effective lifetimes are τ B eff = (0.15 ÷ 0.24) ns and τ C2 eff = (3.1 ÷ 3.5) ns [34].The g FNS,SPS are factors denoting the ratio of the bandhead intensity and the integral band intensity of the respective FNS and SPS bands, which are used to account for the different spectrometric representations of FNS(0,0) and SPS (2,5) bands.These factors depend on the T rot and instrumental function parameter ∆ 1/2 .In our case ∆ 1/2 = 0.19 nm and T rot = 300 K, which, following the procedure suggested in [34], leads to the value gSPS gFNS = 0.75.The FNS(0,0) and SPS(2,5) ratio at the time of the FNS(0,0) maximum determined using kinetic series is 0.185.Therefore, the value of 0.185 × 0.75 = 0.139 was used to scale the normed PMT signals, leading to the The PMT waveforms of FNS(0,0)/SPS(2,5) blend and SPS(2,4) acquired with different bandpass filters.Blue dashed curve demonstrates the subtraction of the SPS(2,4) signal from the FNS(0,0)/SPS (2,5) blend to obtain pure FNS(0,0) signal.(c) Reconstruction of the electric field based on FNS(0,0) and SPS(2,5) intensity ratio.(d) Contributions of different components to the formula (14).Time t = 0 on the x-axis coincides with the onset of the superimposed HV pulse and timescales in (a)-(d) are identical as in figure 9. absolute values presented in figure 13(c).The E/N with its uncertainty band was determined based on the equation ( 14), where we used parameters a lower = 33.358,b lower = 100.1 and a upper = 13.0,b upper = 109.9,for lower and upper boundaries of the calibration range R(E/N).Note that the uncertainty springs from the uncertainty of the rate constants defining the R(E/N) relation [34].The computed E/N reaches the value of 380±100 Td in its peak, and the peak precedes the FNS(0,0) peak by approximately 2 ns.Note, that the instant of E/N peak does not match with the FNS(0,0) peak due to the nature of equation (14).Simply put, the increase in the electric field causes the increase of electron density; therefore, the peaks of FNS(0,0) and SPS signals follow with some delay after the peak of E/N. Figure 13(d) shows the contributions of the different components of the formula (14).It is apparent that the numerator of the formula ( 14) is defined mainly by the I FNS , since the short effective lifetime of N + 2 (B 2 Σ + u ) state causes instant equilibrium of the FNS upper state emission and population.Thus, the FNS derivative part in the formula ( 14) can be neglected (τ FNS eff dI FNS dt << I FNS ).In the denominator (SPS part), during the increase of the emission, the derivative component (τ SPS eff dI SPS dt ) dominates over the absolute value (I SPS ).Later, about 13.5 ns after the emission onset, the (I SPS ) prevails the derivative component, and the contribution of the derivative component is minor.Owing to the quite noisy nature of the I FNS for t > 15 ns, we restricted the reduced electric field determination just to t ∈ (8, 15) ns.

Implications for kinetic modeling
The kinetic model used is similar to the one in detail described in [12]; moreover, the electron processes from [11] were added to describe the build-up of the electron density.The model relies on the reduced electric field obtained as a fit of the E/N data determined in the previous section.The fit, labeled by the black dashed line, shown in figure 13(c) was performed using the sigmoid function suggested for streamer pulses in [11] and defined by the following formula: The parameters A, t 0 , t c , w 1 , w 3 were obtained by the fit and reached the values: 216.8 kV m −1 , 7.74 ns, 8.38 ns, 0.235 ns and 4.09 ns, respectively.The non-zero initial densities N 2 (A 3 Σ + u , v = 0), N 2 (A 3 Σ + u , v = 1) and N( 4 S) are set according to the experimental results of [16,18] to be 10 13 , 10 13 and 7•10 14 cm −3 , respectively, since, a certain density of  [16].Time t = 0 on the x-axis coincides with the onset of the superimposed HV pulse.these species is already built up within two breakdowns preceding the third breakdown during which measurements were performed [16,18].The initial electron density was defined to achieve a maximum electron density between 10 13 and 10 14 cm −3 .The electron density is related to the electric field pulse, increasing exponentially at the onset of the streamer pulse due to electron-impact ionization of N 2 .The loss of electrons in the model is due to electron-ion recombination, which is more efficient at lower E/N.Note that the 0D modeling cannot capture all of the experimentally observed features since some are also related to spatial characteristics during streamer propagation, where both electric field and electron density are strongly spatially dependent.
Figure 14 shows the comparison of the N 2 (A 3 Σ + u , w) density profiles obtained from the experiment [16] and model, using E/N(t) from figure 13(c) as input.Despite all the approximations involved in the calculations, a good agreement on the microsecond timescale is found.The non-monotonic behavior observed in figure 14 for N 2 (A 3 Σ + u , v = 5-9) has already been reported in [16] and is probably the result of two opposing processes: quenching due to vibrational relaxation and excitation by pooling process (16) followed by radiative quenching via SPS.Simply put, once the population of levels v = 0 and 1 of the N 2 (A 3 Σ + u ) metastable state reaches certain values due to vibrational relaxation, the pooling reaction refills the N 2 (B 3 Π g ) and N 2 (C 3 Π u ) states and their relaxation (collisionradiative) causes a local increase in the population of v = 5-9 levels of the N 2 (A 3 Σ + u ) metastable.Also, the time of the N 2 (C 3 Π u ,v = 0-4) peaks obtained by the model (solid lines in figure 15(a)) show good agreement with the experimental data (points in figure 9(b)).The decay of the N 2 (C 3 Π u ,v = 0-4) densities for t ∈ (10, 100) ns is much faster than in the experiment, which reports that the excited state source term (k C im n e n N2 ) in the model diminishes faster than in the experiment.However, to account for this effect, it would require the development of a high-dimensional fluid model, as in [38,39].For t > 100 ns, the 0D model predicts the increase of N 2 (C 3 Π u , v = 0-4) densities and their peak at t ≈ 6 µs due to the pooling from the N 2 (A 3 Σ + u , v = 0) and N 2 (A 3 Σ + u , v = 1) metastable states, see solid lines in figure 15(a).The disagreement of the model and experiment for t ∈ (60, 1000) ns is substantial since the experiment predicts continuous decay of the N 2 (C 3 Π u , v = 0-4) densities, which did not change even for t > 1000 ns, where we obtain spectra at very low signal-to-noise ratio.This is a sign of the possible role of the pooling reactions from the higher vibrational levels of the metastable state N 2 (A 3 Σ + u ).The pooling reactions can be written generally, for pooling from v and w vibrational levels of N 2 (A 3 Σ + u ), as: However, only rate constants for (v, w) combinations (0,0); (1,0) and (1,1) were determined [35,40] and are the only pooling processes included in the 0D kinetic model.The remaining (v, w) combinations were never inspected since it is very challenging to match the metastable levels participating in the pooling reactions with the specific N 2 (C 3 Π u , v ′ = 0-4) products.Moreover, the procedures based on reactants' and products' energy balance also fail; thus, it is difficult to estimate these rate constants [35].In our model, where we consider 22 vibrational levels of N 2 (A 3 Σ + u ), we could potentially arrive at n(n+1) 2 = 253 combinations of (v, w) for n = 22.This can be simplified if we realize that the source term of pooling reactions is proportional to the product, which can be analyzed for the different (v, w) combinations in our specific case of streamer afterglow, supposing similar pooling rate constant magnitude for all combinations.Figure 15 ] products, having the maximum higher than 3•10 27 cm −6 .It is apparent that the (v, w) combinations (1,1); (0,0) and (1,0) dominate for t > 2 µs, while the combinations (3,3); (3,2); (4,3) and (5,3) for t ∈ (0, 2) µs.To illustrate the potential influence of the pooling from the (v >1, w >1) level combinations, we have added to the model the important pooling process within (0, 2) µs time interval, (v, w) = (3,3): its pooling reaction rates k v ′ p reach values (1.0, 2.0, 1.0, 0.7, 0.6) × 10 −11 cm 3 s −1 for respective N 2 (C 3 Π u , v ′ = 0-4) vibrational levels, chosen to be the same as the value for (v, w) = (1,1) combinations suggested by Piper [35]. Figure 15(a) now compares the N 2 (C 3 Π u , v ′ = 0-4) densities from the model, where we do not consider (solid line) and consider (dashed line) the process (17).The agreement of the model with the experiment improves, which gives us the confirmation that the inclusion of the pooling from (v >1, w >1) level combinations is important for t ∈ (0, 2) µs in our case of N 2 at 200 Torr and could be important even for longer times for reduced N 2 pressures.

Discussion
We obtained fundamental information about the onset of luminosity and subsequent Townsend-to-streamer transition.
Thanks to the 4-channel ICCD framing camera, we were able to observe µs formation of diffuse phase covering the anode dielectric surface from which a single cathode directed streamer is born.The propagation of the streamer according to the velocity measurements takes approximately 3 ns and it is in agreement with the E/N peak, which has a rise time shorter than 2 ns according to the FNS/SPS intensity ratio measurements.After the streamer passes through the gap, a surface streamer develops on the dielectric cathode.Moreover, the voltage and current measurements performed in the section 3.1 inform us about the charge leak created by the discharge through the flow of the current to the conductive parts of the discharge chamber.The 14 ns delay between the onset of the current flowing to dielectrics (I C1 ) and current flowing to conductive parts of the chamber (I leakage ) is consistent with the results observed by the 4-channel ICCD in section 3.2, since the bottom row of the figure 6 for t > T 0 + 14 ns shows the light emission also from distant parts of the discharge chamber.This can be a sign of the charge transfer to the grounded chamber walls through the surface streamer developing in the vicinity of both electrodes.
Our experimental results can be compared with the modeling works of Yurgelenas and Wagner [38] and Papageorghiou et al [39].Yurgelenas and Wagner [38] developed a model of DBD in air at atmospheric pressure in the plane-to-plane geometry with the gas gap of 1.2 mm.In their model, they found two different modes during the discharge initiation: (i) the diffuse discharge and (ii) the ionization wave (streamer).Moreover, they observed a different velocities in these two modes.The streamer mode has approximately twice the speed of the diffuse discharge which is similar as in our case, where we determined velocities 9.8×10 5 m s −1 and 2.4×10 6 m s −1 in the anode and the cathode parts of the discharge gap, respectively.According to Yurgelenas and Wagner [38] the inclusion of the ion-electron emission from the cathode dielectric is sufficient to explain the formation of this type of the discharge and the change of the initial conditions in their model (residual surface charges, secondary emission coefficient from cathode) affects the prevalence of the diffuse or ionization wave character and shifts the boundary between these two modes.
In the work of Papageorghiou et al [39] the interaction of a volume streamer with a dielectric cathode is studied in N 2 at atmospheric pressure using a two-dimensional model.They consider plane-to-plane geometry with the gas gap of 1 mm and observe that the electric field increases from 120 kV cm −1 to 300 kV cm −1 during the streamer arrival on the cathode.These electric fields correspond to the values of 450-1100 Td, higher than the values determined in the present work (380 ± 100 Td).Papageorghiou et al [39] model the discharge just by incorporating the photoionization and photoemission from the cathode and they show that the photoionization is a key process also affecting the propagation of surface streamer on the cathode.They explain in detail several effects also observed in our case.The observed glow near the cathode is caused by the volume streamer impact and subsequent creation of surface streamer on the cathode dielectric.They determined the surface streamer velocity to be 4.5×10 5 m s −1 and the velocity was found to be decreasing with radial distance from the cathode spot.This value is close to the surface streamer velocity determined in our case, which is (2±1)×10 5 m s −1 .After the volume streamer contact with the cathode, a local minimum of the axial E/N magnitude is predicted behind the surface streamer layer.This phenomenon is also present in our experiment since we observe the emission extinction at the boundary of the volume and the surface streamers in figures 4 and 6.The volume streamer velocity determined in [39] is ∼7×10 5 m s −1 , close to the value determined in section 3.3.Papageorghiou et al [39] also observe that after the contact of the streamer with the cathode, the electron density in the middle of the gap still doubles.Such behavior is not present in our 0D model, but this could actually explain the slower relaxation of N 2 (C 3 Π u ) observed by the experiment compared to results of 0D model.

Conclusions
In this work, which is the first part of a two-part series, we have performed a time-resolved study of the PIE produced by a streamer monofilament in pure nitrogen at 200 Torr in the wavelength range between 200-500 nm.We combined a four-channel ICCD camera with an ICCD imaging spectrometer (used in both in zero-and first-order of diffraction) to obtain specific signatures of the evolving streamer filament from the very beginning up to the first microsecond.Using the identical point-plane DBD discharge geometry as in our previous works, the present results directly complement recent studies of the post-discharge evolution of N 2 (A 3 Σ + u , v = 0-10) metastable states [16] and N( 4 S) atoms [18].Here, we focused in detail on the evolution of the N 2 (C 3 Π u ) and N + 2 (B 2 Σ + u ) electronic states to reveal the dynamics of streamer monofilament formation and estimate the characteristic E/N responsible for the production of reactive species.Obtained results may be summarized as follows.
• Single-event images acquired with a four-channel ICCD camera revealed diffuse brightness (generated around the temporary anode on a microsecond time scale) before the onset of the streamer.After crossing the gap, the streamer further develops as a surface wave propagating on the dielectric surface of the temporary cathode from the streamercathode impact point followed with some delay by an analogous wave on the surface of the anode dielectric.• Statistical analysis of a set of single-event ICCD images obtained using imaging spectrometer in zero order of diffraction synchronously with PIE PMT waveforms allowed obtaining two characteristic velocities of the luminous wavefront (streamer head) of 9.8×10 5 and 2.4×10 6 m s −1 .• A detailed analysis of the SPS and FNS emissions revealed characteristic VDFs of the N 2 (C 3 Π u ) and N + 2 (B 2 Σ + u ) states in the middle of the DBD gap with nanosecond temporal resolution.In the case of the VDF of N 2 (C 3 Π u ) state, we clearly distinguish distributions driven by streamer-head electrons and later by metastable N 2 (A 3 Σ + u ) species.The initiation of the metastable-driven phase was also demonstrated using emission spectra of NO-γ and Goldstein-Kaplan emission.
• The value of VDF of N + 2 (B 2 Σ + u ) is lower than the FCFlike values (mainly for v = 1 level of the state).Therefore, a detailed vibrational analysis of FNS at different discharge pressures will be aim of our future study.
• The implementation of the measured E/N to the 0D model led to the finding that the inclusion pooling processes from N 2 (A 3 Σ + u , v > 1) vibrational levels leads to better agreement of the model with the experiment.

Figure 1 .
Figure 1.Experimental setup: simplified sketches of (a) DBD point-plane electrode geometry (b) high-voltage waveform and (c) complete stainless-steel reactor with associated optical diagnostics.

Figure 2 .
Figure 2. Block diagram of the experimental set-up with associated optoelectrical diagnostics.

Figure 3 .
Figure 3. (a) Equivalent electric circuit of the discharge cell consisting of lumped elements (the parameters of the elements are: R i = 60 Ω, L i = 2 nH, Cp = 1.6 pF, Rp = 100 MΩ, C 1 = 0.7 pF, C 2 = 0.03 pF, L 2 = 5 µH, T 1 and T 2 are characterized with a delay of 0.38 ns and the electrical impedance Z 0 = 70 Ω.Waveforms of measured driving voltage (red), measured current supplied from the voltage source (blue) and simulated current (green) without discharge (b) and with discharge (c).(d) shows filament current (blue), parasitic leakage current (magenta), and capacitor C 1 current.(e) Represents the time dependence of estimated filament (blue) and leakage (magenta) resistances.Time t = 0 on the x-axis coincides with the onset of the superimposed HV pulse.

Figure 4 .
Figure 4. XXRapidFrame image sequence obtained by integrating UV-vis discharge emission.The images show characteristic luminosity during pre-streamer (Ch1 and Ch2), streamer (Ch3) and post-streamer (Ch4) phases and acquired during a single discharge event by using specific MCP gates delayed appropriately with respect to the HV pulse.The MCP gates were triggered with respect to the onset of the HV pulse at times t = −7 µs, −2.05 µs, 0 ns, and 1 µs for Ch1, Ch2, Ch3, and Ch4, respectively.

Figure 5 .
Figure 5. XXRapidFrame image sequences obtained by integrating UV-vis discharge emission.Horizontal bars in (a) labeled as Ch1-Ch4 illustrate the timing of the Ch1-Ch4 MCP gates with respect to the onset (t = 0 ns) of the HV pulse.Images in panel (b) show the evolution of three selected single events during the HV pulse acquired by using equal Ch1-Ch4 MCP gates of 1 ns.

Figure 6 .
Figure 6.XXRapidFrame image sequence obtained by integrating UV-vis discharge emission.Images in sequence obtained as an average of five selected equivalent events acquired by using MCP gate (Ch2) of 1 ns to visualize evolution of the discharge morphology during the HV pulse.Image labeled as T 0 captures impact of the streamer on the cathode surface.

Figure 7 .
Figure 7. (a) Example of an ICCD image with a gate of 2 ns.(b) Streamer light intensity as a function of both the vertical position and the corrected gate time t gate corr .The time t gate corr is shifted to be null at the beginning of propagation.The color bar indicates the light intensity values.The positions of the top and bottom electrodes are shown using dashed lines.Furthermore, the streamer front positions are linearly fitted as a function of time in three different parts.The following functions were obtained: SF 1 (t) = −0.979t + 5.05 [mm], SF 2 (t) = −2.40t + 8.08 [mm], SF 3 (t) = −0.447t + 2.30 [mm].

Figure 8 .
Figure 8.(a) Characteristic UV spectra revealing bands of the SPS and NO-γ emissions, acquired in kinetic series acquisition mode with the MCP gate/step of 20 ns with respect to the onset (t = 0 ns) of the HV pulse, obtained using the 150 G mm −1 grating and the input slit width of 10 µm.(b) Characteristic PIE in the 410-480 nm spectral interval obtained using the 1200 G mm −1 grating and the input slit width of 50 µm.Emision was acquired during the post-discharge period (upper trace) and residual spectrum (lower trace) obtained after subtracting SPS bands belonging to the ∆v = −4, −5 and −6 sequences (indicated by red triangles); the residual spectrum reveals three Goldstein-Kaplan bands and two NO-γ (2nd diffraction order) bands.Both (a), (b) were measured with statistics of 5000 per spectrum.

Figure 9 .
Figure 9. (a) Characteristic UV spectra revealing bands of the SPS and FNS emissions, acquired in kinetic series acquisition mode with the MCP gate/step of 2 ns with respect to the onset (t = 0 ns) of the HV pulse.(a) Was obtained using the 1200 G mm −1 grating and the input slit width of 50 µm with statistics of 5000 per spectrum.(b) Time evolution of PMT and ICCD measurements of SPS and FNS bands.(c) Time evolution of N 2 (C 3 Πu, v = 0-4) states.The ICCD gate of 2 ns was used.The timescale coincides with figure 3.

Figure 9 (
Figure 9(b)  shows the SPS(0,0), SPS(2,4) and FNS(0,0) PMT profiles as well as the band intensities of SPS(0,3), SPS(1,4) and SPS(2,5) obtained by the spectrometer using the technique of kinetic series.These measurements can be converted to the time evolution of N 2 (C 3 Π u , v = 0-4) profiles as shown in figure9(c).Obviously, the onset of the ICCD intensities (1-8 ns) is more gradual than the onset measured using the PMT, which is due to the 2 ns gate width of ICCD used to collect spectra.Complementary to the preceding, figure10shows the time evolution of the VDF within 950 ns after the discharge onset.

Figure 11 .
Figure 11.Synthetic FNS (upper trace) and SPS (lower trace) emission in the 380-480 nm spectral interval, synthesized assuming the rotational temperature of 300 K, triangular instrumental function (∆ 1/2 = 0.2 nm) and FCF-like vibrational distributions of the N 2 and N + 2 emitting states.The emission intensities of bands in the 450-480 nm interval are multiplied by 25 and 20 for FNS and SPS, respectively.

Figure 12 .
Figure 12.De-coupling of FNS and SPS emissions in 405-435 nm and 435-480 nm intervals.Blends (a), (b) of FNS and SPS bands were obtained using a 1200 G mm −1 grating, an entrance slit width of 50 µm and statistics of 5000, ICCD gate of 2 ns delayed by 4 ns from the onset.The red triangles represent the FNS emission peaks, which are represented also in the upper trace of figure 11.Residuals (c), (d) were obtained by subtracting the ∆v = −4 and −5 SPS bands from (a) and (b), respectively.

Figure 13 .
Figure 13.(a) Peak intensities of FNS and SPS signals obtained from kinetic series with the gate of 2 ns and step of 2 ns.The studied time interval was from 0 to 20 ns.Similarly, as for the PMT signal, we observe that the FNS maximum precedes the maximum of the SPS.(b)The PMT waveforms of FNS(0,0)/SPS(2,5) blend and SPS(2,4) acquired with different bandpass filters.Blue dashed curve demonstrates the subtraction of the SPS(2,4) signal from the FNS(0,0)/SPS(2,5) blend to obtain pure FNS(0,0) signal.(c) Reconstruction of the electric field based on FNS(0,0) and SPS(2,5) intensity ratio.(d) Contributions of different components to the formula(14).Time t = 0 on the x-axis coincides with the onset of the superimposed HV pulse and timescales in (a)-(d) are identical as in figure9.

Figure 14 .
Figure 14.Time evolution of density profiles of N 2 (A 3 Σ + u , v) states obtained by the experiment (symbols) and the model (full lines) for even (a) and odd (b) v. Experimental points are reproduced from[16].Time t = 0 on the x-axis coincides with the onset of the superimposed HV pulse.