Search

The decay mechanism of S 0 → S 1 excited cytosine (Cyt) and the effect of substitution are studied combining jet-cooled spectroscopy (nanosecond resonant two-photon ionization (R2PI) and picosecond lifetime measurements) with CASPT2//CASSCF computations for eight derivatives. For Cyt and five derivatives substituted at N1, C5, and C6, rapid internal conversion sets in at 250-1200 cm -1 above the 0 0 0 bands. The break-off in the spectra correlates with the calculated barriers toward the "C5-C6 twist" conical intersection, which unambiguously establishes the decay mechanism at low S 1 state vibrational energies. The barriers increase with substituents that stabilize the charge shifts at C4, C5, and C6 following ( 1 ππ *) excitation. The R2PI spectra of the clamped derivatives 5,6-trimethyleneCyt (TMCyt) and 1-methyl-TMCyt (1M-TMCyt), which decay along an N3 out-of-plane coordinate, extend up to +3500 and +4500 cm -1 .

The gas-phase ground-state dissociation energy D 0 ( S 0 ) of an isolated and cold bimolecular complex is a fundamental measure of the intermolecular interaction strength between its constituents. Accurate D 0 values are important for the understanding of intermolecular bonding, for benchmarking high-level theoretical calculations, and for the parameterization of dispersion-corrected density functionals or force-field models that are used in fields ranging from crystallography to biochemistry. We review experimental measurements of the gas-phase D 0 ( S 0 ) and D 0 ( S 1 ) values of 55 different M⋅S complexes, where M is a (hetero)aromatic molecule and S is a closed-shell solvent atom or molecule. The experiments employ the triply resonant SEP-R2PI laser method, which involves M-centered ( S 0 → S 1 ) electronic excitation, followed by S 1 → S 0 stimulated emission spanning a range of S 0 state vibrational levels. At sufficiently high vibrational energy, vibrational predissociation of the M⋅S complex occurs. A total of 49 dissociation energies were bracketed to within ≤1.0 kJ/mol, providing a large experimental database of accurate noncovalent interactions.

The ground-state intermolecular dissociation energies D 0 (S 0 ) of supersonic-jet cooled intermolecular complexes of 1-naphthol (1NpOH) with the bi- and tricycloalkanes trans-decalin, cis-decalin, and adamantane were measured using the stimulated-emission-pumping/resonant two-photon ionization (SEP-R2PI) method. Using UV/UV holeburning, we identified two isomers (A and B) of the adamantane and trans-decalin complexes and four isomers (A-D) of the cis-decalin complex. For 1NpOH·adamantane A and B, the D 0 (S 0 ) values are 21.6 ± 0.15 kJ/mol and 21.2 ± 0.32 kJ/mol, those of 1NpOH·trans-decalin A and B are 28.7 ± 0.3 kJ/mol and 28.1 ± 0.9 kJ/mol, and those of 1NpOH·cis-decalin A and B are 28.9 ± 0.15 kJ/mol and 28.7 ± 0.3 kJ/mol. Upon S 0 → S 1 electronic excitation of the 1NpOH moiety, the dissociation energies of adamantane, trans-decalin, and the cis-decalin isomer C change by <1% and those of cis-decalin isomers A, B, and D increase only slightly (1%-3%). This implies that the hydrocarbons are dispersively adsorbed to a naphthalene "face." Calculations using the dispersion-corrected density functional theory methods B97-D3 and B3LYP-D3 indeed predict that the stable structures have face geometries. The B97-D3 calculated D 0 (S 0 ) values are within 1 kJ/mol of the experiment, while B3LYP-D3 predicts D 0 values that are 1.4-3.3 kJ/mol larger. Although adamantane has been recommended as a "dispersion-energy donor," the binding energies of the trans- and cis-decalin adducts to 1NpOH are 30% larger than that of adamantane. In fact, the D 0 value of 1NpOH·adamantane is close to that of 1NpOH·cyclohexane, reflecting the nearly identical contact layer between the two molecules.

The S 1 excited-state lifetime of jet-cooled 1-ethylcytosine (1ECyt) is ∼1 ns, one of the longest lifetimes for cytosine derivatives to date. Here, we analyze its S 0 → S 1 vibronic spectrum using two-color resonant two-photon ionization and UV/UV holeburning spectroscopy. Compared to cytosine and 1-methylcytosine, the S 0 → S 1 spectrum of 1ECyt shows a progression in the out-of-plane "butterfly" mode ν 1 ' , identified by spin-component scaled-second-order coupled-cluster method ab initio calculations. We also report time-resolved S 1 state nonradiative dynamics at ∼20 ps resolution by the pump/delayed ionization technique. The S 1 lifetime increases with the number of ν 1 ' quanta from τ = 930 ps at v 1 ' =0 to 1030 ps at v 1 ' =2, decreasing to 14 ps at 710 cm -1 vibrational energy. We measured the rate constants for S 1 ⇝ S 0 internal conversion and S 1 ⇝ T 1 intersystem crossing (ISC): At the v' = 0 level, k IC is 8 × 10 8 s -1 or three times smaller than 1-methylcytosine. The ISC rate constant from v' = 0 to the T 1 ( 3 ππ * ) state is k ISC = 2.4 × 10 8 s -1 , 10 times smaller than the ISC rate constants of cytosine, but similar to that of 1-methylcytosine. Based on the calculated S 1 ( 1 ππ * ) state radiative lifetime τ rad = 12 ns, the fluorescence quantum yield of 1ECyt is Φ fl ∼ 7% and the intersystem crossing yield is Φ ISC ∼ 20%. We measured the adiabatic ionization energy of 1-ethylcytosine via excitation of the S 1 state as 8.353 ± 0.008 eV, which is 0.38 eV lower than that of amino-keto cytosine. Measurement of the ionization energy of the long-lived T 1 (ππ * ) state formed via ISC reveals that it lies 3.2-3.4 eV above the S 0 ground state.

The stimulated-emission-pumping/resonant 2-photon ionization (SEP-R2PI) method was used to determine the intermolecular dissociation energies D 0 of jet-cooled 1-naphthol(1NpOH)·S complexes, where S is a linear molecule (N 2 , CO, CO 2 , OCS, N 2 O, and ethyne) or symmetric-top molecule (2-butyne) that contains double or triple bonds. The dissociation energies D 0 (S 0 ) are bracketed as follows: 6.68 ± 0.08 kJ/mol for S=N 2 , 7.7 ± 0.8 kJ/mol for CO, 12.07 ± 0.10 kJ/mol for CO 2 , 13.03 ± 0.01 kJ/mol for N 2 O, 14.34 ± 0.08 kJ/mol for ethyne, 15.0 ± 1.35 kJ/mol for OCS, and 29.6 ± 2.4 kJ/mol for 2-butyne. The minimum-energy structures, vibrational wavenumbers, and zero-point vibrational energies were calculated using the dispersion-corrected density functional theory methods such as B97-D3 and B3LYP-D3 with the def2-QZVPP basis set. These predict that N 2 and CO are dispersively bound Face complexes (S bound to a naphthalene Face), while CO 2 , N 2 O, and OCS adsorb into the "Notch" between the naphthyl and OH groups; these are denoted as Notch complexes. Ethyne and 2-butyne form Edge complexes involving H-bonds from the -OH group of 1NpOH to the center of the molecule. The presence of a double or triple bond or an aromatic C=C bond within S does not lead to a specific calculated geometry (Face, Notch or Edge). However, a correlation exists between the structure and the sign of the quadrupole moment component Θ zz of S: negative Θ zz correlates with Face or Notch, while positive Θ zz correlates with Edge geometries.

Mass-selective two-color resonant two-photon ionization (2C-R2PI), UV/UV hole-burning, and infrared (IR) depletion spectra of supersonic jet-cooled 2-pyridone·(H 2 O) n clusters with n = 1-4 have been measured to investigate the local hydration patterns around 2-pyridone (2PY) as a function of cluster size. As shown by others, the IR frequencies of the OH and NH stretches of the n = 1, 2 clusters are characteristic of water wires stretching from the NH to the C═O group of 2PY. We identify two isomers (3A and 3B) of the n = 3 cluster in the 2C-R2PI spectrum and separate them by IR/UV and UV/UV hole-burning techniques. Isomer 3A exhibits a three-membered water wire, extending the n = 1, 2 structural motif. Isomer 3B exhibits bifurcated water wires with the first H 2 O donating to two waters that form H-bonds to the C═O group. This increases the H-bond strength between the NH group of 2PY and the proximal H 2 O molecule, lowering the NH stretch to ∼2800 cm -1 . The n = 4 cluster is also bifurcated with two water wires between the bifurcating H 2 O and the C═O group. The cluster-selective IR spectra are complemented with density-functional calculations using the PW91, B3LYP, B97-D, and M06-2X functionals, where the latter two include long-range dispersive interactions, and with the ab initio correlated SCS-CC2 method. The calculated IR spectra provide firm assignments of the structures of the n = 1-4 cluster structures and allow us to understand the evolution of individual H-bond strengths with increasing cluster size.

