Introduction

Sulfur is the tenth most abundant element in the universe, and sulfur chemistry plays vital importance not only in the biological systems and atmosphere on the Earth but also in interstellar medium (ISM). According to astronomical observations, the two sulfur oxides SO2 and SO have been found to be abundant in molecular clouds1,2,3,4, which are mainly formed through condensation of gas-phase molecules at the surface of dust grains (mostly amorphous silicates) with an onion-like structure at the temperatures of about 10–20 K5,6. The inner layer of the grains mainly consists of hydrogenated ice (H2O) with low concentrations of other H-containing species such as CH3OH, NH3, and CH4. The outer layer is made up of dehydrogenated ices with dominant compositions of CO, CO2, N2, O2, and SO2, and low concentrations of H2O may also be present in the outer layer of the icy mantle. The icy mantle at the surface of cosmic dust grains are the most important carriers of prebiotic molecules, and the composition of the mantles are largely affected by the exchanges between solid ice and gas-phase and also the photochemistry promoted by cosmic irradiations, including UV and X-ray photons from nearby stars. Therefore, the study about the chemical composition of the icy grains and the complex reaction networks is crucial for understanding the evolution of the molecular clouds7,8,9.

Carbon monoxide (CO) is the most abundant composition in the outer layer of icy grains in interstellar medium, and CO-abundant ices have also been found at the surface of many cold interstellar bodies, including comets, icy moons, and planets in the outer solar system. Therefore, the chemistry of CO through successive hydrogen atom addition reactions in the CO-rich outer layer of the interstellar icy grains may play a key role for the formation of complex organic molecules (COMs), which are probably building blocks for the origin of life10,11,12. Recently, it has been shown that simple radicals bearing elements C, N, O, P, or S play vital importance in prebiotic synthesis13. Reaction networks of these radicals in interstellar ice grains and the corresponding geochemical scenarios may help in unveiling the chemical evolution and origins of life. For instance, formyl radical (HCO•) and hydroxycarbonyl radical (HOCO•) are important intermediates in atmospheric and combustion chemistry14, and they are also versatile building blocks in the interstellar formation of biorelevant COMs such as formic acid (HC(O)OH)15, glyoxylic acid (HC(O)C(O)OH)16, and pyruvic acid (CH3C(O)C(O)OH)17 in low temperature (<30 K) interstellar CO and CO2 ices doped with H-containing species CH4 and H2O through barrierless radical-radical association reactions, in which the reactive acyl radicals can be generated through the hydrogenation of carbon oxides with the H-containing molecules embedded in the same ice layer at cosmic radiations18,19. As a simple organic species, HCO• has been observed in many interstellar clouds such as DR 21, Sgr B2, and NGC 202420,21, and its radical recombination reaction with •CH2OH in producing COMs during the phase transition of interstellar CO ices doped with CH3OH and H2O at a typical dense cloud temperature of about 10 K has been recently disclosed22,23,24. The cationic form of HOCO• has been also identified in star-forming regions such as SgrB2(OH) and low‑mass protostar IRAS 04368 + 2557 in L152725,26.

In sharp contrast to the extensively explored mechanisms for the formation of COMs through the photoreactions of H-containing species (e.g., CH3OH, NH3, and CH4) via the intermediacy of organic radicals such as HCO•, HOCO•, and CH3O• in interstellar icy grain mantles, the ice-grain photochemistry of the typical dehydrogenated molecules such as SO, SO2, and the derived sulfur-containing radicals HOS• and HOSO• in astronomical CO and CO2 ices remains barely investigated. On the other hand, the photochemistry of SO and SO2 also has great impact on the sulfur cycle in planetary atmospheres due to the formation and evolution of hazes and clouds in the upper atmospheres of Solar system planets such as Earth27, Venus28, Jupiter29, and the moon Io30. Among these sulfur oxides, SO2 is one of the most common pollutant in the Earth’s atmosphere. Therefore, the SO2 photochemistry has been the focus of enormous attention due to the important role in sulfur cycle by forming sulfuric acid and sulfate aerosols. According to the recent modeling studies31,32, the tropospheric photochemistry of SO2 in the presence of water proceeds mainly through the formation of HOSO• and hydroxyl radical (•OH) after absorption of light in the near UV‒vis range (250–340 nm). In contrast, the UV photolysis (190–220 nm) of SO2 in the gas phase leads to fragmentation by yielding SO33, which is an interstellar species that has been detected in the atmospheres of Venus34 and Io35. Chemically, SO is more reactive and it dimerizes easily to yield elemental sulfur and SO2 via the intermediacy of OSSO, and the dimer has been recognized as a candidate species that contributes to the mysterious near‑UV absorption (320–400 nm) in the yellowish atmosphere of Venus36,37,38.

