Agniva
Banerjee
,
Julien
Mandon
,
Frans
Harren
and
David H.
Parker
*
Radboud University Nijmegen, Institute for Molecules and Materials, Heyendaalseweg 135, Nijmegen, 6525 AJ, The Netherlands. E-mail: parker@science.ru.nl; Tel: +31 243653423
First published on 17th December 2018
Collision-induced absorption between O2 and CO2 molecules associated with the a1Δg (v = 1) ← X3Σ−g (v = 0) band of oxygen around 1060 nm was measured using cavity ring-down spectroscopy. The lineshape for this transition is measured for the first time, and the integrated cross-section is found to be smaller than the only previous report. For pure oxygen, we find an integrated absorption value of (2.10 ± 0.31) × 10−4 cm−2 amg−2 which is in good agreement with the previous reported values. For O2–CO2 collisions we report an integrated value of (6.37 ± 1.09) × 10−5 cm−2 amg−2 which is small but still significant and not accounted for by theory.
Carbon dioxide, CO2, is the most abundant greenhouse gas on Earth2 and industrial progress has resulted in a significant rise in CO2 concentration3 which causes increased surface warming, affecting Earth's climate. Reports of local yearly average CO2 concentrations at Mauna Loa Observatory, Hawaii, for example, indicated values of 407.98 ppm in January 2018, 406.13 ppm in January 2017 and 402.52 ppm in January 2016.4 A better understanding of the sources and sinks of carbon cycles is thus important. High precision satellite measurements of the local concentration of CO2 by ESA's Environmental Satellite (ENVISAT), and NASA's Orbiting Carbon Observatory (OCO-2) provide important new data for this quest.5–7 However, clouds and aerosols complicate accurate determination of CO2 concentration,8 and for these satellites, calibration methods are also needed.9,10
Molecular oxygen, O2, is often used for calibration as it is abundant and has a well-known spectrum.11–13 O2 is unusual in having many bound electronic states below the first dissociation limit, but all of the transitions from the ground electronic state in the infrared, visible, and ultraviolet region below the allowed VUV Schumann–Runge transition are optically forbidden. Very weak magnetic dipole and/or electric quadrupole transitions from the ground state to the lower excited electronic states are possible, along with transitions induced by weak spin–orbit coupling, which can increase significantly during a collision. Under Earth's atmospheric conditions, light absorption is often governed by collision-induced absorption. When one oxygen molecule collides with other atmospheric species such as O2, N2, or CO2, a transition moment can be induced, therefore lifting up the symmetry restrictions slightly,14 the efficiency of this phenomenon depends on the quantum character of the collision partner. The probability of CIA is also strongly dependent on the vibrational level of the O2 excited electronic state reached by CIA, where vertical excitation (Δv = 0) should be favoured for the lower vibrational states of oxygen.
In this paper, we study the O2 transition from X3Σ−g (v = 0) to a1Δg (v = 1), which lies around 1060 nm. Mixing by spin–orbit coupling of the X3Σ−g ground state with the higher lying d1Πg and/or C3Πg Rydberg state gives this transition magnetic dipole allowed character.15 While a strong decrease in the magnetic dipole transition for vibrational levels a (v > 0) due to unfavourable Frank–Condon factors is predicted,16 the collision induced absorption coefficient was found to be surprisingly large.17,18 Literature values for integrated O2–O2 CIA cross sections in cm−2 amg−2 for the near infrared bands of oxygen are listed in Table 1.19–22
O2 | O2–O2 | O2–CO2 | ||
---|---|---|---|---|
X3Σ−g (v = 0) → | a1Δg (v = 0) | a1Δg (v = 1) | a1Δg (v = 0) | a1Δg (v = 1) |
a For a given pressure P (atm) and temperature T (K), the density D in amagat (amg) is expressed as ![]() |
||||
Chagas et al.17 | 3.34 × 10−4 | 2.27 × 10−4 | ||
Tran et al.19 | 2.74 × 10−4 | |||
Tabisz et al.18 | 3.00 × 10−4 | 2.30 × 10−4 | ||
Hermans et al.20 | 2.97 × 10−4 | 2.14 × 10−4 | ||
Fraser et al.21 | 2.13 × 10−4 | |||
Vangvichith et al.22 | 7.80 × 10−5 | 1.85 × 10−4 |
In a recent paper, Karman et al.23 reported ab initio calculations and experimental measurements of CIA between O2–O2 and O2–N2 for the X3Σ−g (v = 0) → a1Δg (v′ = 0, 1, 2) and b1Σ+g (v = 0, 1) transitions in details. Two mechanisms were found to be important: spin–orbit coupling which enhances only the Δv = 0 transition, and for O2–O2 (two paramagnetic molecules) an exchange mechanism involving the total electron spin of the collision complex in which total spin remains conserved while one O2 monomer changes its spin state.24 The spin exchange mechanism was found to also enhance Δv > 0 transitions. Near quantitative agreement was found between ab initio theory and experiment for these two collision pairs.
From fundamental interest, but also due to the importance of CO2 in exoplanet atmospheres,25,26 it is useful to compare CIA of O2–O2 and O2–N2 with that of O2–CO2. Vangvichith et al.22 reported an integrated CIA strength for O2–CO2 for the X (v = 0) → a (v = 0) and a (v = 1) transition. Fraser21 obtained a value of 2.13 × 10−4 cm−2 amg−2, significantly higher than the value from Vangvichith for the (v = 0) transition. We also note that Hidemori and co-workers27 reported a 2–3 times increase in emission intensity with the addition of a heavy atom/molecule for the a (v = 0) to X (v = 0) O2–O2 collision-induced emission transition, which they suggested as being due to an increase in the spin–orbit mixing of the O2 molecule. The only measurement of the CIA for the O2–CO2 X3Σ−g (v = 0) to a1Δg (v = 1) transition, ref. 22, reported only the integrated cross section, i.e., without line-shape information. The limited availability of data, together with the discrepancy in cross section values motivated the present study of the X (v = 0) → a (v = 1) transition for O2–CO2.
![]() | (1) |
The constant value C was obtained from the empty cavity ring-down time at the measured laser wavelength. The Rayleigh scattering values for O2 were taken from Bates32 and for CO2 were taken from Bideau-Mehu33 (which produce values within the error bar) and subtracted from our measured extinction values. The resultant absorption values were fitted with a quadratic equation in order to extract the CIA coefficients at different wavenumbers. The fitting procedure with resultant absorption values for different gas mixtures is shown below in Fig. 2.
This whole procedure was repeated for both pure O2 and O2/CO2 (1:
1 mixture) for all wavelengths. The CIA coefficient values obtained for pure O2 were lowered with a factor of (0.502) or 0.25 because in the mixture only 50% O2 is present and CIA is quadratic in nature. The resultant scaled coefficient values originating from O2–O2 collisions were subtracted from the CIA coefficient values obtained from a (1
:
1) O2–CO2 mixture, which contains all contributions according to their relative partial pressures, so that this correction results in a CIA coefficient originating only due to O2–CO2 collisions.
In order to extract the CIA contribution arising from pure CO2–CO2 collisions (if this is measurable), the same procedure was carried out for pure CO2 gas. Line absorption of CO2 was taken from the HITRAN database,34 corrected for the corresponding pressure values and subtracted from the measured extinction values. The remaining CO2–CO2 absorption signal was very low as was the magnitude of its CIA coefficient when fit with a quadratic equation. Our system may not be sensitive enough to quantify such very low CIA values, as discussed later in this article. In the following analysis, we neglected the contribution arising from CO2–CO2 collisions.
The HITRAN database allows us to identify the most prominent H2O line absorption between 9150–9350 cm−1. To avoid these absorptions, experimental wavelengths were chosen which fall between water line absorptions. Gases with very high purity and heating of the tubes ensured a water absorption value <10−8 cm−1. From Fig. 3, it is clear that line absorptions of CO2 are much more pronounced than those of O2. Although the strongest absorptions for CO2 occur around 9500 cm−1, our diode could measure only until 9480 cm−1. The CO2 lines were so dense, it was very difficult to set the laser at a wavelength without line absorption. In addition, when the molecules collide, due to inelastic collisions, there will be a transfer of intensity. Our line mixing program41 takes into account molecular inelastic collisions towards the line shape which results in the intensity transfer between two adjacent rotational lines following an expression by Fano for this effect.35 Later Tonkov36 introduced one more parameter previously defined by Fano, to allow line mixing between transitions between two different branches, called the ABC model. This model was applied for O2, CO, CO2, CH3F, and CH4 successfully.36–40 Our line mixing model uses a Voigt convolution of line mixing spectra and Gaussian contributions arising from Doppler effect as described by Spiering and co-workers.41 The computational program determines the line mixing correction on the sum of Voigt line shape spectrums based on the collisional parameters, which includes pressure shifts, for each transition value reported in the HITRAN database.34 Due to molecular collisions, the Lorentzian lines need to be modified with a Gaussian function, which resulted in a Voigt shape profile. The corrected line absorption values for both O2 and CO2 were subtracted from the measured extinction. Resulting absorption values when fitted with the square of density, yield the final CIA coefficient values, which are presented in Fig. 4.
![]() | ||
Fig. 3 Line absorptions from the HITRAN database.34 The region of interested in this study is marked in the red box. Note the different scales for CO2versus O2 absorptivity. |
![]() | ||
Fig. 4 Collision-induced absorption measured spectrum; black squares: CIA total (using 50–50% O2–CO2 mixture); red circles: O2–O2 CIA (considering 50% contribution). |
It is clear from Fig. 4 that the total collision-induced absorption coefficient values from a 1:
1 O2–CO2 mixture is higher than the O2–O2 component (considering its 50% contribution), with a maximum value around 9380 cm−1. Since pure CO2 showed a very weak and noisy absorption profile, another approach was used to evaluate any CO2–CO2 contribution which could underlie the O2–CO2 collisional effect. Different O2–CO2 mixture ratios, as shown in Fig. 2, were used to re-calculate collision coefficients of the different possible interactions included in eqn (1). The wavelength of 9380.685 cm-1 (near the O2–O2 CIA peak) was studied, considering the fact that both O2–CO2 and CO2–CO2 (if contributing at all) components contribute the most at this photon energy. For the five different mixtures used in this study, eqn (1) was transformed into five different expressions and by solving these equations, the normalized values obtained for O2–O2, O2–CO2 and CO2–CO2 were (9.65 ± 0.18) × 10−7 cm−1 amg−2, (3.83 ± 0.21) × 10−7 cm−1 amg−2, (6.50 ± 0.95) × 10−9 cm−1 amg−2 respectively. CIA values for O2–O2 and O2–CO2 obtained from (1
:
1) O2
:
CO2 mixture around 9380 cm−1, plotted in Fig. 5, match quite well with the one obtained in this process performed with different mixtures of O2–CO2 at 9380.685 cm−1 and these lie within the error limit. Owing to its low value, we neglected pure CO2–CO2 collision effects for the rest of our calculation. We assume similar behaviour from CO2 for other wavelengths for this band as well. Fig. 5 shows both the normalized O2–O2 and O2–CO2 CIA spectrum. The integrated CIA cross section value obtained from an extrapolated Lorentzian fit with FWHM of (97 ± 2) cm−1 was (6.37 ± 1.09) × 10−5 cm−2 amg−2. The error in the integrated absorption strength originates from the choice of fitting method, and from statistical fluctuations of the measured data points and also from the extrapolation. The peak was observed at (9384 ± 2) cm−1. Uncertainties in temperature and pressure are smaller than statistical error.
No sharp features were observed in the CIA spectrum. The O2–CO2 CIA spectrum is much narrower than the O2–O2 CIA spectrum, as is clearly visible in Fig. 5, and its integrated intensity is about 3.3 times lower. The efficiency, which is a measure of the probability for a transition to occur based purely on Frank–Condon (FC) factors, can be expressed with the following expression
![]() | (2) |
O2–O2 | O2–CO2 | |
---|---|---|
η | η | |
v a = 0 | 1 | 1 |
v a = 1 | 0.76 ± 0.10 | 0.85 ± 0.50 |
The large error associated with ηva = 1 for O2–CO2 is due to uncertainty (46%) associated with the determination of the integrated CIA values. Our cross section for the X (v = 0) → a (v = 1) transition is a factor of almost three smaller than the only previous literature value for this process.22 A possible reason for this discrepancy could arise from the extrapolation method and the accuracy of our line-mixing model, which excludes all the line-absorptions arising due to H2O, O2 and CO2. Among these, the CO2 monomer makes the most significant contribution to the light absorption process at high pressure.
The relatively large O2–CO2 CIA cross section we find for the X (v = 0) → a (v = 1) compared to X (v = 0) → a (v = 0) is not consistent with a recent model of CIA by Karman et al.23 which predicts a very low cross section since CO2 (as well as N2 as discussed in ref. 23) is diamagnetic and cannot follow the favoured total spin exchange mechanism. This disagreement is not surprising considering that the isotropic potentials used in ref. 23 will be less applicable to CO2. Further analysis is necessary to better understand this effect.
The efficiency values for the optical transition are comparable between O2–CO2 & O2–O2 for the above-mentioned transition, considering the error bar. The efficiencies based on the FC factor for this transition is 0.013.17 A change in the internuclear distance in the excited a-state from 1.2156 Å to 1.1152 Å is necessary to modify the FC factors in a way that it agrees with the integrated intensity values, within the error bar, for O2–O2 collisions.30 Comparing the measured efficiency values we can suggest that CO2 deforms the excited a-state to some extent. Despite having a ground state configuration of 1Σ+g, a large electric quadrupole moment and polarizability anisotropy of CO2 could cause more penetration during a collision with the O2 molecule, producing a compression in the excited state bond length of O2.22 Being a larger molecule, CO2 could further increase spin–orbit coupling between the X3Σ−g and higher electronic d1Πg state compared to O2.27 We note that this is an oversimplification of the complex processes occur during collisions.
Chagas et al.17 have reported the effect of CIA on the 1060 nm band in the total absorption in our atmosphere. For a mid-latitude summer atmospheric profile with a solar zenith angle of 60%, a surface albedo of 0.1, and a total incoming solar irradiation of 1368 W m−2, they found that the CIA of O2–O2 at 1060 nm band contributed for 0.30 W m−2. This value was recalculated for our CIA contributions of O2–O2 and O2–CO2 collisions and we see that the CIA at 1060 nm band adds 0.22% to the short wavelength absorption in our atmosphere. About 0.000118% of this is caused by collisions between CO2 and O2 (Table 3).
CIA contributors | O2–O2 | O2–CO2 | Total |
---|---|---|---|
CIA contribution (10−5 cm−2 amg−2) | 1.01 ± 0.04 | (5.30 ± 0.90) × 10−4 | ≈1.01 ± 0.04 |
Sunlight absorption (W m−2) | 0.30 ± 0.03 | (1.58 ± 0.29) × 10−4 | ≈0.30 ± 0.03 |
Short-wave absorption (%) | 0.22 ± 0.05 | (1.18 ± 0.14) × 10−4 | ≈0.22 ± 0.05 |
At sunrise or sunset when longest path lengths are present the absorption process should become most significant. Due to the very low absolute concentration of CO2 gas present in our atmosphere, the O2–CO2 CIA coefficients are very small and it is probably not possible to measure this induced absorption process in Earth's atmosphere.
This journal is © the Owner Societies 2019 |