We have measured the intermolecular dissociation energies D 0 of supersonically cooled 1-naphthol (1NpOH) complexes with solvents S = furan, thiophene, 2,5-dimethylfuran, and tetrahydrofuran. The naphthol OH forms non-classical H-bonds with the aromatic π -electrons of furan, thiophene, and 2,5-dimethylfuran and a classical H-bond with the tetrahydrofuran O atom. Using the stimulated-emission pumping resonant two-photon ionization method, the ground-state D 0 ( S 0 ) values were bracketed as 21.8 ± 0.3 kJ/mol for furan, 26.6 ± 0.6 kJ/mol for thiophene, 36.5 ± 2.3 kJ/mol for 2,5-dimethylfuran, and 37.6 ± 1.3 kJ/mol for tetrahydrofuran. The dispersion-corrected density functional theory methods B97-D3, B3LYP-D3 (using the def2-TZVPP basis set), and ω B97X-D [using the 6-311++G(d,p) basis set] predict that the H-bonded (edge) isomers are more stable than the face isomers bound by dispersion; experimentally, we only observe edge isomers. We compare the calculated and experimental D 0 values and extend the comparison to the previously measured 1NpOH complexes with cyclopropane, benzene, water, alcohols, and cyclic ethers. The dissociation energies of the nonclassically H-bonded complexes increase roughly linearly with the average polarizability of the solvent, α ¯ (S). By contrast, the D 0 values of the classically H-bonded complexes are larger, increase more rapidly at low α ¯ (S), but saturate for large α ¯ (S). The calculated D 0 ( S 0 ) values for the cyclopropane, benzene, furan, and tetrahydrofuran complexes agree with experiment to within 1 kJ/mol and those of thiophene and 2,5-dimethylfuran are ∼3 kJ/mol smaller than experiment. The B3LYP-D3 calculated D 0 values exhibit the lowest mean absolute deviation (MAD) relative to experiment (MAD = 1.7 kJ/mol), and the B97-D3 and ω B97X-D MADs are 2.2 and 2.6 kJ/mol, respectively.

Using the stimulated-emission-pumping/resonant 2-photon ionization (SEP-R2PI) method, we have determined accurate intermolecular dissociation energies D 0 of supersonic jet-cooled intermolecular complexes of 1-naphthol (1NpOH) with alkanes, 1NpOH·S, with S = methane, ethane, propane, and n-butane. Experimentally, the smaller alkanes form a single minimum-energy structure, while 1-naphthol·n-butane forms three different isomers. The ground-state dissociation energies D 0 (S 0 ) for the complexes with propane and n-butane (isomers A and B) were bracketed within ±0.5%, being 16.71 ± 0.08 kJ/mol for S = propane and 20.5 ± 0.1 kJ/mol for isomer A and 20.2 ± 0.1 kJ/mol for isomer B of n-butane. All 1NpOH·S complexes measured previously exhibit a clear dissociation threshold in their hot-band detected SEP-R2PI spectra, but weak SEP-R2PI bands are observed above the putative dissociation onset for the methane and ethane complexes. We attribute these bands to long-lived complexes that retain energy in rotation-type intermolecular vibrations, which couple only weakly to the dissociation coordinates. Accounting for this, we find dissociation energies of D 0 (S 0 ) = 7.98 ± 0.55 kJ/mol (±7%) for S = methane and 14.5 ± 0.28 kJ/mol (±2%) for S = ethane. The D 0 values increase by only 1% upon S 0 → S 1 excitation of 1-naphthol. The dispersion-corrected density functional theory methods B97-D3, B3LYP-D3, and ωB97X-D predict that the n-alkanes bind dispersively to the naphthalene "Face." The assignment of the complexes to Face structures is supported by the small spectral shifts of the S 0 → S 1 electronic origins, which range from +0.5 to -15 cm -1 . Agreement with the calculated dissociation energies D 0 (S 0 ) is quite uneven, the B97-D3 values agree within 5% for propane and n-butane, but differ by up to 20% for methane and ethane. The ωB97X-D method shows good agreement for methane and ethane but overestimates the D 0 (S 0 ) values for the larger n-alkanes by up to 20%. The agreement of the B3LYP-D3 D 0 values is intermediate between the other two methods.

We measured accurate intermolecular dissociation energies D 0 of the supersonic jet-cooled complexes of 1-naphthol (1NpOH) with the noble gases Ne, Ar, Kr, and Xe and with N 2 , using the stimulated-emission pumping resonant two-photon ionization method. The ground-state values D 0 (S 0 ) for the 1NpOH⋅S complexes with S= Ar, Kr, Xe, and N 2 were bracketed to be within ±3.5%; they are 5.67 ± 0.05 kJ/mol for S = Ar, 7.34 ± 0.07 kJ/mol for S = Kr, 10.8 ± 0.28 kJ/mol for S = Xe, 6.67 ± 0.08 kJ/mol for isomer 1 of the 1NpOH⋅N 2 complex, and 6.62 ± 0.22 kJ/mol for the corresponding isomer 2. For S = Ne, the upper limit is D 0 < 3.36 kJ/mol. The dissociation energies increase by 1%-5% upon S 0 → S 1 excitation of the complexes. Three dispersion-corrected density functional theory (DFT-D) methods (B97-D3, B3LYP-D3, and ωB97X-D) predict that the most stable form of these complexes involves dispersive binding to the naphthalene "face." A more weakly bound edge isomer is predicted in which the S moiety is H-bonded to the OH group of 1NpOH; however, no edge isomers were observed experimentally. The B97-D3 calculated dissociation energies D 0 (S 0 ) of the face complexes with Ar, Kr, and N 2 agree with the experimental values within <5%, but the D 0 (S 0 ) for Xe is 12% too low. The B3LYP-D3 and ωB97X-D calculated D 0 (S 0 ) values exhibit larger deviations to both larger and smaller dissociation energies. For comparison to 1-naphthol, we calculated the D 0 (S 0 ) of the carbazole complexes with S = Ne, Ar, Kr, Xe, and N 2 using the same DFT-D methods. The respective experimental values have been previously determined to be within <2%. Again, the B97-D3 results are in the best overall agreement with experiment.

p-Dioxane is non-polar, hence its rotational constants cannot be determined by microwave rotational coherence spectroscopy (RCS). We perform high-resolution gas-phase rotational spectroscopy of para-dioxane-h 8 and -d 8 using femtosecond time-resolved Raman RCS in a gas cell at T = 293 K and in a pulsed supersonic jet at T∼130 K. The inertial tensor of p-dioxane-h 8 is strongly asymmetric, leading to a large number of asymmetry transients in its RCS spectrum. In contrast, the d 8 -isotopomer is a near-oblate symmetric top that exhibits a much more regular RCS spectrum with few asymmetry transients. Fitting the fs Raman RCS transients of p-dioxane-h 8 to an asymmetric-top model yields the ground-state rotational constants A 0 = 5084.4(5) MHz, B 0 = 4684(1) MHz, C 0 = 2744.7(8) MHz, and (A 0 + B 0 )/2 = 4884.5(7) MHz (±1σ). The analogous values for p-dioxane-d 8 are A 0 = 4083(2) MHz, B 0 = 3925(4) MHz, C 0 = 2347.1(6) MHz, and (A 0 + B 0 )/2 = 4002.4(6) MHz. We determine the molecular structure with a semi-experimental approach involving the highly correlated coupled-cluster singles, doubles and iterated triples method and the cc-pCVXZ basis set series from double- to quadruple-zeta (X = D, T, Q). Combining the calculated vibrationally averaged rotational constants A 0 calc (X),B 0 calc (X),C 0 calc (X) for increasing basis-set size X with non-linear extrapolation to the experimental constants A 0 exp ,B 0 exp ,C 0 exp allows to determine the equilibrium ground state structure of p-dioxane. For instance, the equilibrium C-C and C-O bond lengths are r e (CC) = 1.5135(3) Å and r e (CO) = 1.4168(4) Å, and the four axial C-H bond lengths are 0.008 Å longer than the four equatorial ones. The latter is ascribed to the trans-effect (anomeric effect), i.e., the partial delocalization of the electron lone-pairs on the O atoms that are oriented trans, relative to the axial CH bonds.