Herein, we report an experimental study on the photochemistry of the astrochemically relevant sulfur-containing species HOSO•, HOS•, SO2, SO, and OSSO in solid CO and CO2 ices at 16 K (Fig. 1). In addition to the molecular complexes formed between HOSO• and the carbon oxides, the photo-induced hydrogen atom transfer (HAT) to form new complexes consisting of acyl radicals (HCO• and HOCO•) and sulfur oxides (SO2 and SO) has been observed. It is noteworthy that weakly bonded molecular complexes consisting of interstellar species have been considered as potent contributors to the rich chemistry in low-temperature giant molecular clouds39.

Fig. 1: Photochemistry of HOSO• in solid CO and CO2 ices.
figure 1

Proposed pathways for the photo-induced decomposition, isomerization, and hydrogen atom transfer reactions of HOSO• in solid CO and CO2 at 16 K.

Results and Discussion

Isolation of HOSO• in CO and CO2 ices

Thanks to the strong “cage effect” of the solid host matrix materials (e.g., Ne and Ar) at low temperatures (<30 K), the matrix isolation technique has been broadly applied in trapping highly unstable intermediates, including weakly bonded molecule‑radical complexes such as •OH···CO18, •C6H5···H2O40, •OC6H5···H2O41, •OH···H2O42, and •NH2···H2O43. Recently, it has been shown that HOSO• can be efficiently generated in the gas phase through high‑vacuum flash pyrolysis (HVFP, ca. 700 °C) of CHF2S(O)OH (Fig. 2a), and photolysis of HOSO• (Fig. 2b) in solid Ar‑matrix at 10 K yields isomeric HSO2•, fragments SO2/H• together with the caged radical complex •OH···OS44. The absence of free fragments •OH and OS in the matrix suggests that they can hardly escape from the rigid matrix cages. When the pyrolytic generation of HOSO• was performed in the presence of CO by using a 1:20:1000 mixture of CHF2S(O)OH/CO/Ar, the IR spectrum of the isolated species (Supplementary Fig. 1) shows the appearance of new IR bands in the range of 3500–3350 cm−1 for O–H stretching vibrations (ν(OH)), implying complex formation between HOSO• and CO.

Fig. 2: Infrared (IR) spectrum of CHF2S(O)OH and HOSO• in matrixes.
figure 2

a IR spectrum of CHF2S(O)OH in Ar-matrix at 10 K. b IR spectrum for the high‑vacuum flash pyrolysis (HVFP, ca. 700 °C) products of CHF2S(O)OH in Ar-matrix at 10 K. c IR spectrum for the HVFP (ca. 700 °C) products of CHF2S(O)OH in CO-matrix at 16 K.

When using neat CO as the matrix host material, the IR spectrum of the isolated pyrolysis products at 16 K (Fig. 2c) clearly shows the absence of all the IR bands for free HOSO• (Fig. 2b), while the bands for the complex CO···HOSO• become dominant. The ν(OH) mode in HOSO• shifts from 3545.3 cm−1 in Ar‑matrix to 3396.8 cm−1 in CO‑matrix, corresponding to a red‑shift (Δν) of –148.5 cm−1. It is comparable with the shift of the ν(OH) mode in CO‑matrix isolated HOCO• (Δν = –146.9 cm−1) comparing to its IR spectrum in Ar‑matrix45. The assignment of CO···HOSO• is supported by the good agreement with the theoretically calculated shift of –140 cm−1 at the B3LYP‑GD3(BJ)/def2‑TZVP level (Table 1). In line with a stable hydrogen-bonded structure through OC···H bond in the complex, the stretching mode for the terminal S = O moiety (ν(S = O)) is less perturbed than the S–O stretching mode (ν(S–O)) as indicated by the shifts of –3.6 and +25.3 cm−1, respectively. In contrast, the weak deformation mode δ(SOH) exhibits a large blue‑shift of +44.2 cm−1, which is in agreement with the calculated shift of +53 cm−1 at the CCSD(T)/aug‑cc‑pV(T + D)Z level. Note there are several weak satellite bands around the main peaks for the fundamental modes of the complex, they probably arise from the less stable matrix sites or the less abundant complexes of HOSO• consisting two or more CO molecules.

