Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

A Cost Effective Scheme for the Highly Accurate Description of Intermolecular Binding in Large Complexes

There has been a growing interest in quantitative predictions of the intermolecular binding energy of large complexes. One of the most important quantum chemical techniques capable of such predictions is the domain-based local pair natural orbital (DLPNO) scheme for the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)], whose results are extrapolated to the complete basis set (CBS) limit. Here, the DLPNO-based focal-point method is devised with the aim of obtaining CBS-extrapolated values that are very close to their canonical CCSD(T)/CBS counterparts, and thus may serve for routinely checking a performance of less expensive computational methods, for example, those based on the density-functional theory (DFT). The efficacy of this method is demonstrated for several sets of noncovalent complexes with varying amounts of the electrostatics, induction, and dispersion contributions to binding (as revealed by accurate DFT-based symmetry-adapted perturbation theory (SAPT) calculations). It is shown that when applied to dimeric models of poly(3-hydroxybutyrate) chains in its two polymorphic forms, the DLPNO-CCSD(T) and DFT-SAPT computational schemes agree to within about 2 kJ/mol of an absolute value of the interaction energy. These computational schemes thus should be useful for a reliable description of factors leading to the enthalpic stabilization of extended systems.

Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance

Hierarchical, convergent ab initio benchmark computations were performed followed by a systematic analysis of DFT performance for five pericyclic reactions comprising Diels-Alder, 1,3-dipolar cycloaddition, electrocyclic rearrangement, sigmatropic rearrangement, and double group transfer prototypes. Focal point analyses (FPA) extrapolating to the ab initio limit were executed via explicit quantum chemical computations with electron correlation treatments through CCSDT(Q) and correlation-consistent Gaussian basis sets up to aug'-cc-pV5Z. Optimized geometric structures and vibrational frequencies of all stationary points were obtained at the CCSD(T)/cc-pVTZ level of theory. The FPA reaction barriers and energies exhibit convergence to within a few tenths of a kcal mol-1. The FPA benchmarks were used to evaluate the performance of 60 density functionals (eight dispersion-corrected), covering the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. The meta-hybrid M06-2X functional provided the best overall performance [mean absolute error (MAE) of 1.1 kcal mol-1] followed closely by the double-hybrids B2K-PLYP, mPW2K-PLYP, and revDSD-PBEP86 [MAE of 1.4-1.5 kcal mol-1]. The regularly used GGA functional BP86 gave a higher MAE of 5.8 kcal mol-1, but it qualitatively described the trends in reaction barriers and energies. Importantly, we established that accurate yet efficient meta-hybrid or double-hybrid DFT potential energy surfaces can be acquired based on geometries from the computationally efficient and robust BP86/DZP level.

Weakly Bound Complex Formation between HCN and CH3Cl: A Matrix-Isolation and Computational Study

The matrix-isolated infrared spectrum of a hydrogen cyanide-methyl chloride complex was investigated in a solid argon matrix. HCN and CH₃Cl were co-condensed onto a substrate held at 10 K with an excess of argon gas, and the infrared spectrum was measured using Fourier-transform infrared spectroscopy. Quantum chemical geometry optimization, harmonic frequency, and natural bonding orbital calculations indicate stabilized hydrogen- and halogen-bonded structures. The two resulting weakly bound complexes are both composed of one CH₃Cl molecule bound to a (HCN)₃ subunit, where the three HCN molecules are bound head-to-tail in a ring formation. Our study suggests that─in the presence of CH₃Cl─the formation of (HCN)₃ is promoted through complexation. Since HCN aggregates are an important precursor to prebiotic monomers (amino acids and nucleobases) and other life-bearing polymers, this study has astrophysical implications toward the search for life in space.

High-Resolution Infrared Synchrotron Investigation of (HCN)2 and a Semi-Experimental Determination of the Dissociation Energy D0

The high‐resolution infrared absorption spectrum of the donor bending fundamental band ν61 of the homodimer (HCN)₂ has been collected by long‐path static gas‐phase Fourier transform spectroscopy at 207 K employing the highly brilliant 2.75 GeV electron storage ring source at Synchrotron SOLEIL. The rovibrational structure of the ν61 transition has the typical appearance of a perpendicular type band associated with a Σ–Π transition for a linear polyatomic molecule. The total number of 100 assigned transitions are fitted employing a standard semi‐rigid linear molecule Hamiltonian, providing the band origin ν₀ of 779.05182(50) cm⁻¹ together with spectroscopic parameters for the degenerate excited state. This band origin, blue‐shifted by 67.15 cm⁻¹ relative to the HCN monomer, provides the final significant contribution to the change of intra‐molecular vibrational zero‐point energy upon HCN dimerization. The combination with the vibrational zero‐point energy contribution determined recently for the class of large‐amplitude inter‐molecular fundamental transitions then enables a complete determination of the total change of vibrational zero‐point energy of 3.35±0.30 kJ mol⁻¹. The new spectroscopic findings together with previously reported benchmark CCSDT(Q)/CBS electronic energies [Hoobler et al. ChemPhysChem. 19, 3257–3265 (2018)] provide the best semi‐experimental estimate of 16.48±0.30 kJ mol⁻¹ for the dissociation energy D₀ of this prototypical homodimer.