# Reaction of F Atoms with Carbonyl Sulfide:
Experimental and Theoretical Study

Yuri Bedjanian, Manolis Romanias

<figure><img src="image_1.png" /></figure>

► To cite this version:

Yuri Bedjanian, Manolis Romanias. Reaction of F Atoms with Carbonyl Sulfide: Experimental and Theoretical Study. Molecules, 2024, 29 (21), pp.5027. 10.3390/molecules29215027. hal-04759621

HAL Id: hal-04759621

https://hal.science/hal-04759621v1

Submitted on 19 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# Reaction of F Atoms with Carbonyl Sulfide: Experimental and Theoretical Study

Yuri Bedjanian<fnref n="1" />,* and Manolis N. Romanias<fnref n="2" />

**Abstract:** The gas phase reaction of fluorine atoms with carbonyl sulfide, F + OCS → SF + CO (1), has been investigated over a wide temperature range, T = 220–960 K, in a discharge-flow system combined with mass spectrometry for the analysis of the reactive mixture. The reaction rate coefficient was determined as a function of temperature using both absolute (under pseudo-first-order conditions using an excess of F atoms) and relative rate methods: k₁ = (8.95 ± 0.52) × 10⁻¹¹ exp((102 ± 20)/T) cm³ molecule⁻<fnref n="1" /> s⁻<fnref n="1" /> (with an estimated total uncertainty of 15% independent of temperature). The SF radical was observed as the main primary reaction product. The experimental results are supported by theoretical calculations, which evidence that the reaction is exothermic with a weak negative temperature dependence.

**Keywords:** fluorine atom; carbonyl sulfide; rate coefficient; temperature dependence; negative activation energy

## 1. Introduction

Reactions of atomic fluorine with closed-shell molecules are generally very rapid and lead to the generation of diverse active species. This specificity of the reactions of fluorine atoms determines, in addition to theoretical interest, practical interest in these reactions associated with laboratory studies of gas-phase kinetics, where the reactions of F atoms can potentially be used as a source of various atoms and radicals, as well as titration reactions for detecting and measuring the absolute concentrations of the active substances formed. The existing database on the reactions of F atoms is very limited and uncertain: the rate constants of many studied reactions of fluorine atoms are determined with insufficient precision even at room temperature, let alone the temperature dependence [1]. The main reason for this is the rapidity of the reactions (rate constants are often close to the frequency of bimolecular collisions) and the presence of highly reactive species in the reaction products, which imposes certain requirements on detection sensitivity to limit the potential impact of rapid secondary chemistry.

In the present work, we report the results of the combined experimental and theoretical study of the reaction of F atoms with carbonyl sulfide.

$$F + OCS \rightarrow \text{products} \quad (R1)$$

Reaction (R1) is a potential source of SF radicals. It can be also of potential interest for the chemistry of hot near-source volcanic plumes, where OCS and fluorinated compounds are simultaneously emitted [2]. Note that for reaction (1) only one previous study at room temperature is available [3].

The aim of this work was to experimentally determine the rate constant of the title reaction in an extended temperature range T = 220–960 K and to obtain information on

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Bedjanian, Y.; Romanias, M.N.

**Reaction of F Atoms with Carbonyl Sulfide:** Experimental and Theoretical Study. *Molecules* **2024**, *29*, 5027.

https://doi.org/10.3390/molecules29215027

Academic Editor: Enrico Bodo

Received: 16 September 2024

Revised: 17 October 2024

Accepted: 17 October 2024

Published: 24 October 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Institut de Combustion, Aérothermique, Réactivité et Environnement (ICARE), CNRS, 45071 Orléans, France</footnote>

<footnote n="2">² IMT Nord Europe, Institut Mines-Télécom, Centre for Energy and Environment, University Lille, 59000 Lille, France; emmanouil.romanias@imt-nord-europe.fr</footnote>