Table 1 Calculated and observed Infrared (IR) spectra of HOSO• in different matrixes

By analogy, deposition of the pyrolysis products of a 1:50:1000 mixture of CHF2S(O)OH/CO2/Ar at 16 K leads to the formation of the complex CO2···HOSO• (Supplementary Fig. 2). Consistent with the CCSD(T) calculated red-shift of –73 cm−1 (B3LYP: –84 cm−1) for the ν(OH) mode in the complex, the band for HOSO• shifts from 3545.3 cm−1 in Ar‑matrix to 3478.7 cm−1 in CO2‑matrix (Table 1), corresponding to a shift of –66.6 cm−1. Concomitantly, the ν(S–O) mode undergoes a blue-shift by +22.4 cm−1 (CCSD(T): + 21 cm−1; B3LYP: + 24 cm−1).

Weak hydrogen‑bonding interactions of HOSO• with carbon oxides also affect its UV-vis absorption. Recently, a broad absorption centered at 270 nm (λmax) has been observed for HOSO• in Ar-matrix, corresponding to the transition from the ground state (X2A) to the C2A/D2A excited states46. As shown in Fig. 3, the major absorption of HOSO• in CO‑matrix appears in the range of 350–240 nm as a weak band, and its assignment is ascertained with the photochemistry that HOSO• can be efficiently depleted by UV‑light irradiation at 266 nm47. This characteristic absorption for HOSO• is also observable in CO2‑matrix. In sharp contrast to the appearance of the absorption of HSO2• in the range of 320 − 500 nm after the irradiation (266 nm) of HOSO• in Ar‑matrix, same photolysis in CO‑ and CO2‑matrixes results in the occurrence of weaker absorptions in the range of 300 − 400 nm, implying the formation of new species arising from the photoreactions of HOSO• with the carbon oxides.

Fig. 3: Ultraviolet-visible (UV‒vis) spectra of the acid precursor CHF2S(O)OH and its high-vacuum flash pyrolysis (HVFP) products isolated in solid CO and CO2 ices at 16 K.
figure 3

The absorptions for HOSO• and its photolytic reaction products in the matrixes are rendered with orange and green colors, respectively.

Photochemistry of HOSO• in CO and CO2 ices

To unravel the photochemistry of HOSO• in CO‑ and CO2‑matrixes, the IR spectra for the 266 nm laser photolysis products were recorded and the resulting difference spectra reflecting the reactions of HOSO• are depicted in Fig. 4. In contrast to the photodissociation of HOSO• to H•/SO2 and the caged complex •OH···SO in Ar‑matrix (Fig. 4a), its photolysis in solid CO (Fig. 4b) yields CO2, OCS, HCO•, HOCO•, H2CO, HOS•, and SO2. It is noteworthy that the IR frequencies for all these species shift slightly in comparison to those observed for the corresponding species in CO‑matrixes, indicating weak interactions between the neighboring counterpart species formed after the bimolecular reaction of HOSO• with CO inside the rigid CO‑matrix cages. For instance, the two characteristic IR bands of t‑HOCO• for the ν(O–H) and ν(C=O) modes at 3456 and 1833 cm−1 in CO‑matrix shift to 3311.0 and 1831.7 cm−1 due to complexation with SO (Supplementary Table 1). Concomitantly, the IR band of the counterpart SO at 1139.5 cm−1 (CO‑matrix) undergoes blue‑shift to 1140.0 cm−1, and the deformation mode δ(COH) exhibits a larger blue‑shift of +12.8 cm−1. The presence of the less stable conformer c‑HOCO• is evidenced by the band at 1795.2 cm−1, and it is also slightly red‑shifted comparing to the band at 1797 cm−1 for c‑HOCO• in CO‑matrix18. Conformational conversion of c‑HOCO• to the lower‑energy t‑HOCO• happens upon subsequent irradiation at 532 nm. However, the previously reported45 spontaneous transformation of c-HOCO• → t‑HOCO• in N2‑matrix (4.5 K) via quantum mechanical tunneling was not observed in CO-matrix (16 K), which is consistent with the frequently observed environmental effects on the tunneling processes in low‑temperature matrixes48.

Fig. 4: Reactions of HOSO• in matrixes upon irradiation.
figure 4

