Shurik Yatom methane decomposition in a short glow discharge at Princeton, NJ; H atoms peak at cathode; polymerization zone farther

Shurik Yatom methane decomposition in a short glow discharge at Princeton, NJ; H atoms peak at cathode; polymerization zone farther

1

Spatially-resolved methane decomposition in a short glow discharge: insights into suprathermal hydrogen and radical chemistry Shurik Yatom Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ 08543, United States of America Abstract This work presents the spatially resolved mapping of methane decomposition chemistry in a short direct- current glow discharge, revealing segregated zones for dissociation and polymerization in an Ar-CH4 mixture at 400 mTorr. Laser-induced fluorescence (LIF), two-photon absorption LIF (TALIF), and optical emission spectroscopy (OES), are used to map the absolute number densities of atomic hydrogen (H), methylidyne (CH), and dicarbon (C2) radicals. The results reveal a segregated chemical environment driven by the non-local electron kinetics. The primary dissociation of methane is confined to the cathode sheath (< 4 mm), where the H atom density peaks at a value (~1.8×1021 m-3). Analysis of the Hα line profile identifies two distinct suprathermal H atom populations with peak kinetic energies of ~110 eV and ~17.5 eV, attributed to ion-surface reflection and electron-impact dissociation, respectively. The densities of secondary radicals CH and C2 peak further from the cathode at the sheath-negative glow boundary (y≈ 4-6 mm), identifying this region as the primary zone for polymerization, with the peak C2 density (~5.6×1017 m-3) higher than that of CH (~1.5×1015 m-3) by a factor of ~370. The discharge operates under strong non-equilibrium conditions, confirmed by the disparity between the low bulk gas temperature (~570 K) and the high vibrational temperature (3300-4700 K) of the emitting CH(A) radicals, the latter serving as a signature of a specific high-energy dissociative excitation pathway. This work demonstrates how the structured energy landscape of a short glow discharge spatially separates methane dissociation from subsequent polymerization. Kinetic estimates suggest that suprathermal H atoms may play a role in bulk plasma chemistry via abstraction reactions, with predicted CH4 activation rates exceeding those of the thermal H population.

1. Introduction Methane (CH4), the main component of natural gas, is a key molecule in energy and chemical research due to its abundance and environmental impact. It is a vast hydrocarbon resource, making its efficient conversion into value-added products a paramount scientific and technological challenge1. Converting methane into hydrogen (H2) and higher-value hydrocarbons (e.g., ethane, ethylene, acetylene) or carbon nanomaterials is a key research focus2–7. Traditional conversion methods, such as thermal pyrolysis or steam reforming, are energy-intensive8, typically requiring elevated temperatures (often exceeding 1000 °C) and complex catalytic systems to overcome the high stability of the C-H bond. Non-equilibrium, or “cold,” plasmas offer an alternative approach. In these weakly ionized gases, an external electric field energizes light electrons, creating a high-energy electron population while the bulk gas (ions and neutrals) remains near ambient temperature. This deviation from thermodynamic equilibrium enables electron-impact collisions to activate stable molecules like CH49, bypassing high gas temperatures and enabling unique chemical pathways unavailable in thermal processes10.

2

Glow discharge plasmas, a key type of non-equilibrium plasma, are used in processing11, lasers12, analytical chemistry13,14, water treatment15, gas conversion16 and surface cleaning17, including the surfaces on nuclear fusion reactors18. Methane-containing glow discharges are particularly relevant for depositing diamond-like carbon (DLC)19–22. This study focuses on a direct-current (DC) “short” glow discharge, where a small inter-electrode gap prevents the formation of a positive column (PC), creating a unique plasma environment with steep gradients in the electric field and electron energy. The DC glow discharge, a self-sustaining plasma generated between two electrodes in a low-pressure gas, provides a convenient system for studying fundamental plasma phenomena. Its operation is predicated on a cascade of events initiated by the acceleration of free electrons, leading to ionizing collisions. The resulting positive ions are accelerated toward the cathode, where their bombardment induces the secondary emission of new electrons, which in turn sustain the discharge. A fully developed glow discharge exhibits a complex, stratified structure of luminous and dark regions along the inter-electrode axis. These regions include the cathode glow (CG), cathode dark space (CDS, also known as the Crookes or Hittorf dark space), the intensely bright negative glow (NG), the Faraday dark space (FDS), and, if the electrode gap is sufficiently long, a quasi-uniform PC.23,24 This work uses a configuration known as a “short” glow discharge, defined by an inter-electrode distance that is insufficient to support the formation of a PC25,26. The use of large, parallel-plate electrodes allows for the generation of a discharge that is laterally homogeneous and provides clear, unobstructed optical access for diagnostic techniques, a critical requirement for spatially resolved measurements. Finally, short discharges operating at low pressures (e.g., on the left-hand branch of the Paschen curve) typically exhibit a rising current-voltage (I-V) characteristic, which ensures stable, homogeneous operation that is free from the instabilities that can complicate experiments in other discharge regimes. Unlike PC or dielectric barrier discharges, where volume kinetics are often averaged, the short glow discharge regime creates a unique kinetic environment. It enables the spatial isolation of high-energy surface interactions (in the sheath) from volumetric ballistic transport (in the negative glow), allowing for the distinct identification of chemical drivers that are typically obscured in well-mixed reactor geometries. This paper presents an experimental investigation into the spatially resolved chemistry of methane decomposition within a low-pressure (400 mTorr), large-area, short DC glow discharge sustained in an argon-methane mixture. This study employs optical emission spectroscopy to characterize the plasma. Laser spectroscopy, i.e., planar Laser-Induced Fluorescence (LIF) and Two-Photon Absorption LIF (TALIF), is employed to map the axial distributions of key reactive species: atomic hydrogen (H), methylidyne (CH), and dicarbon (C2) radicals. The primary objective is to leverage the unique, structured environment of the short glow discharge to deconstruct the complex reaction network of methane. The choice of these specific species is strategic for untangling the spatially dependent chemistry. • Atomic hydrogen (H): As a direct product of the initial C-H bond breaking in methane, the spatial distribution of H atoms, measured via TALIF, serves as a direct probe of where the primary decomposition of CH₄ occurs. This allows us to distinguish between regions dominated by high-energy electron impact (expected in the cathode sheath) and those driven by energy transfer from argon metastables.

