Large-amplitude motion in 1,4-cyclohexadiene and 1,4-dioxin: theoretical background for joint treatment of spectroscopic, electron diffraction and ab initio data

A simplistic computational procedure for tunneling splittings caused by proton transfer

In this manuscript, we present an approach for computing tunneling splittings for large amplitude motions. The core of the approach is a solution of an effective one-dimensional Schrödinger equation with an effective mass and an effective potential energy surface composed of electronic and harmonic zero-point vibrational energies of small amplitude motions in the molecule. The method has been shown to work in cases of three model motions: nitrogen inversion in ammonia, single proton transfer in malonaldehyde, and double proton transfer in the formic acid dimer. In the current work, we also investigate the performance of different DFT and post-Hartree–Fock methods for prediction of the proton transfer tunneling splittings, quality of the effective Schrödinger equation parameters upon the isotopic substitution, and possibility of a complete basis set (CBS) extrapolation for the resulting tunneling splittings.

Equilibrium molecular structure and spectra of 6-methyl-1,5-diazabicyclo[3.1.0]hexane: joint analysis of gas phase electron diffraction, quantum chemistry, and spectroscopic data

The equilibrium geometry of the boat conformation (Cs point group symmetry) of the 6-methyl-1,5-diazabicyclo[3.1.0]hexane (MDABH) molecule, absolutely dominating under normal conditions, was studied by the gas-phase electron diffraction (GED) method at 20 °C with the involvement of NMR, IR, and Raman spectroscopic data and quantum chemical calculations. The potential function of ring-puckering deformation for the MDABH bicyclic system was calculated at the MP2/aug-cc-pVTZ and B3LYP/cc-pVTZ levels. It was found by MP2 calculation that the total energy of the boat conformation is 3.52 kcal mol-1 lower than that of the chair conformation. For the first time, we recorded the IR and Raman spectra for liquid samples of MDABH and assigned their peculiarities only to boat conformation vibrations using the Pulay technique of scaling quantum chemical force fields. In the case of the chair form, transferability of the refined scale factors was used for reliable prediction of the location of its fundamental frequencies. According to the joint structural analysis of the above data, the most important equilibrium geometric re-parameters for the boat conformation of the MDABH molecule were determined to be (bond lengths in Å; angles in degrees, Cs symmetry): C2N1 = 1.466(2), C2C3 = 1.523(2), N1N5 = 1.512(2), C6N1 = 1.440(2), C6C7 = 1.487(2), ∠C2N1N5 = 106.1(2), ∠N1C2C3) = 110.2(4), ∠C2C3C4 = 99.9(4), ∠N1N5C6 = 58.3(1), ∠N1C6N5 = 63.3(1), ∠N1C6C7 = 114.9(6), ∠C6N1C2 = 111.8(1), ∠N5N1C2C3 = 17.3(1), ∠N1C2C3C4 = -26.8(2), θ = C2C3C4/C2N1N5C4 = -26.2(3), φ = N1C6N5/C2N1N5C4 = 74.0(1). Comparison of these and earlier results showed that the NN bond length in the diaziridine ring is very weakly dependent on the cis- or trans-arrangement of substituents at the nitrogen atoms.

Equilibrium structures of the tetramezine diastereomers and their ratio: joint analysis of gas phase electron diffraction, quantum chemistry, and spectroscopic data

For the first time, we applied a gas-phase electron diffraction (GED) method together with vibrational spectroscopy and quantum chemical calculations to investigate the equilibrium geometries of achiral meso and enantiomeric diastereomers of tetramezine [1,2-bis-(3,3-dimethyldiaziridin-1-yl)ethane] and their ratio in the mixture. In the joint structural analysis of these data, a new approach based on PES parameters is used in the framework of a static molecular model (small amplitude motion approximation). The agreement between the theoretical and experimental molecular intensities is characterized by a divergence factor Rf of 5.9%. The experimental re-parameters of tetramezine diastereomers agreed with our B3LYP/cc-pVTZ and MP2/cc-pVTZ calculations, which predicted the total energy of the meso form (Ci point group symmetry) to be lower than that of the enantiomeric form (C2 point group symmetry), by 6.4 and 4.7 kJ mol-1, respectively. The experimentally measured percentages of the meso and both enantiomeric diastereoisomers at 360 K were 70% and 30%, respectively. We characterized the meso form using 2D NMR spectra. Our GED data are in good agreement with the X-ray diffraction analysis of the meso form. This result reflects the weak effect of intermolecular interactions in the crystal. We assigned the IR spectrum bands of the crystalline meso form using the Pulay technique of scaling quantum chemical force fields. In the case of the enantiomeric form calculated at the same level, transferability of the refined scale factors was used for more reliable prediction of the mutual location and interpretation of its fundamental frequencies.

Application of classical simulations for the computation of vibrational properties of free molecules

In this study, we investigate the ability of classical molecular dynamics (MD) and Monte-Carlo (MC) simulations for modeling the intramolecular vibrational motion. These simulations were used to compute thermally-averaged geometrical structures and infrared vibrational intensities for a benchmark set previously studied by gas electron diffraction (GED): CS₂, benzene, chloromethylthiocyanate, pyrazinamide and 9,12-I₂-1,2-closo-C₂B₁₀H₁₀. The MD sampling of NVT ensembles was performed using chains of Nose-Hoover thermostats (NH) as well as the generalized Langevin equation thermostat (GLE). The performance of the theoretical models based on the classical MD and MC simulations was compared with the experimental data and also with the alternative computational techniques: a conventional approach based on the Taylor expansion of potential energy surface, path-integral MD and MD with quantum-thermal bath (QTB) based on the generalized Langevin equation (GLE). A straightforward application of the classical simulations resulted, as expected, in poor accuracy of the calculated observables due to the complete neglect of quantum effects. However, the introduction of a posteriori quantum corrections significantly improved the situation. The application of these corrections for MD simulations of the systems with large-amplitude motions was demonstrated for chloromethylthiocyanate. The comparison of the theoretical vibrational spectra has revealed that the GLE thermostat used in this work is not applicable for this purpose. On the other hand, the NH chains yielded reasonably good results.