Open Access Article
A. Dorn
*a,
H. Mutafb,
R. Orhanb,
W. Wolff
c,
T. Pfeifer
a and
H. Ahmedovb
aMax Planck Institute for Nuclear Physics, Heidelberg 69117, Germany. E-mail: dornalex@mpi-hd.mpg.de
bTÜBİTAK National Metrology Institute (UME), Kocaeli, Gebze, Turkey
cPhysics Institute, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
First published on 7th April 2026
The influence of extraterrestrial particles like cosmic radiation (CR) on the chemistry and ozone density in the Earth's stratosphere is not well investigated and normally neglected in stratospheric chemistry models. Here we present the commissioning of a lab-based apparatus which aims at simulating conditions in the stratosphere in order to get better insight into the reactions induced by the secondary-particle showers from high-energetic CR which can reach low altitudes. Admixtures of ozone and the halocarbon CHClF2 (R22, chlorodifluoromethane) to atmospheric gases (N2, O2, Ar) were exposed to a glow discharge in the total pressure regime of a few hPa. According to the mass spectrometric analysis of the gas composition the discharge initiates significant ozone depletion by a factor four in the absence of R22. This depletion is strongly enhanced to two orders of magnitude in the presence of R22. The possible underlying reactions are discussed.
Environmental significanceAtmospheric gas mixtures at pressures in the low hPa range are exposed to glow discharges and analyzed by mass spectroscopy. The aim is to simulate stratospheric pressure conditions and to study how charged particle showers initiated by cosmic radiation can influence, e.g. the ozone concentration. A typical halocarbon R22 is introduced and its ozone depleting impact is observed. Presently the influence of extraterrestrial energetic particles on the chemistry and ozone density in the Earth's stratosphere is considered to be negligible and normally neglected in stratospheric chemistry models. On the other hand observations by Lu and Sanche [Q.-B. Lu and L. Sanche, Phys. Rev. Lett. 87, 078501 (2001)] suggest that ozone varies with the cosmic radiation flux. Respective research is relevant and timely. |
One question that is presently under debate is the influence of the energetic particle streams which originate from the Sun or from sources outside the solar system on the troposphere and biosphere. There is the solar particle stream (solar wind) from protons and alpha particles reaching energies up to roughly 500 MeV. Further, there is cosmic radiation (CR) which can have higher energies but has lower flux. Most primary and secondary particles are stopped in the ionosphere, which concerns essentially all solar particles. Only particle showers initiated by the more energetic CR (>1 GeV) can reach lower atmospheric layers like the stratosphere. Here ionization rates are up to 100 s−1 cm−3 and the ion density is around a few thousand cm−3. There are observations that CR influences atmospheric properties such as cloud formation and ozone in the stratosphere.1 In particular Lu demonstrated that stratospheric ozone concentration varies with the 11-year periodicity of the CR flux.2 The significance of these observations and the underlying mechanisms are under debate. Lu suggested a reaction model where anthropogenic chlorofluorocarbons (CFC) reaching the stratosphere are efficiently dissociated in reactions with secondary electrons from CR which are trapped on surfaces of stratospheric ice crystals.3
On the modelling side there are simulations of the atmospheric ionization density as, e.g., the AtRIS code4 which is a GEANT4 based simulation of CR induced atmospheric particle showers providing ionization densities as function of the incoming radiation energy and the altitude. This data can be the basis for modelling the atmospheric ion chemistry and the influence of CR showers on the relevant quantities like the ozone column density.
Lab based experimental studies on ionization of atmospheric gas mixtures are scarce. Cacace et al.5 have performed mass-spectroscopic ionization studies on mixtures of ozone and hydrochlorofluorocarbons (HCFCs) diluted in atmospheric gases. They found particular compound ions formed from ozone and HCFC which will be discussed below in relation with the present findings.
The abundant secondary electrons produced by the CR particle showers can be of importance for atmospheric ion chemistry. We aim at providing an experimental means to study the effects of electrons on ozone and the influence of HCFCs under conditions similar to the stratosphere. Therefore, we have set-up a laboratory apparatus consisting of a vacuum chamber that can be filled with controlled mixtures of gases at pressures between 1 and 4 hPa. A glow discharge (GDC) is used to provide a hot electron gas to induce energy transfer collisions like electron attachment, excitation and ionization processes and, thereby, initiate dissociation, radical formation and subsequent chemical processes modifying the gas concentrations.
Inelastic collisions and secondary electron production by CR are induced by high energetic charged particles within the particle showers. Their energy transfer in inelastic collisions with atmospheric molecules most likely is small. The secondary electron energy spectrum is maximal close to zero energy and falling off quickly for increasing energy6–8 such that for 10 eV secondary electron energy the cross section is about half the cross section for energies below 1 eV. The glow discharge electron energy spectrum also contains mostly low energy electrons in the few eV region which can be described as a non-Maxwellian distribution with a tail to higher energies beyond the ionization energies of the gas.9,10 While the CR induced electron energy distributions in the stratosphere and in a GDC are certainly not identical both are dominated by low energies below 10 eV. Therefore, the main processes induced in molecular gases like excitation, ionization and electron attachment all with subsequent dissociation and production of reactive species are similar.
In our experiment the evolving gas composition is analyzed by means of mass-spectroscopy. First studies were performed with the atmospheric gases (N2, O2) and admixtures of ozone and chlorodifluoromethane (CHClF2) which is the main CFC component in the stratosphere.11 We observe significant modifications of the gas concentrations and in particular a strong reduction of ozone.
An ozone generator (Fischer 500) is utilized to produce a small ozone admixture to the oxygen flow. During experiments the chamber pressure could be varied between 0.1 hPa and 4 hPa by adjusting the gas pumping speed. At the lower pressure range the chamber is pumped by a turbo-molecular pump with an oil-free scroll backing pump. An upstream adjustable throttle valve is used to regulate the pumping speed of the turbo pump and, therefore, the total chamber pressure. In the higher pressure range the chamber is pumped directly by the scroll pump. To achieve the typical experimental conditions with 2 hPa chamber pressure and gas flows of 100–200 sccm a pumping speed of 0.8–1.6 l s−1 is required. Therefore, it takes around 40–80 seconds for a complete turnover of the chamber volume of 63 l.
To induce electron-molecule collisions a glow discharge (GDC) is initiated between a cathode and an anode, both with 5 cm diameter and being 7 cm apart. Typical electrical voltages and currents at 2 hPa pressure are 700 V and 10 mA.
The plasma properties are determined by a commercial Langmuir probe (Impedans plasma measurement). Here a voltage biased tungsten wire-tip is introduced into the plasma. By analyzing measured current–voltage (I–V) characteristic curves generated as different voltages are applied to the probe the electron density, the electron temperature and the plasma potential are determined (see, e.g., 12 and 13). The electron gas parameters showed typical electron temperatures between 5 eV and 10 eV and electron densities in the range of (5–10)·107 cm−3.
In the present study we performed full mass scans from 5 to 84 amu for some specific gas mixtures and we measured time lapses tracking the concentrations of the main gas components while varying the GDC and R22 admission. In the latter case only particular fragment ion masses characteristic for the parent molecules of interest were tracked. These are listed in Table 1. E.g. R22 (CHClF2) shows various ionic fragments in the mass spectrum with the most abundant one being CHF2+ (m/z = 51 amu) with about 90% branching ratio (see NIST reference mass spectrum14). Ozone was identified according to the line for m/z = 48 amu corresponding to the O3+ ion which is produced with well-known cross section and branching ratio.15 For mixtures of ozone and R22 the contribution of the R22-fragment ion CH35Cl+ (m = 48 amu) with relative abundance of 0.7% of the CHF2+-intensity was corrected for. In order to obtain the pure O3+ ion signal this contribution was subtracted. Of particular interest in this study were plasma reactions involving R22 and ozone which can form compounds containing oxygen and chlorine. Stable molecules like Cl2O or Cl2O2 which reach the MS ion source were identified by the characteristic 37ClO+ ionic fragment (m/z = 53 amu). The 35ClO+ ion with m/z = 51 amu is superimposed and dominated by the CHF2+ fragment from R22.
| Gas species | Ion species used for identification | m/z (amu) | Fraction of fragment channel | Other ions with the same m/z |
|---|---|---|---|---|
| N2 | N2+ | 28 | 91% | — |
| O2 | O2+ | 32 | 91% | — |
| O3 | O3+ | 48 | 50% | CH35Cl+ |
| CHClF2 (R22) | CHF2+ | 51 | 90% | 35ClO+ |
| Ar | Ar+ | 40 | 95% | — |
| 37ClO | 37ClO+ | 53 | — | — |
Fig. 2 shows the mass spectrum for R22 at low pressure of 0.05 hPa with GDC on and off. The spectra for both cases are rather similar and can be mostly understood from dissociative ionization of CHClF2 in the ion source of the MS plus some residual gas lines from N2, O2 and H2O. Also heterogenous reactions on surfaces, e.g., with adsorbed molecules like water or the admitted gas species can evolve. Significant changes for GDC on are the increase of hydrogen chloride HCl (m/z = 36, 38) and HF (m/z = 20) which are direct dissociation products of R22. The main ionization channel results in CHF2+ (m/z = 51) and neutral Cl. At sufficient chamber pressure Cl can form Cl2 in three-body collisions Cl + Cl + M → Cl2 + M. Here M is a third collision partner taking energy to stabilize Cl2. Despite the low pressure of 0.05 hPa the respective ions are found at m/z being 70, 72 and 74 with their characteristic ratios of isotopic line intensities I(35Cl35Cl+)/I(37Cl35Cl+)/I(37Cl37Cl+) = 9/6/1. Also doubly charged 35,37Cl2+. Show up at mass over charge ratios 17.5 and 18.5.
In Fig. 3 a mass spectrum for 400 sccm flux of O2, 7 sccm O3 and 11 sccm of R22 resulting in a chamber pressure of 2.4 hPa is shown. Again m/z spectra with and without GDC are shown and for comparison also the m/z-spectrum without R22 and GDC off. While the flux of R22 is the same as for Fig. 2 the Cl2 production becomes more significant due to increased probability of three-body collisions at the higher pressure of 2.4 hPa. Furthermore, reaction products from R22 dissociation products and oxygen/ozone are formed. This is indicated by 37ClO+ at m/z = 53 (35ClO is merged with CHF2+ at m/z = 51) as well as F2O+ (m/z = 54)16 and possibly CFO+ at m/z = 47. Due to their reactivity these radicals will barely reach the MS but react quickly in the plasma. E.g., ClO can bind to another ClO or to a chlorine radical and form the stable molecules Cl2O2 or Cl2O, respectively. When these molecules reach the MS ion source the observed ClO+ ions and other ionic species are formed (for Cl2O, see the NIST reference mass spectrum17). As in Fig. 2 the GDC increases the HCl+ and HF+ intensities. Ozone as seen in the intensity for O3+ at m/z = 48 is significantly decreasing with R22 flux and with GDC on. Without GDC there is a very small reduction in this line intensity for R22 switched off which is due to a small contribution from the R22 fragment CH35Cl+ (m/z = 48) produced in the MS ion source.
In the following we focus on the behavior of the main gas components in time lapses where we adjust the different gas fluxes. In Fig. 4 an exemplary experimental sequence is presented. The ion yield currents taken with the MS are shown as a function of time. The initial gas flows were N2: 200 sccm, O2: 100 sccm and O3: 1.5 sccm. The shown sequence is an excerpt of a longer measurement where the R22 flux was switched of at time t = −80 s. Therefore, the CHF2+ current initially still decreases with time. At t1 the GDC is started (voltage 900 V, current 10 mA). The electron gas temperature was around 6.5 eV and the electron density about 5·107 cm−3. A significant reduction of O3 concentration by a factor of four is observed. In the same time a small increase of the O2+ signal of 1.6% is found. After switching off the GDC at t2 the O3+ ion current recovers to the original value. Subsequently at t3 the CHClF2 gas is admitted with a low flux of 11 sccm which is 3.7% of the combined N2 and O2 fluxes. At t4 the glow discharge is switched on and this results in a strong decrease of ozone by more than two orders of magnitude. Again, there is a small increase of the O2+ signal by 1.1%. In the same time also the CHClF2 concentration is reduced by about 15%. Then at t5 the R22 gas flux is stopped allowing the ozone ion yield to recover to the original value for GDC operation which was also found between t1 and t2. Finally, when the GDC was stopped at t6 the ozone concentration settles at the starting value. Additionally, an important observation is the course of the ion with M = 53 amu. In the stratosphere the relevant mechanism for ozone depletion by CFCs is the release of chlorine initiating the odd chlorine catalytic cycle which involves ClO as the first product (see discussion below). Therefore, we examined mass 53 corresponding to 37ClO (25% abundance) while 35ClO (75% abundance) is hidden under the dominant peak with M = 51 from CHF2+. Indeed, there is significant increase in 53 amu yield under plasma operation if both O3 and R22 flux are on. This leads us to the tentative conclusion that M = 53 amu is a reaction product of both molecules where at least one fragment from an electron induced dissociation is involved.
In the stratosphere important ozone depleting reactions involve the odd nitrogen family NOx (defined as the sum of N, NO, and NO2).18 Potentially such species could be produced also in the present N2/O2 mixture as a result of the GDC. In order to test if reactions involving ionization or dissociation of N2 play a role in the observed ozone depletion another experiment was conducted with N2 replaced by argon (Fig. 5). Initially Ar, O2, O3 and R22 gas fluxes were on and at t1 the discharge was started giving rise to a decrease of ozone by two orders of magnitude. At t2 R22 was shut off resulting in ozone increase. At t4 and t5 the GDC is switched on and off again, respectively, showing ozone decrease by a factor of 3.5. The dependence of the ozone concentration on R22 flux is demonstrated during discharge operation t7–t10 where R22 is briefly increased between t8 and t9. This results in a further decline of ozone.
These sequences demonstrate the strong influence of electron collision induced processes on the ozone concentration.
![]() | (1) |
| Ozone destruction O3 + hv (λ < 320 nm) → O + O2 ⇒ O3 + O → 2O2 | (2) |
If R22 is introduced a number of new fragments are produced from collisions involving free electrons. The possible reactions are dissociative electron attachment (DEA),20 dissociative excitation21 and dissociative ionization.22 In all these processes the most likely product is the chlorine atom or the negative chlorine ion together with the molecular fragment CHF2(+). For ionization we have measured the absolute partial ionization cross sections for various CFCs including R22 from low to high impact energies.23 The dominant fragmentation results in CHF2+ + Cl which is also the channel used for detection of R22 according to Table 1. Chlorine in the stratosphere is very efficient in ozone depletion which evolves in catalytic cycles. In the odd chlorine cycle the following reactions take place.
| Cl + O3 → ClO + O2 | (3) |
| ClO + O → O2 + Cl. | (4) |
Another ozone-depleting reaction in the stratosphere is the odd nitrogen cycle involving NO and NO2:
| NO + O3 → NO2 + O2 | (5) |
| O + NO2 → NO + O2 | (6) |
Another ozone depleting mechanism was found in a previous study by Cacace et al.5 using advanced analytical and mass spectroscopic methods. The authors have found reactions between ionized ozone and CFCs which in case of R22 are:
| O3+ + CHClF2 → CHF2O3+ + Cl | (7) |
| CHF2O3+ → HF2O2+ + CO | (8) |
CO with M = 28 has the same mass as N2 such that we cannot separate it with the present MS. Nevertheless, in the mass scan of Fig. 3 for a mixture of 400 sccm O2, 7 sccm O3 and 11 sccm R22 we find that the initial M = 28 amu peak intensity due to residual N2 gas increases for GDC operation by around 50%. The M = 28 amu peak is unchanged with/without discharge if R22 is not introduced. Therefore, a certain quantity of CO is produced from reactions like eqn (7) and (8). In the full mass scans we found indications of more reaction products from ozone and R22 which are initiated by the discharge like F2O+ (M = 54 amu). This supports the findings of Cacace et al. that there is a rich chemistry between ions and neutrals of oxygen, ozone and F22.
In order to conclude on processes in the Earth's stratosphere several aspects have to be considered. In the stratosphere the ionization density by CR is much lower than in the glow discharge, namely around 10–100 cm−3 s−1 and the resulting electron-ion pair density is around 103–104 cm−3. In comparison the electron-ion density in the present experiment is above 107 cm−3. Therefore, the present experiment strongly enhances the collisional reactions in the gas mixture. On the other hand, the catalytic ozone depleting cycles are running on longer time scales and cannot contribute significantly to modify the ozone concentrations due to the fast gas exchange rate in the present experiment. From these observations possible future studies should aim at higher chamber pressure of 10 to 20 hPa and lower gas turnover rates to better mimic stratospheric conditions. As result of these conditions significantly lower halocarbon concentration would be required to influence ozone concentration.
Possible future studies should analyze the pressure dependence of the observations. E.g. at higher pressure the discharge-induced ozone depletion can be compensated by ozone production in three-body processes O + O2 + M → O3 + M. Also, the gas exchange rate can be reduced in order to enhance the role of the repeating catalytic cycles with respect to the initial collisional processes. In this way the natural stratospheric conditions can be better approached and extrapolations to lower halomethane concentrations and lower electron densities can be performed. In addition, modelling of the atmospheric chemical processes starting from the various collisional reactions initiated by low energetic electrons would be desirable.
| This journal is © The Royal Society of Chemistry 2026 |