<footnote>* Correspondence: yuri-bedjanian@cnrs-orleans.fr</footnote>
the reaction mechanism through direct detection of the reaction products and theoretical thermochemical calculations.

**2. Results and Discussion**

**2.1. Rate Constant of Reaction (R1): Absolute Measurements**

The configuration of the flow reactor used in the absolute measurements of the rate constant of Reaction (R1) is shown in Figure 1.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Diagram of the low temperature flow reactor: configuration used in the absolute measurements of $k_1$.</figcaption></figure>

The measurements of $k_1$ consisted of monitoring the kinetics of the OCS (initial concentration, [OCS]₀, of (1 – 2) × 10¹¹ molecule cm⁻³) consumption under an excess concentration of F atoms ([F] = (0.35 – 6.20) × 10¹² molecule cm⁻³). The advantages of working with excess of F atoms over OCS were as follows: (i) better detection sensitivity for OCS, allowing the use of low initial concentrations of OCS, and (ii) non-involvement of OCS in possible secondary reactions, unlike the fluorine atom. Examples of the exponential decays of OCS, [OCS] = [OCS]₀ × exp(−$k_1'$ × t), where $k_1' = k_1 × [F]$ is the pseudo-first-order rate constant and [F] represents the average F-atom concentration along the reaction zone, are shown in Figure 2.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Plots of OCS concentration against reaction time for various concentrations of fluorine atom (T = 300 K).</figcaption></figure>

The consumption of F atoms due their heterogeneous loss ($k_W \le 10 \text{ s}^{-1}$) and reaction with OCS was insignificant, reaching a maximum of 20% in a few kinetic runs. The diffusion corrections applied to $k_1'$ values, determined from the data like those in Figure 2,
were generally <10% and up to 20% in a few kinetic runs. Examples of second-order plots measured at different temperatures are shown in Figure 3. The rate constants of the Reaction (R1) obtained from a linear through origin least-squares fit of the $k_1'$ data as a function of [F] are given in Table 1. The combined uncertainty of the absolute measurements of the rate constants was estimated to be ≤15%, including (in quadrature) statistical error and those of the measurements of the flows, pressure, temperature and absolute concentration of F atoms.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Pseudo-first-order rate constant $k_1' = k_1 \times [F]$ as a function of F-atom concentration at different temperatures.</figcaption></figure>

<table><caption>**Table 1.** Summary of the measurements of the rate constant of Reaction (R1).</caption><thead><tr><th><i>T</i> (K)</th><th><i>k</i><sub>1</sub><sup>a</sup></th><th>Method<sup>b</sup></th><th>Reactor Surface<sup>c</sup></th></tr></thead><tbody><tr><td>220</td><td>1.45</td><td>AM</td><td>HW</td></tr><tr><td>235</td><td>1.42</td><td>AM</td><td>HW</td></tr><tr><td>255</td><td>1.39</td><td>AM</td><td>HW</td></tr><tr><td>275</td><td>1.24</td><td>AM</td><td>HW</td></tr><tr><td>300</td><td>1.23</td><td>AM</td><td>HW</td></tr><tr><td>310</td><td>1.16</td><td>RRM2</td><td>Q</td></tr><tr><td>320</td><td>1.17</td><td>AM</td><td>Q</td></tr><tr><td>330</td><td>1.25</td><td>AM</td><td>HW</td></tr><tr><td>340</td><td>1.12</td><td>RRM1</td><td>Q</td></tr><tr><td>370</td><td>1.13</td><td>AM</td><td>Q</td></tr><tr><td>375</td><td>1.12</td><td>RRM2</td><td>Q</td></tr><tr><td>420</td><td>1.15</td><td>AM</td><td>Q</td></tr><tr><td>420</td><td>1.08</td><td>RRM1</td><td>Q</td></tr><tr><td>475</td><td>1.09</td><td>RRM2</td><td>Q</td></tr><tr><td>490</td><td>1.06</td><td>AM</td><td>Q</td></tr><tr><td>555</td><td>1.09</td><td>RRM1</td><td>Q</td></tr><tr><td>585</td><td>1.08</td><td>AM</td><td>Q</td></tr><tr><td>650</td><td>0.94</td><td>RRM2</td><td>Q</td></tr><tr><td>730</td><td>1.02</td><td>AM</td><td>Q</td></tr><tr><td>790</td><td>0.97</td><td>RRM1</td><td>Q</td></tr><tr><td>960</td><td>1.05</td><td>AM</td><td>Q</td></tr></tbody></table>

<footnote><sup>a</sup> units of $10^{-10}$ cm$^3$ molecule$^{-1}$ s$^{-1}$, statistical 2σ uncertainty is (2–3)%, total estimated uncertainty on $k_1$ is 15% and 20% for absolute and relative measurements, respectively; <sup>b</sup> AM: absolute measurements; RRM1 and RRM2: two relative rate methods (see text); <sup>c</sup> HW: halocarbon wax; Q: quartz.</footnote>
## 2.2. Rate Constant of Reaction (R1): Relative Rate Measurements

Figure 4 displays the configuration of the flow reactor used in the relative measurements of $k_1$.

<figure><img src="image_8.png" /><figcaption>Figure 4. Diagram of the high temperature flow reactor: configuration used in the relative measurements of $k_1$.</figcaption></figure>

In the first series of relative rate measurements, the rate constant of Reaction (R1) was determined using a relative rate method (RRM1), comparing the decay rates of OCS and reference compound, $Br_2$, simultaneously present in the reactor, in reactions with F atoms (Figure 4).

$$ \ln \frac{[\text{OCS}]_0}{[\text{OCS}]} = \frac{k_1}{k_2} \times \ln \frac{[\text{Br}_2]_0}{[\text{Br}_2]} \quad (1) $$

where the expressions under the logarithm are the ratios of the compound concentration in the absence to that in the presence of F atoms at a given reaction time. Examples of the experimental data are shown in Figure 5.

<figure><img src="image_9.png" /><figcaption>Figure 5. Dependence of $\ln([\text{OCS}]_0/[\text{OCS}])$ on $\ln(\text{Br}_2)_0/[\text{Br}_2]$ observed at T = 340 and 790 K.</figcaption></figure>

The final values of $k_1$ (Table 1) were determined using $k_1/k_2$, derived from the slopes of the straight lines in Figure 5, and the rate constant of the reference reaction.

$$ F + Br_2 \rightarrow Br + FBr \quad (R2) $$

where $k_2 = (1.28 \pm 0.20) \times 10^{-10} \text{ cm}^3 \text{ molecule}^{-1} \text{ s}^{-1}$ (T = 299–940 K) [4].

Although the results of these relative measurements are in good agreement with the absolute data (Table 1), it should nevertheless be noted that the lines in Figure 5 do not
pass through the origin, as one would expect from Equation (1). This seems to indicate additional consumption of Br₂ in the side processes. In search of potential reasons for the observed effect, we conducted additional experiments at T = 420 K by varying the reaction time and initial concentrations of the reagents. The results are presented in Figure 6.

<figure><img src="image_10.png" /><figcaption>Figure 6. Dependence of ln([OCS]₀/[OCS]) on ln(Br₂]₀/[Br₂]) observed at T = 420 K under different experimental conditions.</figcaption></figure>

