Solvation Effects on Quantum Tunneling Reactions

A decisive factor for obtaining high yields and selectivities in organic synthesis is the choice of the proper solvent. Solvent selection is often guided by the intuitive understanding of transition state-solvent interactions. However, quantum-mechanical tunneling can significantly contribute to chemical reactions, circumventing the transition state and thus depriving chemists of their intuitive handle on the reaction kinetics. In this Account, we aim to provide rationales for the effects of solvation on tunneling reactions derived from experiments performed in cryogenic matrices.The tunneling reactions analyzed here cover a broad range of prototypical organic transformations that are subject to strong solvation effects. Examples are the hydrogen tunneling probability for the cis-trans isomerization of formic acid which is strongly reduced upon formation of hydrogen-bonded complexes and the [1,2]H-shift in methylhydroxycarbene where a change in product selectivity is predicted upon interaction with hydrogen bond acceptors.Not only hydrogen but also heavy atom tunneling can exhibit strong solvent effects. The direction of the nearly degenerate valence tautomerization between benzene oxide and oxepin was found to reverse upon formation of a halogen or hydrogen bond with ICF₃ or H₂O. But even in the absence of strong noncovalent interactions such as hydrogen or halogen bonding, solvation can have a decisive effect on tunneling as evidenced by the Cope rearrangement of semibullvalenes via heavy-atom tunneling. Can quantum tunneling be catalyzed? The acceleration of the ring expansion of 1H-bicyclo[3.1.0.]-hexa-3,5-dien-2-one by complexation with Lewis acids provides a proof-of-concept for tunneling catalysis.Two concepts are central for the explanation and prediction of solvation effects on tunneling phenomena: a simple approach expands the Born-Oppenheimer approximation by separating nuclear degrees of freedom into intra- and intermolecular degrees. Intermolecular movements represent the slowest motions within molecular aggregates, thus effectively freezing the position of the solvent in relation to the reactant during the tunneling process. Another useful approach is to treat reactants and products by separate single-well potentials, where the intersection represents the transition state. Thus, stabilization of the reactants via solvation should result in an increase in barrier heights and widths which in turn lowers tunneling probabilities. These simple models can predict trends in tunneling kinetics and provide a rational basis for controlling tunneling reactions via solvation.

Optical van-der-Waals forces in molecules: from electronic Bethe-Salpeter calculations to the many-body dispersion model

Molecular forces induced by optical excitations are connected to a wide range of phenomena, from chemical bond dissociation to intricate biological processes that underpin vision. Commonly, the description of optical excitations requires the solution of computationally demanding electronic Bethe-Salpeter equation (BSE). However, when studying non-covalent interactions in large-scale systems, more efficient methods are desirable. Here we introduce an effective approach based on coupled quantum Drude oscillators (cQDO) as represented by the many-body dispersion model. We find that the cQDO Hamiltonian yields semi-quantitative agreement with BSE calculations and that both attractive and repulsive optical van der Waals (vdW) forces can be induced by light. These optical-vdW interactions dominate over vdW dispersion in the long-distance regime, showing a complexity that grows with system size. Evidence of highly non-local forces in the human formaldehyde dehydrogenase 1MC5 protein suggests the ability to selectively activate collective molecular vibrations by photoabsorption, in agreement with recent experiments.

CC-stretched formic acid: isomerisation, dimerisation, and carboxylic acid complexation

The cis-trans-isomerism of the propiolic acid monomer (HC[triple bond, length as m-dash]C-COOH) is examined with linear Raman jet spectroscopy, yielding the first environment-free vibrational band centres of a higher-energy cis-rotamer beyond formic acid (HCOOH) in addition to all fundamentals and a large number of hot and combination/overtone bands of the trans-conformer. Two near-isoenergetic trans-fundamentals of different symmetry (CC[double bond, length as m-dash]O bend and OH torsion) prove to be a sensitive benchmarking target, as their energetic order is susceptible to the choice of electronic structure method, basis set size, and inclusion of vibrational anharmonicity. For the infrared- and Raman-active C[double bond, length as m-dash]O stretching fundamentals of the cyclic (C_2h) trans-propiolic acid dimer, resonance couplings are found that in part extend to the C_s-symmetric heterodimer of trans-propiolic and trans-formic acid. Exploratory vibrational perturbation theory (VPT2) calculations show that all perturbing states involve displacements of the OH moieties located on the doubly hydrogen bonded ring. The comparison of the infrared spectra of the propiolic acid dimer and its heterodimer with formic acid to that of several other carboxylic acid dimers from the literature reveals a notable similarity regarding a non-fundamental dimer band around 1800 cm^-1, which in most cases is so far unassigned. VPT2 calculations and a simple harmonic model suggest an assignment to a combination vibration of the symmetric and antisymmetric OH torsion.

