Exact geometries from quantum chemical calculations

Reply to a Comment on "The Nature of Chalcogen-Bonding-Type Tellurium-Nitrogen Interactions"

We reply to the comment by J.‐M. Mewes, A. Hansen and S. Grimme (MHG), who challenged the accuracy of our r_e value for the N⋅⋅⋅Te distance in (C₆F₅)Te(CH₂)₃NMe₂ determined by gas electron diffraction (GED). We conclusively demonstrate that MHG′s quoted reference calculations are less accurate than they claim for solid state and gas phase. We show by higher level calculations, that we did not miss substantial contributions from open‐chain conformers. Refinements on simulated scattering data show that such contributions would have had only an almost negligible effect on r_e(N⋅⋅⋅Te). MHG suggested the use of a H0‐tuned GFN method for calculating vibrational corrections r_a−r_e, but this did not change these values substantially. Alternative amplitude calculations using higher level analytic harmonic and numeric cubic force fields (PBE0‐D3BJ/def2‐TZVP) yield a GED value for r_e(N⋅⋅⋅Te)=2.852(25) Å that is well within the experimental error of our original value 2.918(31) Å but far from the 2.67(8) Å predicted by MHG. A now improved error estimation accounts for inaccuracies in the calculated auxiliary values. The gas/solid difference of the weak N⋅⋅⋅Te interaction is in a realistic range compared to other systems involving weak chemical interactions.

We reply to the comment which challenged the accuracy of our r_e value for the N⋅⋅⋅Te distance in (C₆F₅)Te(CH₂)₃NMe₂ determined by gas electron diffraction. We show by higher level calculations, that we did not miss substantial contributions from open‐chain conformers.

Investigation of the Abstraction and Dissociation Mechanism in the Nitrogen Trifluoride Channels: Combined Post-Hartree-Fock and Transition State Theory Approaches

The present paper concludes our series of kinetics studies on the reactions involved in the complex mechanism of nitrogen trifluoride decomposition. Two other related reactions that, along with this mechanism, take part in an efficient boron nitride growth process are also investigated. We report results concerning two abstraction reactions, namely NF2 + N ⇄ 2NF and NF3 + NF ⇄ 2NF2, and two dissociations, N2F4 ⇄ 2NF2 and N2F3 ⇄ NF2 + NF. State-of-the-art electronic structure calculations at the CCSD(T)/cc-pVTZ level of theory were considered to determine geometries and frequencies of reactants, products, and transition states. Extrapolation of the energies to the complete basis set limit was used to obtain energies of all the species. We applied transition state theory to compute thermal rate constants including Wigner, Eckart, Bell, and deformed theory corrections in order to take tunneling effects into account. The obtained results are in good agreement with the experimental data available in the literature and are expected to provide a better phenomenological understanding of the NF3 decomposition role in the boron nitride growth for a wide range of temperature values.

Vibrational Properties of the Isotopomers of the Water Dimer Derived from Experiment and Computations

The water dimer and its 11 deuterated isotopomers are investigated utilizing coupled cluster theory and experimental data as input for a perturbational determination of the isotopomer frequencies. Deuterium substitution reduces the H-bond stretching frequency by maximally 12 cm–1 from 143 to 131 cm–1, which makes a spectroscopic differentiation of H- and D-bonds difficult. However, utilizing the 132 frequencies obtained in this work, the identification of all isotopomers is straightforward. The CCSD(T)/CBS value of the binding energy De is 5.00 kcal mol–1. The binding energy D0 of the water dimer increases upon deuterium substitution from 3.28 to maximally 3.71 kcal mol–1 reflecting a decrease in the zero point energy contribution. The entropy values of the D-isotopomers increase from 73 to 77 entropy units in line with the general observation that a mass increase leads to larger entropies. All 12 isotopomers possess positive free binding energies at 80 K and a reduced pressure of 110 Pa, which means that they can be spectroscopically observed under these conditions.

Identification of the strongest bonds in chemistry

Increasing the effective electronegativity of two atoms forming a triple bond can increase the strength of the latter. The strongest bonds found in chemistry involve protonated species of hydrogen cyanide, carbon monoxide, and dinitrogen. CCSD(T)/CBS (complete basis set) and G4 calculations reveal that bond dissociation energies are misleading strength descriptors. The strength of the bond is assessed via the local stretching force constants, which suggest relative bond strength orders (RBSO) between 2.9 and 3.4 for heavy atom bonding (relative to the CO bond strength in methanol (RBSO = 1) and formaldehyde (RBSO = 2)) in [HCNH](+)((1)Σ(+)), [HCO](+)((1)Σ(+)), [HNN](+)((1)Σ(+)), and [HNNH](2+)((1)Σg(+)). The increase in strength is caused by protonation, which increases the electronegativity of the heavy atom and thereby decreases the energy of the bonding AB orbitals (A, B: C, N, O). A similar effect can be achieved by ionization of a nonbonding or antibonding electron in CO or NO. The strongest bond with a RBSO value of 3.38 is found for [HNNH](2+) using scaled CCSD(T)/CBS frequencies determined for CCSD(T)/CBS geometries. Less strong is the NN bond in [FNNH](2+) and [FNNF](2+).

The gas-phase acidities of the elemental hydrides are functions of electronegativity and bond length