One can note that variation in the reaction time (by a factor of 2) and initial concentrations of OCS and Br₂ (by a factor of 3 and 2.5, respectively) do not affect either the slope of the straight line (k₁/k₂) or its intercept. This observation seems to indicate that the additional consumption of Br₂ in secondary reactions is unlikely. The observed non-zero intercept in Figures 5 and 6 could potentially be due to the presence of O atoms, which are known to be generated in microwave discharge of F₂. Moreover, oxygen atoms react much faster with Br₂ than with OCS.

$$O + Br_2 \rightarrow Br + BrO \quad (R3)$$

$$k_3 = 9.85 \times 10^{-16} \text{ T}^{1.41} \exp(543/T) \text{ cm}^3 \text{ molecule}^{-1} \text{ s}^{-1} (T = 220-950 \text{ K}) \text{ [5]}$$

$$O + \text{OCS} \rightarrow \text{products} \quad (R4)$$

where k₄ = 1.92 × 10⁻¹² (T/298)².⁰⁸ exp(-1524/T) cm³ molecule⁻¹ s⁻¹ (T = 220-960 K) [6]. In addition, the contribution of small amounts of H atoms and OH radicals from the microwave discharge to the consumption of Br₂ cannot be ruled out, considering that H and OH are more reactive towards Br₂ than toward OCS.

