Tunneling splittings in vibrational spectra of non-rigid molecules

Theoretical Study of Proton Tunneling in the Imidazole-Imidazolium Complex

Proton tunneling in the hydrogen-bonded imidazole–imidazolium complex ion has been studied theoretically. Ab initio CASSCF/6-311++G(d,p) calculations concerning geometry optimization and vibrational frequencies have been carried out for equilibrium and transition state structures of the system. Two-dimensional double-well model potentials were constructed on the basis of ab initio results and used to analyze the proton dynamics in the hydrogen bond and the influence of the excitation of low-frequency hydrogen-bond vibrations on the proton tunneling splittings. The energy of tunneling-split vibrational sublevels of the high-frequency tunneling mode have been calculated for its ground and first excited vibrational state for the series of excitations of the coupled low-frequency intramolecular hydrogen-bond modes. The promoting and suppressing effect of the low-frequency modes on the proton splittings was shown in the ground and first excited vibrational state of the tunneling mode. The vibrational sublevels form the two separate semicontinuous bands between which the absorption transitions may occur. This mechanism explains the experimentally observed splitting and doublet-component broadening of the high-frequency N–H stretching infrared (IR) absorption band.

Density Functional Theory and Atoms-in-Molecules Investigation of Intramolecular Hydrogen Bonding in Derivatives of Malonaldehyde and Implications for Resonance-Assisted Hydrogen Bonding

A density functional theory (DFT) and atoms-in-molecules (AIM) analysis has been applied to the intramolecular hydrogen bonding in the enol conformers of malonaldehyde and its fluoro-, chloro-, cyano-, and nitro-substituted derivatives. With the B3LYP/6-311++G(2d,p) method, good agreement between the DFT geometries and published experimental structures has been found. The donor-acceptor distance was also varied in a series of constrained optimizations in order to determine if energetic, structural, and topological trends associated with intermolecular hydrogen bonding remain valid in the intramolecular case. At very short donor-acceptor distances (<2.24 A), the hydrogen is symmetrically located between donor and acceptor; at distances longer than this, the hydrogen bonding is no longer symmetric. The AIM methodology has been applied to explore the topology of the electron density in the intramolecular hydrogen bonds of the chosen model systems. Most AIM properties for intramolecular hydrogen bond distances longer than 2.24 A show smooth trends, consistent with intermolecular hydrogen bonds. Integrated AIM properties have also been used to explore the phenomenon of resonance-assisted hydrogen bonding (RAHB). It is shown that as the donor-acceptor distance is varied, pi-electron density is redistributed among the carbon atoms in the intramolecular hydrogen bond ring; however, contrary to prior studies, the integrated atomic charges on the donor-acceptor atoms were found to be insensitive to variation of hydrogen-bonding distance.

Molecular dynamics integration and molecular vibrational theory. II. Simulation of nonlinear molecules

A series of molecular dynamics (MD) simulations of nonlinear molecules has been performed to test the efficiency of newly introduced semianalytical second-order symplectic time-reversible MD integrators that combine MD and the standard theory of molecular vibrations. The simulation results indicate that for the same level of accuracy, the new algorithms allow significantly longer integration time steps than the standard second-order symplectic leap-frog Verlet method. Since the computation cost per integration step using new MD integrators with longer time steps is approximately the same as for the standard method, a significant speed-up in MD simulation is achieved.

An intrinsic reaction coordinate calculation of the torsion-internal rotation potential of hydrogen peroxide and its isotopomers

Intrinsic reaction coordinate (IRC) calculations of the internal rotation (torsional) potentials for H(2)O(2) and its isotopomers HDO(2) and D(2)O(2) were carried out at the CCSD(T)/CBS//aug-cc-pVDZ level. Two extrapolation methods were used to obtain energies in the complete basis set (CBS) limit. The full IRC potential was constructed from scans from the C(2v) (cis) and C(2h) (trans) transition states to the equilibrium C(2) (gauche) structure. The IRC potential for H(2)O(2) was fit to a five-term Fourier function; coefficients were compared with values obtained from spectroscopic data. The twofold IRC torsional potentials were used to obtain torsional eigenvalues, which yielded values of the transitions between various ntau states. These results compare favorably with Raman and near-infrared data. Our calculations provide values of the cis and trans barriers of 2495 and 364 cm(-1), respectively, which are in good agreement with both previously calculated and experimentally derived values. It appears that coupling between torsional motion and other degrees of freedom is not significant in these molecules.