Cuihong
Zhang‡
abc,
Xiaoxiao
Lin‡
a,
Xiaofeng
Tang
*a,
Christa
Fittschen
c,
Sebastian
Hartweg
d,
Gustavo A.
Garcia
*d,
Bo
Long
e,
Weijun
Zhang
a and
Laurent
Nahon
d
aLaboratory of Atmospheric Physico-Chemistry, Anhui Institute of Optics and Fine Mechanics, HFIPS, Chinese Academy of Sciences, Hefei, 230031 Anhui, China. E-mail: tangxf@aiofm.ac.cn
bGraduate School, University of Science and Technology of China, Hefei, 230026 Anhui, China
cUniversity Lille, CNRS, UMR 8522, PC2A – Physicochimie des Processus de Combustion et de l’Atmosphère, F-59000 Lille, France
dSynchrotron SOLEIL, L’Orme des Merisiers, St. Aubin BP 48, 91192 Gif sur Yvette, France. E-mail: gustavo.garcia@synchrotron-soleil.fr
eSchool of Materials Science and Engineering, Guizhou Minzu University, Guiyang, 550025 Guizhou, China
First published on 21st December 2021
We present a vacuum ultraviolet (VUV) photoionization study of the gas-phase sulfuric acid (H2SO4) molecule in the 11–14 eV energy range by using the method of synchrotron radiation-based double imaging photoelectron photoion coincidence (i2PEPICO) spectroscopy complemented with accurate theoretical calculations. The slow photoelectron spectrum (SPES) of H2SO4 has been acquired and the three electronic states of H2SO4+, X2A, A2A and B2A have been populated and assigned. The adiabatic ionization energy of the H2SO4 molecule towards the X2A cationic ground state is measured at 11.684 ± 0.006 eV, in accordance with high-level calculated findings. With increasing photon energies, the H2SO4+ cation dissociates into HSO3+ and OH fragments and their adiabatic appearance energy is measured at 13.498 ± 0.007 eV. Then, the enthalpies of formation for the species involved in the photoionization and dissociative photoionization have been determined through a thermochemical cycle.
There have been several theoretical studies on H2SO4 and sulfuric acid–water clusters, in which their optimized structures and spectra have been calculated using various levels of theory.9–12 For example, the theoretical results show that H2SO4 has two conformers depending on the relative position of the two OH groups: the most stable one is the trans form having a C2 symmetry, while the cis form has a C1 symmetry and an energy of +4.6 kJ mol−1 with respect to the trans conformer.9 The energy barrier of isomerization was calculated to be 12.6 kJ mol−1 relative to the trans conformer.9 Theoretical studies on the uptake of water molecules, leading to the hydration of H2SO4, have yielded the structure of these complexes, and it has been predicted that progressive hydration increases the probability of proton transfer from the acid to water, with the resulting ion-pair structures corresponding to the global energy minima.9,11
Experimentally, spectroscopic techniques such as microwave rotational spectroscopy, infrared (IR) spectroscopy and Raman spectroscopy have been employed to probe H2SO4 and hydrated complexes.13–17 The structures, vibrational spectra and frequencies of H2SO4 and the H2SO4·H2O cluster were obtained. Mass spectrometry, as a convenient analytical method, has also been applied in the past. For instance, the mass spectrum of H2SO4 vapor was acquired by Snow and Thomas using electron impact ionization mass spectrometry (EIMS) and its ionization energy was roughly measured at 12.4 eV, together with the appearance energies of its major fragment ions, H2O+, SO3+ and HSO3+.18 The positive and negative ion distributions of the sulfuric acid–water clusters as well as their equilibrium constants were measured by Froyd and Lovejoy using EIMS.19,20 Recently, the formation of mixed sulfuric acid–water clusters with different concentrations of water and H2SO4 has been investigated by electron attachment and negative ion mass spectrometry.21 However, the abovementioned mass spectrometry methods lack the energy resolution to distinguish contributions from all possible isomers. Moreover, much of the thermochemical data of H2SO4 and the sulfuric acid–water clusters are not precisely known.22
The purpose of this study is to investigate the vacuum ultraviolet (VUV) photodynamics of H2SO4. Photoionization and dissociative photoionization of H2SO4 in the energy range of 11–14 eV are analyzed by using the method of double imaging photoelectron photoion coincidence (i2PEPICO) spectroscopy complemented with high-level theoretical calculations. Three cationic electronic states of H2SO4+ have been prepared and assigned in the slow photoelectron spectrum (SPES). In particular, the adiabatic ionization energies (AIEs) of the trans and cis conformers of H2SO4, towards the X2A cationic ground state, have been measured precisely for the first time, as well as the adiabatic appearance energy (AE0K) of the HSO3+ fragment ion from the dissociation of H2SO4+. The enthalpies of formation for the species involved in the photoionization and dissociative photoionization are also determined through a thermochemical cycle.
The i2PEPICO spectrometer is composed of a modified Wiley–McLaren ion time-of-flight (TOF) 3D-momentum imaging device and an electron velocity map imaging (VMI) analyzer, equipped with two position-sensitive detectors (PSDs) for the analysis of ions and electrons, respectively.24 A mixture of helium carrier gas (1.5 bars), water vapor (from the vaporization of liquid water at room temperature) and H2SO4 vapor, generated from the vaporization of liquid H2SO4 (95% purity) in an in-vacuum oven heated at 180 °C, was expanded through a 70 μm diameter nozzle to form a continuous molecular beam which was skimmed twice before reaching the center of the spectrometer.25,26 The multiplex coincidence scheme between the ion and electron yields mass-selected photoelectron spectra (PES) as a function of photon energy, which are then reduced to SPES, as previously detailed.27,28 Meanwhile, the coincidence scheme can provide spectral fingerprints to identify and separate isomers, and accurately measure thermochemical data such as fragment appearance energies and bond energies.29–31
We have performed theoretical calculations to assign the experimental spectra. The AIE and the vertical ionization energy (VIE) of the X2A ground electronic state of H2SO4+ were calculated at the CCSD(T)-F12A/aug-cc-pVTZ level of theory as implemented in MOLPRO software.32 The VIEs of the first (A2A) and the second (B2A) excited electronic states of H2SO4+ were computed using time-dependent density functional theory (TD-DFT), TD-M062X/6-311++G(d,p).33 The Franck–Condon factors in the photoionization towards the X2A cationic ground state, as well as for the required geometry optimizations and the vibrational frequency computations with a harmonic approximation, were calculated at the M062X/6-311++G(d,p) level of theory with the Gaussian package using the time-independent adiabatic Hessian Franck–Condon model.33
Fig. 1 (a) Synchrotron photoionization mass spectra matrix in the 11–14 eV energy range and (b) their integrated mass spectrum, together with 5 times magnified data in red. |
SO3 is formed from the equilibrium reaction of H2SO4(g) ↔ SO3(g) + H2O(g) inside the oven, and SO2 is produced from the reduction reaction of H2SO4 with the metallic parts of the oven. The appearance energies of SO3 and SO2 can be obtained from the mass spectra matrix of Fig. 1(a) and both of them agree with their individual ionization energies,35 confirming their assignments as neutral nascent species formed in the oven. Note that SO3+ and SO2+ were also observed in a previous study based on EIMS, but both were assigned as fragment ions from the dissociation of H2SO4+ by Snow and Thomas.18
Normally, in Fig. 1, the mass peak at m/z = 81 should be assigned as the 33SO3 isotopic species in the mass spectra. But, its relative intensity (∼28% of 32SO3) measured in the mass spectrum is much higher than the natural abundance of 33SO3 (∼0.91%),34 and its peak width is much wider than that of the SO3 or H2SO4 mass peak. In addition, as shown in the mass spectra matrix, the appearance energy of the m/z = 81 species is not the same as and much higher than that of SO3 (m/z = 80). Therefore, the m/z = 81 peak in the mass spectra is not contributed from the isotope 33SO3 and corresponds to the HSO3+ fragment ion from the dissociative ionization of H2SO4, as inferred from its broad shape due to the kinetic energy release in the fragmentation. Note that the contribution from the isotope 33SO3 (∼0.91%)34 to the m/z = 81 channel is considered negligible and has not been subtracted here. Similarly, due to its large width, the more intense peak at m/z = 99 is attributed to the H2SO4H+ fragment ion produced from the dissociative photoionization of the H2SO4·H2O complexes formed in the molecular beam, with again a negligible contribution from the H233SO4+ isotopic species. Note that a complete study of water complexes and clusters of sulfuric acid will be reported in a forthcoming work.
The mass-selected SPES corresponding to the H2SO4+ parent ion and the HSO3+ fragment ion are obtained from the photoelectron matrices and presented in Fig. 2(b and d), respectively, as well as their photoionization yields (PIYs).27,28 In Fig. 2(b), the SPES of H2SO4 is rather unstructured although an onset can be observed at hν = 11.684 ± 0.006 eV in the SPES. The uncertainty of the photon energy calibration, the influence of the scanning step size (5 meV) and the signal-to-noise ratio of the SPES have been taken into account and included in the given error bars.36
As mentioned in the introduction, the gas-phase H2SO4 molecule has two conformers, trans and cis, and their optimized structures as the neutrals as well as their cations in the X2A cationic ground state have been calculated and are presented in Fig. 3, agreeing well with previous results.9,37 The trans conformer has a C2 symmetry structure and is the global minimum with an energy of ∼4.6 kJ mol−1 below the cis conformer with a C1 symmetry.9 The temperature of the molecules in the molecular beam is estimated at 150 K (see Section 3.3 below), which would mean that the trans conformer is dominant, with only ∼2% for the cis conformer. However, the flash cooling in the adiabatic expansion might lead to a different conformer population than the one expected assuming thermal equilibrium.
Fig. 3 Optimized structures of H2SO4 (trans and cis conformers) and HSO3, as well as their individual cations, in the ground electronic states. |
Upon photoionization, an electron is removed from the outer valence orbital of the neutral H2SO4 molecule leading to the formation of the H2SO4+ cation. The structure of the trans conformer relaxes in the cation and its symmetry changes from C2 to C1 during the photoionization, as shown in Fig. 3, while the cis conformer geometry undergoes smaller changes. The AIE of the X2A ground electronic state of the trans-H2SO4+ has been calculated at 11.688 eV using the CCSD(T)-F12A/aug-cc-pVTZ method with zero-point corrections (ZPEs).32 Note that no scaling factor has been used for the determination of harmonic frequencies and ZPEs. The calculated values are consistent with the onset of the experimental curve and their differences are within 10 meV. Although some structures are seen in the region where the AIEs are predicted in our experimental spectra, the unfavorable Franck–Condon factors and consequent lack of structure do not allow a definite assignment of the experimental AIE. Nevertheless, in view of the good agreement between the calculated AIE and the feature at 11.684 eV, we tentatively offer this value as the AIE of the trans conformer. The very weak features at 11.643 and 11.609 eV (marked with two asterisks in the 10 times magnified data of Fig. 2b) might be ascribed to the cis conformer, with the former being assigned to the AIE of the X2A ground state, for which the value is calculated at 11.635 eV, and the latter originating from a hot band in the photoionization, i.e., from the v2 vibrational mode of the cis conformer whose frequency is calculated at 289.32 cm−1 using the M062X/6-311++G(d,p) method. However, we note that purely thermal arguments would not allow for the presence of the cis conformer, nor a hot band, at the presently estimated temperature of 150 K. In a previous study based on EIMS, restricted by the electron energy resolution, the ionization energy of H2SO4 was measured at 12.4 ± 0.05 eV (see Table 1), which would actually correspond to the VIE measured in this work.18
State | IE(H2SO4) | AE(HSO3+) | Method | Ref. | |
---|---|---|---|---|---|
a AIE. b VIE. c CCSD(T)-F12A/aug-cc-pVTZ. d TD-M062X/6-311+ +G(d,p). e This work. | |||||
Expt. | 12.4 ± 0.05 | EIMS | 18 | ||
13.9 ± 0.1 | EIMS | 18 | |||
13.498 ± 0.007 | SPES | ||||
X2A (trans) | 11.684 ± 0.006a | SPES | |||
X2A (cis) | 11.643 ± 0.006a | SPES | |||
X2A | 12.53 ± 0.05b | SPES | |||
A2A | 12.80 ± 0.05b | SPES | |||
B2A | 13.20 ± 0.05b | SPES | |||
Calc. | X2A (trans) | 11.688a | CCSD(T)c | ||
X2A (cis) | 11.635a | CCSD(T)c | |||
X2A | 12.54b | CCSD(T)c | |||
A2A | 12.84b | TD-M062Xd | |||
B2A | 13.10b | TD-M062Xd |
The Franck–Condon factors for the photoionization of H2SO4 towards the cationic X2A ground electron state of H2SO4+ have been calculated at the M062X/6-311++G(d,p) level of theory, using the time-independent adiabatic Hessian Franck–Condon model.33 Note that the influences of inconsequential internal rotations are not included in this model and thus, the differences of the two conformers are not considered in the Franck–Condon factors. The simulated PES is subsequently generated by convolving the stick spectrum with a Gaussian function (FWHM = 100 cm−1, full width at half maximum) and is shown in Fig. S1 (ESI†), with an energy offset of 50 meV to fit the experimental result. Although the simulated PES fits the experimental SPES very well at the beginning of the X2A state, an apparent difference can be observed between them at the high energy part, indicating that other factors need to be considered, such as the appearance of cationic excited electronic states. Indeed, the VIEs of the A2A excited electronic state and the B2A electronic state of H2SO4+ have been calculated at 12.84 and 13.10 eV, respectively, at the TD-M062X/6-311++G(d,p) level of theory, and have been assigned in the experimental SPES too. In the SPES of Fig. 2(b), the VIEs of the A2A and B2A cationic states are estimated at 12.80 ± 0.05 and 13.20 ± 0.05 eV, respectively, as listed in Table 1, agreeing well with the calculated results.
Based on the AIE of the X2A ground state of H2SO4+ and the AE0K of the HSO3+ fragment ion, the bond energy of H2SO4+ associated with the HSO3+ and OH fragments is calculated at 1.814 ± 0.008 eV for the trans conformer. Through a thermochemical cycle, as shown in Fig. 5, the enthalpies of formation for the cations of H2SO4+ and HSO3+ are obtained at ΔfH0K(H2SO4+) = 406.2 ± 8.5 kJ mol−1 and ΔfH0K(HSO3+) = 543.9 ± 8.5 kJ mol−1, with the aid of the available thermochemical data in the literature, i.e., ΔfH0K(H2SO4) = −721.2 ± 8.4 kJ mol−139 and ΔfH0K(OH) = 37.252 ± 0.026 kJ mol−1,40 and assuming barrierless fragmentations. In addition, the AIE of the HSO3 radical is calculated at 9.533 eV with the CCSD(T)-F12A/aug-cc-pVTZ high-level of theory,32 so that we report the enthalpy of formation for HSO3 at ΔfH0K(HSO3) = −375.9 ± 8.5 kJ mol−1.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cp05237c |
‡ These authors contributed equally to this work. |
This journal is © the Owner Societies 2022 |