Philipp Ottiger

Strong N-H pi Hydrogen Bonding in Amide-Benzene Interactions

Among the weak intermolecular interactions found in proteins, the amide N--H...pi interaction has been widely observed but remains poorly characterized as an individual interaction. We have investigated the isolated supersonic-jet-cooled dimer of the cis-amide and nucleobase analogue 2-pyridone (2PY) with benzene and benzene-d6. Both MP2 and SCS-MP2 geometry optimizations yield a T-shaped structure with a N--H...pi hydrogen bond to the benzene ring and the C=O group above, but far from the C--H bonds of benzene. The CCSD(T) calculated binding energy at the optimum geometry is De = 25.2 kJ/mol (dissociation energy D0 = 21.6 kJ/mol), corresponding to the H-bond strength of the water dimer or of N--H...O hydrogen bonds. The T-shaped geometry is supported by the infrared-ultraviolet depletion spectra of 2PY x benzene: The N--H stretch vibrational frequency is lowered by 56 cm(-1), and the C=O stretch vibration is lowered by 10 cm(-1), relative to those of bare 2PY, indicating a strong N--H...pi interaction and a weak interaction of the C=O group. The benzene C--H infrared stretches exhibit very small shifts (approximately 2 cm(-1)) relative to benzene, signaling the absence of interactions with the benzene C--H groups. The infrared spectral shifts are consistent with a strong nonconventional pi hydrogen bond and a T-shaped structure for 2PY x benzene. Symmetry-adapted perturbation theory calculations show that the N--H...pi interaction is by far the dominant stabilization factor.

Scope and limitations of the SCS-MP2 method for stacking and hydrogen bonding interactions

Fluorobenzenes are pi-acceptor synthons that form pi-stacked structures in molecular crystals as well as in artificial DNAs. We investigate the competition between hydrogen bonding and pi-stacking in dimers consisting of the nucleobase mimic 2-pyridone (2PY) and all fluorobenzenes from 1-fluorobenzene to hexafluorobenzene (n-FB, with n = 1-6). We contrast the results of high level ab initio calculations with those obtained using ultraviolet (UV) and infrared (IR) laser spectroscopy of isolated and supersonically cooled dimers. The 2PY.n-FB complexes with n = 1-5 prefer double hydrogen bonding over pi-stacking, as diagnosed from the UV absorption and IR laser depletion spectra, which both show features characteristic of doubly H-bonded complexes. The 2-pyridone.hexafluorobenzene dimer is the only pi-stacked dimer, exhibiting a homogeneously broadened UV spectrum and no IR bands characteristic for H-bonded species. MP2 (second-order Møller-Plesset perturbation theory) calculations overestimate the pi-stacked dimer binding energies by about 10 kJ/mol and disagree with the experimental observations. In contrast, the MP2 treatment of the H-bonded dimers appears to be quite accurate. Grimmes spin-component-scaled MP2 approach (SCS-MP2) is an improvement over MP2 for the pi-stacked dimers, reducing the binding energy by approximately 10 kJ/mol. When applied to explicitly correlated MP2 theory (SCS-MP2-R12 approach), agreement with the corresponding coupled-cluster binding energies [at the CCSD(T) level] is very good for the pi-stacked dimers, within +/- 1 kJ/mol for the 2PY complexes with 1-fluorobenzene, 1,2-difluorobenzene, 1,2,4,5-tetrafluorobenzene, pentafluorobenzene and hexafluorobenzene. Unfortunately, the SCS-MP2 approach also reduces the binding energy of the H-bonded species, leading to disagreement with both coupled-cluster theory and experiment. The SCS-MP2-R12 binding energies follow the SCS-MP2 binding energies closely, being about 0.5 and 0.7 kJ/mol larger for the H-bonded and pi-stacked forms, respectively, in an augmented correlation-consistent polarized valence quadruple-zeta basis. It seems that the SCS-MP2 and SCS-MP2-R12 methods cannot provide sufficient accuracy to replace the CCSD(T) method for intermolecular interactions where H-bonding and pi-stacking are competitive.

Vibrational quenching of excitonic splittings in H-bonded molecular dimers: The electronic Davydov splittings cannot match experiment

