More Like this

1. The HOX⋯SO ₂ (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry

   https://doi.org/10.1002/cphc.202000746  

   Kitzmiller, Nathaniel_L; Wolf, Mark_E; Turney, Justin_M; Schaefer, III, Henry_F (December 2020, ChemPhysChem)

   Abstract Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug‐cc‐pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug‐cc‐pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug‐cc‐pwCVDZ‐PP pseudopotential. 27 HOX⋯SO₂complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂complexes in this study ranged from 1.35 to 3.81 kcal mol⁻¹. The single‐interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol⁻¹, while the single‐interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol⁻¹, indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂which may guide further research of related systems. 

   more » « less
2. Semi-experimental equilibrium ( r e SE) and theoretical structures of hydrazoic acid (HN3)

   https://doi.org/10.1063/5.0101064  

   Owen, Andrew N.; Sahoo, Nitai P.; Esselman, Brian J.; Stanton, John F.; Woods, R. Claude; McMahon, Robert J. (July 2022, The Journal of Chemical Physics)

   Hydrazoic acid (HN3) is used as a case study for investigating the accuracy and precision by which a molecular structure—specifically, a semi-experimental equilibrium structure (reSE)—may be determined using current state-of-the-art methodology. The influence of the theoretical corrections for effects of vibration–rotation coupling and electron-mass distribution that are employed in the analysis is explored in detail. The small size of HN3 allowed us to deploy considerable computational resources to probe the basis-set dependence of these corrections using a series of coupled-cluster single, double, perturbative triple [CCSD(T)] calculations with cc-pCVXZ (X = D, T, Q, 5) basis sets. We extrapolated the resulting corrections to the complete basis set (CBS) limit to obtain CCSD(T)/CBS corrections, which were used in a subsequent reSE structure determination. The reSE parameters obtained using the CCSD(T)/cc-pCV5Z corrections are nearly identical to those obtained using the CCSD(T)/CBS corrections, with uncertainties in the bond distances and angles of less than 0.0006 Å and 0.08°, respectively. The previously obtained reSE structure using CCSD(T)/ANO2 agrees with that using CCSD(T)/cc-pCV5Z to within 0.000 08 Å and 0.016° for bond distances and angles, respectively, and with only 25% larger uncertainties, validating the idea that reSE structure determinations can be carried out with significantly smaller basis sets than those needed for similarly accurate, strictly ab initio determinations. Although the purely computational re structural parameters [CCSD(T)/cc-pCV6Z] fall outside of the statistical uncertainties (2σ) of the corresponding reSE structural parameters, the discrepancy is rectified by applying corrections to address the theoretical limitations of the CCSD(T)/cc-pCV6Z geometry with respect to basis set, electron correlation, relativity, and the Born–Oppenheimer approximation, thereby supporting the contention that the semi-experimental approach is both an accurate and vastly more efficient method for structure determinations than is brute-force computation. 

   more » « less
3. Benchmark coupled-cluster lattice energy of crystalline benzene and assessment of multi-level approximations in the many-body expansion

   https://doi.org/10.1063/5.0159410  

   Borca, Carlos H.; Glick, Zachary L.; Metcalf, Derek P.; Burns, Lori A.; Sherrill, C. David (June 2023, The Journal of Chemical Physics)

   The many-body expansion (MBE) is promising for the efficient, parallel computation of lattice energies in organic crystals. Very high accuracy should be achievable by employing coupled-cluster singles, doubles, and perturbative triples at the complete basis set limit [CCSD(T)/CBS] for the dimers, trimers, and potentially tetramers resulting from the MBE, but such a brute-force approach seems impractical for crystals of all but the smallest molecules. Here, we investigate hybrid or multi-level approaches that employ CCSD(T)/CBS only for the closest dimers and trimers and utilize much faster methods like Møller–Plesset perturbation theory (MP2) for more distant dimers and trimers. For trimers, MP2 is supplemented with the Axilrod–Teller–Muto (ATM) model of three-body dispersion. MP2(+ATM) is shown to be a very effective replacement for CCSD(T)/CBS for all but the closest dimers and trimers. A limited investigation of tetramers using CCSD(T)/CBS suggests that the four-body contribution is entirely negligible. The large set of CCSD(T)/CBS dimer and trimer data should be valuable in benchmarking approximate methods for molecular crystals and allows us to see that a literature estimate of the core-valence contribution of the closest dimers to the lattice energy using just MP2 was overbinding by 0.5 kJ mol−1, and an estimate of the three-body contribution from the closest trimers using the T0 approximation in local CCSD(T) was underbinding by 0.7 kJ mol−1. Our CCSD(T)/CBS best estimate of the 0 K lattice energy is −54.01 kJ mol−1, compared to an estimated experimental value of −55.3 ± 2.2 kJ mol−1. 

   more » « less
4. Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory

   https://doi.org/10.3390/ijms242417480  

   Rock, Carly A.; Tschumper, Gregory S. (December 2023, International Journal of Molecular Sciences)

   This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three–five water molecules (Ar(H2O)n=3–5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body–many-body (2b:Mb) and three-body–many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol−1 and by nearly 6 kJ mol−1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm−1 for some hydrogen-bonded OH stretching frequencies. 

   more » « less
5. Conformational comparison of urea and thiourea near the CCSD ( T ) complete basis set limit

   https://doi.org/10.1002/qua.27075  

   Barlow, Kayleigh R.; Tschumper, Gregory S. (December 2022, International Journal of Quantum Chemistry)

   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 C_2vsecond‐order saddle point of urea lies at least 1.5 kcal mol⁻¹above the C₂global 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 C_slocal 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 NH₂groups 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