To be more confident in the relative measurement of k₁, we carried out another series of experiments in which the rate constant of Reaction (R1) was measured relative to that of the reaction of F atoms with Br₂, but within the framework of a different approach (RRM2). The experiments consisted of rapid titration of F atoms in the reaction with a mixture of OCS and Br₂ and measuring the yield of FBr as a function of the [OCS]/[Br₂] ratio. The fraction of the initial concentration of F atoms converted into FBr in Reaction (R2) is equal to the following:

$$[\text{FBr}] = \frac{k_2[\text{Br}_2]}{k_2[\text{Br}_2] + k_1[\text{OCS}] + k_w} \times [\text{F}]_0$$
Rearranging this equation gives the following equation:

$$ \frac{[F]_0}{[FBr]} - 1 = \frac{k_1[OCS]}{k_2[Br_2]} + \frac{k_w}{k_2[Br_2]} \quad (2) $$

At a constant Br₂ concentration, the second term in Equation (2) is constant and k₁/k₂ can be determined as the slope of the linear relationship ([F]₀/[FBr] - 1) versus [OCS]/[Br₂] ratio. The concentration of FBr was monitored in the presence of Br₂ and in the absence of OCS, which corresponds to [FBr]₀ = [F]₀, and in the presence of OCS in the reactor, which corresponds to the fraction of [F]₀ that reacted with Br₂. Examples of the observed experimental data are shown in Figure 7 and final values of k₁ from these experiments are given in Table 1.

<figure><img src="image_11.png" /><figcaption>Figure 7. Yield of FBr upon titration of F atoms in reaction with OCS/Br₂ mixtures at different temperatures.</figcaption></figure>

### 2.3. Rate Constant of Reaction (R1): Temperature Dependence

All the data available for the rate constant of Reaction (R1) are displayed in Figure 8. Concerning the data from the present study, we can note a very good agreement between the relative rate method, which does not require measurements of the absolute concentrations of the species, and the absolute measurements of the reaction rate constant. Negative temperature dependence observed for k₁ can be described by the following Arrhenius equation:

$$ k_1 = (8.95 \pm 0.52) \times 10^{-11} \exp((102 \pm 20)/T) \text{ cm}^3 \text{ molecule}^{-1} \text{ s}^{-1} $$

in the temperature range from 220 to 960 K and with 2σ uncertainties representing the precision of the fit.

Previously, Reaction (R1) has been investigated in only one study [3]. Brunning and Clyne [3], using a discharge-flow system combined with mass spectrometric detection of the species involved, reported k₁ = (1.81 ± 0.16) × 10⁻¹¹ cm³ molecule⁻¹ s⁻¹ at T = 295 K and 1 Torr total pressure. This value is lower by a factor of 7 than that determined in the present work. It is difficult for us to explain such a large discrepancy, especially since both studies were carried out in similar flow reactors and the reaction rate constant was determined using the similar experimental protocol, from the kinetics of OCS consumption in excess of F atoms. It seems that the main reason for the disagreement lies in the determination of the
absolute concentrations of F atoms. In the study of Brunning and Clyne [3], F-atom concentrations were determined using the titration reaction $F + Cl_2 \rightarrow Cl + FCl$. Unfortunately, the article does not provide details of the measurements and does not specify whether only the initial concentration of F atoms was measured or their kinetics along the reaction zone was controlled. In the present work, a similar procedure, the titration of F atoms with excess $Br_2$, was applied for measurements of their absolute concentrations and [F] in the reactor was monitored in situ as a function of the reaction time. In addition, relative measurements of the rate constant, independent of the absolute concentration measurements, support the absolute data for $k_1$.

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Summary of the measurements of the rate constant of Reaction (R1). Error bars on the present data correspond to 15% and 20% uncertainty of the absolute and relative measurements of $k_1$, respectively. Red circle: Brunning and Clyne [3].</figcaption></figure>

