H.
Fleurbaey
,
S.
Kassi
and
A.
Campargue
*
Univ. Grenoble Alpes, CNRS, LIPhy, Grenoble, France. E-mail: Alain.Campargue@univ-grenoble-alpes.fr; Fax: +33 476635495; Tel: +33 476514319
First published on 24th July 2024
The hydrogen dimer, (H2)2, is among the most weakly bound van der Waals complexes and a prototype species for first principles ab initio studies. The detection of the (H2)2 infrared absorption spectrum was reported more than sixty years ago at a temperature of 20 K. Due to the sharp decrease of the (H2)2 abundance with temperature, detection at room temperature was generally considered hardly achievable. Here we report the first room temperature detection of partly resolved rotational structures of (H2)2 by cavity ring down spectroscopy at sub-atmospheric pressures, in the region of the first overtone band of H2 near 1.2 μm. The quantitative analysis of the absorption features observed around ten allowed or forbidden transition frequencies of the monomer provides insight on the structure of this elusive species and a benchmark for future theoretical studies.
With a bond dissociation energy of about 3 cm−1 (∼4.5 K), the hydrogen dimer (H2)2 is indeed among the most weakly bound van der Waals complexes, and is thus expected to have an extremely low concentration at room temperature. In order to theoretically predict its absorption spectrum, the multi-dimensional potential energy surface (PES) of the (H2)2 dimer must be calculated with an accuracy significantly better than 1 cm−1, which represents a real challenge for first principles quantum chemical calculations. Started in the 1940s and 1950s,2–4 the calculations of the interaction potential between two H2 molecules remain nowadays an active field.5–13 Beyond its academic interest, the H2 dimer is of special interest in various domains. As hydrogen is more and more considered as a versatile energy carrier with important needs for storage and transportation at standard temperature conditions, the equilibrium thermodynamics of the H2 dimers in high pressure tanks of hydrogen is particularly relevant in the present period. In planetary science, the (H2)2 spectral signature has been evidenced in the far-infrared (FIR) spectra of Jupiter and Saturn recorded in the Voyager IRIS mission,14–16 and more recently in the FIR spectra of Uranus and Neptune obtained by the Spitzer space telescope.11 The assignment to (H2)2 of the faint narrow structure observed in planetary spectra near the S(0) and S(1) pure rotational lines was confirmed by McKellar from laboratory absorption spectra of H2 at 77 K.15,17,18
Although Coulomb explosion imaging has very recently been used to probe the geometry of the dimer (there is no preferred orientation) and determine the intermolecular separation,19 absorption spectroscopy appears to be a complementary method to provide observations and validation tests for theoretical calculations. Due to the weakness of the binding energy, all previous spectroscopic detections of the H2 dimer were reported at low temperature. The H2 dimer contribution to the absorption spectrum of hydrogen at 20 K was first demonstrated in 1964 by Watanabe and Welsh in the region of the (1–0) fundamental band of H2 centered near 4150 cm−1.20 Later, McKellar and Welsh recorded the same band at a similarly low temperature with a grating spectrograph and an improved spectral resolution, and partly assigned the few weak sharp dimeric features surrounding the electric quadrupole H2 lines.21 A further improvement of the quality of the (H2)2 infrared spectrum in terms of spectral resolution and sensitivity was achieved in ref. 22 by Fourier transform spectroscopy (FTS) with a spectral resolution in the range of 0.04–0.16 cm−1. The FTS recordings were performed at 20 K and pressure values of a few tens Torr and extended up to the first overtone region, (2–0), for H2, D2 and HD, and their mixtures. In summary, the review of the literature indicates that the previous spectra of the H2 dimer were recorded at low temperature and are limited to rotational transitions at 77 K in the FIR and spectra at 20 K in the fundamental band. To the best of our knowledge, no experimental spectra of the H2 dimer were reported since 1991.
Revisiting the (H2)2 spectrum with modern highly sensitive laser-based techniques appears to be timely. The sensitivity nowadays achieved by cavity-enhanced techniques allows for the quantitative detection of trace gases with extremely low concentration (see e.g. ref. 23). In the present work, we adopted the cavity ring down spectroscopy (CRDS) technique and investigated the high sensitivity room temperature absorption spectrum of hydrogen in the 8000–8500 cm−1 region of the (2–0) first overtone band. In his study of the (1–0) band at 20 K (and 25 Torr),22 McKellar observed a few sharp lines assigned to transitions between bound states of the dimer around the H2 quadrupole lines. Our aim is to investigate to which extent similar sharp lines are detectable and preserved at room temperature in the (2–0) band, taking into account the expected sharp decrease of the dimer concentration and the possible change of the appearance of the dimeric rotational structure.
![]() | ||
Fig. 1 The quadrupole lines and CIA in the 8000–8500 cm−1 region of the first overtone band of H2. The binary coefficients (left-hand scale) as recommended by the HITRAN database (black line) or measured by CRDS26 (yellow diamonds) are superimposed to the stick spectrum of the (2–0) electric quadrupolar transitions (right-hand scale). |
Position (cm−1)24 | Intensity (cm per molecule)24 | Pressure (Torr) | |||||||
---|---|---|---|---|---|---|---|---|---|
50 | 100 | 200 | 350 | 500 | 750 | 1000 | |||
Q 2(3) | 8017.183 | 2.36 × 10−28 | × | ||||||
Q 2(2) | 8051.988 | 3.29 × 10−28 | × | × | × | × | × | ||
Q 2(1) | 8075.307 | 2.60 × 10−27 | × | × | × | × | × | × | × |
Q 2(0) | 8087.003 | — | × | × | × | × | × | × | |
Q 1(1) + Q1(3) | 8281.126 | — | × | × | × | ||||
Q 1(1) + Q1(2) | 8298.719 | — | × | × | × | × | |||
Q 1(1) + Q1(1) | 8310.507 | — | × | × | × | × | |||
Q 1(1) + Q1(0) | 8316.420 | — | × | × | × | × | |||
S 2(0) | 8406.361 | 1.56 × 10−27 | × | × | |||||
Q 2(1) + S0(0) | 8429.681 | — | × | × | |||||
Q 2(0) + S0(0) | 8441.376 | — | × |
The present experimental study follows two previous contributions by cavity ring down spectroscopy (CRDS) dedicated to the H2 absorption spectrum in the region. In ref. 25 the frequency of the electric quadrupole transition lines was determined with sub-MHz accuracy. The zero-pressure transition frequencies were obtained by applying a multi-spectrum fit procedure with various beyond-Voigt profile models to a series of spectra recorded at different pressures over a spectral interval limited to less than 0.5 cm−1 around the line center. As illustrated below, this spectral amplitude is insufficient to reveal the dimer spectral signature which extends over several cm−1 around the center of the H2 electric quadrupole lines. A second contribution was devoted to CIA measurements performed during pressure ramps up to 1 atm, at 28 spectral points selected far from absorption lines.26 A purely quadratic pressure dependence was obtained for the CRDS absorption coefficient at each measurement point and the CIA binary coefficients derived with a 1.5% accuracy were found significantly higher than the HITRAN values (see Fig. 1) (note that HITRAN CIA values are calculated values from ref. 27).
The CRDS setup used for the present recordings is similar to that used in ref. 25. Briefly, an external cavity diode laser (ECDL) (Toptica fiber-connected DL pro, 1200 nm) tunable over the region of interest (8000–8650 cm−1) is used as light source and injected in a 1.4-meter-long high-finesse cavity. For the evacuated cavity, the ring-down time τ0 varied from about 160 μs to 340 μs depending on the wavenumber. For each frequency step, the average central emission frequency of the ECDL was actively stabilized using a software based proportional–integral loop, with 100 Hz bandpass, acting on the laser current. For each frequency point, about 40 ring-down events were averaged leading to a minimum detectable absorption coefficient between 2 × 10−12 and 8 × 10−12 cm−1 for a single scan. Following the frequency comb referenced CRDS method (CR-CRDS), each ring-down event was associated to its laser frequency. A self-referenced frequency comb (Model FC 1500–250 WG from Menlo Systems) was used for the frequency calibration of the spectra. The laser frequency is measured “on the fly” by recording the beat note between a fraction of the ECDL light and the nearest tooth of the frequency comb, the tooth number being determined using a Fizeau type wavemeter (HighFinesse WSU7-IR, 5 MHz resolution, 20 MHz accuracy over 10 hours) (see ref. 28 and 29 for details).
The temperature measured with an analog temperature sensor (TSic 501, IST-AG, ±0.1 K accuracy) fixed on the cell surface varied between 292.6 and 294.5 K according to the recordings.
The CR-CRDS absorption spectra of natural dihydrogen (Alphagaz2, 99.9999% chemical purity) were recorded in flow regime in order to reduce spectral interferences with very strong lines of water vapor desorbing from the CRDS cell or from the injection tubes. The hydrogen pressure in the CRDS cell was actively regulated to set values through a needle valve connecting the cell to a turbo pump group and a voltage-controlled valve connecting the cell to the hydrogen cylinder, using a computer based proportional/integral controller. The gas pressure in the cell was continuously monitored by a capacitance gauge (MKS Baratron, 1000 mbar full range).
The spectra were recorded for series of pressure values ranging between 10 and 1000 Torr over about 10 cm−1 around the Q2(1), Q2(2), Q2(3), and S2(0) electric quadrupole lines - see Table 1. Additional recordings were performed around the calculated frequencies of the Q2(0) transition, which is forbidden for symmetry reasons, and several collision-induced transitions involving a double excitation of two H2 molecules either both in the (1–0) fundamental band: Q1(1) + Q1(0–3) or the pure rotational transition S0(0) combined with overtone transitions: S0(0) + Q2(0–1). The corresponding frequency values obtained from ref. 24 are listed in Table 1.
A detailed view of the spectrum around the Q2(1) line before and after removal of the water absorption is presented in Fig. 3 (upper panel). The absorption features surrounding the H2 line are unambiguously assigned to H2 dimer by comparison to the dimer features observed by McKellar in his FTS spectra in the region of the Q1(1) line of the (1–0) fundamental band.22 At first sight, the high similarity between the dimeric structures observed at 20 K and 300 K (and in distinct vibrational bands) is surprising. It will be discussed below.
![]() | ||
Fig. 3 Comparison of the (H2)2 absorption spectrum near the Q2(1) line of the first overtone of H2 and near the Q1(1) line of the (1–0) fundamental band (lower panel). The frequency axes of the two spectra are centered on the H2 transition frequencies and the width of the displayed spectral interval is identical. Upper panel: room temperature CR-CRDS spectrum near the Q2(1) line recorded at 200 Torr. The dimer spectrum (black) is obtained after subtraction of the water vapor contribution (light red), Lower panel: FTS spectrum at 20 K recorded by McKellar in the region of the Q1(1) line. The H2 pressure was 25 Torr and the absorption pathlength was 112 m.22 |
The overview of the different dimeric spectral features observed in the entire region is presented in Fig. 4. Spectra recorded at 350 and 750 Torr are superimposed. Overall, ten absorption structures are assigned to the dimer. Four of them are surrounding an electric quadrupole line of the (2–0) band: Q2(1–3) and S2(0). The remaining six are located near the sum of the frequencies of two H2 transitions either both in the (1–0) fundamental band, Q1(1) + Q1(0–3), or the S0(0) pure rotational excitation combined with an overtone transition, S0(0) + Q2(0), and S0(0) + Q2(1). The calculated frequencies of the H2 allowed and forbidden transitions (given in Table 1) correspond to the blue vertical lines displayed in Fig. 4. It is worth noting that both the relative spectral scale and the intensity scale are the same on the different panels. A roughly estimated baseline due to the CIA was subtracted from each spectrum. Around the four H2 quadrupole lines, the intensity of the dimeric structures appears to roughly scale with the intensity of the line.
For comparison purposes, we give the orders of magnitude of the various absorptions involved. For instance, in the Q2(1) region at 1 atm, the highest absorption coefficient of the dimer structure is on the order of 10−9 cm−1 to be compared with a peak absorption coefficient of about 10−6 cm−1 for the Q2(1) line (see Fig. 2 of ref. 30) and a value of about 3.5 × 10−8 cm−1 for the CIA absorption level (see Fig. 1). The noise equivalent absorption of the recorded CR-CRDS spectra varies from 2 to 8 × 10−12 cm−1, depending on the spectral region.
The pressure dependence of the (H2)2 absorption structure is displayed in Fig. 5 for the Q2(1) region for which the largest series of recordings were performed (see Table 1). A clear broadening of the sharp dimer lines is observed (see also Fig. 4).
Only a partial assignment is available for the different absorption features (see Fig. 5). Following ref. 7 and 22, the different sharp structures observed around Q2(1), Q1(1) + Q1(1), and Q1(1) + Q1(2) are assigned to transitions involving the l = 0 and l = 1 bound levels, thus 1 → 0 and 0 → 1 on the low and high energy side of the monomer transition, respectively. The observation of resolved doublets for both 1 → 0 and 0 → 1, near the Q2(1) and Q1(1) + Q1(1) lines, has been interpreted by McKellar as due to the splitting of the total angular momentum, j, which is the vector sum of the l end-over-end rotational angular momentum of the dimer with j rotational momenta of the two monomers. This splitting reflects the angular anisotropy of the intermolecular potential.22 (Note that, as in the fundamental region, a few additional sharp lines significantly weaker than the main doublets are observed – see Fig. 3 and 6). Such narrow features are not observed around the S2(0) line; in the fundamental band, they are predicted to be much weaker near S lines due to the different type of dipole inducing mechanisms.21,32 Transitions involving higher l quasi-bound states give rise to much broader resonances.
Although the dimer absorption structures appear as a complex superposition of features with distinct shapes and amplitudes, it is worth trying to model them empirically to account for the observations and provide quantitative information for future analysis and comparison to theoretical results. The sharp l = 0 ↔ l = 1 lines could be modeled as a sum of absorption lines with a Voigt line-shape, the Gaussian component being fixed the Doppler broadening of (H2)2 at room temperature. An example of such a reproduction of the sharp lines surrounding the Q1(1) + Q1(1) simultaneous transition is displayed in Fig. 6. The retrieved parameters of the observed sharp absorption features are listed in Table 2. They include the absolute line position and the position value relative to the monomer transition frequency, the Lorentzian width and the normalized integrated absorption coefficient obtained from the ratio of the line area (A) by the pressure squared.
Monomer transition | Position (cm−1) | Relative positiona (cm−1) | Broadening coefficient (cm−1 atm−1) |
A/P
2![]() |
---|---|---|---|---|
Notes.a Difference between the measured position and the nearest transition frequency of the monomer (see Table 1).b Normalized integrated absorption coefficient determined from the quadratic pressure fit of the integrated absorption coefficients in Fig. 7 and equal to ratio of the line area by the pressure squared.c This dimer line is partly obscured by a strong water line which hindered the retrieval of reliable line parameters. | ||||
Q 2(1) | 8072.966 | −2.342 | 0.19 | 2.84 |
8073.266 | −2.042 | 0.19 | 1.14 | |
8073.716 | −1.592 | 0.19 | 0.36 | |
8075.649 | 0.341 | 0.21 | 0.28 | |
8076.007 | 0.699 | 0.21 | 1.95 | |
8076.229 | 0.921 | 0.21 | 4.48 | |
8076.429 | 1.121 | 0.21 | 0.56 | |
Q 1(1) + Q1(1) | 8308.100 | −2.407 | 0.19 | 4.88 |
8308.220 | −2.287 | 0.19 | 0.61 | |
8308.401 | −2.107 | 0.19 | 2.44 | |
8308.950 | −1.557 | 0.19 | 0.49 | |
8310.501 | −0.006 | 0.19 | 0.25 | |
8311.129 | 0.622 | 0.19 | 2.76 | |
8311.351 | 0.844 | 0.19 | 6.07 | |
8311.601 | 1.094 | 0.19 | 0.25 | |
Q 1(1) + Q1(2) | 8296.25 | −2.47 | ||
8299.792 | 1.073 | 0.20 | 1.37 |
The broader, l = 1 ↔ l = 2 structures, however, have distinctively asymmetric line shapes resembling Fano lines. This is particularly striking in the structure surrounding the S2(0) line (Fig. 4). A more detailed analysis of these features is beyond the scope of this paper.
The pressure dependence of the integrated absorption coefficient and of the Lorentzian width is shown for five of the sharp lines in Fig. 7. Let us note that the observed Lorentzian broadening of all the sharp lines (left panel) is similar and proportional to pressure with no significant broadening at zero-pressure. This confirms the absence of predissociation effects for the l = 0 and 1 states involved in the considered transitions. Note that the obtained pressure broadening coefficient of (H2)2 by H2 [about γD = 0.20(1) cm−1 atm−1(HWHM)] is significantly larger than the self-broadening coefficient of the H2 monomer (about 0.05 cm−1 atm−1 (ref. 1)). It is interesting to compare γD to γc = 1/(2πτc) where τc is the average time between two collisions. If each collision of the dimer with an H2 molecule leads to the dissociation of the dimer, γc is expected to be close to γD. The average time between two collisions is the ratio of the mean free path to the average molecular velocity: where d = r1 + r2 is the effective kinetic diameter (r1 and r2 are the kinetic radii of (H2)2 and H2, respectively) and
where μ is the reduced mass of (H2)2 and H2 and n is the volume density of H2 molecules. The kinetic radius of H2 is 1.445 Å. If we assume that the kinetic radius of the dimer coincides with the Lennard-Jones radius of about 3.4 Å,19 we obtain a broadening coefficient of 0.21 cm−1 atm−1 close to the measured value and thus confirming that each collision dissociates the fragile dimer structure.
A clear negative position shift is observed for the apparent center of the dimer structure compared to the position of the monomer line. The average position value of the four sharp lines near Q2(1) and Q1(1) + Q1(1) is shifted by −0.6906 and −0.7624 cm−1 compared to the corresponding monomer line center, respectively. This value can be compared to a value of −0.44 cm−1 obtained for the Q1(1) transition of the fundamental band22 - see Fig. 3. The negative value of the shift indicates that the dimer potential well is deeper when the monomer is vibrationally excited.
As concerns the separation of the average value of the 1 → 0 and 0 → 1 doublets, we obtain a value of 3.00 cm−1, which would roughly correspond to 4B in a rigid molecule (B is the dimer rotational constant). From this B value, we derive an average distance between the two H2 molecules of about 4.5 Å which is consistent with the distribution of H2–H2 intermolecular distances recently measured by Coulomb explosion imaging.19 Interestingly, from the (H2)2 spectra in the fundamental band at 20 K which led Watanabe & Welsh20 to identify bound states for the first time, sixty years ago, a 4.2–4.6 Å range of intermolecular distances was already proposed.
The measured value of the integrated absorption coefficient is the product of the dimer concentration (in cm−3) by the line intensity (in cm per molecule). As the latter is unknown, our quantitative measurements do not allow for an evaluation of the dimer concentration. Recently, starting from the equation of state (EOS), Halpern evaluated the equilibrium constant of the H2 dimerization, 2H2 ↔ (H2)2, in the 25–45 K temperature range.13 Statistical thermodynamics (ST) calculations were found consistent with the EOS results. Under our request, Halpern extended his ST calculations up to 300 K and obtained a value of the KD constant of 3.7 × 10−6 atm−1 indicating that, at 1 atm, the relative concentration of the H2 dimer is 3.7 ppm. Note that the room temperature KD constant is about four orders of magnitude smaller than its 20 K value. As an energy vector, hydrogen gas is currently pressurized in tanks at room temperature and pressures of several hundred atmospheres. Based on the above KD constant, the dimer relative concentration in these tanks reaches values higher than 0.1%.
Regarding spectral congestion, with its very sparse absorption spectrum, the hydrogen molecule is a more favorable case than water vapor for the detection of its dimer but the very small binding energy of about 3 cm−1 (to be compared to a value of 1105 cm−1 (ref. 36) in the case of the water dimer) leads to much lower relative abundances in the case of (H2)2 (the room temperature equilibrium constant is about four orders of magnitude larger in the case of water dimer37).
The comparison of the contribution of the (H2)2 and (H2O)2 spectrum to the hydrogen and water vapor absorption, respectively, is interesting. The water absorption is modeled as a sum of two contributions: the rovibrational lines and the absorption (self)-continuum. The semi-empirical MT_CKD model38,39 which is the standard model of the water continuum implemented in climate and weather prediction models is a far-wing line shape model of the monomer. The basis of the MT_CKD model is thus not related to the water dimer even though it is now established that an important part of the water continuum is due to the water dimer. Conversely, for the (non-polar) hydrogen molecule, the absorption also includes both absorption lines and a continuum (see Fig. 1) but the monomer lines are very weak electric-quadrupole lines, and the continuum is in this case mainly due to transient electric dipole moments induced by collisions between H2 molecules (commonly known as CIA). The dimeric absorption is inseparably linked to the CIA since both arising from the same induced dipole mechanism. The hydrogen CIA values recommended by the HITRAN database are calculated values obtained by neglecting the dimer contribution.27 At room temperature, broad absorption features of (H2)2 may extend over several tens cm−1 around the H2 transition frequencies and thus contribute to the H2 continuum. The present measurements confirm that the dimer contribution is significantly weaker than the CIA but the most intense sharp lines of (H2)2 may represent up to 3% of the CIA level (Note that both have a quadratic pressure dependence).
In spite of its weakness, it is necessary to examine to which extent the (H2)2 spectrum may bias metrological measurements of the H2 quadrupole lines. In recent years, highly accurate measurements of the H2 (and HD) transition frequencies have been used to challenge the most accurate theoretical calculations, including relativistic and quantum electrodynamics (QED) effects (e.g. ref. 40). Experimentally, transition frequencies were determined with a few tens kHz error bars for a few transitions recorded either in saturation regime (Lamb dips) (e.g. ref. 41) or in the Doppler-limited regime (e.g. ref. 25). While Lamb dip measurements are unaffected by dimers, a possible bias on the H2 line centers determined in Doppler regime has to be considered. Although a detailed error budget remains to be performed, from the present observations, it appears that the spectral interference between the (H2)2 and H2 absorption is too limited to significantly bias the derived H2 transition frequencies.
More than 30 years after the last report on the absorption spectrum of the (H2)2 dimer and while all previous studies were performed at low temperature (e.g. 20 K in the fundamental band), the high sensitivity of the CRDS technique has allowed us to measure the first (H2)2 spectrum at room temperature (and sub-atmospheric pressures). As only a few low-l states give rise to clear absorption features, the observed (H2)2 spectrum appears to be mostly independent of the temperature.
Empirical parameters were retrieved for the dimeric features observed around ten allowed and forbidden H2 transition frequencies in the region of the first overtone band, near 8000 cm−1. We hope that the obtained absolute position and absorption values together with the provided widths of the sharp and broad features will stimulate theoretical efforts and serve as a benchmark for future validation tests. Indeed, in spite of its apparent simplicity, to the best of our knowledge, no first-principles calculations of the (H2)2 spectrum are available for a direct comparison with our infrared measurements. Due to the very small binding energy of the hydrogen van der Waals dimer and the absence of a preferred geometry, a highly accurate six-dimensional potential energy surface is required. The evaluation of the strengths and dissociation rates of the observed (H2)2 transitions is an additional challenge to predict the spectrum. The first quantitative intensity information on the (H2)2 absorption obtained in this work will be valuable to check the consistency between the 2H2 ↔ (H2)2 dimerization constant and the line strengths provided by theory. In the present context of an increasing use of molecular hydrogen as an energy vector, the knowledge of the (H2)2 abundance in highly pressurized tanks used for transportation at room temperature might have some industrial impact.
As a final conclusion, the present work opens new perspectives for direct room temperature detection of dimers or complexes of the major molecular species, with possible impact in the study of the Earth and planetary atmospheres.
This journal is © the Owner Societies 2024 |