Locating Cytosine Conical Intersections by Laser Experiments and Ab Initio Calculations.

The decay mechanism of S0 → S1 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 000 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 S1 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.

Intermolecular dissociation energies of 1-naphthol complexes with large dispersion-energy donors: Decalins and adamantane.

The ground-state intermolecular dissociation energies D0(S0) 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 D0(S0) 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 S0 → S1 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 D0(S0) values are within 1 kJ/mol of the experiment, while B3LYP-D3 predicts D0 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 D0 value of 1NpOH·adamantane is close to that of 1NpOH·cyclohexane, reflecting the nearly identical contact layer between the two molecules.

Benchmark Experimental Gas-Phase Intermolecular Dissociation Energies by the SEP-R2PI Method.

The gas-phase ground-state dissociation energy D0(S0) of an isolated and cold bimolecular complex is a fundamental measure of the intermolecular interaction strength between its constituents. Accurate D0 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 D0(S0) and D0(S1) 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 (S0 → S1) electronic excitation, followed by S1 → S0 stimulated emission spanning a range of S0 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.

Excited-state vibrations, lifetimes, and nonradiative dynamics of jet-cooled 1-ethylcytosine.

The S1 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 S0 → S1 vibronic spectrum using two-color resonant two-photon ionization and UV/UV holeburning spectroscopy. Compared to cytosine and 1-methylcytosine, the S0 → S1 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 S1 state nonradiative dynamics at ∼20 ps resolution by the pump/delayed ionization technique. The S1 lifetime increases with the number of ν1 ' quanta from τ = 930 ps at v1 '=0 to 1030 ps at v1 '=2, decreasing to 14 ps at 710 cm-1 vibrational energy. We measured the rate constants for S1 ⇝ S0 internal conversion and S1 ⇝ T1 intersystem crossing (ISC): At the v' = 0 level, kIC is 8 × 108 s-1 or three times smaller than 1-methylcytosine. The ISC rate constant from v' = 0 to the T1(3ππ*) state is kISC = 2.4 × 108 s-1, 10 times smaller than the ISC rate constants of cytosine, but similar to that of 1-methylcytosine. Based on the calculated S1(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 S1 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 T1(ππ*) state formed via ISC reveals that it lies 3.2-3.4 eV above the S0 ground state.

Face, Notch, or Edge? Intermolecular dissociation energies of 1-naphthol complexes with linear molecules.

The stimulated-emission-pumping/resonant 2-photon ionization (SEP-R2PI) method was used to determine the intermolecular dissociation energies D0 of jet-cooled 1-naphthol(1NpOH)·S complexes, where S is a linear molecule (N2, CO, CO2, OCS, N2O, and ethyne) or symmetric-top molecule (2-butyne) that contains double or triple bonds. The dissociation energies D0(S0) are bracketed as follows: 6.68 ± 0.08 kJ/mol for S=N2, 7.7 ± 0.8 kJ/mol for CO, 12.07 ± 0.10 kJ/mol for CO2, 13.03 ± 0.01 kJ/mol for N2O, 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 N2 and CO are dispersively bound Face complexes (S bound to a naphthalene Face), while CO2, N2O, 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.

Intermolecular dissociation energies of hydrogen-bonded 1-naphthol complexes.

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.

Transition from Water Wires to Bifurcated H-Bond Networks in 2-Pyridone·(H2O) n, n = 1-4 Clusters.

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·(H2O) 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 H2O 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 H2O molecule, lowering the NH stretch to ∼2800 cm-1. The n = 4 cluster is also bifurcated with two water wires between the bifurcating H2O 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.

Intermolecular dissociation energies of 1-naphthol·n-alkane complexes.

