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
Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as đȘ(đ9) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the đȘ(đ7) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the đȘ(đ9) CCSD(TQ). Further, we examine the performance of Goodsonâs continued fraction method in the estimation of CCSDT(Q)Î contributions to reaction energies as well as a ânewâ method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Î/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal molâ1 and an L2D of 0.52 kcal molâ1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably âballparkâ how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.
more »
« less
- Award ID(s):
- 2143725
- PAR ID:
- 10581541
- Publisher / Repository:
- ACS
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry A
- Volume:
- 128
- Issue:
- 36
- ISSN:
- 1089-5639
- Page Range / eLocation ID:
- 7720 to 7732
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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âŻSO2binary 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âŻSO2complexes 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âŻSO2complexes in this study ranged from 1.35 to 3.81â kcalâmolâ1. The singleâinteraction hydrogen bonded complexes spanned a range of 2.62 to 3.07â kcalâmolâ1, while the singleâinteraction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06â kcalâmolâ1, indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO2which may guide further research of related systems.more » « less
-
Abstract While limited coupled cluster theory isformallynonvariational, it is not broadly appreciated whether this is a major issuein practice. We carried out a detailed comparison withde factofull CI energies for a relatively large and diverse set of molecules near equilibrium geometries. Fully iterative limited CC methods such as CCSDT, CCSDTQ, and CCSDTQ5 do represent practical upper bounds to the FCI energy, as do CCSDT-3, CCSDTQ-3, and CCSD(cT). While quasiperturbative approaches such as CCSD(T) and especially CCSDT(Q) may significantly over-correlate molecules if there is significant static correlation, this is much less of an issue with Lambda approaches such as CCSDT(Q)Î.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 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
-
We have developed a new composite model chemistry method called WMS (WuhanâMinnesota scaling method) with three characteristics: (1) a composite scheme to approximate the complete configuration interaction valence energy with the affordability condition of requiring no calculation more expensive than CCSD(T)/jul-cc-pV(T+d)Z, (2) low-cost methods for the inner-shell correlation contribution and scalar relativistic correction, and (3) accuracy comparable to methods with post-CCSD(T) components. The new method is shown to be accurate for the W4-17 database of 200 atomization energies with an average mean unsigned error (averaged with equal weight over strongly correlated and weakly correlated subsets of the data) of 0.45 kcal mol â1 , and the performance/cost ratio of these results compares very favorably to previously available methods. We also assess the WMS method against the DBH24-W4 database of diverse barrier heights and the energetics of the reactions of three strongly correlated Criegee intermediates with water. These results demonstrate that higher-order correlation contributions necessary to obtain high accuracy for molecular thermochemistry may be successfully extrapolated from the lower-order components of CCSD(T) calculations, and chemical accuracy can now be obtained for larger and more complex molecules and reactions.more » « less
An official website of the United States government