The S(1)/S(2) state exciton splittings of symmetric doubly hydrogen-bonded gas-phase dimers provide spectroscopic benchmarks for the excited-state electronic couplings between UV chromophores. These have important implications for electronic energy transfer in multichromophoric systems ranging from photosynthetic light-harvesting antennae to photosynthetic reaction centers, conjugated polymers, molecular crystals, and nucleic acids. We provide laser spectroscopic data on the S(1)/S(2) excitonic splitting Δ(exp) of the doubly H-bonded o-cyanophenol (oCP) dimer and compare to the splittings of the dimers of (2-aminopyridine)(2), [(2AP)(2)], (2-pyridone)(2), [(2PY)(2)], (benzoic acid)(2), [(BZA)(2)], and (benzonitrile)(2), [(BN)(2)]. The experimental S(1)/S(2) excitonic splittings are Δ(exp) = 16.4 cm(-1) for (oCP)(2), 11.5 cm(-1) for (2AP)(2), 43.5 cm(-1) for (2PY)(2), and <1 cm(-1) for (BZA)(2). In contrast, the vertical S(1)/S(2) energy gaps Δ(calc) calculated by the approximate second-order coupled cluster (CC2) method for the same dimers are 10-40 times larger than the Δ(exp) values. The qualitative failure of this and other ab initio methods to reproduce the exciton splitting Δ(exp) arises from the Born-Oppenheimer (BO) approximation, which implicitly assumes the strong-coupling case and cannot be employed to evaluate excitonic splittings of systems that are in the weak-coupling limit. Given typical H-bond distances and oscillator strengths, the majority of H-bonded dimers lie in the weak-coupling limit. In this case, the monomer electronic-vibrational coupling upon electronic excitation must be accounted for; the excitonic splittings arise between the vibronic (and not the electronic) transitions. The discrepancy between the BO-based splittings Δ(calc) and the much smaller experimental Δ(exp) values is resolved by taking into account the quenching of the BO splitting by the intramolecular vibronic coupling in the monomer S(1) ← S(0) excitation. The vibrational quenching factors Γ for the five dimers (oCP)(2), (2AP)(2), (2AP)(2), (BN)(2), and (BZA)(2) lie in the range Γ = 0.03-0.2. The quenched excitonic splittings Γ[middle dot]Δ(calc) are found to be in very good agreement with the observed splittings Δ(exp). The vibrational quenching approach predicts reliable Δ(exp) values for the investigated dimers, confirms the importance of vibrational quenching of the electronic Davydov splittings, and provides a sound basis for predicting realistic exciton splittings in multichromophoric systems.

Jet-Cooled 2-Aminopyridine Dimer: Conformers and Infrared Vibrational Spectra

The 2-aminopyridine dimer, (2AP)(2), is linked by two N-H...N hydrogen bonds, providing a model for the Watson-Crick configurations of the adenine or cytosine self-dimers. Structure optimization of (2AP)(2) at the MP2 level with the aug-cc-pVQZ basis set establishes the existence of two nearly degenerate conformers with C(i) and C(2) symmetry. Adding complete basis set extrapolation and DeltaCCSD(T) corrections gives binding energies D(e) = 10.70 and 10.72 kcal/mol, respectively. Both isomers are chiral, each giving rise to a pair of enantiomers. The potential energy surface of (2AP)(2) is calculated along the 2AP amino flip coordinates, revealing a 4-fold minimum low-energy region with a planar C(2h) symmetric and four asymmetric transition structures. The mass-selective resonant two-photon ionization (R2PI) spectra of supersonically cooled (2AP)(2) were remeasured. Three different species (A-C) were separated and characterized by UV/UV depletion spectroscopy and by infrared (IR) depletion spectroscopy in the 2600-3800 cm(-1) range. The R2PI and IR spectra of species A and B are very similar, in agreement with the prediction of two conformers of (2AP)(2). The IR bands are assigned to the H-bonded N-H(b) stretch, the N-H(2) bend overtone, and the free N-H(f) stretch of (2AP)(2), based on the calculated IR spectra, thereby extending and correcting previous assignments. Conformer A is tentatively assigned as the C(2) conformer. The UV spectrum of species C is very different from those of A and B, its IR spectrum exhibiting additional O-H stretching bands. C is assigned to the (2AP)(2).H(2)O cluster, based on the agreement of its IR spectrum with calculated IR spectra. Complete dissociation into the (2AP)(2)(+) ion occurs upon ionization.

Vibrational quenching of excitonic splittings in H-bonded molecular dimers: adiabatic description and effective mode approximation.

