An official website of the United States government
Despite the interest in sulfur monoxide (SO) among astrochemists, spectroscopists, inorganic chemists, and organic chemists, its interaction with water remains largely unexplored. We report the first high level theoretical geometries for the two minimum energy complexes formed by sulfur monoxide and water, and we report energies using basis sets as large as aug-cc-pV(Q+d)Z and correlation effects through perturbative quadruple excitations. One structure of SO⋯H 2 O is hydrogen bonded and the other chalcogen bonded. The hydrogen bonded complex has an electronic energy of −2.71 kcal mol −1 and a zero kelvin enthalpy of −1.67 kcal mol −1 , while the chalcogen bonded complex has an electronic energy of −2.64 kcal mol −1 and a zero kelvin enthalpy of −2.00 kcal mol −1 . We also report the transition state between the two structures, which lies below the SO⋯H 2 O dissociation limit, with an electronic energy of −1.26 kcal mol −1 and an enthalpy of −0.81 kcal mol −1 . These features are much sharper than for the isovalent complex of O 2 and H 2 O, which only possesses one weakly bound minimum, so we further analyze the structures with open-shell SAPT0. We find that the interactions between O 2 and H 2 O are uniformly weak, but the SO⋯H 2 O complex surface is governed by the superior polarity and polarizability of SO, as well as the diffuse electron density provided by sulfur's extra valence shell.
more »
« less
Examination of the structures, energetics, and vibrational frequencies of small sulfur‐containing prototypical dimers, (H 2 S) 2 and H 2 O/H 2 S
The optimized geometries, vibrational frequencies, and dissociation energies from MP2 and CCSD(T) computations with large correlation consistent basis sets are reported for (H2S)2and H2O/H2S. Anharmonic vibrational frequencies have also been computed with second‐order vibrational perturbation theory (VPT2). As such, the fundamental frequencies, overtones, and combination bands reported in this study should also provide a useful road map for future spectroscopic studies of the simple but important heterogeneous H2O/H2S dimer in which the hydrogen bond donor and acceptor can interchange, leading to two unique minima with very similar energies. Near the CCSD(T) complete basis set limit, the HOH⋯SH2configuration (H2O donor) lies only 0.2 kcal mol−1below the HSH⋯OH2structure (H2S donor). When the zero‐point vibrational energy is included, however, the latter configuration becomes slightly lower in energy than the former by <0.1 kcal mol−1. © 2018 Wiley Periodicals, Inc.
more »
« less
- Award ID(s):
- 1664998
- PAR ID:
- 10078167
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 40
- Issue:
- 1
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- p. 229-236
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The global minima of urea and thiourea were characterized along with other low‐lying stationary points. Each structure was optimized with the CCSD(T) method and triple‐ζcorrelation consistent basis sets followed by harmonic vibrational frequency computations. Relative energies evaluated near the complete basis set limit with both canonical and explicitly correlated CCSD(T) techniques reveal several subtle but important details about both systems. These computations resolve a discrepancy by demonstrating that the electronic energy of the C2vsecond‐order saddle point of urea lies at least 1.5 kcal mol−1above the C2global minimum regardless of whether the structures were optimized with MP2, CCSD, or CCSD(T). Additionally, urea effectively has one minimum instead of two because the electronic barrier for inversion at one amino group in the Cslocal minimum vanishes at the CCSD(T) CBS limit. Characterization of both systems with the same ab initio methods and large basis sets conclusively establishes that the electronic barriers to inversion at one or both NH2groups in thiourea are appreciably smaller than in urea. CCSDT(Q)/cc‐pVTZ computations show higher‐order electron correlation effects have little impact on the relative energies and are consistently offset by core correlation effects of opposite sign and comparable magnitude.more » « less
-
The spectroscopic identification of Bi 4 has been very elusive. Two constitutional Bi 4 isomers of T d and C 2v symmetry are investigated and each is found to be a local energetic minimum. The optimized geometries and vibrational frequencies of these two isomers are obtained at the CCSD(T)/cc-pVQZ-PP level of theory, utilizing the Stoll, Metz, and Dolg 60-electron effective core potential. The fundamental frequencies of the T d isomer are obtained at the same level of theory. The focal point analysis method, from a maximum basis set of cc-pV5Z-PP, and proceeding to a maximum correlation method of CCSDTQ, was employed to determine the dissociation energy of Bi 4 ( T d ) into two Bi 2 and the adiabatic energy difference between the C 2v and T d isomers of Bi 4 . These quantities are predicted to be +65 kcal mol −1 and +39 kcal mol −1 , respectively. Two electron vertical excitation energies between the T d and C 2v electronic configurations are computed to be 156 kcal mol −1 for the T d isomer and 9 kcal mol −1 for the C 2v isomer. The most probable approach to laboratory spectroscopic identification of Bi 4 is via an infrared spectrum. The predicted fundamentals (cm −1 ) with harmonic IR intensities in parentheses (km mol −1 ) are 94(0), 123(0.23), and 167(0) for the T d isomer. The moderate IR intensity for the only allowed fundamental may explain why Bi 4 has yet to be observed. Through natural bond orbital analysis, the C 2v isomer of Bi 4 was discovered to exhibit “long-bonding” between the furthest apart ‘wing’ atoms. This long-bonding is postulated to be facilitated by the σ-bonding orbital between the ‘spine’ atoms of the C 2v isomer.more » « less
-
This study systematically characterizes the four homogeneous and six heterogeneous hydrogen-bonded dimers formed by hydrogen halide pairs (HX/HY where X, Y = F, Cl, Br, and I). The notation HX⋯HY indicates the direction of the hydrogen bond from the HY donor to the HX acceptor. All stationary points reported for these ten dimer systems are fully optimized utilizing the MP2 and CCSD(T) ab initio methods in conjunction with quadruple-ζ correlation-consistent basis sets augmented with diffuse functions, and their nature is verified by harmonic vibrational frequency computations. The electronic dissociation energies (De) for all ten global minima are evaluated near the CCSD(T) complete basis set (CBS) limit via extrapolation schemes. These values are 19.11, 8.32, 7.38, and 6.22 kJ mol−1 for the homogeneous dimers of HF, HCl, HBr, and HI, respectively. For the heterogeneous pairs, the lighter hydrogen halide is consistently the donor in the global minimum configuration, with De ranging from 12.23 kJ mol−1 for HCl⋯HF to 7.22 kJ mol−1 for HI⋯HBr near the CCSD(T) CBS limit. Interestingly, not all heterodimer donor/acceptor permutations correspond to minima. For example, the HCl⋯HBr configuration is identified as a local minimum at all levels of theory employed in this investigation, whereas the in-plane barrier for donor/acceptor exchange vanishes for HCl⋯HI and HBr⋯HI when larger quadruple-ζ basis sets are utilized. For the seven dimer systems containing Br and/or I, the structures, energetics, and vibrational frequencies computed using conventional valence-only electron correlation procedures are similar to those obtained using an expanded valence treatment that includes the (n − 1)d subvalence electrons associated with Br and I.more » « less
-
Abstract We have performed a series of highly accurate calculations between CO2and the 20 naturally occurring amino acids for the investigation of the attractive noncovalent interactions. Different nucleophilic groups present in the amino acid structures were considered (α‐NH2, COOH, side groups), and the stronger binding sites were identified. A database of accurate reference interactions energies was compiled as computed by explicitly‐correlated coupled‐cluster singles‐and‐doubles, together with perturbative triples extrapolated to the complete‐basis‐set limit. The CCSD(F12)(T)/CBS reference values were used for comparing a variety of popular density functionals with different basis sets. Our results show that most density functionals with the triple‐zeta basis set def2‐TZVPP align with the CCSD(F12)(T)/CBS reference values, but errors range from 0.1 kcal/mol up to 1.0 kcal/mol.more » « less
