A simplistic computational procedure for tunneling splittings caused by proton transfer

Hydrogen-Atom Tunneling in a Homochiral Environment

The role-exchanging concerted torsional motion of two hydrogen atoms in the homochiral dimer of trans-1,2-cyclohexanediol was characterized through a combination of broadband rotational spectroscopy and theoretical modeling. The results reveal that the concerted tunneling motion of the hydrogen atoms leads to the inversion of the sign of the dipole moment components along the a and b principal axes, due to the interchange motion that cooperatively breaks and reforms one intermolecular hydrogen bond. This motion is also coupled with two acceptor switching motions. The energy difference between the two ground vibrational states arising from this tunneling motion was determined to be 29.003(2) MHz. The corresponding wavefunctions suggest that the two hydrogen atoms are evenly delocalized on two equivalent potential wells, which differs from the heterochiral case where the hydrogen atoms are confined in separate wells, as the permutation-inversion symmetry breaks down. This intriguing contrast in hydrogen-atom behavior between homochiral and heterochiral environments could further illuminate our understanding of the role of chirality in intermolecular interactions and dynamics.

Study on depolymerization kinetics of formic acid dimers in binary mixture

In this study, polarization Raman spectra were collected for binary mixtures of formic acid/methanol and formic acid/acetonitrile with different volume fractions. The broad band of formic acid in the CO vibration region was divided into four vibration peaks, corresponding to CO symmetric and anti-symmetric stretching vibration from cyclic dimer, CO stretching from open dimer, and CO stretching from the free monomer. The experiments showed that as the volume fraction of formic acid in the binary mixture decreased, the cyclic dimer gradually converted to the open dimer, and at a volume fraction of 0.1, fully depolymerized into monomer form (free monomer, solvated monomer, and hydrogen bonding monomer clusters with solvent). The contribution percentage of the total CO stretching intensity of each structure at different concentrations was quantitatively calculated using high resolution infrared spectroscopy, and the results were consistent with the conclusions predicted by polarization Raman spectroscopy. Concentration-triggered 2D-COS synchronous and asynchronous spectra also confirmed the kinetics of formic acid diluted in acetonitrile. This work provides a spectroscopic method for studying the structure of organic compounds in solution and concentration-triggering kinetics in mixtures.

Metadynamics simulations with Bohmian-style bias potential

Here, we present a parametrization of the metadynamics simulations for reactions involving breaking the chemical bonds along a single collective variable coordinate. The parameterization is based on the similarity between the bias potential in metadynamics and the quantum potential in the de Broglie-Bohm formalism. We derive the method and test it on two prototypical reaction types: proton transfer and breaking of the cyclohexene cycle (reversed Diels-Alder reaction).

Rotational Spectroscopy of 2-Furoic Acid and Its Dimer: Conformational Distribution and Double Proton Tunneling

Structural and tunneling properties of the 2-furoic acid (FA) monomer and dimer were investigated using rotational spectroscopy and DFT calculations. CREST, a conformational ensemble space exploration tool, was used to identify all possible low energy conformations of the FA monomer and dimer, followed by the DFT geometry optimization and harmonic frequency calculations. Broadband rotational spectra in the 2-6 GHz and in the 8-12 GHz regions were recorded in a supersonic jet expansion. The monomeric FA was found to exist dominantly as three different conformers: I , II , and III in a jet, with I and II taking on the trans -COOH configuration while III having the cis -COOH configuration. For the FA dimer, only the I - II conformer was observed experimentally, whereas the symmetric I - I and II - II conformers were not observed because of their zero dipole moments. The analysis of the splittings in the rotational transitions of I - II allowed one to extract the tunneling splitting to be 1056.0(12) MHz. The barrier height was determined to be ~442 cm -1 using the scaled potential energy scans at several different levels of theory.

Internal rotation arena: Program performances on the low barrier problem of 4-methylacetophenone

In the rotational spectroscopy community, several popular codes are available to treat multiple internal rotors in a molecule. In terms of the pros and cons of each code, it is often a difficult task to decide which program to apply to a specific internal rotation problem. We faced this issue when dealing with the spectroscopic fingerprint of 4-methylacetophenone (4MAP), recently investigated in the microwave region, which we here extended into the millimeterwave region. The methyl group attached to the phenyl ring in 4MAP undergoes internal rotation with a very low barrier of only 22 cm^-1. The acetyl methyl group features a much higher barrier of about 580 cm^-1. The performances of a program using the so-called "local" approach in terms of Herschbach's perturbative treatment, SPFIT, as well as three programs XIAM, ERHAM, and ntop, representing "global" fits, were tested. The results aim at helping spectroscopists in the decision on how to tackle their own internal rotation problems.