---
title: "Shurik Yatom methane decomposition in a short glow discharge at Princeton, NJ; H atoms peak at cathode; polymerization zone farther"
sdDatePublished: "2026-09-16T04:09:00Z"
source: "https://arxiv.org/pdf/2609.17345"
topics:
  - name: "science and technology"
    identifier: "medtop:13000000"
locations:
  - "Trenton"
---


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