a Infrared (IR) difference spectrum reflecting the change of the Ar‑matrix isolated high‑vacuum flash pyrolysis (HVFP) products of CHF2S(O)OH upon irradiation at 266 nm (80 min, 10 K). b IR difference spectrum reflecting the change of the CO‑matrix isolated HVFP products of CHF2S(O)OH upon irradiation at 266 nm (7 min, 16 K). c IR difference spectrum reflecting the change of the CO2‑doped Ar-matrix (50: 1000) isolated HVFP products of CHF2S(O)OH upon irradiation at 266 nm (30 min, 16 K). d IR difference spectrum reflecting the change of the Ar‑matrix isolated HVFP products of 18O‑labeled CHF2S(O)OH upon irradiation at 266 nm (32 min, 10 K). e IR difference spectrum reflecting the change of the CO‑matrix isolated HVFP products of 18O‑labeled CHF2S(O)OH upon irradiation at 266 nm (22 min, 16 K). The symbol “16/18O” refers to a 1:1 mixture of the species containing 16O and 18O.

Consistent with the photodecomposition of HOSO• (→ H• + SO2) in Ar‑matrix, its photolysis in solid CO also causes H–O bond fragmentation followed by CO‑trapping of the mobile hydrogen atoms to afford HCO•, which acts as a hydrogen donor through weak interaction with the counterpart SO2 in the same CO-matrix cage by forming complex HCO•···SO2. In this complex, the ν(C‑H) mode shifts to 2493.8 cm−1 in comparison to the same mode at 2488 and 2483 cm−1 for HCO• in CO‑ and Ar‑matrixes49. The ν(C = O) and δ(COH) modes in HCO• and the two stretching modes of SO2 display small red-shifts (Supplementary Table 2) in the complex. The changes of the two stretching modes of SO2 in HCO•···SO2 (Δνasym = –8.2 cm−1; Δνsym = –0.3 cm−1) are smaller than those observed in other SO2‑contaning complexes such as H2O2···SO2 (Δνasym = –12.9 cm−1; Δνsym = –2.1 cm−1)50. Further combination of HCO• with mobile hydrogen atoms in the matrix during the laser photolysis affords H2CO, and it is also perturbed by the neighboring molecule (SO2) as evidenced by the appearance of two IR bands at 1736.7 and 1734.9 cm−1 for the ν(C = O) mode in H2CO. Note, formation of HCO• and H2CO has been previously observed in a solid CO/H2 ice mixture at 8 K after irradiation with ultra‑high vacuum UV (~160 nm) light19.

The mechanism for the formation of HOS• during the photolysis of HOSO• in CO-matrix is intriguing. A plausible pathway is the HAT in HOCO•···SO ( → CO2···HOS•). In the CO2···HOS• complex, the ν(O–H) mode shifts by –154.3 cm−1 in comparison with the same mode observed in HOS• in solid para‑H251, while the ν(S–O) mode displays a blue-shift of +8.9 cm−1 (Supplementary Table 3). Accordingly, the two IR bands for CO2 at 2346.7 and 659.3 cm−1 exhibit shoulders that can be assigned to the complex (Fig. 5). Particularly, the bending mode δ(CO2) in the complex appears as a broad band in the range of 665–650 cm−1 due to removal of vibrational degeneracy upon complexation with HOS•. Similar splitting for the nondegenerate δ(CO2) modes has been observed in other CO2‑complexes such as HKrCCH···CO252.

Fig. 5: Sections of Infrared (IR) spectra showing the two bands for the two vibrational modes of CO2 in CO‑matrix isolated CO2 and CO2···HOS•.
figure 5

The lower trace corresponds to the IR spectrum of CO2 in CO-matrix at 16 K. The upper trace corresponds to the IR spectrum of CO2···HOS• in CO-matrix at 16 K.

The spectroscopic identification of the photoproducts of HOSO• in CO‑matrix is further supported by 18O‑isotope labeling experiments. Using 18O‑labeled sulfinic acid as the precursor, a 1:1:1:1 mixture of HOSO•, H18OSO•, HOS18O•, and H18OS18O• can be generated, and their distinction can be assured with the revolved IR bands at 1049.7, 1045.2, 1043.4, and 1039.5 cm−1 for the δ(SOH) mode in the four isotopologues, however, the remaining IR fundamental modes appear as doublets due to closeness of the isotopic shifts (Fig. 4d). The photochemistry of these isotopologues in CO-matrix (Fig. 4e) provides useful information for probing the reaction mechanism. The sole formation of HCO• (without HC18O•) and a 1:2:1 mixture of SO2, OS18O, and S18O2 among the photoproducts confirms the route for the straightforward HAT from HOSO• to CO. The absence of the IR bands for HOC18O• and H18OC18O• rules out the possibility for the formation of HOCO• via direct hydrogenation of CO2, since the singly 18O‑enriched CO2 is present in the same matrix. In contrast, the association of the photolytically generated •OH/•18OH with CO yields the experimentally observed HOCO• and H18OCO• in 1:1 ratio.

