More Like this

1. A comparison between hydrogen and halogen bonding: the hypohalous acid–water dimers, HOX⋯H ₂ O (X = F, Cl, Br)

   https://doi.org/10.1039/C9CP00422J  

   Wolf, Mark E.; Zhang, Boyi; Turney, Justin M.; Schaefer, Henry F. (March 2019, Physical Chemistry Chemical Physics)

   Hypohalous acids (HOX) are a class of molecules that play a key role in the atmospheric seasonal depletion of ozone and have the ability to form both hydrogen and halogen bonds. The interactions between the HOX monomers (X = F, Cl, Br) and water have been studied at the CCSD(T)/aug-cc-pVTZ level of theory with the spin free X2C-1e method to account for scalar relativistic effects. Focal point analysis was used to determine CCSDT(Q)/CBS dissociation energies. The anti hydrogen bonded dimers were found with interaction energies of −5.62 kcal mol −1 , −5.56 kcal mol −1 , and −4.97 kcal mol −1 for X = F, Cl, and Br, respectively. The weaker halogen bonded dimers were found to have interaction energies of −1.71 kcal mol −1 and −3.03 kcal mol −1 for X = Cl and Br, respectively. Natural bond orbital analysis and symmetry adapted perturbation theory were used to discern the nature of the halogen and hydrogen bonds and trends due to halogen substitution. The halogen bonds were determined to be weaker than the analogous hydrogen bonds in all cases but close enough in energy to be relevant, significantly more so with increasing halogen size. 

   more » « less
2. A quantitative assessment of deformation energy in intermolecular interactions: How important is it?

   https://doi.org/10.1063/5.0155895  

   Sargent, Caroline T.; Kasera, Raina; Glick, Zachary L.; Sherrill, C. David; Cheney, Daniel L. (June 2023, The Journal of Chemical Physics)

   Dimer interaction energies have been well studied in computational chemistry, but they can offer an incomplete understanding of molecular binding depending on the system. In the current study, we present a dataset of focal-point coupled-cluster interaction and deformation energies (summing to binding energies, De) of 28 organic molecular dimers. We use these highly accurate energies to evaluate ten density functional approximations for their accuracy. The best performing method (with a double-ζ basis set), B97M-D3BJ, is then used to calculate the binding energies of 104 organic dimers, and we analyze the influence of the nature and strength of interaction on deformation energies. Deformation energies can be as large as 50% of the dimer interaction energy, especially when hydrogen bonding is present. In most cases, two or more hydrogen bonds present in a dimer correspond to an interaction energy of −10 to −25 kcal mol−1, allowing a deformation energy above 1 kcal mol−1 (and up to 9.5 kcal mol−1). A lack of hydrogen bonding usually restricts the deformation energy to below 1 kcal mol−1 due to the weaker interaction energy. 

   more » « less
3. The non-covalently bound SO⋯H ₂ O system, including an interpretation of the differences between SO⋯H ₂ O and O ₂ ⋯H ₂ O

   https://doi.org/10.1039/C8CP05749D  

   Misiewicz, Jonathon P.; Noonan, Julia A.; Turney, Justin M.; Schaefer, Henry F. (November 2018, Physical Chemistry Chemical Physics)

   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
4. A Highly Correlated, Multireference Study of the Lowest Lying Singlet and Triplet States of the Four Thiophene Diradicals

   https://doi.org/10.1002/jcc.70044  

   Pandian, Joshua; Vu, Khanh; Muya, Jules Tshishimbi; Parker, Anna; Ancajas, Christine_Mae F; Saldana‐Greco, Diomedes; Yewer, Tabitha; Parish, Carol (January 2025, Journal of Computational Chemistry)

   ABSTRACT The energies and geometries of the lowest lying singlet and triplet states of the four diradicals formed by removing two H atoms from thiophene have been characterized. We utilized the highly correlated, multireference methods configuration interaction with single and double excitations with and without the Pople correction for size‐extensivity (MR‐CISD+Q and MR‐CISD) and averaged quadratic coupled cluster theory (MR‐AQCC). CAS (8,7) and CAS (10,8) active spaces involving σ, σ*, π, and π* orbitals were employed along with the cc‐pVDZ and cc‐pVTZ basis sets. The larger active space included the two electrons in the nonbonding sp²hybrid orbital on sulfur. We find that all didehydro isomers exist as planar, stable ground state singlets. The singlet‐triplet (S‐T) adiabatic gaps range from 15 to 25 kcal/mol while the vertical splittings are 21–35 kcal/mol. The 2,3 isomer has the lowest absolute ground state singlet energy and the largest adiabatic and vertical S‐T splitting. The ground states of the 2,3‐, and 2,5‐didehydrothiophene isomers are predicted to exhibit the smallest and largest diradical character, respectively, based on their electronic structures, spin densities and bonding analysis. To our knowledge, no experimental excitation energies of any of the didehydrothiophene isomers are available, and our computed MR‐AQCC/cc‐pVTZ data are believed to be among the most accurate computed results. This extensive study shows a competitive performance between MR‐AQCC and MR‐CISD+Q. 

   more » « less
5. Large transition state stabilization from a weak hydrogen bond

   https://doi.org/10.1039/D0SC02806A  

   Vik, Erik C.; Li, Ping; Maier, Josef M.; Madukwe, Daniel O.; Rassolov, Vitaly A.; Pellechia, Perry J.; Masson, Eric; Shimizu, Ken D. (July 2020, Chemical Science)

   A series of molecular rotors was designed to study and measure the rate accelerating effects of an intramolecular hydrogen bond. The rotors form a weak neutral O–H⋯OC hydrogen bond in the planar transition state (TS) of the bond rotation process. The rotational barrier of the hydrogen bonding rotors was dramatically lower (9.9 kcal mol −1 ) than control rotors which could not form hydrogen bonds. The magnitude of the stabilization was significantly larger than predicted based on the independently measured strength of a similar O–H⋯OC hydrogen bond (1.5 kcal mol −1 ). The origins of the large transition state stabilization were studied via experimental substituent effect and computational perturbation analyses. Energy decomposition analysis of the hydrogen bonding interaction revealed a significant reduction in the repulsive component of the hydrogen bonding interaction. The rigid framework of the molecular rotors positions and preorganizes the interacting groups in the transition state. This study demonstrates that with proper design a single hydrogen bond can lead to a TS stabilization that is greater than the intrinsic interaction energy, which has applications in catalyst design and in the study of enzyme mechanisms. 

   more » « less