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 timeresolved 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 mechanismsThis research has been supported by the Schweiz. Nationalfonds (Grant Nos. 121993 and 132540), the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR) from Catalonia (Spain) (Grant No. 2014SGR1202), the Ministerio de Economía y Competividad (MINECO) from Spain (Grant No. CTQ2015-69363-P), and the National Natural Science Foundation of China (Grant No. 21303007

Low-lying excited states and nonradiative processes of 9-methyl-2-aminopurine

he 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 0A′′1 (l = 0) and 1E ″ (l = ±1) torsional levels. These and further torsional bands that appear up to 000+230 cm−1 allow to fit the threefold (V 3) barriers of the torsional potentials as ∣∣V′′3∣∣=50 cm−1 in the S 0 and ∣∣V′3∣∣=126 cm−1 in the S 1 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 V′′3=20 cm−1 and V′3=115 cm−1. The 000 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 000 contour of 9M-2AP, the 1ππ* lifetime is τ ⩾ 120 ps, reflecting a rapid nonradiative transition

Excited-State Structure, Vibrations, and Nonradiative Relaxation of Jet-Cooled 5-Fluorocytosine

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 000 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 (knr ≥ 1012 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

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

NH···π hydrogen bonds occur frequently between the amino acid side groups in proteins and peptides. Data-mining studies of protein crystals find that ~80% of the T-shaped histidine···aromatic contacts are CH···π, and only ~20% are NH···π interactions. We investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet-cooled imidazole·benzene (Im·Bz) complex as a model for the NH···π interaction between histidine and phenylalanine. Ground- and excited-state dispersion-corrected density functional calculations and correlated methods (SCS-MP2 and SCS-CC2) predict that Im·Bz has a Cs-symmetric T-shaped minimum-energy structure with an NH···π hydrogen bond to the Bz ring; the NH bond is tilted 12° away from the Bz C₆ axis. IR depletion spectra support the T-shaped geometry: The NH stretch vibrational fundamental is red shifted by −73 cm⁻¹ relative to that of bare imidazole at 3518 cm⁻¹, indicating a moderately strong NH···π interaction. While the Sₒ(A1g) → S₁(B₂u) origin of benzene at 38 086 cm⁻¹ is forbidden in the gas phase, Im·Bz exhibits a moderately intense Sₒ → S₁ origin, which appears via the D₆h → 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⁻¹). 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)₂ cluster are observed via fragmentation into the Im·Bz+ mass channel. The spectra of Im·Bz and Im·Bz₂ are cleanly separable by IR hole burning. The UV spectrum of Im·Bz₂ exhibits two 000 bands corresponding to the Sₒ → S₁ excitations of the two inequivalent benzenes, which are symmetrically shifted by −86/+88 cm⁻¹ relative to the 000 band of benzene

Accurate computations of the structures and binding energies of the imidazole ... benzene and pyrrole ... benzene complexes

Using explicitly-correlated coupled-cluster theory with single and double excitations, the intermolecular distances and interaction energies of the T-shaped imidazole⋯⋯benzene and pyrrole⋯⋯benzene complexes have been computed in a large augmented correlation-consistent quadruple-zeta basis set, adding also corrections for connected triple excitations and remaining basis-set-superposition errors. The results of these computations are used to assess other methods such as Møller–Plesset perturbation theory (MP2), spin-component-scaled MP2 theory, dispersion-weighted MP2 theory, interference-corrected explicitly-correlated MP2 theory, dispersion-corrected double-hybrid density-functional theory (DFT), DFT-based symmetry-adapted perturbation theory, the random-phase approximation, explicitly-correlated ring-coupled-cluster-doubles theory, and double-hybrid DFT with a correlation energy computed in the random-phase approximation