## 2.4. Theoretical Findings

There are three different pathways for the title reaction, which are shown below:

$$OCS + F \rightarrow SF + CO \qquad (R1a)$$

$$\rightarrow FO + CS \qquad (R1b)$$

$$\rightarrow OC(F)S (+M) \qquad (R1c)$$

Reaction (R1b) was calculated to be highly endothermic with energy of 461 kJ mol⁻¹ using the G3//MP2/aug-cc-pVDZ. It is highly improbable that this channel will occur and thus it is not further discussed in the following section of the manuscript. On the contrary, pathways (R1a) and (R1c) are exothermic. The optimized geometries of reactants, products, TS, and intermediate species formed along the reaction coordinate with the G3B3 method are shown in Figure 9. The relative energies with the corresponding zero-point energies (ZPEs) for each conformation with the two different computational methods are given in Table 2. As displayed in Table 2, there is a relatively good agreement between the relative energies calculated with the two different methods. The reaction of F with OCS initially proceeds with the formation of a stable adduct (named pre-adduct in the following) where F atoms approach the OCS from the S site. The conformation of this adduct locked with G3B3 is displayed in Figure 9.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Profile of the potential energy surface for the reaction of F atoms with OCS. Structure optimizations and single-point energy calculations (in kJ mol⁻¹), carried out using the quantum composite G3B3 method, are also presented. The red and black lines correspond to the S abstraction (R1a) and F-atom addition (R1c) pathways, respectively. The bond lengths are in Å.</figcaption></figure>

<table><caption>**Table 2.** Calculated energies (kJ mol⁻¹) for the title reaction using the G3B3 and G3//MP2/aug-cc-pVDZ methods. Spin coupling corrections for F and SF were applied.</caption><thead><tr><th>Method</th><th>Pre-Adduct</th><th>TS</th><th>Post-Adduct</th><th>Final Products</th></tr></thead><tbody><tr><td>G3B3</td><td>-71.3</td><td>OCS +F → SF + CO (R1a)</td><td>-42.2</td><td>-32.2</td></tr><tr><td>G3//MP2/aug-cc-pVDZ</td><td>-67.7</td><td>-6.9</td><td>-2.0</td><td>-27.8</td></tr><tr><td>G3B3</td><td>-71.3</td><td>OCS +F → OC(F)S (R1c)</td><td>-33.3</td><td>-149.6</td></tr><tr><td>G3//MP2/aug-cc-pVDZ</td><td>-67.7</td><td>-3.6</td><td>-36.4</td><td>-146.6</td></tr></tbody></table>