3

• Methylidyne radical (CH): This radical is a secondary or tertiary product of methane decomposition. Comparing its spatial profile, obtained via LIF, to that of atomic hydrogen provides crucial information about the subsequent steps in the chemical reaction chain and the transport of reactive species from their point of origin. • Dicarbon radical (C2): The C2 radical is a key precursor to the formation of higher-order hydrocarbons and soot27,27. Its spatial distribution, measured by LIF, acts as an indicator for the onset of polymerization chemistry28, revealing where in the discharge the pathways shift from simple fragmentation to C-C bond formation. By correlating the distinct spatial profiles of these primary (H), intermediate (CH), and polymer-precursor (C2) species with the non-local properties of the discharge, we aim to provide a fundamental understanding of how the spatially variant Electron Energy Distribution Function dictates the reaction pathways. Furthermore, the detailed spatial maps of these radicals will provide critical data to benchmark and validate kinetic models. The results will also guide future research by pinpointing regions of specific chemical activity. 2. Experimental setup and methods

Figure 1. Experimental setup The experiments were conducted in a stainless steel, 6-way cross chamber (8-inch diameter) with three fused silica windows for optical access and a fast-entry door with a window (Figure 1). The discharge was sustained between two parallel, coaxial electrodes: a 20 mm diameter copper anode, mounted on the top flange via a Wilson seal for vertical adjustment, and a 40 mm diameter brass cathode, fixed 17 mm below the anode and grounded. The anode, isolated by a quartz sleeve to prevent spurious discharges, was connected to a DC power supply (Glassman, 1 kV, 30 mA) through a 1 MΩ ballast resistor. The discharge current was 0.736 mA, calculated from the voltage drop across the resistor. The chamber was evacuated to a base pressure of 30 mTorr using a rotary pump, then filled with a continuous flow of argon (Ar) and methane (CH4) mixture (85% Ar, 15% CH4) at 30 sccm and 5 sccm, respectively, regulated by mass flow controllers (Alicat Scientific). To minimize background gas contamination, the chamber was purged with the working gas mixture for at least 30 min. prior to data acquisition. During operation, the continuous flow of 35 sccm maintains a steady-state working pressure

4

of 400 mTorr, continuously displacing residual species. The impact of the base pressure on discharge chemistry was assessed via optical emission spectroscopy: no emission lines attributable to N₂, NO, or OH were detected in the OES spectrum under discharge conditions. Furthermore, no Hα emission or atomic hydrogen signal was detected when operating the discharge in a pure Ar flow under otherwise identical conditions, confirming that any hydrogen-containing residual gas has a negligible contribution to the observed H atom production and chemistry. The working pressure was monitored by a Baratron capacitance manometer. The reactor’s operation deposits thin carbonaceous films on the electrodes, requiring frequent cleaning with fine sandpaper and ethanol. Two optical emission diagnostics were applied: • Filtered imaging: An iCCD camera (PI-MAX 4, Princeton Instruments) with a 656.3 nm filter (1 nm FWHM) captured 2D Hα emission to estimate the discharge’s axial current cross-section. • Spectral imaging: The plasma was imaged onto a monochromator’s entrance slit (HRS Spectra Pro 750, Princeton Instruments), coupled to an iCCD camera, producing spectra with axial spatial resolution29. CH(A-X) emission spectra were fitted using LIFBASE software to calculate rotational (Trot) and vibrational (Tvib) temperatures. Calibration used a Hg-Ar lamp (Oriel 6035). LIF and TALIF measured species densities: • CH and C2 LIF: A 20 Hz laser system (frequency-tripled Nd:YAG laser pumping an optical parametric oscillator at 355 nm) was used, as described in prior H-TALIF30 and C2-LIF studies31. • H-Atom TALIF: A 100 Hz diode-pumped laser (NT230, Ekspla) enabled faster signal collection for planar H-TALIF. The laser beam was expanded (×3.75) using a two-lens telescope (f = 40 mm, f = 150 mm) and focused into a planar sheet (height: ~19.5 mm, waist: ~521 µm) using a cylindrical lens (f = 250 mm). Beam parameters and spectroscopic transitions are listed in Table 1. Fluorescence was imaged using an iCCD camera with a 3rd Gen HBf intensifier, an objective lens, and bandpass filters to isolate signals. Laser fluence (2.5-2.8 mJ/cm²) was kept low to avoid power bro