In the IR difference spectrum for the 18O-labeling experiments (Fig. 4e), each band for the ν(O–H) and ν(S–O) modes in HOS• splits into doublet due to the additional presence of H18OS•, the corresponding 16/18O-isotopic shifts –10.8 and –31.5 cm−1 show good agreement with the calculated values –8 and –31 cm−1, respectively. Assuming a bimolecular reaction of HOCO• and SO for the formation of HOS• and CO2, the observation of more CO2 than HOS• (based on the experimental and calculated IR band intensities) in the photochemistry indicates that there is an alternative pathway for producing CO2. Given the generation of HO• and SO in the photolysis of HOSO•, the photochemistry of SO in CO ice was also studied. Co-condensation of gaseous SO in the presence of CO (ca. 1:1000) at 16 K also yields OSSO, which can be completely destroyed by UV‑light irradiation at 365 nm (Supplementary Fig. 3). The formation of SO2 and OCS coincides with the photodissociation of OSSO to SO2 and sulfur atom with subsequent CO-trapping reaction (S + CO → OCS). When changing the irradiation source to a 193 nm laser, depletion of SO occurs by forming CO2 and OCS (Supplementary Fig. 3), and the mechanism can be reasonably explained that monomeric SO dissociates to sulfur and oxygen atoms followed by association reactions with CO. Therefore, formation of CO2 and OCS in the photochemistry of HOSO• in solid CO is very likely caused by the photofragmentation of the initially generated SO. This mechanism is also consistent with the sole observation of none-isotopically labeled OCS and 1:1 mixture of CO2/OC18O in the 18O‑lalebing experiments (Fig. 4e). Additionally, traces of HCO• form during laser irradiation of the impurity H2O in solid CO.

The photolytic depletion of HOSO• in CO‑matrix at 266 nm is extremely fast, and the depletion becomes much less efficient under subsequent UV‑light irradiation at 365 nm (Supplementary Fig. 4 and 5). Prolonged 266 nm laser irradiation (36 min) leads to complete depletion of HOSO•. Further irradiation at 365 nm promotes reverse HAT from HOCO• to SO2 by reforming HOSO• and CO2. The photosensitivity of HOCO• towards the UV‑light (365 nm) strongly indicates that the aforementioned weak absorptions in the UV‑vis spectra in the range of 300 − 400 nm (Fig. 3) belongs to this carbonyl radical. It is also consistent with the observation in the IR spectra that HOCO• forms after the photolysis of HOSO• in CO‑matrix (Fig. 4). The absorption for HCO• above 400 nm53 was not observed in the UV‑vis spectra due to low intensity.

The photo‑induced HAT of HOSO• also occurs in CO2-doped Ar-matrix at 16 K (Fig. 4c), yielding a new molecular complex HOCO•···SO2 (Supplementary Table 1). The strongest band for the ν(C = O) mode in the complex locates at 1818.9 cm−1, and it is lower than the same mode in HOCO•···SO and HOCO• at 1831.7 and 1833 cm−1, respectively. The very weak δ(COH) mode undergoes small blue‑shift to 1278.7 cm−1 comparing to the frequencies of 1261 and 1265.8 cm−1 for HOCO• in solid CO‑matrix and in HOCO•···SO. In contrast, the two SO2 stretching modes at 1350.9 and 1146.4 cm−1 are red-shifted in comparison to free SO2 at 1355.0 and 1152.2 cm−1. Unlike the generation of two conformers of HOCO• in the photochemistry of HOSO• in CO‑matrix (Fig. 4b), only one conformer was generated in the photolysis of HOSO• in CO2‑doped Ar‑matrix (Supplementary Fig. 2). Additional formation of OCS among the photolysis products implies the trapping reaction of sulfur atoms (SO → S + O) by CO2, however, the elusive intermediate OSCO54 was not observed. Formation of OCS and HCO• was also observed during the 193 nm laser irradiation of a mixture of H2O and SO2 in CO‑matrix (Supplementary Fig. 6).

Calculated structures of molecule-radical complexes