We measure the S 0  → S 1 spectrum and time-resolved S 1 state nonradiative dynamics of the "clamped" cytosine derivative 5,6-trimethylenecytosine (TMCyt) in a supersonic jet, using two-color resonant two-photon ionization (R2PI), UV/UV holeburning, and ns time-resolved pump/delayed ionization. The experiments are complemented with spin-component scaled second-order approximate coupled cluster (SCS-CC2), time-dependent density functional theory, and multi-state second-order perturbation-theory (MS-CASPT2) ab initio calculations. While the R2PI spectrum of cytosine breaks off ∼500 cm -1 above its 0 0 0 band, that of TMCyt extends up to +4400 cm -1 higher, with over a hundred resolved vibronic bands. Thus, clamping the cytosine C 5 -C 6 bond allows us to explore the S 1 state vibrations and S 0  → S 1 geometry changes in detail. The TMCyt S 1 state out-of-plane vibrations ν 1 ' , ν 3 ' , and ν 5 ' lie below 420 cm -1 , and the in-plane ν 11 ' , ν 12 ' , and ν 23 ' vibrational fundamentals appear at 450, 470, and 944 cm -1 . S 0   →  S 1 vibronic simulations based on SCS-CC2 calculations agree well with experiment if the calculated ν 1 ' , ν 3 ' , and ν 5 ' frequencies are reduced by a factor of 2-3. MS-CASPT2 calculations predict that the ethylene-type S 1  ⇝ S 0 conical intersection (CI) increases from +366 cm -1 in cytosine to >6000 cm -1 in TMCyt, explaining the long lifetime and extended S 0  → S 1 spectrum. The lowest-energy S 1  ⇝ S 0 CI of TMCyt is the "amino out-of-plane" (OP X ) intersection, calculated at +4190 cm -1 . The experimental S 1  ⇝ S 0 internal conversion rate constant at the S 1 (v ' =0) level is k IC =0.98-2.2⋅10 8 s -1 , which is ∼10 times smaller than in 1-methylcytosine and cytosine. The S 1 (v ' =0) level relaxes into the T 1 ( 3 ππ*) state by intersystem crossing with k ISC =0.41-1.6⋅10 8 s -1 . The T 1 state energy is measured to lie 24 580±560 cm -1 above the S 0 state. The S 1 (v ' =0) lifetime is τ=2.9 ns, resulting in an estimated fluorescence quantum yield of Φ fl =24%. Intense two-color R2PI spectra of the TMCyt amino-enol tautomers appear above 36 000 cm -1 . A sharp S 1 ionization threshold is observed for amino-keto TMCyt, yielding an adiabatic ionization energy of 8.114±0.002 eV.

The ground-state dissociation energy, D0(S0), of isolated intermolecular complexes in the gas phase is a fundamental measure of the interaction strength between the molecules. We have developed a three-laser, triply resonant pump-dump-probe technique to measure dissociation energies of jet-cooled M•S complexes, where M is an aromatic chromophore and S is a closed-shell 'solvent' molecule. Stimulated emission pumping (SEP) via the S0→S1 electronic transition is used to precisely 'warm' the complex by populating high vibrational levels v" of the S0 state. If the deposited energy E(v") is less than D0(S0), the complex remains intact, and is then mass- and isomer-selectively detected by resonant two-photon ionization (R2PI) with a third (probe) laser. If the pumped level is above D0(S0), the hot complex dissociates and the probe signal disappears. Combining the fluorescence or SEP spectrum of the cold complex with the SEP breakoff of the hot complex brackets D0(S0). The UV chromophores 1-naphthol and carbazole were employed; these bind either dispersively via the aromatic rings, or form a hydrogen bond via the -OH or -NH group. Dissociation energies have been measured for dispersively bound complexes with noble gases (Ne, Kr, Ar, Xe), diatomics (N2, CO), alkanes (methane to n-butane), cycloalkanes (cyclopropane to cycloheptane), and unsaturated compounds (ethene, benzene). Hydrogen-bond dissociation energies have been measured for H2O, D2O, methanol, ethanol, ethers (oxirane, oxetane), NH3 and ND3.

The S 1 /S 2 splitting of the m-cyanophenol dimer, (mCP) 2 and the delocalization of its excitonically coupled S 1 /S 2 states are investigated by mass-selective two-color resonant two-photon ionization and dispersed fluorescence spectroscopy, complemented by a theoretical vibronic coupling analysis based on correlated ab initio calculations at the approximate coupled cluster CC2 and SCS-CC2 levels. The calculations predict three close-lying ground-state minima of (mCP) 2 : The lowest is slightly Z-shaped (C i -symmetric); the second-lowest is <5 cm -1 higher and planar (C 2h ). The vibrational ground state is probably delocalized over both minima. The S 0 → S 1 transition of (mCP) 2 is electric-dipole allowed (A g → A u ), while the S 0 → S 2 transition is forbidden (A g → A g ). Breaking the inversion symmetry by 12 C/ 13 C- or H/D-substitution renders the S 0 → S 2 transition partially allowed; the excitonic contribution to the S 1 /S 2 splitting is Δ exc = 7.3 cm -1 . Additional isotope-dependent contributions arise from the changes of the m-cyanophenol zero-point vibrational energy upon electronic excitation, which are Δ iso ( 12 C/ 13 C) = 3.3 cm -1 and Δ iso (H/D) = 6.8 cm -1 . Only partial localization of the exciton occurs in the 12 C/ 13 C and H/D substituted heterodimers. The SCS-CC2 calculated excitonic splitting is Δ el = 179 cm -1 ; when multiplying this with the vibronic quenching factor Γ vibron exp = 0.043, we obtain an exciton splitting Δ vibron exp = 7.7 cm -1 , which agrees very well with the experimental Δ exc = 7.3 cm -1 . The semiclassical exciton hopping times range from 3.2 ps in (mCP) 2 to 5.7 ps in the heterodimer (mCP-h)·(mCP-d). A multimode vibronic coupling analysis is performed encompassing all the vibronic levels of the coupled S 1 /S 2 states from the v = 0 level to 600 cm -1 above. Both linear and quadratic vibronic coupling schemes were investigated to simulate the S 0 → S 1 /S 2 vibronic spectra; those calculated with the latter scheme agree better with experiment.

Intermolecular dissociation energies D 0 (S 0 ) of the supersonic jet-cooled complexes of 1-naphthol (1NpOH) with cyclopentane, cyclohexane, and cycloheptane were determined to within <0.5% using the stimulated-emission pumping resonant two-photon ionization method. The ground state D 0 (S 0 ) values are bracketed as 20.23±0.07 kJ/mol for 1NpOH⋅cyclopentane, 20.34±0.04 kJ/mol for 1NpOH⋅cyclohexane, and 22.07±0.10 kJ/mol for two isomers of 1NpOH⋅cycloheptane. Upon S 0 →S 1 excitation of the 1-naphthol chromophore, the dissociation energies of the 1NpOH⋅cycloalkane complexes increase from 0.1% to 3%. Three dispersion-corrected density functional theory (DFT) methods predict that the cycloalkane moieties are dispersively bound to the naphthol face via London-type interactions, similar to the "face" isomer of the 1-naphthol⋅cyclopropane complex [S. Maity et al., J. Chem. Phys. 145, 164304 (2016)]. The experimental and calculated D 0 (S 0 ) values of the cyclohexane and cyclopentane complexes are practically identical, although the polarizability of cyclohexane is ∼20% larger than that of cyclopentane. Investigation of the calculated pairwise atomic contributions to the D2 dispersion energy reveals that this is due to subtle details of the binding geometries of the cycloalkanes relative to the 1-naphthol ring. The B97-D3 DFT method predicts dissociation energies within about ±1% of experiment, including the cyclopropane face complex. The B3LYP-D3 and ωB97X-D calculated dissociation energies are 7-9 and 13-20% higher than the experimental D 0 (S 0 ) values. Without dispersion correction, all the complexes are calculated to be unbound.

The 1-naphthol⋅cyclopropane intermolecular complex is formed in a supersonic jet and investigated by resonant two-photon ionization (R2PI) spectroscopy, UV holeburning, and stimulated emission pumping (SEP)-R2PI spectroscopy. Two very different structure types are inferred from the vibronic spectra and calculations. In the "edge" isomer, the OH group of 1-naphthol is directed towards a C-C bond of cyclopropane, the two ring planes are perpendicular. In the "face" isomer, the cyclopropane is adsorbed on one of the π-aromatic faces of the 1-naphthol moiety, the ring planes are nearly parallel. Accurate ground-state intermolecular dissociation energies D 0 were measured with the SEP-R2PI technique. The D 0 (S 0 ) of the edge isomer is bracketed as 15.35 ± 0.03 kJ/mol, while that of the face isomer is 16.96 ± 0.12 kJ/mol. The corresponding excited-state dissociation energies D 0 (S 1 ) were evaluated using the respective electronic spectral shifts. Despite the D 0 (S 0 ) difference of 1.6 kJ/mol, both isomers are observed in the jet in similar concentrations, so they must be separated by substantial potential energy barriers. Intermolecular binding energies, D e , and dissociation energies, D 0 , calculated with correlated wave function methods and two dispersion-corrected density-functional methods are evaluated in the context of these results. The density functional calculations suggest that the face isomer is bound solely by dispersion interactions. Binding of the edge isomer is also dominated by dispersion, which makes up two thirds of the total binding energy.