The relative energies calculated for the pre-adduct with G3B3 and G3//MP2/aug-cc-pVDZ were -71.3 kJ mol⁻¹ and -67.7 kJ mol⁻¹, respectively, corresponding to a relative stable conformation. We locked two different TS conformations corresponding to the abstraction (R1a) and addition (R1c) pathways, respectively. In both cases, submerged reaction barriers were calculated, indicating a possible inverse temperature dependence of the corresponding rate coefficient, in agreement with the experimental observations. Concerning the (R1a), channel the calculated energies were -6.9 kJ mol⁻¹ and -2.0 kJ mol⁻¹ with G3B3 and G3//MP2/aug-cc-pVDZ, respectively. The relative energy of the TS calculated for (R1c) was significantly lower, in the range of -33 kJ mol⁻¹. A post-adduct was locked for the (R1a) channel with the G3B3 theory (but not with G3//MP2/aug-cc-pVDZ) with a calculated energy of -42.2 kJ mol⁻¹ before the final products formation (i.e., CO and SF). Regarding channel (R1c), the final product was the OC(F)S radical, with a calculated energy of -149.6 and -146.6 kJ mol⁻¹ with G3B3 and G3//MP2/aug-cc-pVDZ, respectively. It should be noted that a small basis set was used for the ab-initio theoretical calculations. While the use of more advanced basis sets would yield more accurate energy estimations, the objective of these calculations was to support the experimental data and provide a relative comparison of the two reaction pathways.

The intermediate energized adduct can dissociate back into reactants or to reaction products (Reaction (R1a)), or be collisionally stabilized (Reaction (R1c)). For the addition
channel to take place (adduct stabilization), relatively high pressure is usually required to remove excess energy of the excited adduct via collisions with a third body (the bath gas in our experiments). It is clear that the competition between pathways (R1a) and (R1c) will be pressure and temperature dependent with a possible significant contribution of addition channel (R1c) at high pressures. However, under the experimental conditions of our measurements, at a total pressure of 2 Torr, a noticeable contribution to $k_1$ (~$10^{-10}$ cm³ molecule⁻¹ s⁻¹) of the addition channel seems very unlikely. For this reason, we did not consider it necessary to measure the reaction rate constant as a function of pressure. Using mass spectrometry, we detected the formation of the SF radical in Reaction (R1) directly at *m/z* = 51 (SF⁺) and upon its conversion to BrSF (*m/z* = 130/132, BrSF⁺) in reaction with Br₂. Although the absolute concentrations of SF were not measured, judging by the relative mass spectrometric signals, the SF formation channel (1a) is the main, if not the only, reaction channel under the experimental conditions of the study. Reaction (R1) turns out to be a convenient and effective source of SF radicals, which opens up prospects and opportunities for detailed laboratory studies of elementary reactions involving SF.

# 3. Materials and Methods

## 3.1. *Experimental Setup*

The experimental setup, consisting of a discharge-flow reactor combined with a mod-
ulated molecular beam mass spectrometer with electron impact ionization, has been used
extensively in the past [4–8]. Therefore, only a brief description will be presented herein. It
contains the gas introduction vacuum lines, the flow reactor, and the differentially pumped
stainless steel high-vacuum chamber, which houses the quadrupole mass spectrometer
(Balzers, QMG 420, Balzers Aktiengesellschaft, Liechtenstein). The chemical composition
of the reactive system was monitored by sampling the gas-phase molecules from the flow
reactor and detecting them using the mass spectrometer. The gas mixture sampled from the
flow reactor was modulated by a tuning-fork chopper (35 Hz), ionized through the impact
with high kinetic energy electrons (30 eV) emitted by the ion source of the mass spectrome-
ter and detected using an electron multiplier. Subsequently, the mass spectrometric signals
were filtered and amplified with a lock-in amplifier and recorded for further analysis.

The reaction time in the flow reactor was determined by the position of the tip of the movable injector relative to the sampling cone of the mass spectrometer (Figures 1 and 4) and by the linear flow velocity in the reactor, which was in the range (1720–2460) cm s⁻¹.

All experiments were carried out at a total pressure of 2 Torr, using He as the carrier gas. To provide an extended temperature range for kinetic measurements, two flow reactors were used. The first one used at low temperatures (220–330 K) consisted of a Pyrex tube (45 cm length and 2.4 cm i.d.); temperature regulation was achieved by circulating thermostated ethanol (Figure 1). The second flow reactor (T = 310–960 K) consisted of an electrically heated quartz tube (45 cm length and 2.5 cm i.d.) with water-cooled attachments (Figure 4) [8].