Weakly bonded complexes consisting of simple radicals (e.g., HO• and HOO•) and small molecules (e.g., CO, CO2, and H2O) play important roles in gas‑phase chemistry, as they may serve as key intermediates in clusters formation in atmospheric and astrochemical processes at the low‑temperature surface of dust and ice grains39. Therefore, the structures and reactivity of these molecule‑radical complexes have been the focus of comprehensive experimental and computational studies. Despite the importance of HOSO•, HOCO•, and HCO• in gas-phase chemistry has been increasingly recognized, their complexes remain scarcely investigated. Using the UCCSD(T)/aug‑cc‑pV(T + D)Z method, the structures, energies, and bonding properties for the new molecule-radical complexes involved in the photochemistry of HOSO• in CO and CO2 ices were calculated (Fig. 6).

Fig. 6: The CCSD(T)/aug-cc-pV(T + D)Z calculated molecular structures and noncovalent interaction analysis for all the observed complexes.
figure 6

The bond lengths (Å) for the monomers are given in square brackets, and the lengths for the intermolecular hydrogen bonds are shown in blue. The calculated gradient isosurfaces (s = 0.6 au) are colored on a blue‑green‑red scale according to values of sign(λ2)ρ, ranging from –0.035 to 0.02 au. Blue indicates strong attractive interactions, and red indicates strong repulsive interactions.

In sharp contrast to a favorable cis‑planar structure for free HOSO•, the HOSO moiety in CO···HOSO• and CO2···HOSO• are nonplanar. The global minimum of CO···HOSO• prefers OH···CO linkage and the hydrogen atom is dramatically tilted out of the OSO plane by 33.9°. The corresponding hydrogen bond length is 2.168 Å, which is shorter than the hydrogen bond in •OH···CO (2.341 Å)55. The binding energy (De = 3.9 kcal mol−1) is higher than those in •OH···CO (De = 2.3 kcal mol−1)56 and •OH···OS (De = 3.2 kcal mol−1)57. The isomer bearing a OH···OC linkage between HOSO• and CO is less stable by ca. 5 kcal mol−1. The hydrogen atom in CO2···HOSO• is also tilted from the OSO plane by a dihedral angle of 30.5°, and the CO2 molecule is slightly bent by 1.8°. The shorter hydrogen bond in CO2···HOSO• (1.990 Å) than that in CO···HOSO• (2.168 Å) coincides with the larger stabilizing interaction (De = 5.3 kcal mol−1). The noncovalent interaction (NCI) analysis suggests strong attractive hydrogen bond interactions and simultaneous weak repulsive interactions of the terminal S = O moiety with CO and CO2 in these complexes (Fig. 6).

Unlike the distortion of the complexed cis‑HOSO•, the HOCO moiety in the hydrogen-bonded HOCO•···SO and HOCO•···SO2 complexes keeps a favorable trans‑planar configuration. The shorter hydrogen bond in the HOCO•···SO (1.893 Å) than that in the latter (1.909 Å) is also consistent with a higher stabilizing energy of 6.1 kcal mol−1 (5.4 kcal mol−1 in HOCO•···CO2). The binding energy in HOCO•···SO is larger than that in the hydrogen-bonded H2O···SO complex (De = 3.1 kcal mol−1)58, while the latter remains yet experimentally unobserved. The weak interaction between HCO• and SO2 in HCO•···SO2 (De = 3.5 kcal mol−1) is facilitated by forming a five‑membered ring through concerted contacts of hydrogen bond CH–OS (2.590 Å) and chalcogen bond CO–SO (2.978 Å). Similar five‑membered ring structure has also been predicted for the detectable HO2•···SO2 complex (De = 4.6 kcal mol−1)59. The hydrogen bond interaction in the planar CO2···HOS• complex (De = 4.5 kcal mol−1) is stronger than the intermolecular O···C contact in the T‑shaped van der Waals complex H2O···CO2 (De = 2.8 kcal mol−1)60.

Energetically, HCO•···SO2 and HOCO•···SO are higher than CO···HOSO• by 37.5 and 27.1 kcal mol−1, respectively, whereas the secondary HAT in HOCO•···SO to form CO2···HOS• is highly exothermic by releasing –51.5 kcal mol−1. Therefore, the overall process for the oxidation of CO to CO2 by reaction with HOSO• is exothermic by –14.0 kcal mol−1, which is comparable with the energy (−20.0 kcal mol−1) for the •OH radical promoted oxidation14. The HAT process in CO2···HOSO• to form HOCO•···SO2 is endothermic by 38.5 kcal mol−1, which is lower than the H–O bond dissociation energy in HOSO• (44.1 kcal mol−1)57.

