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
-
-
more » « less
Abstract The structures of zinc carbene ZnCH2and zinc carbyne HZnCH, and the conversion transition states between them are optimized at B3LYP/aug‐cc‐pVTZ, MP2/aug‐cc‐pVTZ, and CCSD/aug‐cc‐pVTZ levels of theory. The thermodynamic energies with CCSD(T) method are further extrapolated to basis set limit through a series of basis sets of aug‐cc‐pVXZ (X=D, T, Q, 5). The Zn−C bonding characteristics are interpreted by molecular plots, Laplacian of density plots, the integrated delocalization indices, net atomic charges, and derived atomic hardness. On the one hand, the studies demonstrated the efficiency of DFT method in structure optimizations and the accuracy of CBS method in obtaining thermodynamic energies; On the other hand, the density analysis of CCSD/aug‐cc‐pVDZ density demonstrates that both the sharing interaction and ionic interaction are important in ZnCH2ad HZnCH. The3B1state of ZnCH2is the global minimum and formed in visible light, but its small bond dissociation energy (47.0 kcal/mol) cannot keep the complex intact under UV light (79.4–102.1 kcal/mol). However, the3Σ−state of HZnCH can survive the UV light due to the greater Zn−C dissociation energy (100.7 kcal/mol). The delocalization indices of Zn…C in both3B1of ZnCH2(0.777) and the3Σ−state of HZnCH (0.785) are close to the delocalization index of the single C−C bond of ethane (0.841), i. e. the nomenclature of Zinc carbene and Zinc carbyne is incorrect. The stronger Zn−C bond in the3Σ−state of HZnCH than in the3B1state of ZnCH2can be attributed to the larger charge separation in the former. It was found that the Zn−C bonds in related Zinc organic compounds were also single bonds no matter whether the organic groups are CR, CR2, or CR3. The ionic interactions were discussed in terms of the atomic hardness that were in turn related to ionization energy and electron affinity. The unique combination of covalent and ionic characteristics in the Zn−C bonds of organic Zinc compounds could be the origin of many interesting applications of organic Zinc reagents.
-
more » « less
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 C2v second‐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. -
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
-
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.