Intricate Conformational Tunneling in Carbonic Acid Monomethyl Ester

Disentangling internal and external effects is a key requirement for understanding conformational tunneling processes. Here we report the s- trans/ s- cis tunneling rotamerization of carbonic acid monomethyl ester (1) under matrix isolation conditions and make comparisons to its parent carbonic acid (3). The observed tunneling rate of 1 is temperature-independent in the 3-20 K range and accelerates when using argon instead of neon as the matrix material. The methyl group increases the effective half life (τ_eff) of the energetically disfavored s- trans-conformer from 3-5 h for 3 to 11-13 h for 1. Methyl group deuteration slows the rotamerization further (τ_eff ≈ 35 h). CCSD(T)/cc-pVQZ//MP2/aug-cc-pVTZ computations of the tunneling probability suggest that the rate should be almost unaffected by methyl substitution or its deuteration. Thus the observed relative rates are puzzling, and they disagree with previous explanations involving fast vibrational relaxation after the tunneling event facilitated by the alkyl rotor.

Nuclear dynamics and phase polymorphism in solid formic acid

We apply a unique sequence of structural and dynamical neutron-scattering techniques, augmented with density-functional electronic-structure calculations, to establish the degree of polymorphism in an archetypal hydrogen-bonded system - crystalline formic acid. Using this combination of experimental and theoretical techniques, the hypothesis by Zelsmann on the coexistence of the β₁ and β₂ phases above 220 K is tested. Contrary to the postulated scenario of proton-transfer-driven phase coexistence, the emerging picture is one of a quantitatively different structural change over this temperature range, whereby the loosening of crystal packing promotes temperature-induced shearing of the hydrogen-bonded chains. The presented work, therefore, solves a fifty-year-old puzzle and provides a suitable framework for the use neutron-Compton-scattering techniques in the exploration of phase polymorphism in condensed matter.

Acetic acid-water complex: The first observation of structures containing the higher-energy acetic acid conformer

Non-covalent interaction of acetic acid (AA) and water is studied experimentally by IR spectroscopy in a nitrogen matrix and theoretically at the MP2 and coupled-cluster with single and double and perturbative triple excitations [CCSD(T)]/6-311++G(2d,2p) levels of theory. This work is focused on the first preparation and characterization of complexes of higher-energy (cis) conformer of AA with water. The calculations show three 1:1 structures for the trans-AA⋯H2O complexes and three 1:1 structures for the cis-AA⋯H2O complexes. Two trans-AA⋯H2O and two cis-AA⋯H2O complexes are found and structurally assigned in the experiments. The two cis-AA⋯ ⋅ H2O complexes are obtained by annealing of a matrix containing water and cis-AA molecules prepared by selective vibrational excitation of the ground-state trans form. The less stable trans-AA⋯H2O complex is obtained by vibrational excitation of the less stable cis-AA⋯H2O complex. In addition, the 1:2 complexes of trans-AA and cis-AA with water molecules are studied computationally and the most stable forms of the 1:2 complexes are experimentally identified.

Theoretical Study of Vibrationally Averaged Dipole Moments for the Ground and Excited C═O Stretching States of trans-Formic Acid

Recent experimental studies of trans-formic acid (FA) in solid para-hydrogen (pH2) highlighted the importance of vibrationally averaged dipole moments for the interpretation of the high-resolution infrared (IR) spectra, in particular for the C═O stretch (ν3) mode. In this report, dipole moments for the ν3 ground (v = 0) and excited (v = 1, 2, 3, and 4) anharmonic vibrational states in trans-FA are investigated using two different approaches: a single mode approximation, where the vibrational states are obtained from the solution of the one-dimensional Schrödinger equation for the harmonic normal coordinate, and a limited vibrational configuration interaction (VCI) approximation. Density functional theory (B3LYP, BPW91) and correlated ab initio (MP2 and CCSD(T)) electronic methods were employed with a number of double- and triple-ζ and correlation consistent basis sets. Both single mode and VCI approaches show comparable agreement with experimental data, which is more dependent on the level of theory used. In particular, the BPW91/cc-pVDZ level appears to perform remarkably well. Effects of solvation of FA in solid state Ar and pH2 matrices were simulated at the BPW91/cc-pVDZ level using a conductor-like polarized continuum model (CPCM). The Ar and pH2 solid-state matrices cause quite a substantial increase in the FA dipole moments. Compared to gas-phase calculations, the CPCM model for pH2 better reproduces the experimental FA spectral shifts caused by interaction with traces of ortho-hydrogen (oH2) species in solid pH2. The validity of the single mode approach is tested against the multidimensional VCI results, suggesting that the isolated (noninteracting) mode approximation is valid up to the third vibrationally excited state (v = 3). Finally, the contribution of the ground anharmonic vibrational states of the remaining modes to the resulting ν3 single mode dipole moments is examined and discussed.