Using the stimulated-emission-pumping/resonant 2-photon ionization (SEP-R2PI) method, we have determined accurate intermolecular dissociation energies D0 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 D0(S0) 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 D0(S0) = 7.98 ± 0.55 kJ/mol (±7%) for S = methane and 14.5 ± 0.28 kJ/mol (±2%) for S = ethane. The D0 values increase by only 1% upon S0 → S1 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 S0 → S1 electronic origins, which range from +0.5 to -15 cm-1. Agreement with the calculated dissociation energies D0(S0) 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 D0(S0) values for the larger n-alkanes by up to 20%. The agreement of the B3LYP-D3 D0 values is intermediate between the other two methods.

Intermolecular dissociation energies of dispersively bound complexes of aromatics with noble gases and nitrogen.

We measured accurate intermolecular dissociation energies D0 of the supersonic jet-cooled complexes of 1-naphthol (1NpOH) with the noble gases Ne, Ar, Kr, and Xe and with N2, using the stimulated-emission pumping resonant two-photon ionization method. The ground-state values D0(S0) for the 1NpOH⋅S complexes with S= Ar, Kr, Xe, and N2 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⋅N2 complex, and 6.62 ± 0.22 kJ/mol for the corresponding isomer 2. For S = Ne, the upper limit is D0 < 3.36 kJ/mol. The dissociation energies increase by 1%-5% upon S0 → S1 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 D0(S0) of the face complexes with Ar, Kr, and N2 agree with the experimental values within <5%, but the D0(S0) for Xe is 12% too low. The B3LYP-D3 and ωB97X-D calculated D0(S0) values exhibit larger deviations to both larger and smaller dissociation energies. For comparison to 1-naphthol, we calculated the D0(S0) of the carbazole complexes with S = Ne, Ar, Kr, Xe, and N2 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.