Implications in interstellar sulfur chemistry

Our experimental results demonstrate that organic radicals HCO• and HOCO• can be produced by UV-irradiation of HOSO• in astronomical CO and CO2 ices via the hydrogen atom transfer (HAT) reactions of the initially formed molecular complexes at 16 K, and the acyl radicals also form stable molecular clusters with sulfur oxides through strong hydrogen bonding interactions. Considering the facile formation of HOSO• from the photoreactions of SO2 with H2O, the generation of the two important building blocks HCO• and HOCO• from the photoreactions mimics the chemical evolution network of the dehydrogenated sulfur-containing molecules SO and SO2 in the outer layer of the CO-dominant interstellar icy grains in molecular clouds at a typical temperature of about 10 K. Alternatively, other radicals (e.g., •CH3, •CH2OH, and •CN) derived from the photoreactions of the interstellar carbon- or nitrogen-containing molecules (e.g., CH4, CH3OH, and HCN) in the cryogenic astronomical ices may also be present and react further with the acyl radicals through barrierless radical-radical association reactions to form more complex organic molecules (COMs). Additionally, the formation of OCS in the photochemistry of HOSO•, SO and SO2 in CO and CO2 ices may also contribute to understanding the interstellar sulfur chemistry, since OCS not only serves as a prebiotic activating agent for amino acid polymerization in forming peptides under mild conditions in aqueous solution61,62,63, also it involves in the reduction of CO264 and acts as a condensing agent in phosphate chemistry65.

Clearly, the HAT processes of HOSO• may serve as the link connecting the chemistry of SO2 and the chemistry of carbon oxides (CO and CO2) in interstellar ices, although the ubiquity of HAT in chemistry, biology, and industry has been well recognized66. The uncovered chemical network for the formation of the astrochemically relevant organic radicals HCO• and HOCO• in the simple systems containing the primordial substrates (e.g., H2O, SO, SO2, CO, and CO2) might aid in disclosing the intriguing mechanism for the chemical evolution of biomolecules such as organic acids in dense molecular clouds, where barrierless radical-radical reactions at low temperatures (10–20 K) are assumed to happen spontaneously for the formation of COMs67,68. Furthermore, the spectroscopic identification and photochemistry of the new complexes consisting of the astrochemically relevant radicals HOSO•, HCO•, HOCO•, and HOS• will help in understanding the chemical composition and abundances in the interstellar medium, since molecular complexes are known to contribute to the formation of interstellar media and nucleation of aerosols in diverse planetary atmospheres, such as the CO‑rich interstellar comet 2I/Borisov69 and SO2‑rich atmosphere of Venus70.

Aside from the role in ice-grain photochemistry, the planetary atmospheric chemistry of sulfur-containing species also attracts enormous interest due to importance in astrochemical reactions and planetary geology. Particularly, the photochemistry of SO2 plays a fundamental role in the sulfur cycle in the Venusian atmosphere. In addition to the primary contribution from the volcanic eruptions71, photolysis of H2SO4 vapor is another source of SO2 as the planet is completely enshrouded by the acid droplets clouds. As demonstrated by previous field measurements72 and modeling34, photolysis of H2SO4 vapor yields SO3 and H2O, followed by further photodecomposition of SO3 to SO2 and SO. This simple photochemistry can enhance the abundances of SO2 (66 ± 5 ppb) and SO (31 ± 4 ppb) in the cold (ca. –80 °C) mesosphere of Venus at altitudes of 85–100 km73, where the CO abundances (ca. 30 ppm) vary with altitude partially due to reactions with sulfur compounds (e.g., SO2, SO, and S2)74. Hence, the disclosed photochemical reactions between sulfur oxides and carbon oxides in the presence of H2O (an average abundance of 30 ppm at 30–45 km altitude)75 may affect the composition of Venusian atmosphere. Recently, Limaye et al.76 proposed that the lower cloud layer of Venus (50–70 km) is an important target for study, since biorelevant organic acids might be generated through the iron-catalyzed metabolic redox reactions of the abundant CO2, CO, H2O, and SO2 under favorable chemical and physical conditions. The observed photoreactions of sulfur-containing species with CO2 and CO at low temperatures suggests that the hydrogenation of carbon oxides to the organic radicals HCO• and HOCO• for the formation of organic acids might be facilitated by HAT reactions with the derived radicals HOSO• and HSO•.