The quenching of the excitonic splitting in hydrogen-bonded molecular dimers has been explained recently in terms of exciton coupling theory, involving Försters degenerate perturbation theoretical approach [P. Ottiger, S. Leutwyler, and H. Köppel, J. Chem. Phys. 136, 174308 (2012)]. Here we provide an alternative explanation based on the properties of the adiabatic potential energy surfaces. In the proper limit, the lower of these surfaces exhibits a double-minimum shape, with an asymmetric distortion that destroys the geometric equivalence of the excitonically coupled monomers. An effective mode is introduced that exactly reproduces the energy gain and amount of distortion that occurs in a multi-dimensional normal coordinate space. This allows to describe the quenched exciton splitting as the energy difference of the two (S(1) and S(2)) vibronic band origins in a one-dimensional (rather than multi-dimensional) vibronic calculation. The agreement with the earlier result (based on Förster theory) is excellent for all five relevant cases studied. A simple rationale for the quenched exciton splitting as nonadiabatic tunneling splitting on the lower double-minimum potential energy surface is given.

Large-amplitude vibrations of an N-H center dot center dot center dot pi hydrogen bonded cis-amide-benzene complex

The ground-state N-H...pi interaction of 2-pyridone.benzene (2PY.Bz) has been studied by infrared-UV depletion spectroscopy of the supersonic-jet cooled complex [P. Ottiger et al., J. Phys. Chem. B (2009) 113, 2937]. Here, we investigate the large-amplitude vibrations of 2PY.Bz and its d(1)-2PY and benzene-d(6) isotopologues in the S(1) state, using two-color resonant two-photon ionization and UV-holeburning spectroscopies, complemented by RI-CC2 and SCS-RI-CC2 calculations of the S(1) state. The latter predict a tilted T-shaped structure with an N-H...pi hydrogen bond to the benzene ring, similar to the S(0) state. The binding energy is predicted to increase by 1.5 kJ mol(-1) upon S(1)<--S(0) excitation, in close agreement with the experimental value of 1.2 kJ mol(-1). The vibronic band structure up to 60 cm(-1) above the 0 band is dominated by large-amplitude delta tilting excitations, reflecting a change in the tilt angle of the T-shaped complex. The S(0) and S(1) state delta potentials were fitted to experiment, yielding a single minimum in the S(0) state and a double-minimum S(1) potential with delta(min) = +/-13 degrees. The second large-amplitude vibration is the theta twisting or benzene internal-rotation mode. Due to the C(6) symmetry of the benzene moiety the S(0) and S(1) state theta potentials are sixfold symmetric. Analysis of the theta band structure reveals that the S(0) and S(1)theta potentials are mutually aligned and that the internal rotation barriers are V(6)(S(0)) < 0.2 kJ mol(-1) and V(6)(S(1)) = 0.10(1) kJ mol(-1), in close agreement with the calculations. Weaker excitations of the totally symmetric intermolecular vibrations chi (shear), omega (bend) and sigma (stretch) vibrations are also observed. The 2PY intramolecular nu(1) overtone, corresponding to an 2PY amide out-of-plane twist distortion, lies approximately 30% higher than in bare 2PY, reflecting the hindrance of this motion by the strong N-H...pi interaction.

Excitonic splitting and coherent electronic energy transfer in the gas-phase benzoic acid dimer

