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
A comparison between hydrogen and halogen bonding: the hypohalous acid–water dimers, HOX⋯H 2 O (X = F, Cl, Br)
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
- Award ID(s):
- 1661604
- PAR ID:
- 10118793
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 21
- Issue:
- 11
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 6160 to 6170
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The ability of two anions to interact with one another is tested in the context of pairs of TrX 4 − homodimers, where Tr represents any of the triel atoms B, Al, Ga, In, or Tl, and X refers to a halogen substituent F, Cl, or Br. None of these pairs engage in a stable complex in the gas phase, but the situation reverses in water where the two monomers are held together by Tr⋯X triel bonds, complemented by stabilizing interactions between X atoms. Some of these bonds are quite strong, notably those involving TrF 4 − , with interaction energies surpassing 30 kcal mol −1 . Others are very much weaker, with scarcely exothermic binding energies. The highly repulsive electrostatic interactions are counteracted by large polarization energies.more » « less
-
A quantitative assessment of deformation energy in intermolecular interactions: How important is it?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
-
This study systematically characterizes the four homogeneous and six heterogeneous hydrogen-bonded dimers formed by hydrogen halide pairs (HX/HY where X, Y = F, Cl, Br, and I). The notation HX⋯HY indicates the direction of the hydrogen bond from the HY donor to the HX acceptor. All stationary points reported for these ten dimer systems are fully optimized utilizing the MP2 and CCSD(T) ab initio methods in conjunction with quadruple-ζ correlation-consistent basis sets augmented with diffuse functions, and their nature is verified by harmonic vibrational frequency computations. The electronic dissociation energies (De) for all ten global minima are evaluated near the CCSD(T) complete basis set (CBS) limit via extrapolation schemes. These values are 19.11, 8.32, 7.38, and 6.22 kJ mol−1 for the homogeneous dimers of HF, HCl, HBr, and HI, respectively. For the heterogeneous pairs, the lighter hydrogen halide is consistently the donor in the global minimum configuration, with De ranging from 12.23 kJ mol−1 for HCl⋯HF to 7.22 kJ mol−1 for HI⋯HBr near the CCSD(T) CBS limit. Interestingly, not all heterodimer donor/acceptor permutations correspond to minima. For example, the HCl⋯HBr configuration is identified as a local minimum at all levels of theory employed in this investigation, whereas the in-plane barrier for donor/acceptor exchange vanishes for HCl⋯HI and HBr⋯HI when larger quadruple-ζ basis sets are utilized. For the seven dimer systems containing Br and/or I, the structures, energetics, and vibrational frequencies computed using conventional valence-only electron correlation procedures are similar to those obtained using an expanded valence treatment that includes the (n − 1)d subvalence electrons associated with Br and I.more » « less
-
Abstract The replacement of a CH group of benzene by a triel (Tr) atom places a positive region of electrostatic potential near the Tr atom in the plane of the aromatic ring. This σ‐hole can interact with an X lone pair of XCCH (X=F, Cl, Br, and I) to form a triel bond (TrB). The interaction energy between C5H5Tr and FCCH lies in the range between 2.2 and 4.4 kcal/mol, in the order Tr=B+cation above the ring pulls density toward itself and thus magnifies the Tr σ‐hole. The TrB to the XCCH nucleophile is thereby magnified as is the strength of the TrB. This positive cooperativity is particularly large for Tr=B.more » « less
An official website of the United States government