Methods

Sample preparation

Difluoromethylsulfinic acid (CHF2S(O)OH) was synthesized by reaction of hydrogen chloride (HCl) with sodium difluromethylsulfinate (CHF2S(O)ONa)44. Specifically, freshly dried HCl (10 mmol) was condensed into a glass vessel containing solid CHF2S(O)ONa (0.14 g, 1 mmol) at –196 °C (liquid nitrogen bath), and the mixture was warmed to −110 °C (cold ethanol bath) and kept for overnight reaction. Then, the reaction mixture was slowly warmed to ca. −60 °C and the volatile part was pumped (10 pa) through the vacuum line consisting of two successive cold U-traps at −60 (cold ethanol bath) and −196 °C, and the acid CHF2S(O)OH (ca. 50 mg, 0.5 mmol) was obtained in the first cold trap. The 18O-enriched CHF2S(O)OH was synthesized through hydrolysis of CHF2S(O)Cl with water (18O, 97%, Eurisotop), from which a 1:1:1:1 mixture of CHF2S(O)OH, CHF2S(18O)OH, CHF2S(O)18OH, and CHF2S(18O)18OH (Supplementary Fig. 7) was obtained according to the IR spectrum of its decomposition product HOSO• and also the IR spectrum of its reaction product (SO2) with CO (Fig. 4e).

Matrix‑isolation spectroscopy

Matrix infrared (IR) spectra are recorded using an FT‑IR spectrometer (Bruker 70 V) in a reflectance mode with a transfer optic. A KBr beam splitter and liquid-nitrogen-cooled mercury cadmium telluride (MCT) detector are used in the mid‑IR region (5000‒450 cm−1). For each spectrum, 200 scans at a resolution of 0.5 cm−1 are co‑added. Matrix ultraviolet-visible (UV‒vis) spectra in the range of 190‒800 nm are recorded using a Perkin Elmer Lambda 850+ spectrometer with a scanning speed of 1 nm s−1.

For the preparation of the matrix, gaseous sample is mixed by passing matrix gas (Ar, CO or CO2) through a cold U‑trap (−10 °C) containing ca. 50 mg of the acid precursor (CHF2S(O)OH). Then, the mixture of acid vapor in matrix gas (a ratio of ca. 1:1000) is passed through an aluminum oxide furnace (2.0 mm, i.d.:1.0 mm), which can be heated over a length of ∼25 mm by a tantalum wire (o.d. 0.4 mm, resistance 0.4 Ω) and immediately deposited (2 mmol h−1) in a high vacuum (∼10−5 Pa) onto a gold‑plated copper block matrix support (10 K for Ar‑matrix, 16 K for CO‑matrix) for IR spectroscopy and onto a CaF2 window (16 K) for UV‒vis spectroscopy using closed‑cycle helium cryostats (Sumitomo Heavy Industries, SRDK‑408D2‑F50H) inside the vacuum chambers. For the preparation of the CO2- or CO-doped Ar-matrix, a premix of Ar with CO2 or CO (a ratio of 1:20) was used as the matrix gas. Temperatures at the second stage of the cold heads are controlled and monitored using East Changing TC 290 (IR spectroscopy) and Lake‑Shore 335 (UV‒vis spectroscopy) digital cryogenic temperature controller silicon diodes (DT‑670). The voltage and current use in the pyrolysis experiments are 7 V and 2.9 A, respectively. Photolysis experiments were performed using the Nd3+:YAG laser (266 nm, MPL‑F‑266, 10 mW) and UV lamp (365 nm, 24 W).

Quantum chemistry calculation

Molecular geometries of stationary points for the monomers and complexes were first calculated at the B3LYP‑GD3(BJ)/def2-TZVP77 level of theory with the Gaussian 16 software package78. The dispersion correction using the D3 version of Grimme’s dispersion with Becke-Johnson damping, GD3(BJ)79, is necessary to obtain reliable structures for the hydrogen-bonded complexes. Then, the structures were further optimized at the CCSD(T)/aug‑cc‑pV(T + D)Z80 level of theory with MOLPRO (ver. 2019.1) software81. The non‑covalent interactions (NCIs) analyses are carried out at the B3LYP‑GD3(BJ)/def2‑TZVP level of theory with Bader’s quantum atoms‑in‑molecules (QTAIM)82 and Johnson’s NCI analyses83 on basis of the CCSD(T)/aug‑cc‑pV(T + D)Z level of theory optimized structures.