The benzoic acid dimer, (BZA)(2), is a paradigmatic symmetric hydrogen bonded dimer with two strong antiparallel hydrogen bonds. The excitonic S(1)/S(2) state splitting and coherent electronic energy transfer within supersonically cooled (BZA)(2) and its (13)C-, d(1) -, d(2) -, and (13)C/d(1) - isotopomers have been investigated by mass-resolved two-color resonant two-photon ionization spectroscopy. The (BZA)(2)-(h - h) and (BZA)(2)-(d - d) dimers are C(2h) symmetric, hence only the S(2) ← S(0) transition can be observed, the S(1) ← S(0) transition being strictly electric-dipole forbidden. A single (12)C/(13)C or H/D isotopic substitution reduces the symmetry of the dimer to C(s), so that the isotopic heterodimers (BZA)(2) - (13)C, (BZA)(2) -(h - d), (BZA)(2) -(h(13)C-d), and (BZA)(2) -(h - d(13)C) show both S(1) ← S(0) and S(2) ← S(0) bands. The S(1)/S(2) exciton splitting inferred is Δ(exc) = 0.94 ± 0.1 cm(-1). This is the smallest splitting observed so far for any H-bonded gas-phase dimer. Additional isotope-dependent contributions to the splittings, Δ(iso), arise from the change of the zero-point vibrational energy upon electronic excitation and range from Δ(iso) = 3.3 cm(-1) upon (12)C/(13)C substitution to 14.8 cm(-1) for carboxy H/D substitution. The degree of excitonic localization/delocalization can be sensitively measured via the relative intensities of the S(1) ← S(0) and S(2) ← S(0) origin bands; near-complete localization is observed even for a single (12)C/(13)C substitution. The S(1)/ S(2) energy gap of (BZA)(2) is Δ(calc) (exc)=11 cm(-1) when calculated by the approximate second-order perturbation theory (CC2) method. Upon correction for vibronic quenching, this decreases to Δ(vibron) (exc)=2.1 cm(-1) [P. Ottiger et al., J. Chem. Phys. 136, 174308 (2012)], in good agreement with the observed Δ(exc) = 0.94 cm(-1). The observed excitonic splittings can be converted to exciton hopping times τ(exc). For the (BZA)(2)-(h - h) homodimer τ(exc) = 18 ps, which is nearly 40 times shorter than the double proton transfer time of (BZA)(2) in its excited state [Kalkman et al., ChemPhysChem 9, 1788 (2008)]. Thus, the electronic energy transfer is much faster than the proton-transfer in (BZA)(2)(∗).

Excitonic Splitting, Delocalization, and Vibronic Quenching in the Benzonitrile Dimer

The excitonic S1/S2 state splitting and the localization/delocalization of the S1 and S2 electronic states are investigated in the benzonitrile dimer (BN)2 and its (13)C and d5 isotopomers by mass-resolved two-color resonant two-photon ionization spectroscopy in a supersonic jet, complemented by calculations. The doubly hydrogen-bonded (BN-h5)2 and (BN-d5)2 dimers are C2h symmetric with equivalent BN moieties. Only the S0 → S2 electronic origin is observed, while the S0 → S1 excitonic component is electric-dipole forbidden. A single (12)C/(13)C or 5-fold h5/d5 isotopic substitution reduce the dimer symmetry to Cs, so that the heteroisotopic dimers (BN)2-(h5 – h5(13)C), (BN)2-(h5 – d5), and (BN)2-(h5 – h5(13)C) exhibit both S0 → S1 and S0 → S2 origins. Isotope-dependent contributions Δiso to the excitonic splittings arise from the changes of the BN monomer zero-point vibrational energies; these range from Δiso((12)C/(13)C) = 3.3 cm(–1) to Δiso(h5/d5) = 155.6 cm(–)1. The analysis of the experimental S1/S2 splittings of six different isotopomeric dimers yields the S1/S2 exciton splitting Δexc = 2.1 ± 0.1 cm(–1). Since Δiso(h5/d5) ≫ Δexc and Δiso((12)C/(13)C) > Δexc, complete and near-complete exciton localization occurs upon (12)C/(13)C and h5/d5 substitutions, respectively, as diagnosed by the relative S0 → S1 and S0 → S2 origin band intensities. The S1/S2 electronic energy gap of (BN)2 calculated by the spin-component scaled approximate second-order coupled-cluster (SCS-CC2) method is Δel(calc) = 10 cm(–1). This electronic splitting is reduced by the vibronic quenching factor Γ. The vibronically quenched exciton splitting Δel(calc)·Γ = Δvibron(calc) = 2.13 cm(–1) is in excellent agreement with the observed splitting Δexc = 2.1 cm(–1). The excitonic splittings can be converted to semiclassical exciton hopping times; the shortest hopping time is 8 ps for the homodimer (BN-h5)2, the longest is 600 ps for the (BN)2(h5 – d5) heterodimer.

Structure and intermolecular vibrations of perylene·trans-1,2-dichloroethene, a weak charge-transfer complex.

