1. Introduction
Radiation effects in solid methane and methane-containing frozen films, so-called ices, attract much attention in diverse fields of science and technology. Being widely present in the Universe [
1,
2,
3,
4,
5,
6] methane plays important role in astrophysics and astrochemistry [
7]. An overview of laboratory experiments on the relationship of methane to prebiotic chemistry is presented in [
8]. Properties of methane, their transformation caused by the continuous radiation of outer space, and physico-chemical processes became the topic of a large body of research. Numerous efforts have been made to simulate radiation-induced processes in the laboratory using a variety of techniques surveyed in Rev. [
9]. Various types of ionizing radiation were used to study radiation behavior of methane and methane-rich ices – ions [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19], electrons [
20,
21,
22,
23,
24,
25,
26,
27,
28,
29] and photons [
26,
30,
31,
32,
33,
34,
35].
Interaction of solid methane with neutrons is closely connected to the technology of cryogenic moderators [
36,
37,
38,
39,
40,
41,
42]. Solid methane is widely used for moderating hot neutrons to cold and ultra-cold ones because its efficiency is approximately 3.5 times that of the commonly used liquid hydrogen-based moderator [
36]. Despite the attractive neutronic properties of solid methane its use in cryogenic moderators faces a problem especially for the high intensity neutron flaxes. After long exposure to neutron radiation, as it was found by Carpenter [
36,
37], solid methane moderator experienced sudden, violent, warm-ups and pressure rise in its vessel in some cases with followed vessel destruction. This phenomenon, called by „burp“ effect, significantly constrains the current use of solid methane moderators. The reason of this phenomenon according [
36,
37] is a large spontaneous release of energy. The moderator material stores a part of energy absorbed from radiation in view of frozen-in products of radiolysis. At high density of reactive species their recombination and the expansion of hydrogen, which builds up in the solid methane during irradiation by neutrons [
36] result in fast energy release which destroys a moderator vessel. The pressing character of this problem stimulated the thorough study the interaction of neutron radiation with solid methane [
37,
38,
40,
41,
42]. According to Carpenter’s model [
37] all recombination processes in solid methane moderators proceed in one stage and are controlled by the same activation energy for defect diffusion. However, commissioning tests of the ISIS TS2 solid methane moderator, in which a „burp“-like effect was observed [
40], and the study of relaxation phenomena in solid methane pre-irradiated with an electron beam [
28] led to the conclusion that the radiolysis defect recombination process happens in two stages, at different temperatures, and is therefore controlled by two different activation energies. Modifying Carpenter’s approach Kirichek and coauthors implemented two different types of radiation defects (H atoms and CH
3 radicals), with different recombination rates and thermal activation energies [
40]. The model suggested satisfactory described results of the commissioning tests the ISIS TS2 solid methane moderator. The possible role of radiolysis product recombination processes in cryo-volcanism on comets was discussed in [
37,
40].
Strong explosive-like delayed ejection of particles was observed in solid methane exposed to an electron beam of the subthreshold energy [
29,
43] that is, under conditions when the impact mechanisms of defect formation and desorption are excluded. The effect was observed at low temperature upon reaching a critical irradiation dose of 100 eV per methane molecule and was accompanied by the sample temperature rise and a flash of the luminescence. The outburst of particles was preceded by oscillations of particles ejection with increasing amplitude. The period of these oscillations depended on the intensity of the irradiation, it decreases with an increase in the beam current density. In [
44] a model has been proposed that qualitatively describes the appearance of two types of self-oscillations upon electronic excitation of solid CH
4, their periods and delay time of burst. Two types of self-oscillations found are a periodic change of temperature and concentration with time of CH
4 decay products – H atoms and CH
3 radicals upon irradiation. The found patterns of self-oscillations was shown to determine the temporal dynamics of the delayed explosive particle ejection upon irradiation with the subthreshold electrons. The discovered effect is a new manifestation of self-oscillations similar to those considered in [
45]. To the best of our knowledge, there were only a few studies in which a delayed in time burst of particles, the so-called “delayed desorption”, from solid methane exposed to other types of ionizing radiation, was observed [
46,
47]. In [
46] solid CD
4 was irradiated with MeV He
+ and H
+. When the beam was turned on, weak desorption of D
2 first appeared. Upon reaching the threshold fluence density, the yield of D
2 rapidly increased by more than an order of magnitude, and then gradually decreased. The threshold fluence values for irradiation with 1.5 MeV He
+ ions and 1.5 MeV H
+ were 3 10
14 ions cm
-2 and 9 10
15 ions cm
-2, respectively. In [
47], solid CH
4 and CD
4 were irradiated with 9.0 MeV α particles and 7.3 MeV protons. Doses up to 145 eV caused pressure shocks rising the chamber pressure by several orders of magnitude. This process released up to 90% of the molecules of the solid target into the gas phase. The authors concluded that the H and CH
3 radicals play a major role in this phenomenon, and the observed explosion resembles a nonequilibrium process. However, in these experiments no regular oscillations were observed. The detection of complex organic molecules (COMs) in the gas phase of cold molecular clouds at such low temperatures of the cryogenic environment that they would have to condense on the surfaces raised the question of the reasons of this phenomenon. Several hypotheses have been proposed [
48,
49,
50,
51,
52] and radical–radical induced explosive desorption of ice-coated interstellar nanoparticles among them [
52]. In this article, the authors present the direct observation by Fourier transform infrared (FTIR) spectroscopy of rapid radical–radical reactions of formyl (HCO•) and methyl (•CH3) radicals and the reaction-induced explosive desorption during the exposure of methane and carbon monoxide ices to superthreshold 5 keV electrons. However, the formation and recombination of the H atom remained beyond the reach of this experimental approach.
To get more insight into the phenomenon the matrix isolation method was introduced in combination with a set of emission spectroscopy methods [
53,
54,
55]. It should be noted that in contrast to pure methane and mixtures of methane with other molecular gases of astrophysical interest there are only a few studies performed in matrices of rare gases [
56,
57,
58]. Photolysis of methane in Kr matrix was studied using synchrotron radiation [
56]. Strong VUV absorption and emission features were attributed to atomic carbon C and the CH radical. Neither CH
2 nor CH
3 were observed in Kr matrix. In contrast to that photolysis of methane in Ne matrix resulted in appearance of wide range of hydrocarbons including CH
3 and carbon clusters C
n with number of atoms up to n=20 [
58].
Below, we briefly summarize the results of our previous studies of radiation effects induced by an electron beam in Ar matrices doped with CH
4 [
53,
54,
55,
59]. In the cathodoluminescence (CL) spectra, the following products of radiation-induced methane transformation were recorded: H, CH and C. It has been shown that because of the small penetration depth of electrons the bulk of the matrix is excited preferentially by photons of the matrix with an energy of 9.8 eV (the most intense emission band of the self-trapped excitons), in other words “internal photolysis” occurs. Analysis of an efficiency of different channels of methane transformation under these photons and neutralization reactions indicated that the CH radical can be considered as a signature of the CH
3 radicals [
55]. Monitoring of the so-called nonstationary luminescence (NsL), viz. luminescence under nonstationary conditions (upon external heating under beam) in combination with the detection of thermally stimulated exoelectron emission (TSEE) and optical emission spectroscopy made it possible to reveal the contribution of the reactions of charged and neutral species. The dynamics of the main products of methane fragmentation were traced and it was found that behavior of CH radicals (therefore, CH
3 radicals) and hydrogen atoms is different and depends on the concentration of methane in the matrix. This means that the assumption about the same concentration of CH
3 radicals and H atoms based only on the reaction: CH
4 + Δ
E → CH
3 + H, accepted in the delayed explosive desorption models [
37,
40,
44], is not entirely correct. Long-term and short-term oscillations of the total particle yields were found and interpreted as thermo-concentration self-oscillations similar to those considered in [
44] for the case of pure methane. It has been shown that they can be initiated either by the spontaneous release of stored energy upon reaching a critical concentration of reactants, or by external heating of the doped matrix to stimulate diffusion and release the stored chemical energy in subsequent recombination reactions [
59]. Preliminary measurements of the dose dependence of the optical emission of the H atom, carried out on a lightly doped Ar matrix (C = 1%) at LHe temperature, found a puzzling fact – an increase in this emission with a pressure burst [
54]. The dynamics of optical emissions from highly doped matrices upon low temperature irradiation remained unexplored.
Here we present an outgrowth of this research with a focus on H atom and CH
3 radical behavior at different concentrations of CH
4 in Ar matrix and processes induced by energy transfer from the matrix to dopant and products of radiolysis. Three channels of energy transfer are considered – charge transfer (by holes), transfer by free excitons, and self-trapped excitons. Special attention is paid to the recombination of H atoms and its monitoring in the matrix. This reaction proceeds nonradiatively with energy release of 218 kJ mole
-1. A hydrogen molecule can be monitored by being transferred to an excited state. The thresholds of excitation of the lowest singlet (B
1 Σ
u+) and triplet (a
3 Σ
g+) excited states [
60] are below the bands of free excitons of Ar – 12.06 eV for
Γ(3/2) and 12.24 eV for
Γ(1/2) [
61]. So, the transfer of energy by free excitons to the hydrogen molecule can lead to its transition to the lowest excited states B
1 Σ
u+ and a
3 Σ
g+. Electric dipole transitions from the coupled triplet state a
3 Σ
g+ into an unstable lower state b
3 Σ
u+ lead to the appearance of a wide emission continuum [
62] followed by dissociation of H
2 molecule. The detection of H atoms proceeds under conditions of competition between exciton self-trapping and energy transfer to the dopant and radiolysis products. A nonlinear concentration behavior of the optical emission of H atoms in doped Ar matrices has been found. Real-time correlated monitoring of optical emissions (H atom and CH
3 radicals), particle ejection, and temperature revealed a nonmonotonic behavior of optical yields with a strong luminescence flash after almost an hour of exposure, which correlated with surge of temperature and particle ejection. The connection of this phenomenon with the processes of energy transfer and radical-radical recombination is discussed.
3. Results and Discussion
As already mentioned, the electronic states of the methane molecule are dissociative, so its luminescence spectrum is determined by the products of its radiolysis. Main channels of CH
4 degradation upon photon flax with the threshold wavelengths [
69] in the solid Ar emission range are:
In accordance with this we detected following emissive products of methane degradation: H, CH, C and H
2, as can be seen from the CL spectrum shown in
Figure 3.
The spectrum contains both the emission bands of the matrix and the emission bands of the radiolysis products. Emissions are observed both from the solid phase and from the gaseous due to the electronically induced desorption. The most intense feature of this spectrum at 127 nm is the well-known emission of the molecular-type self-trapped excitons (STE) – Ar
2*, corresponding to the bound-free transitions from
1,3Σ
u+ states to a repulsive part of the ground state
1Σ
g+ [
61]. At the short-wavelength edge of the Ar
2* band, recorded in the second order, the Lyman-ɑ line, which belongs to the desorbed into the gas phase excited H atoms, was detected. Mechanisms of the excited H atoms formation and their desorption from Ar matrices, doped with CH
4, were discussed recently [
71]. Four scenarios for the desorption of excited H atoms were proposed, including energy transfer by holes and excitons of the Ar matrix. Own atoms of the matrix and molecules also desorb in excited states emitting corresponding lines and band (the feature with unresolved structure in
Figure 3). Their desorption is driven by the exciton mechanism [
72]. A necessary condition of the electronically stimulated desorption is the localization of electronic excitation on the atoms or molecules in the near-surface region with the release of energy exceeding the binding energy per atom ε
b (for atom of Ar matrix ε
b=88.8 meV [
61]). For the case of Ar matrix desorption of excited particles (including H
* atoms) is facilitated by negative electron affinity of Ar χ = -0.4 eV [
61], in other words, excited particles are ejected under the repulsive forces between an excited atom and atoms of surrounding lattice, so-called “cavity-ejection” mechanism [
72]. Interestingly, the desorption of excited argon molecules in the states
1,3Σ
u+(v’) was also detected as a result of self-trapped hole (Ar
2+) neutralization reaction [
73]. The so-called “third continuum” of Ar at 200 nm, assigned to the radiative transitions from the excited (Ar
2+)
* state [
74], was quenched by admixture of methane. Note that the Ar
2*-band recorded in the second order overlaps with a wide emission continuum of molecular hydrogen associated with transitions from the bound triplet state a
3Σ
g+ to the unstable lower state b
3Σ
u+, which makes it difficult to monitor the H
2 molecule by this transition. According to [
75], the maximum of the spectral distribution of the radiative transition probability from the level v'=0 of the term a
3 Σ
g+ is about 260 nm. A weak band at 146 nm observed in the CL spectrum was tentatively assigned to the band 0-6 of the Lyman series B
1Σ
u+ → X
1Σ
g+ of molecular hydrogen. However, it was not possible to use this band for the H
2 monitoring because of its low intensity. In the visible range we detected emission bands of other radiolysis products – CH radical and C atom. CH radicals were registered by the emission bands at 432 nm (the A
2∆ → X
2∏ transition) and 387 nm (the B
2Σ
- → X
2∏ transition). C atoms were recorded by the emission lines at 470 nm (the
1S →
3P transition) and 295 nm (the
5S
0 →
3P transition). The spectrum also contains a broad band at 184 nm unidentified at present and impurity bands of CO (the Cameron system
a 3Π→
X1 Σ
+).
The positions and half-widths of the observed bands remain practically unchanged over the entire concentration range studied. Main changes in the spectra relate to changes in relative intensity of bands. Upon doping of Ar matrices, a strong quenching of the emission band of self-trapped excitons (that is the M-band) was observed as shown in
Figure 4.
This demonstrates quite efficient energy transfer by free excitons to dopant and radiolysis products, since the diffusion length L
dif of thermalized free excitons in solid Ar (L
dif ~100 nm [
76,
77]) exceeds the distance between the dopants already at C=0.1%.
H atoms in the matrices were traced by a CL band centered at 166 nm, which belongs to the excimer Ar
2H
* [
70]. Formation of Ar
2H
* occurs with the participation of the matrix hole Ar
+ at the expense of the large affinity of Ar atom to a proton (proton affinity PA – 369.2 kJ mol
-1 [
64]). The quantum chemical calculations [
78] on the Rg
nH
+ species showed that the configuration with two solvent atoms of the D
∞h configuration with a proton between them is the energetically favorable configuration of the proton solvated in all Rare Gas matrices. The recombination of Ar
2H
+ cations with electrons results in appearance of the excimers Ar
2H
*. It is interesting to follow the change in the relative intensity of this excimer band with the dopant concentration. The intensity of the Ar
2H* emission band, hereinafter referred to as H, measured with respect to the Ar
2 band of self-trapped excitons, denoted M, is shown in
Figure 5.
Attention is drawn to the nonlinear behavior of the relative intensity of the H band which reflects concentration of H atoms. As can be seen from the channels (1) – (5) of CH
4 fragmentation, H atoms are formed in reactions (1), (3) and (4) with the branching ratios BR=0.5: 0.2 and 0.1 respectively for a photon energy of 9.8 eV and BR=0.25; 0.5 and 0.2 for a photon energy of 12.1 eV according to [
79], in which an analysis of energy-dependent branching ratios, the so-called breakdown curves, was presented. Let us analyze the reactions of secondary products of methane transformation CH
3 and CH
2 based on the data [
79]:
and
The emission of self-trapped excitons (9.8 eV) can produce H atoms from CH
3 via channels (6), (7) and (9) with BR=0.28, 0.5 and 0.25 correspondingly. An excitation of CH
2 by these photons induces its dissociation via channels (10) and (11) with BR=0.45 for both channels. In this case 2H atoms are produced in correlated way in channel (11). When excited by free Ar excitons (12.06 eV), the main contribution to the formation of the H atom is made by channels (6), (7) and (9) with BR=0.2, 0.8 and 0.3, respectively, moreover, reaction (7) produces 2H atoms simultaneously. Formation of H atoms from CH
2 occurs most effectively via channel (11) with BR=1. Channel (10) does not operate (BR=0) at this photon energy. It should be noted that H atoms which are formed in these reactions have an excess of kinetic energy, e.g. in channel (1) the experimental mean kinetic energy of the neutral lighter fragment H appeared to be 3.1 eV [
79]. In most other channels this energy exceeds 1 eV and such “fast” H atoms may diffuse for quite a long distance facilitating H + H recombination after thermalization. The mobility of H atoms in rare Gas solids and the stability of trapping sites have been the subject of numerous studies, incl. [
80,
81,
82,
83,
84,
85]. The case of an Ar matrix was considered in [
81,
84]. In Ar matrices the hydrogen atom can occupy two types of sites – interstitial
Oh and substitutional ones. According [
81] the EPR signals due to interstitially trapped hydrogen atoms in octahedral sites disappear near 16 K in solid Ar. The atoms trapped in the substitutional sites remain trapped [
84]. However, their excitation via energy transfer leads to the excimers Ar
2H
* formation and the resumption of the diffusion process. As a result, the processes of H atom production and their recombination compete each other upon irradiation.
Our study of the dose dependences of the emission of H atoms and CH radicals performed at low temperature showed a sharply nonmonotonic behavior of these emissions. As follows from the analysis of BR reactions (4), (7), (8), and (10), CH
3 radicals remain the main source of CH radical formation upon excitation by free excitons of the matrix; so the conclusion that CH is a CH
3 marker [
55] remains valid.
Figure 6 shows the dose dependencies of the H and CH emission measured simultaneously with the pressure in the experimental chamber.
The dose behavior of both optical emissions was quite nonmonotonic. We first observed an increase in both optical emissions and then an exit to saturation after 350 seconds of exposure to an electron beam, followed by a very gradual decrease. Weak bursts on the radiation yield curves were detected – two on the curve for the CH band and one on the curve for the H band. These peaks correlated with small peaks on the temperature curve. The main change of yields occurred after 3250 s of irradiation. This phenomenon is quite similar to the observed in pure solid methane [
29,
43] interpreted as the result of thermo-concentration self-oscillations [
44]. Details of the emission yields together with change of temperature and pressure are shown in
Figure 7a,b.
After prolonged irradiation at about 3250 s, the temperature and intensity of the H band, which is proportional to the number of H atoms in the matrix, begin to increase, then the intensity of the H band decreases, while the temperature and pressure sharply increase. It should be noted that since the temperature sensor (Si diode) was mounted on the back of a substrate cooled by LHe, it was not possible to determine the actual surface temperature of the Ar film. However, its significant increase during explosive emission of particles is indicated by a sharp drop in the TSEE yield, measured after explosive emission and redistribution of the TSEE yield with respect to the TSEE measured after irradiation of the sample with a subcritical dose. The beginning sublimation of the sample allows one to estimate only the lower boundary of the heating pulse — 30 K. In reality, the temperature rise in the pulse exceeds the triple point temperature of Ar – 87.78 K [
64]. The temperature pulse has three maxima and its duration was 100 s. After the first burst of pressure, which lasted about 20 s another maximum of pressure rise was observed but longer and less intense. During this period the change in the intensity of the H band followed the temperature course. Such an unexpected, at first glance, behavior of the H band can be understood if we consider more carefully the processes of energy transfer to the products of ongoing reactions. A high temperature of the sample during the burst is a “fingerprint” of the occurrence of exothermic reactions, in particular, the reaction of H atom recombination. Transfer of energy by free excitons of Ar matrix can populate low electronic states of H
2 molecule: the singlet state B
1 Σ
u+ and the triplet state a
3 Σ
g+ as can be seen from
Figure 8.
Transitions from the bound triplet state a
3 Σ
g+ to the repulsion curve of the lower state b
3 Σ
u+ result in the dissociation of the hydrogen molecule and the appearance of two “hot” H atoms, which, after thermalization, form centers responsible for the emission of the H band. Thus, the processes of recombination of H atoms with the formation of H
2 molecules proceed in competition with the photon-induced dissociation of H
2, which affects the dose dependence of the H band. An increase in the H band intensity during a burst of pressure was also observed at a methane concentration C=1% in the Ar matrix [
54]. However, this fact was misinterpreted as the absence of a contribution from the recombination of H atoms to the delayed desorption, since long before the pressure burst in this experiment, we observed a sharp drop in the intensity of the H band, which was attributed to the recombination of H atoms.
Note that the dose dependences of the H band intensity measured on lightly (C=1%) and highly (C=10%) doped matrices are different. While in the case of a highly doped Ar matrix, the burst was observed after the stage of a slow decrease in the intensity of the H band (see
Figure 6), in the case of a lightly doped matrix, the burst followed its growth as shown in
Figure 9. The decrease in the burst delay time in this experiment is associated with the use of higher current density.
At a methane concentration of C=5%, the dose dependence of the H band is similar to that recorded for 1%. It should be noted that the concentration of the radiolysis products is significantly lower than that of the dopant since many different primary dissociation channels operate simultaneously. Due to this, the conditions for excitation of H
2 by free excitons change with concentration: with increasing C, the efficiency of excitation of H
2 by excitons rises and therefore the contribution of the recombination reaction of H atoms to the H band increases. Accordingly, the flash effect accompanying the ejection of particles appeared to be much more pronounced at 10%, as can be seen from the comparison of
Figure 7a and
Figure 9.
It is interesting to compare the dose dependence of the H band with that detected for the CH band in highly doped matrix (see
Figure 6 and
Figure 7b). An intensity of the CH band, which is a marker of the CH
3 radical, follows the temperature course with an exception of one local event at about 3340 s when the temperature and intensity of the CH band change in the opposite way. Such a behavior became clear if consider the balance and interplay between C
2H
6 and CH
3. As follows from Refs. [
63,
88,
89] the photodissociation (absorption) cross section of ethane is close to that of methane and both free and self-trapped excitons fall into the range of ethane absorption. There are a number of primary dissociation channels which include reactions: C
2H
6 +
hν → CH
3 + CH
3 and C
2H
6 +
hν → CH
3* + CH
3 [
88]. Note that the enthalpy Δ
rH of the second reaction with the excited CH
3* radical formation Δ
rH = 9.54 eV [
88], is close to the photon energy of self-trapped excitons in Ar (9.8 eV). The enthalpy of the first reaction is essentially lower – Δ
rH = 3.81 eV [
88], and both of them can proceed. Unfortunately, at present, to the best of our knowledge, there is no quantitative information on branching ratios of these reactions, however, as it was shown [
88] the reaction with C-C bond fission yielding two CH
3 radicals refers to active channels of the C
2H
6 photodissociation. Based on this conclusion, we can expect similar behavior of the CH and H bands during the burst. The results obtained (see Fig. 7 a and b) suggest that the luminescence flashes detected in these bands can be assign as indicators of the recombination of H atoms and CH
3 radicals with the release of thermal energy, which ensures the explosive ejection of particles.