skip to main content


Title: Conformational comparison of urea and thiourea near the CCSD ( T ) complete basis set limit
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
Award ID(s):
2154403
NSF-PAR ID:
10443357
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal of Quantum Chemistry
Volume:
123
Issue:
8
ISSN:
0020-7608
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    Thes‐homodesmotic method for computing conventional strain energies (CSE) has been extended for the first time to bicyclic systems and to individual rings within these systems. Unique isodesmic, homodesmotic, and hyperhomodesmotic reactions originate from thes‐homodesmotic method. These are used to investigate 12 bicyclic systems comprising cyclopropane and cyclobutane and how the CSE of each system compares to the sum of the individual rings within each. Equilibrium geometries, harmonic vibrational frequencies, and the corresponding electronic energies and zero point vibrational energy corrections are computed for all relevant molecules using second‐order perturbation theory and density functional theory (B3LYP) with the correlation consistent basis sets cc‐pVDZ and cc‐pVTZ. Single‐point CCSD(T) energies are computed at the MP2/cc‐pVTZ optimized geometries to ascertain the importance of higher order correlation effects. Results indicate that CSEs are additive when the two rings are separated by one or two bonds and somewhat additive in other cases.

     
    more » « less
  4. 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
  5. Bond dissociation energies (BDE) are key descriptors for molecules and are among the most sought-after properties in chemistry. Despite their importance, the accurate prediction of BDE’s for transition metal species can be particularly daunting for both experiment and computation. Experimental data has been limited and, when available, often has large error bars, making the critical evaluation and identification of suitable computational methods difficult. However, recent advancements in the experimental determination of BDE’s with techniques such as Velocity Map Imaging and 2 Photon Ionization now provide useful gauges for computational strategies and new methodologies, providing energies with unprecedented accuracies. The vanadium diatomics (VX, X=B, C, N, O, F, Al, Si, P, S, Cl) have been challenging for computational chemistry methods, and, thus, a new experimental gauge enables methods to be reevaluated and developed for these species. Herein, the super-correlation consistent Composite (super-ccCA or s-ccCA), a new thermochemical scheme centered around CCSD(T)/complete basis set (CBS) limit computations with additional contributions that account for scalar-relativistic effects, and coupled cluster contributions beyond CCSD(T) up to quintuple excitations has been considered. The agreement between determinations made by the s-ccCA scheme and by recent experiment is excellent, demonstrating the utility of the new approach in addressing challenging metal systems, even those of multireference nature. In light of recent experimental BDE’s, the longstanding correlation consistent composite approach (ccCA) is also evaluated for the VX species and find that the mean absolute deviation (MAD) is greatly reduced compared to previously used experimental values. 
    more » « less