The gas-phase Brønsted acidities of the group 1, group 2, and main group elemental hydrides (XHn) are shown to be a combined function of the bond length, electronegativity, and position in the periodic table, via a separation of the acidity into coulombic and electronic reorganization enthalpy parts. The Coulombic acidity is defined as the enthalpy to separate unit positive and negative charges from the neutral acid's X—H bond length to infinity; the reorganization enthalpy is the difference between that and the measured acidity, and represents the enthalpy required to reorganize the electrons of the neutral acid, creating an ion pair at the original bond distance. Predictions are made for the gas-phase Brønsted acidities of several elemental hydrides for which this quantity is not known.Key words: acidity, gas phase, coulomb, elements, hydride.

Theoretical UV Circular Dichroism of Aliphatic Cyclic Dipeptides

Four cyclic dipeptides (piperazine-2,5-diones), cyclo(L-Pro-Gly), cyclo(L-Pro-L-Leu), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-L-Ala), were modeled from crystal structure data. Conformations resulting from energy minimization using molecular mechanics were compared with traditional ab initio and density functional theory geometric optimizations for each dipeptide. In all computational cases, the gas phase was assumed. The pi-pi transition feature of the UV circular dichroic (CD) spectra was predicted for each peptide structure via the classical dipole interaction model. The dipole interaction model predicted CD spectra that qualitatively agreed with experiment when MP2 or DFT geometries were used. By coupling MP2 or DFT geometric optimizations with the classical physics method of the dipole interaction model, significantly better CD spectra were calculated than those using geometries obtained by molecular mechanics. Thus, one can couple quantum mechanical geometries with a classical physics model for calculation of circular dichroism.

Gas-Phase Infrared and ab Initio Study of the Unstable CF3CNO Molecule and Its Stable Furoxan Ring Dimer

The unstable trifluoroacetonitrile N-oxide molecule, CF3CNO, has been generated in high yield in the gas phase from CF3BrC=NOH and studied for the first time by gas-phase mid-infrared spectroscopy. Cold trapping of this molecule followed by slow warming forms the stable ring dimer, bis(trifluoromethyl)furoxan, also investigated by gas-phase infrared spectroscopy. The spectroscopy provides an investigation into the vibrational character of the two molecules, the assignments supported by calculations of the harmonic vibrational frequencies using in the case of CF3CNO both ab initio (CCSD(T)) and density functional theory (B3LYP) and B3LYP for the ring dimer. The ground-state structures of both molecules were investigated at the B3LYP level of theory, with CF3CNO further investigated using coupled-cluster. The CCSD(T) method suggests a slightly bent (C(s)) structure for CF3CNO, while the B3LYP method (with basis sets ranging from 6-311G(d) to cc-pVTZ) suggests a close-to-linear or linear CCNO chain. The CCN bending potential in CF3CNO was explored at the CCSD(T)(fc)/cc-pVTZ level, with the results suggesting that CF3CNO exhibits strong quasi-symmetric top behavior with a barrier to linearity of 174 cm(-1). Since both isomerization and dimerization are feasible loss processes for this unstable molecule, the relative stability of CF3CNO with respect to the known cyanate (CF3OCN), isocyanate (CF3NCO), and fulminate (CF3ONC) isomers and the mechanism of the dimerization process to the ring furoxan and other isomers were studied with density functional theory.

Focal-Point Conformational Analysis of Ethanol, Propanol, and Isopropanol

Conformational analysis of three small alcohols--ethanol, propanol, and isopropanol--was carried out by systematically improving the basis set and the level of electron correlation. Correlation energy contributions to conformational energies are strongly basis-set-dependent but accurate energy contributions can be obtained by extrapolation to the basis-set limit. At the basis-set limit, second- and third-order electron correlation effects play a significant role for rotations around the CC-OH, HC-CO, and CC-CO bonds. Specifically, second- and third-order correlation effects strongly stabilize structures in which the hydroxylic hydrogen eclipses with the adjacent carbon; a lesser stabilization is present in structures where the CC-OH moiety is in the gauche form. Fourth-order correlation effects to the CC-OH rotation are small due to a partial cancellation of the singles, doubles, and quadruples contribution by the triples contribution. Electron correlation significantly lowers barriers for methyl-group rotations in ethanol and isopropanol, and in these cases the fourth-order correlation effects are noticeable. The relatively large overall importance of third-order correlation energy contributions raises a concern that the inability to accurately estimate this slowly converging contribution may become a limiting factor when highly accurate conformational energies in larger molecules are sought.

The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.

Complete basis set extrapolated potential energy, dipole, and polarizability surfaces of alkali halide ion-neutral weakly avoided crossings with and without applied electric fields

Complete basis set extrapolations of alkali halide (LiF, LiCl, NaF, NaCl) energy, dipole, and polarizability surfaces are performed with and without applied fields along the internuclear axis using state-averaged multireference configuration interaction. Comparison between properties (equilibrium separation, dissociation energy, crossing distance, diabatic coupling constant, dipole, and polarizability) derived from the extrapolated potential energy (or dipole) surfaces are made with those obtained from direct extrapolation from the basis set trends. The two extrapolation procedures are generally found to agree well for these systems. Crossing distances from this work are compared to those of previous work and values obtained from the Rittner potential. Complete basis set extrapolated crossing distances agree well with those derived from the Rittner potential for LiF, but were significantly larger for LiCl, NaF, and NaCl. The results presented here serve as an important set of benchmark data for the development of new-generation many-body force fields that are able to model charge transfer.