The vibronic spectra of strong charge-transfer complexes are often congested or diffuse and therefore difficult to analyze. We present the spectra of the π-stacked complex perylene trans-1,2-dichloroethene, which is in the limit of weak charge transfer, the electronic excitation remaining largely confined to the perylene moiety. The complex is formed in a supersonic jet, and its S0 ↔ S1 spectra are investigated by two-color resonant two-photon ionization (2C-R2PI) and fluorescence spectroscopies. Under optimized conditions, vibrationally cold (T(vib) ≈ 9 K) and well resolved spectra are obtained. These are dominated by vibrational progressions in the “hindered-rotation” Rc intermolecular vibration with very low frequencies of 11 (S0) and 13 cm(–1) (S1). The intermolecular Tz stretch and the Ra and Rb bend vibrations are also observed. The normally symmetry-forbidden intramolecular 1a(u) “twisting” vibration of perylene also appears, showing that the π- stacking interaction deforms the perylene moiety, lowering its local symmetry from D2h to D2. We calculate the structure and vibrations of this complex using six different density functional theory (DFT) methods (CAM-B3LYP, BH&HLYP, B97-D3, ωB97X-D, M06, and M06-2X) and compare the results to those calculated by correlated wave function methods (SCS-MP2 and SCS-CC2). The structures and vibrational frequencies predicted with the CAM-B3LYP and BH&HLYP methods disagree with the other calculations and with experiment. The other four DFT and the ab initio methods all predict a π-stacked “centered” structure with nearly coplanar perylene and dichloroethene moieties and intermolecular binding energies of D(e) = −20.8 to −26.1 kJ/mol. The 000 band of the S0 → S1 transition is red-shifted by δν = −301 cm(–1) relative to that of perylene, implying that the D(e) increases by 3.6 kJ/mol or 15% upon electronic excitation. The intermolecular vibrational frequencies are assigned to the calculated Rc, Tz, Ra, and Rb vibrations by comparing to the observed/calculated frequencies and S0 ↔ S1 Franck–Condon factors. Of the three TD-DFT methods tested, the hybrid-meta-GGA functional M06-2X shows the best agreement with the experimental electronic transition energies, spectral shifts, and vibronic spectra, closely followed by the ωB97X-D functional, while the M06 functional gives inferior results.

Modeling the Histidine–Phenylalanine Interaction: The NH···π Hydrogen Bond of Imidazole·Benzene

NH···π hydrogen bonds occur frequently between the amino acid side groups in proteins and peptides. Data-mining studies of protein crystals find that ∼80% of the T-shaped histidine···aromatic contacts are CH···π, and only ∼20% are NH···π interactions. We investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet-cooled imidazole·benzene (Im·Bz) complex as a model for the NH···π interaction between histidine and phenylalanine. Ground- and excited-state dispersion-corrected density functional calculations and correlated methods (SCS-MP2 and SCS-CC2) predict that Im·Bz has a Cs-symmetric T-shaped minimum-energy structure with an NH···π hydrogen bond to the Bz ring; the NH bond is tilted 12° away from the Bz C6 axis. IR depletion spectra support the T-shaped geometry: The NH stretch vibrational fundamental is red shifted by -73 cm(-1) relative to that of bare imidazole at 3518 cm(-1), indicating a moderately strong NH···π interaction. While the S0(A1g) → S1(B2u) origin of benzene at 38 086 cm(–1) is forbidden in the gas phase, Im·Bz exhibits a moderately intense S0 → S1 origin, which appears via the D(6h) → Cs symmetry lowering of Bz by its interaction with imidazole. The NH···π ground-state hydrogen bond is strong, De=22.7 kJ/mol (1899 cm–1). The combination of gas-phase UV and IR spectra confirms the theoretical predictions that the optimum Im·Bz geometry is T shaped and NH···π hydrogen bonded. We find no experimental evidence for a CH···π hydrogen-bonded ground-state isomer of Im·Bz. The optimum NH···π geometry of the Im·Bz complex is very different from the majority of the histidine·aromatic contact geometries found in protein database analyses, implying that the CH···π contacts observed in these searches do not arise from favorable binding interactions but merely from protein side-chain folding and crystal-packing constraints. The UV and IR spectra of the imidazole·(benzene)2 cluster are observed via fragmentation into the Im·Bz+ mass channel. The spectra of Im·Bz and Im·Bz2 are cleanly separable by IR hole burning. The UV spectrum of Im·Bz2 exhibits two 000 bands corresponding to the S0 → S1 excitations of the two inequivalent benzenes, which are symmetrically shifted by -86/+88 cm(-1) relative to the 000 band of benzene