Fluorine atoms were generated from the microwave discharge in trace amounts of F₂ in He in a ceramic (Al₂O₃) tube (Figures 1 and 4), which was used to avoid potential reactions of F atoms with the glass surface inside the microwave cavity. F atoms were detected at *m/z* = 98/100 (FBr⁺) after being converted to FBr prior to their sampling by Reaction (R2) with excess Br₂ (Figure 1).

The absolute concentrations of F atoms were determined through their titration with excess Br₂. The measurement procedure consisted of the chemical conversion of F atoms into FBr upon their complete consumption in Reaction (R2), which relates the initial concentration of F atoms to the consumed fraction of [Br₂] and the formed concentration of FBr: [F]₀ = [FBr] = Δ[Br₂]. The absolute concentrations of OCS and Br₂ in the reactor were calculated from their flow rates obtained from pressure drop measurements in their mixtures in He stored in calibrated volumes.

The purities of the gases used were as follows: He > 99.9995% (Alphagaz, Air Liquide France Industrie, Paris, France), passed through liquid nitrogen trap; Br₂ > 99.99% (Aldrich,
St. Louis, MO, USA); F₂, 5% in helium (Alphagaz), OCS > 99.99% (15% in He, Messer, Bad Soden, Germany).

## 3.2. Computational Methodology

All electronic structure calculations were carried out using the Gaussian 16 program suite [9]. The calculations were performed using two different methods. The first one concerned the quantum composite G3B3 method [10], which is a variation of G3 theory [11]. This compound method consists of geometry optimization and vibrational frequency calculations at the B3LYP/6-31G(d) [12] level of theory. Thereafter, a series of single-point energy calculations at four different higher levels of theory QCISD(T)/6-31G(d), MP4/6-31+G(d), MP4/6-31G(2df,p), and MP2/G3Large were performed using the optimized structure. Calculations were also carried out with the G3//MP2/aug-cc-pVDZ method. More precisely, all geometry optimization and vibrational frequency calculations were conducted at the MP2 [13]/aug-cc-pVDZ [14] level of theory. Subsequently, single-point energy calculations were realized for the optimized structures, manually applying the G3 theory and scaling the harmonic frequencies by a factor of 0.959 as adopted from the CCCBDB database [15]. Further improvements in electronic energies were achieved by applying spin-orbit coupling corrections. In addition, the potential energy surfaces were computed by performing intrinsic reaction coordinate (IRC) calculations in order to verify that the located transition state (TS) structures are connected to the proper reactants and products.

## 4. Conclusions

In this study, a discharge-flow reactor combined with an electron impact ionization mass spectrometer was employed to study the reaction of fluorine atoms with carbonyl sulfide. The reaction rate constant was determined at a total pressure of 2 Torr as a function of temperature in a wide temperature range (T = 220–960 K) using an absolute and two different relative rate methods. The rate constant data obtained by the three methods are in good agreement with each other and provide the following Arrhenius expression: $k_1 = (8.95 \pm 0.52) \times 10^{-11} \exp((102 \pm 20)/T) \text{ cm}^3 \text{ molecule}^{-1} \text{ s}^{-1}$ (with an estimated total uncertainty of 15% independent of temperature), with a $k_1$ value that is a factor of 7 higher compared to the only previous measurement at room temperature. The theoretical calculations performed within this study support the negative sign of the activation energy. The SF radical was observed as the reaction product. The experimental observations suggest that the SF-forming channel, F + OCS → SF + CO, is the main, if not the only, reaction channel under the experimental conditions of the study.