Accurate gas-phase structure of para-dioxane by fs Raman rotational coherence spectroscopy and ab initio calculations.

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-h8 and -d8 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-h8 is strongly asymmetric, leading to a large number of asymmetry transients in its RCS spectrum. In contrast, the d8-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-h8 to an asymmetric-top model yields the ground-state rotational constants A0 = 5084.4(5) MHz, B0 = 4684(1) MHz, C0 = 2744.7(8) MHz, and (A0 + B0)/2 = 4884.5(7) MHz (±1σ). The analogous values for p-dioxane-d8 are A0 = 4083(2) MHz, B0 = 3925(4) MHz, C0 = 2347.1(6) MHz, and (A0 + B0)/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 A0calc(X),B0calc(X),C0calc(X) for increasing basis-set size X with non-linear extrapolation to the experimental constants A0exp,B0exp,C0exp allows to determine the equilibrium ground state structure of p-dioxane. For instance, the equilibrium C-C and C-O bond lengths are re(CC) = 1.5135(3) Å and re(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.

Planarizing cytosine: The S1 state structure, vibrations, and nonradiative dynamics of jet-cooled 5,6-trimethylenecytosine.

We measure the S0 → S1 spectrum and time-resolved S1 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 000 band, that of TMCyt extends up to +4400 cm-1 higher, with over a hundred resolved vibronic bands. Thus, clamping the cytosine C5-C6 bond allows us to explore the S1 state vibrations and S0 → S1 geometry changes in detail. The TMCyt S1 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. S0 → S1 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 S1 ⇝ S0 conical intersection (CI) increases from +366 cm-1 in cytosine to >6000 cm-1 in TMCyt, explaining the long lifetime and extended S0 → S1 spectrum. The lowest-energy S1 ⇝ S0 CI of TMCyt is the "amino out-of-plane" (OPX) intersection, calculated at +4190 cm-1. The experimental S1 ⇝ S0 internal conversion rate constant at the S1(v'=0) level is kIC=0.98-2.2⋅108 s-1, which is ∼10 times smaller than in 1-methylcytosine and cytosine. The S1(v'=0) level relaxes into the T1(3ππ*) state by intersystem crossing with kISC=0.41-1.6⋅108 s-1. The T1 state energy is measured to lie 24 580±560 cm-1 above the S0 state. The S1(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 S1 ionization threshold is observed for amino-keto TMCyt, yielding an adiabatic ionization energy of 8.114±0.002 eV.

Measuring Intermolecular Binding Energies by Laser Spectroscopy.

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.

Excitonic Splitting and Vibronic Coupling Analysis of the m-Cyanophenol Dimer.

The S1/S2 splitting of the m-cyanophenol dimer, (mCP)2 and the delocalization of its excitonically coupled S1/S2 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 (Ci-symmetric); the second-lowest is <5 cm-1 higher and planar (C2h). The vibrational ground state is probably delocalized over both minima. The S0 → S1 transition of (mCP)2 is electric-dipole allowed (Ag → Au), while the S0 → S2 transition is forbidden (Ag → Ag). Breaking the inversion symmetry by 12C/13C- or H/D-substitution renders the S0 → S2 transition partially allowed; the excitonic contribution to the S1/S2 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(12C/13C) = 3.3 cm-1 and Δiso(H/D) = 6.8 cm-1. Only partial localization of the exciton occurs in the 12C/13C and H/D substituted heterodimers. The SCS-CC2 calculated excitonic splitting is Δel = 179 cm-1; when multiplying this with the vibronic quenching factor Γvibronexp = 0.043, we obtain an exciton splitting Δvibronexp = 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 S1/S2 states from the v = 0 level to 600 cm-1 above. Both linear and quadratic vibronic coupling schemes were investigated to simulate the S0 → S1/S2 vibronic spectra; those calculated with the latter scheme agree better with experiment.

Accurate dissociation energies of two isomers of the 1-naphthol⋠cyclopropane complex.

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 D0 were measured with the SEP-R2PI technique. The D0(S0) 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 D0(S1) were evaluated using the respective electronic spectral shifts. Despite the D0(S0) 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, De, and dissociation energies, D0, 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.

The excited-state structure, vibrations, lifetimes, and nonradiative dynamics of jet-cooled 1-methylcytosine.

We have investigated the S0 → S1 UV vibronic spectrum and time-resolved S1 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 S1 state. The S0 → S1 vibronic spectrum breaks off at a vibrational excess energy Eexc ∼ 500 cm-1, indicating that a barrier in front of the ethylene-type S1⇝S0 conical intersection is exceeded, which is calculated to lie at Eexc = 366 cm-1. The S1⇝S0 internal conversion rate constant increases from kIC = 2 ⋅ 109 s-1 near the S1(v = 0) level to 1 ⋅ 1011 s-1 at Eexc = 516 cm-1. The 1ππ∗ state of 1MCyt also relaxes into the lower-lying triplet T1 (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 kIC; it increases from kISC = 2 ⋅ 108 s-1 near S1(v = 0) to kISC = 2 ⋅ 109 s-1 at Eexc = 516 cm-1. The T1 state energy is determined from the onset of the time-delayed photoionization efficiency curve as 25 600 ± 500 cm-1. The T2 (3nπ∗) state lies >1500 cm-1 above S1(v = 0), so S1⇝T2 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 S0 → S1 spectrum. The effect of methylation on the radiationless decay to S0 and ISC to T1 is small, as shown by the similar break-off of the spectrum and the similar computed mechanisms.

Do Hydrogen Bonds Influence Excitonic Splittings?

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.

Experimental and Theoretical Determination of Dissociation Energies of Dispersion-Dominated Aromatic Molecular Complexes.

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.

Gas-Phase Cytosine and Cytosine-N1-Derivatives Have 0.1-1 ns Lifetimes Near the S1 State Minimum.

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.

Intermolecular dissociation energies of dispersively bound 1-naphthol⋠cycloalkane complexes.

Intermolecular dissociation energies D0(S0) 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 D0(S0) 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 S0→S1 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 D0(S0) 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 D0(S0) values. Without dispersion correction, all the complexes are calculated to be unbound.

Intersystem crossing rates of S1 state keto-amino cytosine at low excess energy.

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.