We have investigated the S 0 → S 1 UV vibronic spectrum and time-resolved S 1 state dynamics of jet-cooled amino-keto 1-methylcytosine (1MCyt) using two-color resonant two-photon ionization, UV/UV holeburning and depletion spectroscopies, as well as nanosecond and picosecond time-resolved pump/delayed ionization measurements. The experimental study is complemented with spin-component-scaled second-order coupled-cluster and multistate complete active space second order perturbation ab initio calculations. Above the weak electronic origin of 1MCyt at 31 852 cm -1 about 20 intense vibronic bands are observed. These are interpreted as methyl group torsional transitions coupled to out-of-plane ring vibrations, in agreement with the methyl group rotation and out-of-plane distortions upon 1 ππ ∗ excitation predicted by the calculations. The methyl torsion and ν 1 ' (butterfly) vibrations are strongly coupled, in the S 1 state. The S 0 → S 1 vibronic spectrum breaks off at a vibrational excess energy E exc ∼ 500 cm -1 , indicating that a barrier in front of the ethylene-type S 1 ⇝S 0 conical intersection is exceeded, which is calculated to lie at E exc = 366 cm -1 . The S 1 ⇝S 0 internal conversion rate constant increases from k IC = 2 ⋅ 10 9 s -1 near the S 1 (v = 0) level to 1 ⋅ 10 11 s -1 at E exc = 516 cm -1 . The 1 ππ ∗ state of 1MCyt also relaxes into the lower-lying triplet T 1 ( 3 ππ ∗ ) state by intersystem crossing (ISC); the calculated spin-orbit coupling (SOC) value is 2.4 cm -1 . The ISC rate constant is 10-100 times lower than k IC ; it increases from k ISC = 2 ⋅ 10 8 s -1 near S 1 (v = 0) to k ISC = 2 ⋅ 10 9 s -1 at E exc = 516 cm -1 . The T 1 state energy is determined from the onset of the time-delayed photoionization efficiency curve as 25 600 ± 500 cm -1 . The T 2 ( 3 nπ ∗ ) state lies >1500 cm -1 above S 1 (v = 0), so S 1 ⇝T 2 ISC cannot occur, despite the large SOC parameter of 10.6 cm -1 . An upper limit to the adiabatic ionization energy of 1MCyt is determined as 8.41 ± 0.02 eV. Compared to cytosine, methyl substitution at N1 lowers the adiabatic ionization energy by ≥0.32 eV and leads to a much higher density of vibronic bands in the S 0 → S 1 spectrum. The effect of methylation on the radiationless decay to S 0 and ISC to T 1 is small, as shown by the similar break-off of the spectrum and the similar computed mechanisms.

The dissociation energy (D0) of an isolated and cold molecular complex in the gas-phase is a fundamental measure of the strength of the intermolecular interactions between its constituent moieties. Accurate D0 values are important for the understanding of intermolecular bonding, for benchmarking high-level theoretical calculations, and for the parametrization of force-field models used in fields ranging from crystallography to biochemistry. We review experimental and theoretical methods for determining gas-phase D0 values of M·S complexes, where M is a (hetero)aromatic molecule and S is a closed-shell "solvent" atom or molecule. The experimental methods discussed involve M-centered (S0 → S1) electronic excitation, which is often followed by ionization to the M(+)·S ion. The D0 is measured by depositing a defined amount of vibrational energy in the neutral ground state, giving M(‡)·S, the neutral S1 excited state, giving M*·S, or the M(+)·S ion ground state. The experimental methods and their relative advantages and disadvantages are discussed. Based on the electronic structure of M and S, we classify the M·S complexes as Type I, II, or III, and discuss characteristic properties of their respective potential energy surfaces that affect or hinder the determination of D0. Current theoretical approaches are reviewed, which comprise methods based on a Kohn-Sham reference determinant as well as wave function-based methods based on coupled-cluster theory.

Ultraviolet radiative damage to DNA is inefficient because of the ultrafast S1 ⇝ S0 internal conversion of its nucleobases. Using picosecond pump-ionization delay measurements, we find that the S1((1)ππ*) state vibrationless lifetime of gas-phase keto-amino cytosine (Cyt) is τ = 730 ps or ∼ 700 times longer than that measured by femtosecond pump-probe ionization at higher vibrational excess energy, Eexc. N1-Alkylation increases the S1 lifetime up to τ = 1030 ps for N1-ethyl-Cyt but decreases it to 100 ps for N1-isopropyl-Cyt. Increasing the vibrational energy to Eexc = 300-550 cm(-1) decreases the lifetimes to 20-30 ps. The nonradiative dynamics of S1 cytosine is not solely a property of the amino-pyrimidinone chromophore but is strongly influenced by the N1-substituent. Correlated excited-state calculations predict that the gap between the S2((1)nOπ*) and S1((1)ππ*) states decreases along the series of N1-derivatives, thereby influencing the S1 state lifetime.

The excitonic splitting and vibronic quenching of the inversion-symmetric homodimers of benzonitrile, (BN)2, and meta-cyanophenol, (mCP)2, are investigated by two-color resonant two-photon ionization spectroscopy. These systems have very different hydrogen bond strengths: the OH···N≡C bonds in (mCP)2 are ∼10 times stronger than the CH···N≡C hydrogen bonds in (BN)2. In (BN)2 the S0((1)Ag) → S1((1)Ag) transition is electric-dipole forbidden, while the S0((1)Ag) → S2((1)Bu) transition is allowed. The opposite holds for (mCP)2 due to the different transition dipole moment vector alignment. The S0 → S1S2 spectra of the dimers are compared and their excitonic splittings and vibronic quenchings are investigated by measuring the (13)C-substituted heterodimer isotopomers, for which the centrosymmetry is broken and both transitions are allowed. The excitonic splittings are determined as Δexc = 2.1 cm(-1) for (BN)2 and Δexc = 7.3 cm(-1) for (mCP)2. The latter exhibits a much stronger vibronic quenching, as the purely electronic splitting resulting from ab initio calculations is determined to be Δcalc = 179 cm-1, while in (BN)2 the calculated splitting is Δcalc = 10 cm(-1). The monomer site-shifts upon dimerization and comparing certain vibrations that deform the hydrogen bonds confirm that the OH···N≡C hydrogen bond is much stronger than the CH···N≡C bond. We show that the H-bonds have large effects on the spectral shifts, but little or no influence on the excitonic splitting.

The amino-keto tautomer of supersonic jet-cooled cytosine undergoes intersystem crossing (ISC) from the v = 0 and low-lying vibronic levels of its S1((1)ππ(∗)) state. We investigate these ISC rates experimentally and theoretically as a function of S1 state vibrational excess energy Eexc. The S1 vibronic levels are pumped with a ∼5 ns UV laser, the S1 and triplet state ion signals are separated by prompt or delayed ionization with a second UV laser pulse. After correcting the raw ISC yields for the relative S1 and T1 ionization cross sections, we obtain energy dependent ISC quantum yields QISC (corr)=1%-5%. These are combined with previously measured vibronic state-specific decay rates, giving ISC rates kISC = 0.4-1.5 ⋅ 10(9) s(-1), the corresponding S1⇝S0 internal conversion (IC) rates are 30-100 times larger. Theoretical ISC rates are computed using SCS-CC2 methods, which predict rapid ISC from the S1; v = 0 state with kISC = 3 ⋅ 10(9) s(-1) to the T1((3)ππ(∗)) triplet state. The surprisingly high rate of this El Sayed-forbidden transition is caused by a substantial admixture of (1)nOπ(∗) character into the S1((1)ππ(∗)) wave function at its non-planar minimum geometry. The combination of experiment and theory implies that (1) below Eexc = 550 cm(-1) in the S1 state, S1⇝S0 internal conversion dominates the nonradiative decay with kIC ≥ 2 ⋅ 10(10) s(-1), (2) the calculated S1⇝T1 ((1)ππ(∗)⇝(3)ππ(∗)) ISC rate is in good agreement with experiment, (3) being El-Sayed forbidden, the S1⇝T1 ISC is moderately fast (kISC = 3 ⋅ 10(9) s(-1)), and not ultrafast, as claimed by other calculations, and (4) at Eexc ∼ 550 cm(-1) the IC rate increases by ∼50 times, probably by accessing the lowest conical intersection (the C5-twist CI) and thereby effectively switching off the ISC decay channels.