**Author Contributions:** Conceptualization, experimental investigation, data analysis and interpretation, writing, Y.B.; computation, methodology, writing, M.N.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partly funded by ANR through the PIA (Programme d'Investissement d'Avenir, grant number ANR-10-LABX-100-01).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data supporting reported results are available in this article.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Manion, J.A.; Huie, R.E.; Levin, R.D.; Burgess, D.R.; Orkin, V.L.; Tsang, W.; McGivern, W.S.; Hudgens, J.W.; Knyazev, V.D.; Atkinson, D.B.; et al. NIST Chemical Kinetics Database, NIST Standard Reference Database 17, Version 7.0 (Web Version), Release 1.6.8, Data version 2015.12, National Institute of Standards and Technology, Gaithersburg, Maryland, 20899–28320. Available online: http://kinetics.nist.gov/ (accessed on 13 September 2024).

2. Halmer, M.M.; Schmincke, H.U.; Graf, H.F. The annual volcanic gas input into the atmosphere, in particular into the stratosphere: A global data set for the past 100 years. *J. Volcanol. Geotherm. Res.* **2002**, *115*, 511–528. [CrossRef]
3. Brunning, J.; Clyne, M.A.A. Elementary reactions of the SF radical. Part 1.—Rate constants for the reactions F + OCS → SF + CO and SF + SF → SF₂ + S. *J. Chem. Soc. Faraday Trans.* **2** ***1984*, *80*, 1001–1014. [CrossRef]

4. Bedjanian, Y. Kinetics and Products of the Reactions of Fluorine Atoms with CINO and Br₂ from 295 to 950 K. *J. Phys. Chem. A* **2017**, *121*, 8341–8347. [CrossRef] [PubMed]

5. Bedjanian, Y.; Kalyan, C. Rate constants of the reactions of O(³P) atoms with Br₂ and NO₂ over the temperature range 220–950 K. *Int. J. Chem. Kinet.* **2019**, *51*, 476–483. [CrossRef]

6. Bedjanian, Y. Experimental Study of the Reaction of O(³P) with Carbonyl Sulfide between 220 and 960 K. *J. Phys. Chem. A* **2022**, *126*, 4080–4086. [CrossRef] [PubMed]

7. Bedjanian, Y. Rate Coefficients of the Reactions of Fluorine Atoms with H₂S and SH over the Temperature Range 220–960 K. *Molecules* **2022**, *27*, 8365. [CrossRef] [PubMed]

8. Morin, J.; Romanias, M.N.; Bedjanian, Y. Experimental study of the reactions of OH radicals with propane, *n*-pentane, and *n*-heptane over a wide temperature range. *Int. J. Chem. Kinet.* **2015**, *47*, 629–637. [CrossRef]

9. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. *Gaussian 16, Revision C.01*; Gaussian, Inc.: Wallingford CT, USA, 2016.

10. Baboul, A.G.; Curtiss, L.A.; Redfern, P.C.; Raghavachari, K. Gaussian-3 theory using density functional geometries and zero-point energies. *J. Chem. Phys.* **1999**, *110*, 7650–7657. [CrossRef]

11. Curtiss, L.A.; Raghavachari, K.; Redfern, P.C.; Rassolov, V.; Pople, J.A. Gaussian-3 (G3) theory for molecules containing first and second-row atoms. *J. Chem. Phys.* **1998**, *109*, 7764–7776. [CrossRef]

12. Francl, M.M.; Pietro, W.J.; Hehre, W.J.; Binkley, J.S.; Gordon, M.S.; DeFrees, D.J.; Pople, J.A. Self-consistent molecular orbital methods. XXIII. A polarization-type basis set for second-row elements. *J. Chem. Phys.* **1982**, *77*, 3654–3665. [CrossRef]

13. Møller, C.; Plesset, M.S. Note on an Approximation Treatment for Many-Electron Systems. *Phys. Rev.* **1934**, *46*, 618–622. [CrossRef]

14. Kendall, R.A.; Dunning, T.H., Jr.; Harrison, R.J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. *J. Chem. Phys.* **1992**, *96*, 6796–6806. [CrossRef]

15. Johnson, R. Computational Chemistry Comparison and Benchmark Database. 2018. Available online: http://cccbdb.nist.gov/ (accessed on 13 September 2024). [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.