Structural Distortions Accompanying Noncovalent Interactions: Methane–Water, the Simplest C–H Hydrogen Bond

Simple Composite Approach to Efficiently Estimate Basis Set Limit CCSD(T) Harmonic Frequencies and Reaction Thermochemistry

We introduce a simple strategy that combines the G3(MP2) composite method and explicitly correlated coupled cluster CCSD(T)-F12 method to efficiently estimate complete basis set CCSD(T) molecular geometries and harmonic vibrational frequencies at the cost of a double-ζ basis set calculation. Based on a large test set of 61 neutral, ionic, and open-shell molecules, and additionally 31 molecules in the HFREQ2014 data set, we demonstrate that this composite strategy has an average accuracy of 2 cm-1 or better relative to complete basis set CCSD(T) values. Using this approach, we estimated 696 CCSD(T)/CBS reaction energies of small to medium-sized systems containing up to 6 heavy atoms and confirmed existing approximations that use small basis set density functional theory methods [e.g., M06-2X/6-31+G(d)] to calculate thermal contributions to reaction enthalpies and Gibbs free energies that are accurate to within 0.2 kcal mol-1 on average.

The effect of alkylation on the micro-solvation of ethers revealed by highly localized water librational motion

The specific far-infrared spectral signatures associated with highly localized large-amplitude out-of-plane librational motion of water molecules have recently been demonstrated to provide sensitive spectroscopic probes for the micro-solvation of organic molecules [Mihrin et al., Phys. Chem. Chem. Phys. 21(4), 1717 (2019)]. The present work employs this direct far-infrared spectroscopic approach to investigate the non-covalent intermolecular forces involved in the micro-solvation of a selection of seven ether molecules with systematically varied alkyl substituents: dimethyl ether, diethyl ether, diisopropyl ether, ethyl methyl ether, t-butyl methyl ether, and t-butyl ethyl ether. The ranking of the observed out-of-plane water librational band signatures for this selected series of ether–water complexes embedded in inert neon matrices at 4 K reveals information about the interplay of directional intermolecular hydrogen bond motifs and non-directional and long-range dispersion interactions for the micro-solvated structures. These far-infrared observables differentiate minor subtle effects introduced by specific alkyl substituents and serve as rigorous experimental benchmarks for modern quantum chemical methodologies of various levels of scalability, which often fail to accurately predict the structural variations and corresponding vibrational signatures of the closely related systems. The accurate interaction energies of the series of ether–water complexes have been predicted by the domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations, followed by a local energy decomposition analysis of the energy components. In some cases, the secondary dispersion forces are in direct competition with the primary intermolecular hydrogen bonds as witnessed by the specific out-of-plane librational signatures.

Normal-Mode Vibrational Analysis of Weakly Bound Oligomers at Constrained Stationary Points of Arbitrary Order

Following the full realization of the importance of noncovalent interactions, finding and characterizing stationary points (SP), of various order, for weakly bound oligomers have become important tasks for computational chemists. An efficient algorithm and an associated computer code, called oligoCGO, are described, facilitating constrained geometry optimization of oligomers of arbitrary structure and complexity and normal-mode vibrational analysis at nonstationary geometries. To minimize the adverse effects of nonzero forces on harmonic vibrational analyses at constrained stationary points (cSP), two residual gradient correction (RGC) schemes are proposed. RGC₁, for which a rigorous justification is given, is based on ignoring the remaining forces in internal-coordinate space. RGC₂ modifies the geometry of the cSP in a single Newton step and recalculates the Cartesian Hessian at this updated geometry. As demonstrated by 10 examples related to the water-water, water-methane, and methane-methane dimers as well as the methane trimer, without RGC the harmonic analysis of cSPs may result in even qualitatively incorrect results when compared to reference values obtained at the nearby unconstrained SPs (uSP). Both RGC protocols work exceedingly well, and the corrected harmonic wavenumbers of the cSPs are very close to their uSP counterparts.

Hydrogen bonded dimers of ketocoumarin in the solid state and alcohol:water binary solvent: fluorescence spectroscopy, crystal structure and DFT investigation

Fluorescence emission spectrum of ketocoumarin dimers in an alcohol:water binary mixture and the solid state.

Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states

Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H₂O)₂, CH₄·H₂O, and (CH₄)₂.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol^-1 to -0.33 kcal mol^-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of a_e(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero ²J(C_M-H_M···O_F) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

Rovibrational quantum dynamical computations for deuterated isotopologues of the methane–water dimer

Rovibrational states of methane–water isotopologues are computed in a variational procedure and the wave functions are analyzed in terms of the rigid-rotor and coupled-rotors models.

Unravelling hydrogen bonding interactions of tryptamine–water dimer from neutral to cation

The physical origin for the three intermolecular hydrogen bonds in the neutral and cationic forms of the tryptamine–water dimer is explored.