Femtosecond time-resolved Raman rotational coherence spectroscopy (RCS) is employed to determine accurate rotational, vibration–rotation coupling constants, and centrifugal distortion constants of cyclopentane (C5H10). Its lowest-frequency vibration is a pseudorotating ring deformation that interconverts 10 permutationally distinct but energetically degenerate "twist" minima interspersed by 10 "bent" conformers. While the individual twist and bent structures are polar asymmetric tops, the pseudorotation is fast on the time scale of external rotation, rendering cyclopentane a fluxionally nonpolar symmetric top molecule. The pseudorotational level pattern corresponds to a one-dimensional internal rotor with a pseudorotation constant Bps ≈ 2.8 cm(-1). The pseudorotational levels are significantly populated up to l = ± 13 at 298 K; <10% of the molecules are in the l = 0 level. The next-higher vibration is the “radial” ν23 ring deformation mode at 273 cm–1, which is far above the pseudorotational fundamental. Femtosecond Raman RCS measurements were performed in a gas cell at T = 293 K and in a pulsed supersonic jet at T ≈ 90 K. The jet cooling reduces the pseudorotational distribution to l < ±8 and eliminates the population of ν23, allowing one to determine the rotational constant as A0 = B0 = 6484.930(11) MHz. This value is ∼300 times more precise than the previous value. The fit of the RCS transients reveals that the rotation–pseudorotation coupling constant αe,psB = −0.00070(1) MHz is diminutive, implying that excitation of the pseudorotation has virtually no effect on the B0 rotational constant of cyclopentane. The smallness of αe,psB can be realized when comparing to the vibration–rotation coupling constant of the ν23 vibration, αe,23B = -9.547(1) MHz, which is about 104 times larger.

