Torsional vibrations of two thiol groups in the HSOSH, DSOSH, and DSOSD molecules: 2D PES study in CBS limit
Graphical abstract
Introduction
The HSOSH molecule belongs to a group of molecules described by the general formula HXYXH, where X is some atom, and Y is some atom or functional group. These molecules are characterized by the presence of two stable conformers, where tunneling is feasible between these structurally equivalent states. It can also be argued that they all belong to the C2V(M) molecular symmetry group [[1], [2], [3], [4], [5], [6]]. The structure and conformational properties of many members of this group are being actively studied. In particular, in several papers, 2D potential energy surfaces (PES) caused by the internal rotation of hydroxyl (OH) and thiol (SH) groups around single bonds were calculated in molecules: HO(CH2)OH [1,7,8], HOOOH [2,9,10], HOSOH [5], and HSSSH [11]. These studies show that the heights of potential barriers separating the conformers, their stability, and the difference in energies of the ground states of the trans- and cis- conformers differ significantly. Molecules with mixed symmetric combinations of O and S atoms (HOSOH and HSOSH), which also belong to this group, have received much less attention in the literature. Although, it is known that these molecules can be formed during volcanic eruptions [12], as well as in interstellar space [[13], [14], [15]]. In paper [16] the geometric parameters of the equilibrium configurations of the trans- and cis- conformers of the HSOSH molecule at the MP2/6-311G** level of theory were calculated. It was found that the energies of their ground states differ insignificantly. Based on the results in [16], the trans- conformer is energetically moderately preferable than the cis- conformer. At the same time, the barriers to internal rotation of thiol groups in the HSOSH molecule have yet to be calculated. The energies of the torsional states and frequencies of tunneling between the structurally equivalent configurations of the conformers of this molecule have not been found yet.
In this work, the 2D PES due to the rotation of thiol groups in the HSOSH molecule is calculated for the first time. Besides, 2D surfaces of the kinetic coefficients and components of the dipole moment of this molecule were calculated. The energies of the first thirty stationary torsional states were found and determined their symmetry ыз, which further enabled us to compute the torsional IR spectra of the trans- and cis- conformers of the molecule at different temperatures. A comparison of the kinetic coefficients, barriers to internal rotation, as well as the frequencies of torsional vibrations and tunneling frequencies of the HSOSH molecule with the comparable characteristics of the molecules of methanediol [1,3] and hydrogen trioxide [2,4], is performed.
Section snippets
Symmetry properties
Fig. 1 shows the equilibrium configurations of the trans- and cis-conformers of the HSOSH molecule.
Considering Fig. 1, it helps to see the trans- conformer of the HSOSH molecule belongs to the C2 point group, and the cis- conformer belongs to the CS point group. However, as shown in detail in [[1], [2], [3], [4]], the fact that the tunneling is feasible between structurally equivalent configurations in both conformers enables us to assign HSOSH and DOSOD molecules to the C2V(M) molecular
Calculation details
The potential energy of the HSOSH molecule and its other physical characteristics, as functions of the angles of internal rotation , were calculated at the nodes of a two-dimensional uniform grid. The values of the torsional coordinates and on the grid were varied from 4° to 356° in steps of 8°. Thus, molecular parameters were calculated at 45 * 45 = 2025 nodes on the coordinate plane . In each node, when calculating the molecular characteristics, values of the torsional
Discussion of the results
The optimization of the structure of the HSOSH molecule was performed at the MP2/cc-pVTZ and MP2/cc-pVQZ levels of the theory. Then the values of the equilibrium geometric parameters of the trans- and cis-conformers of the molecule were extrapolated to the complete basis set (CBS) limit. Table 1, Table 2 represent the values of the geometric parameters, rotational constants, and kinetic parameters of the two conformers of the HSOSH molecule obtained at the MP2/CBS level of theory.
As follows
Conclusions
This paper is a pioneering analysis of torsional vibrations of HSOSH, DSOSH, and DSOSD molecules in the anharmonic approximation and by constructing 2D PES at the MP2/CBS level of theory, caused by rotation of thiol groups. It was found that in the analyzed molecule the energies of the trans- and cis- conformers are quite close, and the energy of the latter in the ground torsional state is only 35.6 cm−1 higher than the corresponding energy of the trans- conformer.
It was shown that the kinetic
Autors statement
G. Pitsevich - Statement of the problem, theoretical justification of the calculations, analysis of the results
A. Malevich – Development of software codes in the Mathematica shell for processing the results of calculations of 2D PES and kinematic coefficients and numerical solution of the Schrödinger equation
V. Zheutok – 2D PES calculations, visualization of results
A. Khrapunova – 2D PES calculations, contributed to the discussion of materials
U. Sapeshka – 2D PES calculations, analysis of the
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper
References (38)
- et al.
The hydroxyl groups internal rotations in a methanediol molecule
J. Mol. Spectr.
(2019) - et al.
The torsional spectrum of the hydrogen troxide molecule
Chem. Phys.
(2020) - et al.
Quantum aspects of torsional vibrations in the HO3H, DO3H and DO3D molecules
Spectrochim. Acta Part A
(2020) Ab initio studies of the structures and energies of methandiol and 1,1-ethenediol
Chem. Phys.
(1986)- et al.
Multipole approximation of distant pair energies in local MP2 calculations
Chem. Phys. Lett.
(1998) - et al.
Basis-set convergence in correlated calculations on Ne, N2, and H2O
Chem. Phys. Lett.
(1998) - et al.
Exact quantum mechanical kinetic energy operator in valence coordinates for internal motions of a polyatomic molecule
Chem. Phys. Lett.
(1999) - et al.
Hydrogen bonded pyridine N-oxide/trichloroacetic acid complex in polar media: 2D potential energy surface and O–H…..O vibration analysis using exact vibrational Hamiltonian
J. Mol. Struct.
(2014) - et al.
Vibration-rotation intensities in h2o-type molecules application to the 2ν2, ν1 and ν3 bands of H2O
J. Mol. Spectr.
(1975) - et al.
Temperature dependence of the intensity of the vibration-rotational absorption band ν2 of H2O trapped in an argon matrix
Spectrochim. Acta, Part A
(2017)
Barriers to internal rotation and tunnelling splittings of the torsional states in the HO(CH2)OH, DO(CH2)OH and DO(CH2)OD molecules
Mol.Phys.
Torsional states and tunneling probability in HOSOH, DOSOD and DOSOH molecules analyzed at the CBS limit
J. Phys. Chem. A
Molecular Symmetry and Spectroscopy
The Description of molecular structure
Proc. Nat. Acad. Sci. USA
Theoretical determination of molecular structure and conformation. II. Hydrogen trioxide—a model compound for studying the conformational modes of geminal double rotors and five membered rings
J.Chem.Phys.
In search of the next Holy Grail of polyoxide chemistry: Explicitly correlated ab initio full quartic force fields for HOOH, HOOOH, HOOOOH, and their isotopologues
J. Chem. Phys.
Tunneling and parity violation in trisulfane (HSSSH): an almost ideal molecule for detecting parity violation in chiral molecules
Chem. Phys. Chem.
Hydroxysulfinyl radical and sulfinic acid are stable species in the gas phase
J. Am. Chem. Soc.
Sulfur depeletion in dense clouds and circumstellar regions. I. H2S ice abundance and UV-photochemical reactions in the H2O-matrix
Astron. Astrophys.
Cited by (5)
Potential and kinetic interaction of two internal tops in molecules belonging to C<inf>2v</inf>(M) molecular symmetry group
2024, Journal of Molecular SpectroscopyModelling of the torsional IR spectra of the HSSSH, DSSSH, and DSSSD molecules
2023, Computational and Theoretical ChemistryTorsional motions of the free and H-bonded hydroxyl groups of the catechol molecule
2022, Journal of Molecular SpectroscopyCitation Excerpt :We also note that very close results to the data obtained at the MP2/acc-pVTZ level of theory were obtained by the authors of [20], where the calculated frequency of the fundamental torsional vibration of the free hydroxyl group almost coincides with the experimental value, and the calculated frequency of the fundamental torsional vibration of the donor hydroxyl group is 8 cm−1 higher than the experimental value. Note that, as in previous works [26–29], the values of the frequencies of the torsional vibrations of the OH groups, calculated in the harmonic approximation, are overestimated, while the values of these frequencies, calculated in the anharmonic approximation, are significantly underestimated compared to the experimental data (see Table 3). As noted above the authors of [21] were able to determine experimentally the splitting of the first excited torsional state of CTL, which is associated with the excitation of the fundamental torsional vibration of the free OH group.
Torsional IR spectra of three conformers of the resorcinol molecule
2023, Molecular PhysicsExplicitly correlated study of the torsional vibrations of HSOSH molecule. Comparison with MP2/CBS(T,Q) level of theory
2022, Molecular Crystals and Liquid Crystals