After decades of research on molecular excitons, only few molecular dimers are available on which exciton and vibronic coupling theories can be rigorously tested. In centrosymmetric H-bonded dimers consisting of identical (hetero)aromatic chromophores, the monomer electronic transition dipole moment vectors subtract or add, yielding S 0 → S 1 and S 0 → S 2 transitions that are symmetry-forbidden or -allowed, respectively. Symmetry breaking by 12 C/ 13 C or H/D isotopic substitution renders the forbidden transition weakly allowed. The excitonic coupling (Davydov splitting) can then be measured between the S 0 → S 1 and S 0 → S 2 vibrationless bands. We discuss the mass-specific excitonic spectra of five H-bonded dimers that are supersonically cooled to a few K and investigated using two-color resonant two-photon ionization spectroscopy. The excitonic splittings Δ calc predicted by ab initio methods are 5-25 times larger than the experimental excitonic splittings Δ exp . The purely electronic ab initio splittings need to be reduced ("quenched"), reflecting the coupling of the electronic transition to the optically active vibrations of the monomers. The so-called quenching factors Γ < 1 can be determined from experiment ( Γ exp ) and/or calculation ( Γ calc ). The vibronically quenched splittings Γ · Δ calc are found to nicely reproduce the experimental exciton splittings.

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = -451 and -858 cm(-1) relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm(-1) in S0 and 55 cm(-1) in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of -D0 = 28-35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10(3)-10(4) times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ∼2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J -, K -, A -, and C-type transients, as well as a novel transient (S-type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants A0 = 5637.68(20) MHz, B0 = 1428.23(37) MHz, and C0 = 1138.90(48) MHz (1σ uncertainties). Combining the A0, B0, and C0 constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths re(C1-C2) = 1.3849(4) Å, re(C2-C3) = 1.3917(4) Å, re(C-F) = 1.3422(3) Å, and re(C2-H2) = 1.0791(5) Å.

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

We observe the weak S0 → S2 transitions of the T-shaped benzene dimers (Bz)2 and (Bz-d6)2 about 250 cm(-1) and 220 cm(-1) above their respective S0 → S1 electronic origins using two-color resonant two-photon ionization spectroscopy. Spin-component scaled (SCS) second-order approximate coupled-cluster (CC2) calculations predict that for the tipped T-shaped geometry, the S0 → S2 electronic oscillator strength fel(S2) is ∼10 times smaller than fel(S1) and the S2 state lies ∼240 cm(-1) above S1, in excellent agreement with experiment. The S0 → S1 (ππ(∗)) transition is mainly localized on the "stem" benzene, with a minor stem → cap charge-transfer contribution; the S0 → S2 transition is mainly localized on the "cap" benzene. The orbitals, electronic oscillator strengths fel(S1) and fel(S2), and transition frequencies depend strongly on the tipping angle ω between the two Bz moieties. The SCS-CC2 calculated S1 and S2 excitation energies at different T-shaped, stacked-parallel and parallel-displaced stationary points of the (Bz)2 ground-state surface allow to construct approximate S1 and S2 potential energy surfaces and reveal their relation to the "excimer" states at the stacked-parallel geometry. The fel(S1) and fel(S2) transition dipole moments at the C2v-symmetric T-shape, parallel-displaced and stacked-parallel geometries are either zero or ∼10 times smaller than at the tipped T-shaped geometry. This unusual property of the S0 → S1 and S0 → S2 transition-dipole moment surfaces of (Bz)2 restricts its observation by electronic spectroscopy to the tipped and tilted T-shaped geometries; the other ground-state geometries are impossible or extremely difficult to observe. The S0 → S1/S2 spectra of (Bz)2 are compared to those of imidazole ⋅ (Bz)2, which has a rigid triangular structure with a tilted (Bz)2 subunit. The S0 → S1/ S2 transitions of imidazole-(benzene)2 lie at similar energies as those of (Bz)2, confirming our assignment of the (Bz)2 S0 → S2 transition.

The S2←S0 vibronic spectrum of the ortho-cyanophenol dimer (oCP)2 is analyzed in a joint experimental and theoretical investigation. Vibronic excitation energies up to 750 cm(-1) are covered, which extends our previous analysis of the quenching of the excitonic splitting in this and related species [Kopec et al., J. Chem. Phys. 137, 184312 (2012)]. As we demonstrate, this necessitates an extension of the coupling model. Accordingly, we compute the potential energy surfaces of the ortho-cyanophenol dimer (oCP)2 along all relevant normal modes using the approximate second-order coupled cluster method RI-CC2 and extract the corresponding coupling constants using the linear and quadratic vibronic coupling scheme. These serve as the basis to calculate the vibronic spectrum. The theoretical results are found to be in good agreement with the experimental highly resolved resonant two-photon ionization spectrum. This allows to interpret key features of the excitonic and vibronic interactions in terms of nodal patterns of the underlying vibronic wave functions.

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.

The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B0 = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B0 value, RR-RCS measurements in a room temperature gas cell give the rotational constants Bv of the five lowest-lying thermally populated vibrationally excited states ν7/8, ν9, ν11/12, ν13, and ν14/15. Their Bv constants differ from B0 by between -1.02 MHz and +2.23 MHz. Combining the B0 with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys. 111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths re(C-C) = 1.3866(3) Å and re(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ re bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths rg(C-C)=1.3907(3) Å and rg(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction rg bond lengths measured in the 1960s.

2-Aminopurine (2 AP) is a fluorescent isomer of adenine and has a fluorescence lifetime of ~11 ns in water. It is widely used in biochemical settings as a site-specific fluorescent probe of DNA and RNA structure and base-flipping and -folding. These assays assume that 2 AP is intrinsically strongly fluorescent. Here, we show this not to be the case, observing that gas-phase, jet-cooled 2-aminopurine and 9-methyl-2-aminopurine have very short fluorescence lifetimes (156 ps and 210 ps, respectively); they are, to all intents and purposes, non-fluorescent. We find that the lifetime of 2-aminopurine increases dramatically when it is part of a hydrate cluster, 2 AP · (H2O)n, where n = 1-3. Not only does it depend on the presence of water molecules, it also depends on the specific hydrogen-bonding site to which they attach and on the number of H2O molecules at that site. We selectively microhydrate 2-aminopurine at its sugar-edge, cis-amino or trans-amino sites and see that its fluorescence lifetime increases by 4, 50 and 95 times (to 14.5 ns), respectively.

The S0 → S1 vibronic spectrum and S1 state nonradiative relaxation of jet-cooled keto-amino 5-fluorocytosine (5FCyt) are investigated by two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm(–1) resolution. The 0(0)(0) rotational band contour is polarized in-plane, implying that the electronic transition is (1)ππ*. The electronic transition dipole moment orientation and the changes of rotational constants agree closely with the SCS-CC2 calculated values for the (1)ππ* (S1) transition of 5FCyt. The spectral region from 0 to 300 cm(–1) is dominated by overtone and combination bands of the out-of-plane ν1′ (boat), ν2′ (butterfly), and ν3′ (HN–C6H twist) vibrations, implying that the pyrimidinone frame is distorted out-of-plane by the (1)ππ* excitation, in agreement with SCS-CC2 calculations. The number of vibronic bands rises strongly around +350 cm(–1); this is attributed to the (1)ππ* state barrier to planarity that corresponds to the central maximum of the double-minimum out-of-plane vibrational potentials along the ν1′, ν2′, and ν3′ coordinates, which gives rise to a high density of vibronic excitations. At +1200 cm(–1), rapid nonradiative relaxation (k(nr) ≥ 10(12) s(–1)) sets in, which we interpret as the height of the (1)ππ* state barrier in front of the lowest S1/S0 conical intersection. This barrier in 5FCyt is 3 times higher than that in cytosine. The lifetimes of the ν′ = 0, 2ν1′, 2ν2′, 2ν1′ + 2ν2′, 4ν2′, and 2ν1′ + 4ν2′ levels are determined from Lorentzian widths fitted to the rotational band contours and are τ ≥ 75 ps for ν′ = 0, decreasing to τ ≥ 55 ps at the 2ν1′ + 4ν2′ level at +234 cm(–1). These gas-phase lifetimes are twice those of S1 state cytosine and 10–100 times those of the other canonical nucleobases in the gas phase. On the other hand, the 5FCyt gas-phase lifetime is close to the 73 ps lifetime in room-temperature solvents. This lack of dependence on temperature and on the surrounding medium implies that the 5FCyt nonradiative relaxation from its S1 ((1)ππ*) state is essentially controlled by the same ~1200 cm(–1) barrier and conical intersection both in the gas phase and in solution.

The UV spectrum of the adenine analogue 9-methyl-2-aminopurine (9M-2AP) is investigated with one- and two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm(-1) resolution in a supersonic jet. The electronic origin at 32,252 cm(-1) exhibits methyl torsional subbands that originate from the 0A1'' (l = 0) and 1E(″) (l = ±1) torsional levels. These and further torsional bands that appear up to 00 (0)+230 cm(-1) allow to fit the threefold (V3) barriers of the torsional potentials as |V3''|=50 cm(-1) in the S0 and |V3'|=126 cm(-1) in the S1 state. Using the B3LYP density functional and correlated approximate second-order coupled cluster CC2 methods, the methyl orientation is calculated to be symmetric relative to the 2AP plane in both states, with barriers of V3''=20 cm(-1) and V3'=115 cm(-1). The 00 (0) rotational band contour is 75% in-plane (a/b) polarized, characteristic for a dominantly long-axis (1)ππ(*) excitation. The residual 25% c-axis polarization may indicate coupling of the (1)ππ(*) to the close-lying (1)nπ(*) state, calculated at 4.00 and 4.01 eV with the CC2 method. However, the CC2 calculated (1)nπ oscillator strength is only 6% of that of the (1)ππ(*) transition. The (1)ππ(*) vibronic spectrum is very complex, showing about 40 bands within the lowest 500 cm(-1). The methyl torsion and the low-frequency out-of-plane ν1' and ν2' vibrations are strongly coupled in the (1)ππ(*) state. This gives rise to many torsion-vibration combination bands built on out-of-plane fundamentals, which are without precedence in the (1)ππ(*) spectrum of 9H-2-aminopurine [S. Lobsiger, R. K. Sinha, M. Trachsel, and S. Leutwyler, J. Chem. Phys. 134, 114307 (2011)]. From the Lorentzian broadening needed to fit the 00 (0) contour of 9M-2AP, the (1)ππ(*) lifetime is τ ⩾ 120 ps, reflecting a rapid nonradiative transition.

While keto-amino cytosine is the dominant species in aqueous solution, spectroscopic studies in molecular beams and in noble gas matrices show that other cytosine tautomers prevail in apolar environments. Each of these offers two or three H-bonding sites (Watson-Crick, wobble, sugar-edge). The mass- and isomer-specific S1 ← S0 vibronic spectra of cytosine·2-pyridone (Cyt·2PY) and 1-methylcytosine·2PY are measured using UV laser resonant two-photon ionization (R2PI), UV/UV depletion, and IR depletion spectroscopy. The UV spectra of the Watson-Crick and sugar-edge isomers of Cyt·2PY are separated using UV/UV spectral hole-burning. Five different isomers of Cyt·2PY are observed in a supersonic beam. We show that the Watson-Crick and sugar-edge dimers of keto-amino cytosine with 2PY are the most abundant in the beam, although keto-amino-cytosine is only the third most abundant tautomer in the gas phase. We identify the different isomers by combining three different diagnostic tools: (1) methylation of the cytosine N1-H group prevents formation of both the sugar-edge and wobble isomers and gives the Watson-Crick isomer exclusively. (2) The calculated ground state binding and dissociation energies, relative gas-phase abundances, excitation and the ionization energies are in agreement with the assignment of the dominant Cyt·2PY isomers to the Watson-Crick and sugar-edge complexes of keto-amino cytosine. (3) The comparison of calculated ground state vibrational frequencies to the experimental IR spectra in the carbonyl stretch and NH/OH/CH stretch ranges strengthen this identification.

The picosecond (ps) timescale is relevant for the investigation of many molecular dynamical processes such as fluorescence, nonradiative relaxation, intramolecular vibrational relaxation, molecular rotation and intermolecular energy transfer, to name a few. While investigations of ultrafast (femtosecond) processes of biological molecules, e.g. nucleobases and their analogues in the gas phase are available, there are few investigations on the ps time scale. We have constructed a ps pump-ionization setup and a ps streak camera fluorescence apparatus for the determination of lifetimes of supersonic jet-cooled and isolated molecules and clusters. The ps pump-ionization setup was used to determine the lifetimes of the nucleobase analogue 2-aminopurine (2AP) and of two 2AP˙(H2O)n water cluster isomers with n=1 and 2. Their lifetimes lie between 150 ps and 3 ns and are strongly cluster-size dependent. The ps streak camera setup was used to determine accurate fluorescence lifetimes of the uracil analogue 2-pyridone (2PY), its self-dimer (2PY)2, two isomers of its trimer (2PY)3 and its tetramer (2PY)4, which lie in the 7-12 ns range.

2-Aminopurine (2AP) is an adenine analogue with a high fluorescence quantum yield in water solution, which renders it a useful real-time probe of DNA structure. We report the ultraviolet (UV) and infrared (IR) spectra of size-selected and jet-cooled 9H-2AP·(H2O)n clusters with n = 2 and 3. Mass- and species-specific UV/UV holeburning spectroscopy allows to separate the UV spectra of four cluster isomers in the 31,200–33,000 cm(–1) spectral region with electronic band origins at 31339, 31450, 31891, and 32163 cm(–1). Using IR/UV depletion spectroscopy in combination with B3LYP calculated harmonic vibrational frequencies, the H-bonding topologies of two isomers of the n = 2 and of two isomers of the n = 3 cluster are identified. One n = 2 isomer (denoted 2A) forms a water dimer chain between the N9H and N3 atoms at the sugar-edge site, the other isomer (denoted 2D) binds one H2O at the sugar-edge site and the other at the trans-amino site between the N1 atom and the NH2 group. For 2-aminopurine·(H2O)3, one isomer (denoted 3A) forms an H-bonded water wire at the sugar-edge site, while isomer 3B accommodates two H2O molecules at the sugar-edge and one at the trans-amino site. The approximate second-order coupled cluster (CC2) method predicts the adiabatic S1 ← S0 transitions of 9H-2-aminopurine and six water cluster isomers with n = 1–3 in very good agreement with the experimental 0(0)(0) frequencies, with differences of <0.6%. The stabilization of the S1(ππ*) state of 2-aminopurine by water clusters is highly regiospecific: Isomers with one or two H2O molecules H-bonded in the trans-amino position induce large spectra red shifts, corresponding to 1ππ* state stabilization of 10–12 kJ/mol, while water-wire cluster solvation at the sugar-edge leads to much smaller stabilization. The evolution of the IR spectra of the water-wire clusters with n = 1–3 that are H-bonded to the sugar-edge site is discussed. Qualitatively different regions (denoted I to IV) can be attributed to the different free and H-bonded OH, NH, NH2 and OH···OH water-wire stretch vibrations.

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.

Mass- and isomer-selected infrared spectra of 2-pyridone·(NH3)n clusters with n = 1-3 were measured in the NH and CH stretch fundamental region (2400-3700 cm(-1)) using infrared (IR) laser depletion spectroscopy combined with resonant two-photon ionization UV laser detection. The IR depletion spectra reveal three different H-bonding topologies of these clusters: The n = 1 and 2 clusters form ammonia bridges stretching from the N-H to the C═O group of the cis-amide function of 2-pyridone (2PY), giving rise to intense and strongly red-shifted (2PY)NH and ammonia NH stretch bands. For n = 3, two isomers (3X and 3Y) are observed in the IR spectra: The spectrum of 3X is compatible with an ammonia-bridge structure like n = 2, with the third NH3 accepting an H-bond from C(6)-H of 2PY. The IR spectrum of 3Y exhibits a broad IR band in the 2500-3000 cm(-1) range and is characteristic of a bifurcated double-bridged structure in which the first NH3 accepts an H-bond from the (2PY)NH and donates two H-bonds to the other two ammonias, both of which donate to the C═O group of 2PY. This double-donor/double-bridge H-bonding pattern increases the acceptor strength of the first ammonia and dramatically lowers the (2PY)NH stretching frequency to ∼2700 cm(-1). For all clusters the ammonia 2ν4 HNH bend overtones in the 3180-3320 cm(-1) region gain intensity by anharmonic coupling (Fermi resonance) to the hydrogen-bonded ammonia NH stretches, which are red-shifted into the 3250-3350 cm(-1) region. The experimental results are supported by optimized structures, vibrational frequencies, and IR intensities calculated using density-functional theory with the B3LYP and PW91 functionals, as well as with the more recent functionals B97-D and M06-2X, which are designed to include long-range dispersive interactions.

We have measured the mass- and tautomer-specific S0 → S1 vibronic spectra and S1 state lifetimes of the keto–amino tautomer of cytosine cooled in supersonic jets, using two-color resonant two-photon ionization (R2PI) spectroscopy at 0.05 cm(–1) resolution. The rotational contours of the 0(0)(0) band and nine vibronic bands up to +437 cm(–1) are polarized in the pyrimidinone plane, proving that the electronic excitation is to a 1ππ* state. All vibronic excitations up to +437 cm(–1) are overtone and combination bands of the low-frequency out-of-plane ν1′ (butterfly), ν2′ (boat), and ν3′ (H–N–C6–H twist) vibrations. UV vibronic spectrum simulations based on approximate second-order coupled-cluster (CC2) calculations of the ground and 1ππ* states are in good agreement with the experimental R2PI spectrum, but only if the calculated ν1′ and ν2′ frequencies are reduced by a factor of 4 and anharmonicity is included. Together with the high intensity of the ν1′ and ν2′ overtone vibronic excitations, this implies that the 1ππ* potential energy surface is much softer and much more anharmonic in the out-of-plane directions than predicted by the CC2 method. The 1ππ* state lifetime is determined from the Lorentzian line broadening necessary to reproduce the rotational band contours: at the 0(0)(0) band it is τ = 44 ps, remains at τ = 35–45 ps up to +205 cm(–1), and decreases to 25–30 ps up to +437 cm(–1). These lifetimes are 20–40 times longer than the 0.5–1.5 ps lifetimes previously measured with femtosecond pump–probe techniques at higher vibrational energies (1500–3800 cm(–1)). Thus, the nonradiative relaxation rate of keto–amino cytosine close to the 1ππ* state minimum is k(nr) 2.5 × 10(10) s(–1), much smaller than at higher energies. An additional nonradiative decay channel opens at +500 cm(–1) excess energy. Since high overtone bands of ν1′ and ν2′ are observed in the R2PI spectrum but only a single weak 2ν3′ band, we propose that ν3′ is a promoting mode for nonradiative decay, consistent with the observation that the ν3′ normal-mode eigenvector points toward the “C6-puckered” conical intersection geometry.

In the present study a benchmark set of medium-sized and large aromatic organic molecules with 10-78 atoms is presented. For this test set 0-0 transition energies measured in supersonic jets are compared to those calculated with DFT and the B3LYP functional, ADC(2), CC2 and the spin-scaled CC2 variants SOS-CC2 and SCS-CC2. Geometries of the ground and excited states have been optimized with these methods in polarized triple zeta basis sets. Zero-point vibrational corrections have been calculated with the same methods and basis sets. In addition the energies have been corrected by single point calculations with a triple zeta basis augmented with diffuse functions, aug-cc-pVTZ. The deviations of the theoretical results from experimental electronic origins, which have all been measured in the gas phase with high-resolution techniques, were evaluated. The accuracy of SOS-CC2 is comparable to that of unscaled CC2, whereas ADC(2) has slightly larger errors. The lowest errors were found for SCS-CC2. All correlated wave function methods provide significantly better results than DFT with the B3LYP functional. The effects of the energy corrections from the augmented basis set and the method-consistent calculation of the zero-point vibrational corrections are small. With this benchmark set reliable reference data for 0-0 transition energies for larger organic chromophores are available that can be used to benchmark the accuracy of other quantum chemical methods such as new DFT functionals or semi-empirical methods for excitation energies and structures and thereby augments available benchmark sets augments present benchmark sets which include mainly smaller molecules.

Gas-phase cytosine exists in five different tautomer/rotamer forms 1, 2a, 2b, 3a, and 3b. We determine the threshold ionization energy (IE) of the keto-amino tautomer 1 as 8.73 ± 0.02 eV, using resonant two-photon ionization mass spectrometry in a supersonic molecular beam via the (1)ππ* excited state. This is the first IE threshold measurement for the biologically relevant tautomer 1. The IE of the thermal gas-phase mixture of cytosine has been measured as 8.60 ± 0.05 eV by Kostko et al. using single-photon VUV photoionization [Phys. Chem. Chem. Phys., 2010, 12, 2860]. Given the tautomer distribution and ionization energies calculated in that work, our determination of the keto-amino tautomer IE implies that the IE measured by Kostko et al. is dominated by the enol-amino tautomers 2a and 2b. Upon excitation of keto-amino cytosine to its (1)ππ* state, relaxation occurs to a lower-lying long-lived state. The IE threshold measured via this state places its energy about 0.69 eV below the (1)ππ* state, in good agreement with the triplet T1 energy of keto-amino cytosine calculated by several high-level ab initio methods. The identification of keto-amino cytosine T1 is the basis for characterizing the intersystem crossing rates into and the photochemical reactions of this long-lived state.

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)(∗).

The quenching of the excitonic splitting in hydrogen-bonded molecular dimers has been explained recently in terms of exciton coupling theory, involving Förster's 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.

We investigate the UV vibronic spectrum and excited-state nonradiative processes of jet-cooled 5-methylcytosine (5MCyt) using two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm(–1) resolution. In contrast to cytosine, which shows only five bands above its electronic origin, the lowest electronic transition of 5MCyt exhibits about 25 low-frequency vibronic bands that extend to 0(0)(0) + 450 cm(–1), allowing to extract detailed information on the excited-state electronic and nuclear structure. Most bands are overtones and combinations of the out-of-plane vibrations ν'(1), ν'(2), and ν'(3). Their large intensities reflect butterfly-, boat-, and twist-deformations of the 5MCyt framework upon electronic excitation. From the rotational contours of the 0(0)(0), 1(0)(2), 2(0)(2), and 3(0)(2) bands, the transition is found to be polarized along the in-plane a/b axes, characteristic of a (1)ππ* transition. Approximate second-order coupled-cluster (CC2) and time-dependent B3LYP calculations both predict that 5MCyt undergoes an out-of-plane deformation in its (1)ππ* (S(2)) state but both methods overestimate the out-of-plane ν'(1), ν'(2), and ν'(3) vibrational frequencies by a factor of 3–5. The TD-B3LYP (1)ππ* transition dipole moment direction is 10%:90% a:b, in good agreement with experiment. From the Lorentzian line shape contributions needed to fit the rotational contours, a lower limit to the 5MCyt (1)ππ* state lifetime at the 0(0)(0), 1(0)(2), 2(0)(2), and 3(0)(2) bands is determined as τ ≥ 30 ps. These values are in stark contrast to the ultrafast (picosecond) lifetimes measured for jet-cooled cytosine by femtosecond pump–probe techniques. They also confirm the observation from the R2PI spectrum that 5-methylation of cytosine increases its excited-state lifetime. The higher out-of-plane overtone and combination bands disappear from the spectrum by ~460 cm(–1), signaling the onset of lifetimes τ < 0.5 ps, induced by excitation of the ν'(1), ν'(2), and ν'(3) modes. Above 460 cm(–1) the in-plane fundamentals 9(0)(1), 10(0)(1), and 17(0)(1) are weakly observed, implying that the in-plane vibrations couple weakly to the out-of-plane modes and have lifetimes τ > 3 ps up to 720 cm(–1) excess vibrational energy.

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.

The infrared (IR) spectra of the supersonic-jet cooled 9H- and 7H-tautomers of 2-aminopurine (2AP) and of the 9H-2-aminopurine·H(2)O monohydrate clusters have been measured by mass- and species-selective IR-UV double resonance spectroscopy in the 3200-3900 cm(-1) region, covering the N-H and O-H stretching vibrations. The spectra are complemented by density functional (B3LYP and PW91) and by second-order Møller-Plesset (MP2) calculations of the electronic energies and vibrational frequenciesof the respective 2AP tautomers and clusters. The 9H- and 7H-2-aminopurine tautomers were definitively identified by the shifts of their NH and NH(2) symmetric and asymmetric stretching frequencies and by comparison to the B3LYP/TZVP calculated IR spectra. The H-bond topologies of the two previously observed 9H-2-aminopurine·H(2)O isomers (Sinha. R. K.; et al. J. Phys. Chem. A2011, 115, 6208) are definitively identified as the "sugar-edge" isomer A and the "trans-amino-bound" isomer B by comparing their IR spectra to the calculated frequencies and IR intensities of the cluster isomers A, B, C, and D, as well as to the IR spectrum of 9H-2AP. The sugar-edge isomer A involves N9-H···OH(2) and HOH···N3 hydrogen bonds and is predicted to be the most stable form. The amino-bound isomer B involves NH(2)···OH(2) and HOH···N1 hydrogen bonds and is calculated to lie 2.5 kJ/mol above isomer A. The H-bond topology of the "cis-amino-bound" isomer C is symmetrically related to isomer B, with a hydrogen bond to the N3 of the pyrimidine group. However, it is calculated to lie 7 kJ/mol above isomer A and indeed is not observed in the supersonic jet. Isomer D involves a single H-bond to the N7 position, is predicted to be 14 kJ/mol above A and is therefore not observed.

We report the analysis of the S1<--S0 rotational band contours of jet-cooled 5-methyl-2-hydroxypyrimidine (5M2HP), the enol form of deoxythymine. Unlike thymine, which exhibits a structureless spectrum, the vibronic spectrum of 5M2HP is well structured, allowing us to determine the rotational constants and the methyl group torsional barriers in the S0 and S1 states. The 0(0)(0), 6a(0)(1), 6b(0)(1), and 14(0)(1) band contours were measured at 900 MHz (0.03 cm(-1)) resolution using mass-specific two-color resonant two-photon ionization (2C-R2PI) spectroscopy. All four bands are polarized perpendicular to the pyrimidine plane (>90% c type), identifying the S1<--S0 excitation of 5M2HP as a 1nπ* transition. All contours exhibit two methyl rotor subbands that arise from the lowest 5-methyl torsional states 0A" and 1E". The S0 and S1 state torsional barriers were extracted from fits to the torsional subbands. The 3-fold barriers are V3" = 13 cm(-1) and V3' = 51 cm(-1); the 6-fold barrier contributions V6" and V6' are in the range of 2-3 cm(-1) and are positive in both states. The changes of A, B, and C rotational constants upon S1 <--S0 excitation were extracted from the contours and reflect an “anti-quinoidal” distortion. The 0(0)(0) contour can only be simulated if a 3 GHz Lorentzian line shape is included, which implies that the S1(1nπ*) lifetime is ~55 ps. For the 6a(0)(1) and 6b(0)(1) bands, the Lorentzian component increases to 5.5 GHz, reflecting a lifetime decrease to ~30 ps. The short lifetimes are consistent with the absence of fluorescence from the 1nπ* state. Combining these measurements with the previous observation of efficient intersystem crossing (ISC) from the S1 state to a long-lived T1 (3nπ*) state that lies ~2200 cm(-1) below [S. Lobsiger, S. et al. Phys. Chem. Chem. Phys. 2010, 12, 5032] implies that the broadening arises from fast intersystem crossing with k(ISC) ≈ 2 × 10(10) s(-1). In comparison to 5-methylpyrimidine, the ISC rate is enhanced by at least 10 000 by the additional hydroxy group in position 2.

We combine the technique of femtosecond degenerate four-wave mixing (fs-DFWM) with a high repetition-rate pulsed supersonic jet source to obtain the rotational coherence spectrum (RCS) of cold cyclohexane (C(6)H(12)) with high signal/noise ratio. In the jet expansion, the near-parallel flow pattern combined with rapid translational cooling effectively eliminate dephasing collisions, giving near-constant RCS signal intensities over time delays up to 5 ns. The vibrational cooling in the jet eliminates the thermally populated vibrations that complicate the RCS coherences of cyclohexane at room temperature [Brügger, G.; et al. J. Phys. Chem. A 2011, 115, 9567]. The rotational cooling reduces the high-J rotational-state population, yielding the most accurate ground-state rotational constant to date, B(0) = 4305.859(9) MHz. Based on this B(0), a reanalysis of previous room-temperature gas-cell RCS measurements of cyclohexane gives improved vibration-rotation interaction constants for the ν(32), ν(6), ν(16), and ν(24) vibrational states. Combining the experimental B(0)(C(6)H(12)) with CCSD(T) calculations yields a very accurate semiexperimental equilibrium structure of the chair isomer of cyclohexane.

The excitonic splitting between the S(1) and S(2) electronic states of the doubly hydrogen-bonded dimer 2-pyridone[middle dot]6-methyl-2-pyridone (2PY·6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C(2h) symmetric (2-pyridone)(2) homodimer, in which the S(1) ← S(0) transition is symmetry-forbidden but the S(2) ← S(0) transition is allowed, the symmetry-breaking by the additional methyl group in 2PY·6M2PY leads to the appearance of both the S(1) and S(2) origins, which are separated by Δ(exp) = 154 cm(-1). When combined with the separation of the S(1) ← S(0) excitations of 6M2PY and 2PY, which is δ = 102 cm(-1), one obtains an S(1)/S(2) exciton coupling matrix element of V(AB, el) = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S(1)/S(2) ← S(0) spectrum of 2PY·6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a(') vibration and to the intermolecular σ(') stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a(')) = 0.15 and C(σ(')) = 0.05, which places this dimer in the weak-coupling limit. However, the S(1)/S(2) state exciton splittings Δ(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra- and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY·6M2PY the electronic splitting is quenched by a factor of ~30 (i.e., the vibronic quenching factor is Γ(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine·6-methyl-2-pyridone (2HP·6M2PY) are also observed and assigned., (© 2011 American Institute of Physics)

We combine femtosecond time-resolved rotational coherence spectroscopy with high-level ab initio theory to obtain accurate structural information for the nonpolar molecules cyclohexane (C(6)H(12)) and cyclohexane-d(12) (C(6)D(12)). We measured the rotational B(0) and centrifugal distortion constants D(J), D(JK) of the v = 0 states of C(6)H(12) and C(6)D(12) to high accuracy, for example, B(0)(C(6)H(12)) = 4306.08(5) MHz, as well as B(v) for the vibrationally excited states ν(32), ν(6), ν(16) and ν(24) of C(6)H(12) and additionally ν(15) for C(6)D(12). To successfully reproduce the experimental RCS transient, the overtone and combination levels 2ν(32), 3ν(32), ν(32) + ν(6), and ν(32) + ν(16) had to be included in the RCS model calculations. The experimental rotational constants are compared to those obtained at the second-order Møller-Plesset (MP2) level. Combining the experimental and calculated rotational constants with the calculated equilibrium bond lengths and angles allows determination of accurate semiexperimental equilibrium structure parameters, for example, r(e)(C-C) = 1.526 ± 0.001 Å, r(e)(C-H(axial)) = 1.098 ± 0.001 Å, and r(e)(C-H(equatorial)) = 1.093 ± 0.001 Å. The equilibrium C-C bond length of C(6)H(12) is only 0.004 Å longer than that of ethane. The effect of ring strain due to the unfavorable gauche interactions is mainly manifested as small deviations from the C-C-C, C-C-H(axial), and C-C-H(equatorial) angles from the tetrahedral value., (© 2011 American Chemical Society)

The N-H···π hydrogen bond is an important intermolecular interaction in many biological systems. We have investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet cooled complex of pyrrole with benzene and benzene-d(6) (Pyr·Bz, Pyr·Bz-d(6)). DFT-D density functional, SCS-MP2 and SCS-CC2 calculations predict a T-shaped and (almost) C(s) symmetric structure with an N-H···π hydrogen bond to the benzene ring. The pyrrole is tipped by ω(S(0)) = ±13° relative to the surface normal of Bz. The N···ring distance is 3.13 Å. In the S(1) excited state, SCS-CC2 calculations predict an increased tipping angle ω(S(1)) = ±21°. The IR depletion spectra support the T-shaped geometry: The NH stretch is redshifted by -59 cm(-1), relative to the "free" NH stretch of pyrrole at 3531 cm(-1), indicating a moderately strong N-H···π interaction. The interaction is weaker than in the (Pyr)(2) dimer, where the NH donor shift is -87 cm(-1) [Dauster et al., Phys. Chem. Chem. Phys., 2008, 10, 2827]. The IR C-H stretch frequencies and intensities of the Bz subunit are very similar to those of the acceptor in the (Bz)(2) dimer, confirming that Bz acts as the acceptor. While the S(1)←S(0) electronic origin of Bz is forbidden and is not observable in the gas-phase, the UV spectrum of Pyr·Bz in the same region exhibits a weak 0 band that is red-shifted by 58 cm(-1) relative to that of Bz (38 086 cm(-1)). The origin appears due to symmetry-breaking of the π-electron system of Bz by the asymmetric pyrrole NH···π hydrogen bond. This contrasts with (Bz)(2), which does not exhibit a 0 band. The Bz moiety in Pyr·Bz exhibits a 6a band at 0 + 518 cm(-1) that is about 20× more intense than the origin band. The symmetry breaking by the NH···π hydrogen bond splits the degeneracy of the ν(6)(e(2g)) vibration, giving rise to 6a' and 6b' sub-bands that are spaced by ∼6 cm(-1). Both the 0 and 6 bands of Pyr·Bz carry a progression in the low-frequency (10 cm(-1)) excited-state tipping vibration ω', in agreement with the change of the ω tipping angle predicted by SCS-MP2 and SCS-CC2 calculations.